open TELEMAC-MASCARET The mathematically superior suite of solvers

Could someone help me to add wind wave in my present TELEMAC model.

I have TPXO driving on the seaside boundary and water discharge at river boundary.
The present model works very well, but I don't know how to add wind to my present model. I am not familar with FORTRAN at the moment, could someone edit the steering file and Fortran file for me.

Thanks in advance!

Attached is my present steering file, fortran file and discharge file in the version of v6p3. The desired wind file is also attached.
The geometry file and boundary condition file is a bit big, if necessary I can send separately.

Zhang Min

Following is the steering file


/

---

/
FORTRAN FILE : 'princi3d_eastchinasea01_v6p3.f'
GEOMETRY FILE : './geom1958_friction.slf'
BOUNDARY CONDITIONS FILE : './geom1958_friction.cli'
LIQUID BOUNDARIES FILE : './river_flood_ave.txt'
3D RESULT FILE : 'res3d_1958YE_09tide_flood_spring.slf'
2D RESULT FILE : 'res2d_1958YE_09tide_flood_spring.slf'
/Version 6.2
/

---

/
TITLE : '6. TPXO tidal flow. Tracers.'
/

---

/

VARIABLES FOR 2D GRAPHIC PRINTOUTS : 'U,V,S,H'
/VARIABLES FOR 3D GRAPHIC PRINTOUTS : 'Z,U,V,W,TA1,TA2'
VARIABLES FOR 3D GRAPHIC PRINTOUTS : 'Z,U,V,W'
/

---

/
TIME STEP = 50.0
NUMBER OF TIME STEPS = 25056/50112
GRAPHIC PRINTOUT PERIOD = 18/15minutes a result
LISTING PRINTOUT PERIOD = 30
MASS-BALANCE : YES
/

---

/
NUMBER OF HORIZONTAL LEVELS = 5
/

---

/
INITIAL CONDITIONS : 'TPXO SATELLITE ALTIMETRY'
INITIAL ELEVATION = 0.0
/

---

/
/ PRESCRIBED TRACERS VALUES = 0.0; 0.0; 1.0; 2.0
/

---

/
LAW OF BOTTOM FRICTION = 5 /Nikuradse
/FRICTION COEFFICIENT FOR THE BOTTOM = 0.001
HORIZONTAL TURBULENCE MODEL : 1
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES = 0.1
VERTICAL TURBULENCE MODEL : 2
MIXING LENGTH MODEL : 3 /Nezu and Nakagawa
TURBULENCE MODEL FOR THE BOTTOM : 2
DAMPING FUNCTION : 3 /Munk and Anderson
COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES = 1.E-6
/

---

/
CORIOLIS : YES
CORIOLIS COEFFICIENT = 7.66E-5
/

---

/
IMPLICITATION FOR DEPTH = 0.55
IMPLICITATION FOR VELOCITIES = 0.55
MAXIMUM NUMBER OF ITERATIONS FOR DIFFUSION OF VELOCITIES = 100 / default is 60
/

---

/
TIDAL FLATS : YES
OPTION FOR THE TREATMENT OF TIDAL FLATS : 1
TREATMENT OF NEGATIVE DEPTHS : 2
MASS-LUMPING FOR DEPTH = 1.0
MATRIX STORAGE : 3
/

---

/
ORIGINAL DATE OF TIME = 2009;8;10
ORIGINAL HOUR OF TIME = 0;0;0
GEOGRAPHIC SYSTEM : 2
ZONE NUMBER IN GEOGRAPHIC SYSTEM = 51
TIDAL DATA BASE : 2
OPTION FOR TIDAL BOUNDARY CONDITIONS : 1
/

---

/
BINARY DATABASE 1 FOR TIDE : '../tpxo/h_tpxo7.2'
BINARY DATABASE 2 FOR TIDE : '../tpxo/u_tpxo7.2'
/

---

/
PRESCRIBED FLOWRATES = 0.0; 0.0
PRESCRIBED ELEVATIONS = 0.0; 0.0
VELOCITY PROFILES = 4;4
/

---

/
&ETA
&FIN



Following is the wind file



TIME FUAIR FVAIR
s m/s m/s
0. 1. 0.
86400. 2. 0.
172800. 3. 0.
259200. 4. 0.
345600. 5. 0.
432000. 6. 0.
518400. 7. 0.
604800. 5. 0.
691200. 10. 0.
2505600. 10. 0.