Modèle de Boris Dewitte; Page web d' Alexandre Ganachaud.

• Characteristics

• 124E to 80W by 2 degree; 28.75S to 28.75N by 0.5 degree

• Parameters (see Dewitte 2000)

  Ocean	m1	m2	m3
  Phase speed (m/s)	2.7	1.6	1.1
  Proj Coef Pn	0.6	0.6	0.1
  Friction rn (months)=cn^-0.5	30	23	18
  rho_a Cd (kg/m3)	0.0018
  Coupling efficiency	1
  Mixed layer
  Efficiency factor ga(x)	0.5	to	1.0
  Coeff Tsub	func(x)
  Thermal damp alpha_s	100	1/day
  Friction rs	2	1/day
  Atmosphere
  Wind div. feedback	no
  Alpha (m2/s3/deg)	0.043

• • Note that efficiency factor=Wentrainment/Wupwelling
  • Time step=10 days (5/15/25)
  • Mixed layer depth: 50m
  • Parameters in zfsuf.dat

• Thermocline depth: 150m

• Directory contents
• BDUNCOUP:
• DATA: wind forcing and SSTobs (info and infosst)
• ersfsumano.84302x480: (84x30x2x480)=124E to 28? by 2deg; 29S to 29N by 2deg, first table is Tx, second is Ty. Example:
  parameter(nx=84,ny=30,nt=480)
dimension tab(nx,ny,2)
do k=1,nt
read(10) tab
tab(x,x,1)=Tx
tab(x,x,2)=Ty
enddo

• CLIM: (binary) seasonal cycle: rc=Rasmunsen and Carpenter; .zeb=Zebiak; .bd=Dewitte; rcdiv not used.
• CLIMF: (ascii) formatted climatology
• PARAM: parameters: beta, betaf, mtau, msst, stddev, pn
• RUN.UNC/STD:
• Initial Conditions (out1970-2001.std)
• Output files
• same as for uncoupled
• See read_icmop.m
• BDCOUP:
• DATA: normally not used but can read wind for nudging
• CLIM: same a UNCOUP
• CLIMF: formatted climatology
• PARAM: parameters: beta, betaf, mtau, msst, stddev, pn
• RUN.UNC/STD:
• Initial Conditions (out1970-2001.std)
• Output files  (later use ALWAYS i=longitude (1..34); j=latitude (1..30)
• SST(30,34) 29N to 29S (inverted) by 2 deg; 101.25E to 286.?? by 5.625 deg
• H1(30,34)-thermocline for each mode, as SST
• H1m(30,34,nmodes)- same but decomposed into modes
• SL-sea level (30,34) as SST
• HTAU(34,30)-taux/tauy, here 29S to 29N = tau(t+1) calcule par le modele atmospherique

• Procedure
• make creates zfsuf
• zfsuf.dat (see setup.f) and  zfsuf.dat_spinup for spinup
• RUN.UNC/STD/toto: takes ICs
• RUN.UNC/STD/out61-70.std to restart if stopped
• DATA/ersfsumano.84302x480: wind
• PARAM/betaf: beta atmosphere NOT USED
• PARAM/pn844f: wind proj. coef. as function of x. NOT USED sensitivity other wize manually specified in zfsuf.dat
• NSTART =  0 (spinup) = 2 (ICs)
• TFIND=108.5: where to start off in the restart file. time to find in out61_70_std, relative to Jan70
• TZERO=0.5  Means that 15/01/1961 is t0 then TFIND=108.5 (months) is 15/01/1970 (DT=0.33333m)
• TENDD=360.5 when to stop.
• PNN =1 takes values from PARAM/pn844f
• NMSTAT =1: uses statistical atmosphere
• NREWND =  11: rewnd overwrites files (only last saved); =10: save all time steps
• run zfsuf < zfsuf.dat
• Programs
• Sous-Répertoire BDUNCOUP

• akcalc.fcomputes Ak using characteristics at every point
• amach.f cfunctions i1mach, d1mach, r1mach, amach, amtest,amsub1
• atmosta.fstatistic atmosphere
• bndary.fu,v,h at eastern and western boundaries
• constc.f c
• constmv.f computes parameters of considered vertical mode
• erfin.f c
• ermsg.f error message
• fft2c.f fast Fourier transform
• forcfsuf.f force the coupled model with fsusm4yr
• ggubfs.f c
• ierm1.f c
• ierv1.f c
• mloop.f computes HBn, HBn and Vn then H1n, U1n, V1n
cycles through the processes of the model until tendd
• nrdhist.f write output (to modify if want formatted)
• openfl.f open files (to modify if want formatted)
• projwind.f compute projection of the wind on each vertical mode then horizontal
• qqnqf.f normal random deviate generator (subroutines mdnris,merfi, uertst,uspkd,ugetio)
• readwind.f c
• scsevl.f Chebyshev series
• serf.f 1-erf function
• serm1.f c
• serv1.f c
• setup.f gives info about parameters in the model; writes cr.out and creates /COMMON/ZDATOFZATMC/
• sinits.f initializes the orthogonal series
• tridag.f solves the tridiagonal system
• uhcalc.f produces u,h,v for every time step but the first
• uhinit.f produces u,h,v for first time step
• zfsufm.f main program
• ztmfsuf.f subroutines:
• tinit
• uvinit
• zavg
• summodv
• stress
• trid

• Variables
• WM(j,i,2,12): seasonal climatology for winds
• Algorithm

• For each mode:

1. compute parameters corresponding to the considered vertical mode for each mode
2. compute the projection of the wind on each vertical mode then on K and R horizontal modes
3. compute HBn, Vn then H1n, U1n, V1n (n=next)
1. Ak
2. U and H
3. Ak at western boundary
4. compute oceanic field of vertical mode on coarse grid.
5. <OU EST LA COUCHE DE MELANGE/SST??????????>
6. <Procedure de calcul ?? >