Changeset 3959


Ignore:
Timestamp:
2013-07-08T18:46:13+02:00 (7 years ago)
Author:
agn
Message:

Gurvan's FK + separate rn_rho_c_mle

Location:
branches/2013/dev_r3948_NOC_FK/NEMOGCM
Files:
1 added
7 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/GYRE/EXP00/iodef.xml

    r3940 r3959  
    4747     <field field_ref="mldr10_1"     name="somxl010"  /> 
    4848     <field field_ref="mldkz5"       name="somixhgt"  /> 
    49         </file> 
     49           <!-- variables available with MLE --> 
     50          <field field_ref="Lf_NHpf"    name="Lf_NHpf"  long_name="MLE:_Lf=NH/f"     /> 
     51       </file> 
    5052    
    5153   <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 
    5254          <field field_ref="uoce"         name="vozocrtx"  /> 
    5355          <field field_ref="utau"         name="sozotaux"  /> 
     56          <!-- variables available with MLE --> 
     57          <field field_ref="psiu_mle"     name="psiu_mle"    long_name="MLE_streamfunction_along_i-axis" /> 
    5458        </file> 
    5559    
     
    5761          <field field_ref="voce"         name="vomecrty"  />  
    5862          <field field_ref="vtau"         name="sometauy"  />  
     63          <!-- variables available with MLE --> 
     64          <field field_ref="psiv_mle"     name="psiv_mle"    long_name="MLE_streamfunction_along_j-axis" /> 
    5965        </file> 
    6066    
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/GYRE/EXP00/namelist

    r3795 r3959  
    586586   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    587587   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     588/ 
     589!----------------------------------------------------------------------- 
     590&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
     591!----------------------------------------------------------------------- 
     592   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     593   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     594   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     595   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
     596   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
     597   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     598   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
     599   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE 
     600   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
    588601/ 
    589602!---------------------------------------------------------------------------------- 
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef_default.xml

    r3940 r3959  
    6464     <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           /> 
    6565     <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          /> 
     66          <!-- variables available with MLE --> 
     67          <field field_ref="Lf_NHpf"    name="Lf_NHpf"  long_name="MLE:_Lf=NH/f"     /> 
    6668     <field field_ref="taum"         /> 
    6769     <field field_ref="mldkz5"       /> 
     
    7375     <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
    7476     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
     77          <!-- variables available with MLE --> 
     78          <field field_ref="psiu_mle"     name="psiu_mle"    long_name="MLE_streamfunction_along_i-axis" /> 
    7579   </file> 
    7680    
     
    7983     <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
    8084     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
     85          <!-- variables available with MLE --> 
     86          <field field_ref="psiv_mle"     name="psiv_mle"    long_name="MLE_streamfunction_along_j-axis" /> 
    8187   </file> 
    8288    
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist

    r3795 r3959  
    581581   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    582582   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     583/ 
     584 !----------------------------------------------------------------------- 
     585&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
     586!----------------------------------------------------------------------- 
     587   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     588   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     589   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     590   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
     591   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
     592   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     593   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
     594   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE 
     595   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
    583596/ 
    584597!---------------------------------------------------------------------------------- 
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/SHARED/field_def.xml

    r3905 r3959  
    2929         <field id="mldr10_1"     long_name="Mixed Layer Depth 0.01 ref.10m"            unit="m"                        /> 
    3030         <field id="rhop"         long_name="potential density (sigma0)"                unit="kg/m3" grid_ref="grid_T_3D"/> 
     31         <!-- variables available with MLE --> 
     32         <field id="Lf_NHpf"      long_name="MLE: Lf = N H / f"                         unit="m"                        /> 
    3133         <!-- next variables available with key_diahth --> 
    3234         <field id="mlddzt"       long_name="Thermocline Depth (max dT/dz)"             unit="m"                        /> 
     
    140142         <field id="uoce"         long_name="ocean current along i-axis"                  unit="m/s"  grid_ref="grid_U_3D" /> 
    141143         <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"      unit="m3/s" grid_ref="grid_U_3D" /> 
     144         <!-- variables available with MLE --> 
     145         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"             unit="m3/s"  grid_ref="grid_U_3D"  /> 
    142146         <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    143147         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
     
    159163         <field id="voce"         long_name="ocean current along j-axis"                  unit="m/s"  grid_ref="grid_V_3D" /> 
    160164         <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"      unit="m3/s" grid_ref="grid_V_3D" /> 
     165         <!-- variables available with MLE --> 
     166         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"             unit="m3/s"  grid_ref="grid_U_3D"  /> 
    161167         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    162168         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"              unit="m/s"  grid_ref="grid_V_3D" /> 
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r3718 r3959  
    66   !! History :  2.0  !  2005-11  (G. Madec)  Original code 
    77   !!            3.3  !  2010-09  (C. Ethe, G. Madec)  merge TRC-TRA + switch from velocity to transport 
     8   !!            4.0  !  2011-06  (G. Madec)  Addition of Mixed Layer Eddy parameterisation 
    89   !!---------------------------------------------------------------------- 
    910 
     
    2122   USE traadv_qck      ! QUICKEST scheme           (tra_adv_qck    routine) 
    2223   USE traadv_eiv      ! eddy induced velocity     (tra_adv_eiv    routine) 
     24   USE traadv_mle      ! ML eddy induced velocity  (tra_adv_eiv    routine) 
    2325   USE cla             ! cross land advection      (cla_traadv     routine) 
    2426   USE ldftra_oce      ! lateral diffusion coefficient on tracers 
     
    99101         &              CALL tra_adv_eiv( kt, nit000, zun, zvn, zwn, 'TRA' )    ! add the eiv transport (if necessary) 
    100102      ! 
     103      IF( ln_mle    )   CALL tra_adv_mle( kt, nit000, zun, zvn, zwn, 'TRA' )    ! add the mle transport (if necessary) 
    101104      CALL iom_put( "uocetr_eff", zun )                                         ! output effective transport       
    102105      CALL iom_put( "vocetr_eff", zvn ) 
     
    206209      ENDIF 
    207210      ! 
     211      CALL tra_adv_mle_init          ! initialisation of the Mixed Layer Eddy parametrisation (MLE) 
     212      ! 
    208213   END SUBROUTINE tra_adv_init 
    209214 
  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r3294 r3959  
    2424 
    2525   PUBLIC   zdf_mxl       ! called by step.F90 
     26 
     27   REAL(wp), PUBLIC ::   rho_c = 0.01_wp    ! density criterion for mixed layer depth 
     28   REAL(wp), PUBLIC ::   avt_c = 5.e-4_wp   ! Kz criterion for the turbocline depth 
    2629 
    2730   INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   nmln    !: number of level in the mixed layer (used by TOP) 
     
    6366      !! ** Method  :   The mixed layer depth is the shallowest W depth with  
    6467      !!      the density of the corresponding T point (just bellow) bellow a 
    65       !!      given value defined locally as rho(10m) + zrho_c 
     68      !!      given value defined locally as rho(10m) + rho_c 
    6669      !!               The turbocline depth is the depth at which the vertical 
    6770      !!      eddy diffusivity coefficient (resulting from the vertical physics 
     
    7679      INTEGER  ::   iikn, iiki          ! temporary integer within a do loop 
    7780      INTEGER, POINTER, DIMENSION(:,:) ::   imld                ! temporary workspace 
    78       REAL(wp) ::   zrho_c = 0.01_wp    ! density criterion for mixed layer depth 
    79       REAL(wp) ::   zavt_c = 5.e-4_wp   ! Kz criterion for the turbocline depth 
    8081      !!---------------------------------------------------------------------- 
    8182      ! 
     
    9899         DO jj = 1, jpj 
    99100            DO ji = 1, jpi 
    100                IF( rhop(ji,jj,jk) > rhop(ji,jj,nla10) + zrho_c )   nmln(ji,jj) = jk      ! Mixed layer 
    101                IF( avt (ji,jj,jk) < zavt_c                     )   imld(ji,jj) = jk      ! Turbocline  
     101               IF( rhop(ji,jj,jk) > rhop(ji,jj,nla10) + rho_c )   nmln(ji,jj) = jk      ! Mixed layer 
     102               IF( avt (ji,jj,jk) < avt_c                     )   imld(ji,jj) = jk      ! Turbocline  
    102103            END DO 
    103104         END DO 
Note: See TracChangeset for help on using the changeset viewer.