Skip to content. | Skip to navigation

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

Personal tools

This is SunRain Plone Theme

Navigation

You are here: Home / Members / Alexandre Ganachaud / 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
  • 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