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 12912 for NEMO/branches/UKMO/NEMO_4.0.2_FKOSM/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2020-05-12T14:51:37+02:00 (4 years ago)
Author:
cguiavarch
Message:

OSMOSIS and Fox-Kemper changes (merged from NEMO_4.0.1_FKOSM)

Location:
NEMO/branches/UKMO/NEMO_4.0.2_FKOSM/cfgs/SHARED
Files:
3 edited
1 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.2_FKOSM/cfgs/SHARED/domain_def_nemo.xml

    r12276 r12912  
    1010     </domain> 
    1111 
    12       <domain id="1point" domain_ref="grid_T" > 
    13          <zoom_domain ibegin="139" jbegin="119" ni="1" nj="1"/> 
    14       </domain> 
    15  
    16  
    17      <!--   Eq section --> 
     12     <domain id="1point" domain_ref="grid_T" > 
     13       <zoom_domain ibegin="1" jbegin="1" ni="1" nj="1"/> 
     14     </domain> 
     15     <domain id="danger_T" domain_ref="grid_T" > 
     16       <zoom_domain id="danger_T"  ibegin="45" jbegin="196" ni="5" nj="5"/> 
     17     </domain> 
     18     <domain id="danger_U" domain_ref="grid_U" > 
     19       <zoom_domain id="danger_U"  ibegin="45" jbegin="196" ni="5" nj="5"/> 
     20     </domain> 
     21      <domain id="danger_V" domain_ref="grid_V" > 
     22       <zoom_domain id="danger_V"  ibegin="45" jbegin="196" ni="5" nj="5"/> 
     23     </domain> 
     24     <domain id="danger_W" domain_ref="grid_W" > 
     25       <zoom_domain id="danger_W"  ibegin="45" jbegin="196" ni="5" nj="5"/> 
     26    </domain> 
     27    <!--   Eq section --> 
    1828     <domain id="EqT" domain_ref="grid_T" > <zoom_domain id="EqT"/> </domain> 
    1929     <!--   TAO : see example above   --> 
  • NEMO/branches/UKMO/NEMO_4.0.2_FKOSM/cfgs/SHARED/field_def_nemo-oce.xml

    r12288 r12912  
    215215 
    216216      <field_group id="OSMOSIS_T" grid_ref="grid_T_2D"> 
     217        <field id="hml"                 long_name="mixed layr depth"                         unit="m"       /> 
     218        <field id="hbl"                 long_name="boundary layer depth"                     unit="m"       /> 
     219        <field id="dh"                  long_name="Pycnocline thickness"                     unit=" m"      /> 
     220        <field id="ibld"                long_name="index of boundary layer depth"            unit="#"       /> 
     221        <field id="imld"                long_name="index of mixed layer depth"            unit="#"       /> 
     222        <field id="zhbl"                long_name="boundary layer depth -grid"                     unit="m"       /> 
     223        <field id="zhml"                long_name="mixed layer depth - grid"                        unit="m"       /> 
     224        <field id="zdh"                 long_name="Pycnocline  depth - grid"                 unit=" m"      /> 
     225        <field id="zustke"              long_name="magnitude of stokes drift  at T-points"   unit="m/s"     /> 
     226        <field id="us_x"        long_name="i component of active Stokes drift"                      unit="m/s"     /> 
     227        <field id="us_y"        long_name="j component of active Stokes drift"                      unit="m/s"     /> 
     228        <field id="dstokes"             long_name="stokes drift  depth scale"                unit="m"       /> 
    217229        <field id="zwth0"               long_name="surface non-local temperature flux"       unit="deg m/s" /> 
    218230        <field id="zws0"                long_name="surface non-local salinity flux"          unit="psu m/s" /> 
    219         <field id="hbl"                 long_name="boundary layer depth"                     unit="m"       /> 
    220         <field id="hbli"                long_name="initial boundary layer depth"             unit="m"       /> 
    221         <field id="dstokes"             long_name="stokes drift  depth scale"                unit="m"       /> 
    222         <field id="zustke"              long_name="magnitude of stokes drift  at T-points"   unit="m/s"     /> 
    223231        <field id="zwstrc"              long_name="convective velocity scale"                unit="m/s"     /> 
     232        <field id="zustar"              long_name="friction velocity"                        unit="m/s"     /> 
    224233        <field id="zwstrl"              long_name="langmuir velocity scale"                  unit="m/s"     /> 
    225         <field id="zustar"              long_name="friction velocity"                        unit="m/s"     /> 
    226         <field id="zhbl"                long_name="boundary layer depth"                     unit="m"       /> 
    227         <field id="zhml"                long_name="mixed layer depth"                        unit="m"       /> 
     234        <field id="zvstr"               long_name="mixed velocity scale"                     unit="m/s"     /> 
     235        <field id="zla"                 long_name="langmuir number"                          unit="m/s"     /> 
    228236        <field id="wind_wave_abs_power" long_name="\rho |U_s| x u*^2"                        unit="mW"      /> 
    229237        <field id="wind_wave_power"     long_name="U_s \dot  tau"                            unit="mW"      /> 
    230238        <field id="wind_power"          long_name="\rho  u*^3"                               unit="mW"      /> 
    231239 
    232         <!-- extra OSMOSIS diagnostics --> 
     240       <!-- interior BL OSMOSIS diagnostics --> 
    233241        <field id="zwthav"              long_name="av turb flux of T in ml"                  unit="deg m/s" /> 
    234242        <field id="zt_ml"               long_name="av T in ml"                               unit="deg"     /> 
     243        <field id="zhol"                long_name="Hoenekker number"                         unit="#"       /> 
     244        <field id="zws_ent"            long_name="entrainment turb flux of S"                unit="10^-3 m/s" /> 
    235245        <field id="zwth_ent"            long_name="entrainment turb flux of T"               unit="deg m/s" /> 
    236         <field id="zhol"                long_name="Hoenekker number"                         unit="#"       /> 
    237         <field id="zdh"                 long_name="Pycnocline  depth - grid"                 unit=" m"      /> 
    238       </field_group> 
    239  
    240       <field_group id="OSMOSIS_W" grid_ref="grid_W_3D" operation="instant" > 
     246        <field id="zwb_ent"            long_name="entrainment turb flux of buoyancy"         unit="m^2/s^-3" /> 
     247  
     248        <field id="zdt_bl"             long_name="temperature jump at base of BL"                 unit="deg"      /> 
     249        <field id="zds_bl"             long_name="salinity jump at base of BL"                 unit="10^-3"      /> 
     250        <field id="zdb_bl"             long_name="buoyancy jump at base of BL"                 unit="m/s^2"      /> 
     251        <field id="zdu_bl"             long_name="u jump at base of BL"                       unit="m/s"      /> 
     252        <field id="zdv_bl"             long_name="v jump at base of BL"                       unit="m/s"      /> 
     253 
     254        <!-- extra OSMOSIS diagnostics for debugging --> 
     255       <field id="zsc_uw_1_0"       long_name="zsc u-momentum flux on T after Stokes"                       unit="m^2/s^2" /> 
     256        <field id="zsc_uw_1_f"       long_name="zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     257        <field id="zsc_vw_1_f"       long_name="zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     258        <field id="zsc_uw_2_f"       long_name="2nd zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     259        <field id="zsc_vw_2_f"       long_name="2nd zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
     260        <field id="zuw_bse"       long_name="base u-flux T-points"                          unit="m^2/s^2" /> 
     261        <field id="zvw_bse"       long_name="base v-flux T-points"                          unit="m^2/s^2" /> 
     262 
     263       <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     264         <field id="hmle"          long_name="OBL FK-layer thickness"                                     unit="m"        /> 
     265        <field id="mld_prof"              long_name="FK-layer depth index"                  unit="#" /> 
     266        <field id="zmld"          long_name="target FK-layer thickness"                                     unit="m"        /> 
     267        <field id="zwb_fk"          long_name="FK b-flux"                                     unit="m^2 s^-3"        /> 
     268        <field id="zwb_fk_b"          long_name="layer averaged FK b-flux"                 unit="m^2 s^-3"       /> 
     269        <field id="zdiff_mle"          long_name="max FK diffusivity in MLE"       unit=" 10^-4 m^2 s^-1"       /> 
     270        <field id="zvel_mle"          long_name="FK velocity scale in MLE"       unit=" m s^-1"       /> 
     271    </field_group> 
     272 
     273      <field_group id="OSMOSIS_W" grid_ref="grid_W_3D" > 
     274        <field id="zviscos"       long_name="BL viscosity"   unit="m^2/s" /> 
    241275        <field id="ghamt"       long_name="non-local temperature flux"                       unit="deg m/s" /> 
    242276        <field id="ghams"       long_name="non-local salinity flux"                          unit="psu m/s" /> 
    243277        <field id="zdtdz_pyc"   long_name="Pycnocline temperature gradient"                  unit=" deg/m"  /> 
    244       </field_group> 
     278        <field id="zdsdz_pyc"   long_name="Pycnocline salinity gradient"                  unit=" 10^-3/m"  /> 
     279        <field id="zdbdz_pyc"   long_name="Pycnocline buoyancy gradient"                  unit=" s^-2"  /> 
     280        <field id="zdudz_pyc"   long_name="Pycnocline u gradient"                  unit=" s^-2"  /> 
     281        <field id="zdvdz_pyc"   long_name="Pycnocline v gradient"                  unit=" s^-2"  /> 
     282 
     283        <!-- extra OSMOSIS diagnostics for debugging --> 
     284         <field id="ghamu_00"       long_name="initial non-local u-momentum flux"   unit="m^2/s^2" /> 
     285        <field id="ghamv_00"       long_name="initial non-local v-momentum flux"   unit="m^2/s^2" /> 
     286        <field id="ghamu_0"       long_name="after dstokes non-local u-momentum flux"   unit="m^2/s^2" /> 
     287        <field id="ghamu_f"       long_name="after Coriolis non-local u-momentum flux"   unit="m^2/s^2" /> 
     288        <field id="ghamv_f"       long_name="after Coriolis  non-local v-momentum flux"   unit="m^2/s^2" /> 
     289        <field id="ghamu_b"       long_name="after buoyancy added non-local u-momentum flux"   unit="m^2/s^2" /> 
     290        <field id="ghamv_b"       long_name="after buoyancy added  non-local v-momentum flux"  unit="m^2/s^2" /> 
     291        <field id="ghamu_1"       long_name="after entrainment non-local u-momentum flux"   unit="m^2/s^2" /> 
     292        <field id="ghamv_1"       long_name="after entrainment  non-local v-momentum flux"  unit="m^2/s^2" /> 
     293     </field_group> 
    245294 
    246295      <field_group id="OSMOSIS_U" grid_ref="grid_U_2D" > 
    247296        <field id="ghamu"       long_name="non-local u-momentum flux"   grid_ref="grid_U_3D" unit="m^2/s^2" /> 
    248         <field id="us_x"        long_name="i component of Stokes drift"                      unit="m/s"     /> 
    249       </field_group> 
     297       <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     298       <field id="zdtdx"          long_name="FK  T x-gradient"                                     unit=" deg C m^-1"        /> 
     299        <field id="zdsdx"          long_name="FK  S x-gradient"                                     unit=" 10^-3 m^-1"        /> 
     300        <field id="dbdx_mle"          long_name="FK  B x-gradient"                                     unit=" s^-2"        /> 
     301     </field_group> 
    250302 
    251303      <field_group id="OSMOSIS_V" grid_ref="grid_V_2D" > 
    252304        <field id="ghamv"       long_name="non-local v-momentum flux"   grid_ref="grid_V_3D" unit="m^2/s^2" /> 
    253         <field id="us_y"        long_name="j component of Stokes drift"                      unit="m/s"     /> 
     305        <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     306        <field id="zdtdy"          long_name="FK T y-gradient"                                     unit=" deg C m^-1"        /> 
     307        <field id="zdsdy"          long_name="FK S y-gradient"                                     unit=" 10^-3 m^-1"        /> 
     308        <field id="dbdy_mle"          long_name="FK B y-gradient"                                     unit=" s^-2"        /> 
    254309      </field_group> 
    255310       
     
    263318          <field id="emp_oce"      long_name="Evap minus Precip over ocean"         standard_name="evap_minus_precip_over_sea_water"                                     unit="kg/m2/s"   /> 
    264319          <field id="emp_ice"      long_name="Evap minus Precip over ice"           standard_name="evap_minus_precip_over_sea_ice"                                       unit="kg/m2/s"   /> 
    265           <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="1e-3/m2/s" /> 
    266           <field id="fmmflx"       long_name="Water flux due to freezing/melting"                                                                                        unit="kg/m2/s"   /> 
     320          <field id="saltflx"      long_name="Total downward salinity flux"             standard_name="total_downward_salinity_flux"                                   unit="1e-3/m2/s" /> 
     321          <field id="fmmflx"       long_name="Water flux due to freezing/melting"     standard_name="water_flux_due_to_freezing/melting"                  unit="kg/m2/s"   /> 
    267322          <field id="snowpre"      long_name="Snow precipitation"                   standard_name="snowfall_flux"                                                        unit="kg/m2/s"   /> 
    268323          <field id="runoffs"      long_name="River Runoffs"                        standard_name="water_flux_into_sea_water_from_rivers"                                unit="kg/m2/s"   /> 
     
    699754     <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"   unit="1e-3/s" /> 
    700755 
    701      <!-- --> 
     756     <!-- ln_zdfosm=T only (OSMOSIS-OBL) --> 
     757     <field id="ttrd_osm"      long_name="temperature-trend: OSM-OSBL non-local forcing"                             unit="degC/s" /> 
     758     <field id="strd_osm"      long_name="salinity   -trend: OSM-OSBL non-local forcing"                             unit="1e-3/s" /> 
     759 
     760 
     761    <!-- --> 
    702762     <field id="ttrd_dmp"      long_name="temperature-trend: interior restoring"        unit="degC/s" /> 
    703763     <field id="strd_dmp"      long_name="salinity   -trend: interior restoring"        unit="1e-3/s" /> 
     
    735795     <field id="strd_zdfp_e3t"     unit="1e-3/s * m"  >  strd_zdfp * e3t </field> 
    736796 
     797          <!-- ln_zdfosm=T only (OSMOSIS-OBL) --> 
     798     <field id="ttrd_osm_e3t"      long_name="temperature-trend: OSM-OSBL non-local forcing"                             unit="degC/s * m" >  ttrd_osm * e3t </field> 
     799     <field id="strd_osm_e3t"      long_name="salinity   -trend: OSM-OSBL non-local forcing"                             unit="1e-3/s * m" >  strd_osm * e3t </field> 
     800      
    737801     <!-- --> 
    738802     <field id="ttrd_dmp_e3t"      unit="degC/s * m"  >  ttrd_dmp * e3t </field> 
     
    742806     <field id="ttrd_npc_e3t"      unit="degC/s * m"  >  ttrd_npc * e3t </field> 
    743807     <field id="strd_npc_e3t"      unit="1e-3/s * m"  >  strd_npc * e3t </field> 
    744      <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3ts </field> 
    745      <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3ts </field> 
     808     <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3t_surf </field> 
     809     <!-- <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3t_surf </field> --> 
    746810     <field id="ttrd_qsr_e3t"      unit="degC/s * m"  >  ttrd_qsr * e3t </field> 
    747811     <field id="ttrd_bbc_e3t"      unit="degC/s * m"  >  ttrd_bbc * e3t </field> 
    748812 
    749813     <!-- OMIP  layer-integrated trends --> 
     814     <field id="sfd"      long_name="Total downward salt flux"                                                     unit="kg/(m^2 s)" > saltflux * 1026.0 * 0.001 </field> 
     815     <field id="wfo"      long_name="Total downward FW mass flux"                                                   unit="kg/(m^2 s)" > -empmr </field> 
     816 
    750817     <field id="ttrd_totad_li"    long_name="layer integrated heat-trend: total advection"         unit="W/m^2"     > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field> 
    751818     <field id="strd_totad_li"    long_name="layer integrated salt-trend: total advection"         unit="kg/(m^2 s)"    > strd_totad_e3t * 1026.0 * 0.001  </field> 
     819     <field id="ttrd_osm_li"    long_name="layer integrated heat-trend: non-local OSM"         unit="W/m^2"     > ttrd_osm_e3t * 1026.0 * 3991.86795711963 </field> 
     820     <field id="strd_osm_li"    long_name="layer integrated salt-trend: non-local OSM"         unit="kg/(m^2 s)"    > strd_osm_e3t * 1026.0 * 0.001  </field> 
    752821     <field id="ttrd_evd_li"      long_name="layer integrated heat-trend: EVD convection"          unit="W/m^2"    > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field> 
    753822     <field id="strd_evd_li"      long_name="layer integrated salt-trend: EVD convection"          unit="kg/(m^2 s)"  > strd_evd_e3t * 1026.0 * 0.001  </field> 
     823     <field id="ttrd_npc_li"      long_name="layer integrated heat-trend: NPC convection"          unit="W/m^2"    > ttrd_npc_e3t * 1026.0 * 3991.86795711963 </field> 
     824     <field id="strd_npc_li"      long_name="layer integrated salt-trend: NPC convection"          unit="kg/(m^2 s)"  > strd_npc_e3t * 1026.0 * 0.001  </field> 
    754825     <field id="ttrd_iso_li"      long_name="layer integrated heat-trend: isopycnal diffusion"     unit="W/m^2" > ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field> 
    755826     <field id="strd_iso_li"      long_name="layer integrated salt-trend: isopycnal diffusion"     unit="kg/(m^2 s)" > strd_iso_e3t * 1026.0 * 0.001  </field> 
     
    758829     <field id="ttrd_qns_li"      long_name="layer integrated heat-trend: non-solar flux + runoff" unit="W/m^2" grid_ref="grid_T_2D"> ttrd_qns_e3t * 1026.0 * 3991.86795711963 </field> 
    759830     <field id="ttrd_qsr_li"      long_name="layer integrated heat-trend: solar flux"              unit="W/m^2"  grid_ref="grid_T_3D"> ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field> 
     831     <field id="ttrd_bbc_li"      long_name="layer integrated heat-trend: geothermal heating "  unit="W/m^2" > ttrd_bbc_e3t * 1026.0 * 3991.86795711963 </field> 
    760832     <field id="ttrd_bbl_li"      long_name="layer integrated heat-trend: bottom boundary layer "  unit="W/m^2" > ttrd_bbl_e3t * 1026.0 * 3991.86795711963 </field> 
    761833     <field id="strd_bbl_li"      long_name="layer integrated salt-trend: bottom boundary layer "  unit="kg/(m^2 s)" > strd_bbl_e3t * 1026.0 * 0.001  </field> 
  • NEMO/branches/UKMO/NEMO_4.0.2_FKOSM/cfgs/SHARED/namelist_ref

    r12658 r12912  
    257257   ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    258258   ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     259 
     260   ln_humi_dpt    = .false.   !  Supply dewpoint tempearture instead of specific humidity (is true for ERA5) 
     261 
    259262      ! 
    260263      rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
     
    10881091&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T) 
    10891092!----------------------------------------------------------------------- 
    1090    ln_use_osm_la = .false.      !  Use namelist  rn_osm_la 
     1093   ln_use_osm_la = .false.     !  Use   rn_osm_la 
    10911094   rn_osm_la     = 0.3         !  Turbulent Langmuir number 
    10921095   rn_osm_dstokes     = 5.     !  Depth scale of Stokes drift (m) 
     
    11031106      !                        !  = 1: Pierson Moskowitz wave spectrum 
    11041107      !                        !  = 0: Constant La# = 0.3 
     1108   ln_osm_mle = .false.        !  Use integrated FK-OSM model 
     1109/ 
     1110!----------------------------------------------------------------------- 
     1111&namosm_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF) 
     1112!----------------------------------------------------------------------- 
     1113   rn_osm_mle_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     1114   nn_osm_mle          = 0         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     1115   rn_osm_mle_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_osm_mle=0) 
     1116   rn_osm_mle_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_osm_mle=0) 
     1117   rn_osm_mle_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     1118   rn_osm_mle_rho_c =    0.01      ! delta rho criterion used to calculate MLD for FK 
     1119   rn_osm_mle_thresh = 0.0005      ! delta b criterion used for FK MLD criterion 
     1120   rn_osm_mle_tau     = 172800.   ! time scale for FK-OSM (seconds)  (case rn_osm_mle=0) 
    11051121/ 
    11061122!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.