New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4188 – NEMO

Changeset 4188


Ignore:
Timestamp:
2013-11-14T10:40:55+01:00 (10 years ago)
Author:
agn
Message:

Divide zpsi_[uv]w by box width to give units in m2/s, correct grid for psiv_mle, update namelists

Location:
branches/2013/dev_r3948_NOC_FK/NEMOGCM
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/AMM12/EXP00/namelist

    r3948 r4188  
    550550   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    551551   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     552/ 
     553!----------------------------------------------------------------------- 
     554&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
     555!----------------------------------------------------------------------- 
     556   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     557   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     558   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     559   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
     560   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
     561   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     562   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
     563   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE 
     564   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
    552565/ 
    553566!---------------------------------------------------------------------------------- 

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

    r3959 r4188  
    4747     <field field_ref="mldr10_1"     name="somxl010"  /> 
    4848     <field field_ref="mldkz5"       name="somixhgt"  /> 
    49            <!-- variables available with MLE --> 
     49           <!-- variables available with MLE 
    5050          <field field_ref="Lf_NHpf"    name="Lf_NHpf"  long_name="MLE:_Lf=NH/f"     /> 
     51      --> 
    5152       </file> 
    5253    
    … …  
    5455          <field field_ref="uoce"         name="vozocrtx"  /> 
    5556          <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" /> 
     57          <!-- variables available with MLE 
     58          <field field_ref="psiu_mle"     name="psiu_mle"  long_name="MLE_streamfunction_along_i-axis" /> 
     59     --> 
    5860        </file> 
    5961    
    … …  
    6163          <field field_ref="voce"         name="vomecrty"  />  
    6264          <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" /> 
     65          <!-- variables available with MLE 
     66          <field field_ref="psiv_mle"     name="psiv_mle"  long_name="MLE_streamfunction_along_j-axis" /> 
     67     --> 
    6568        </file> 
    6669    

  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist

    r3815 r4188  
    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/GYRE_PISCES/EXP00/iodef.xml

    r3940 r4188  
    4747     <field field_ref="mldr10_1"     name="somxl010"  /> 
    4848     <field field_ref="mldkz5"       name="somixhgt"  /> 
     49           <!-- variables available with MLE 
     50          <field field_ref="Lf_NHpf"    name="Lf_NHpf"  long_name="MLE:_Lf=NH/f"     /> 
     51      --> 
    4952        </file> 
    5053    
    … …  
    5255          <field field_ref="uoce"         name="vozocrtx"  /> 
    5356          <field field_ref="utau"         name="sozotaux"  /> 
     57          <!-- variables available with MLE 
     58          <field field_ref="psiu_mle"     name="psiu_mle"  long_name="MLE_streamfunction_along_i-axis" /> 
     59     --> 
    5460        </file> 
    5561    
    … …  
    5763          <field field_ref="voce"         name="vomecrty"  />  
    5864          <field field_ref="vtau"         name="sometauy"  />  
     65          <!-- variables available with MLE 
     66          <field field_ref="psiv_mle"     name="psiv_mle"  long_name="MLE_streamfunction_along_j-axis" /> 
     67     --> 
    5968        </file> 
    6069    

  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef_default.xml

    r3959 r4188  
    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"     /> 
    6866     <field field_ref="taum"         /> 
    6967     <field field_ref="mldkz5"       /> 
    7068     <field field_ref="mldr10_1"     /> 
     69           <!-- variables available with MLE 
     70          <field field_ref="Lf_NHpf"    name="Lf_NHpf"  long_name="MLE:_Lf=NH/f"     /> 
     71      --> 
    7172   </file> 
    7273    
    … …  
    7576     <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
    7677     <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" /> 
     78          <!-- variables available with MLE 
     79          <field field_ref="psiu_mle"     name="psiu_mle"  long_name="MLE_streamfunction_along_i-axis" /> 
     80     --> 
    7981   </file> 
    8082    
    … …  
    8385     <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
    8486     <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" /> 
     87          <!-- variables available with MLE 
     88          <field field_ref="psiv_mle"     name="psiv_mle"  long_name="MLE_streamfunction_along_j-axis" /> 
     89     --> 
    8790   </file> 
    8891    

  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml

    r3905 r4188  
    6666     <field field_ref="erp"          name="sowafldp"  /> 
    6767     <field field_ref="ist_ipa"      name="soicetem"  /> 
    68      <field field_ref="icealb_cea"   name="soicealb"  />    
     68     <field field_ref="icealb_cea"   name="soicealb"  />  
     69           <!-- variables available with MLE 
     70          <field field_ref="Lf_NHpf"    name="Lf_NHpf"  long_name="MLE:_Lf=NH/f"     /> 
     71      --> 
    6972   </file> 
    7073    
    … …  
    7376     <field field_ref="uoce_eiv"     name="vozoeivu"  /> 
    7477     <field field_ref="utau"         name="sozotaux"  /> 
     78          <!-- variables available with MLE 
     79          <field field_ref="psiu_mle"     name="psiu_mle"  long_name="MLE_streamfunction_along_i-axis" /> 
     80     --> 
    7581   </file> 
    7682    
    … …  
    7985     <field field_ref="voce_eiv"     name="vomeeivv"  /> 
    8086     <field field_ref="vtau"         name="sometauy"  /> 
     87          <!-- variables available with MLE 
     88          <field field_ref="psiv_mle"     name="psiv_mle"  long_name="MLE_streamfunction_along_j-axis" /> 
     89     --> 
    8190   </file> 
    8291    

  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r3905 r4188  
    562562   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    563563   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     564/ 
     565!----------------------------------------------------------------------- 
     566&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
     567!----------------------------------------------------------------------- 
     568   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     569   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     570   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     571   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
     572   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
     573   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     574   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
     575   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE 
     576   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
    564577/ 
    565578!----------------------------------------------------------------------- 

  • branches/2013/dev_r3948_NOC_FK/NEMOGCM/CONFIG/SHARED/field_def.xml

    r3959 r4188  
    164164         <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"      unit="m3/s" grid_ref="grid_V_3D" /> 
    165165         <!-- variables available with MLE --> 
    166          <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"             unit="m3/s"  grid_ref="grid_U_3D"  /> 
     166         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"             unit="m3/s"  grid_ref="grid_V_3D"  /> 
    167167         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    168168         <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_mle.F90

    r3995 r4188  
    245245         CALL iom_put( "Lf_NHpf" , zLf_NH  )    ! Lf = N H / f 
    246246         ! 
    247          zpsim_u(jpi,:) = 0._wp   ;   zpsim_u(:,jpj) = 0._wp      ! set a value where it is undefined 
    248          zpsim_v(jpi,:) = 0._wp   ;   zpsim_v(:,jpj) = 0._wp      ! before the output 
     247         ! divide by cross distance to give streamfunction with dimensions m^2/s 
     248         DO jk = 1, ikmax+1 
     249            zpsi_uw(:,:,jk) = zpsi_uw(:,:,jk)/e2u(:,:) 
     250            zpsi_vw(:,:,jk) = zpsi_vw(:,:,jk)/e1v(:,:) 
     251         END DO 
    249252         CALL iom_put( "psiu_mle", zpsi_uw )    ! i-mle streamfunction 
    250253         CALL iom_put( "psiv_mle", zpsi_vw )    ! j-mle streamfunction 
Note: See TracChangeset for help on using the changeset viewer.