Laboratoire d’Etudes en Géophysique et Océanographie Spatiales

Personal tools

This is SunRain Plone Theme

Navigation

You are here: Home / icmboris

# Instructions ICM

by Alexandre Ganachaud — last modified Feb 05, 2013 11:31 AM

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.f`c``omputes Ak using characteristics at every point`
• amach.f` c``functions i1mach, d1mach, r1mach, amach, amtest,amsub1`
• atmosta.f`statistic atmosphere`
• bndary.f`u,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 ?? >

Document Actions