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 15574 for NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2021-12-03T20:32:50+01:00 (3 years ago)
Author:
techene
Message:

#2605 #2715 trunk merged into dev_r14318_RK3_stage1

Location:
NEMO/branches/2021/dev_r14318_RK3_stage1
Files:
12 edited
1 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14318_RK3_stage1

    • Property svn:externals
      •  

        old new  
        99 
        1010# SETTE 
        11 ^/utils/CI/sette@14244        sette 
         11^/utils/CI/sette@HEAD        sette 
         12 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/axis_def_nemo.xml

    r14072 r15574  
    2929  <axis id="section_ice" n_glo="4" /> 
    3030  <axis id="gau" /> 
     31  <!-- PISCES Sediment vertical axis definition --> 
     32  <axis id="profsed" long_name="Vertical S levels" unit="cm" positive="down"/> 
    3133</axis_definition> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/domain_def_nemo.xml

    r14072 r15574  
    44  <!--   T grid   --> 
    55  <domain id="grid_T" long_name="grid T"/> 
     6  <domain id="grid_T_inner" long_name="grid T inner"/> 
    67 
    78  <!--   My zoom: example of hand defined zoom   --> 
     
    167168  <!--   U grid   --> 
    168169  <domain id="grid_U" long_name="grid U"/> 
     170  <domain id="grid_U_inner" long_name="grid U inner"/> 
    169171  <!--   Eq section   --> 
    170172  <domain id="EqU" domain_ref="grid_U" > <zoom_domain id="EqU"/> </domain> 
     
    173175  <!--   V grid   --> 
    174176  <domain id="grid_V" long_name="grid V"/> 
     177  <domain id="grid_V_inner" long_name="grid V inner"/> 
    175178  <!--   Eq section : no V point on the Equator...   --> 
    176179 
     
    178181  <!--   W grid   --> 
    179182  <domain id="grid_W" long_name="grid W"/> 
     183  <domain id="grid_W_inner" long_name="grid W inner"/> 
    180184  <!--   Eq section   --> 
    181185  <domain id="EqW" domain_ref="grid_W" > <zoom_domain id="EqW"/> </domain> 
     
    183187  <!--   F grid   --> 
    184188  <domain id="grid_F" long_name="grid F"/> 
     189  <domain id="grid_F_inner" long_name="grid F inner"/> 
    185190 
    186191  <!--   zonal mean grid   --> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/field_def_nemo-ice.xml

    r14072 r15574  
    178178      <field id="hfxcndtop"    long_name="Net conductive heat flux at the ice surface (neg = ice cooling)"  standard_name="conductive_heat_flux_at_sea_ice_surface"    unit="W/m2" /> 
    179179      <field id="hfxcndbot"    long_name="Net conductive heat flux at the ice bottom (neg = ice cooling)"   standard_name="conductive_heat_flux_at_sea_ice_bottom"     unit="W/m2" /> 
     180      <field id="hfxmelt"      long_name="Melt heat flux at the ice surface"       unit="W/m2" /> 
     181      <field id="hfxldmelt"    long_name="Heat flux in the lead for ice melting"   unit="W/m2" /> 
     182      <field id="hfxldgrow"    long_name="Heat flux in the lead for ice growth"    unit="W/m2" /> 
    180183 
    181184      <!-- diags --> 
     
    345348      <field id="SH_icearea"      long_name="Sea ice area South"                     standard_name="sea_ice_area_s"                     unit="1e6_km2"  /> 
    346349 
    347       <!-- available with ln_icediaout --> 
     350      <!-- available with ln_icediahsb --> 
     351      <!-- global forcings  --> 
    348352      <field id="ibgfrcvoltop"    long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)"   unit="km3"      /> 
    349353      <field id="ibgfrcvolbot"    long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3"      /> 
     
    354358      <field id="ibgfrchfxbot"    long_name="global mean heat flux below ice (on top of ocean) "                                        unit="W/m2"     /> 
    355359 
     360      <!-- global drifts (conservation checks) --> 
    356361      <field id="ibgvolume"       long_name="drift in ice/snow volume (equivalent ocean volume)"                                        unit="km3"      /> 
    357362      <field id="ibgsaltco"       long_name="drift in ice salt content (equivalent ocean volume)"                                       unit="pss*km3"  /> 
     
    359364      <field id="ibgheatfx"       long_name="drift in ice/snow heat flux"                                                               unit="W/m2"     /> 
    360365 
     366      <!-- global contents --> 
    361367      <field id="ibgvol_tot"      long_name="global mean ice volume"                                                                    unit="km3"      /> 
    362368      <field id="sbgvol_tot"      long_name="global mean snow volume"                                                                   unit="km3"      /> 
    363369      <field id="ibgarea_tot"     long_name="global mean ice area"                                                                      unit="km2"      /> 
    364       <field id="ibgsalt_tot"     long_name="global mean ice salt content"                                                              unit="1e-3*km3" /> 
     370      <field id="ibgsalt_tot"     long_name="global mean ice salt content"                                                              unit="pss*km3" /> 
    365371      <field id="ibgheat_tot"     long_name="global mean ice heat content"                                                              unit="1e20J"    /> 
    366372      <field id="sbgheat_tot"     long_name="global mean snow heat content"                                                             unit="1e20J"    /> 
     373      <field id="ipbgvol_tot"     long_name="global mean ice pond volume"                                                               unit="km3"      /> 
     374      <field id="ilbgvol_tot"     long_name="global mean ice pond lid volume"                                                           unit="km3"      /> 
    367375 
    368376    </field_group> 
     
    496504    </field_group> 
    497505 
     506    <!--============================--> 
     507    <!--  CONSERVATION diagnostics  --> 
     508    <!--============================--> 
     509 
    498510    <field_group id="ICE_globalbudget"  grid_ref="grid_scalar" > 
    499       <!-- global contents --> 
    500511      <field field_ref="ibgvol_tot"       name="ibgvol_tot"   /> 
    501512      <field field_ref="sbgvol_tot"       name="sbgvol_tot"   /> 
     
    504515      <field field_ref="ibgheat_tot"      name="ibgheat_tot"  /> 
    505516      <field field_ref="sbgheat_tot"      name="sbgheat_tot"  /> 
    506  
    507       <!-- global drifts (conservation checks) --> 
    508       <field field_ref="ibgvolume"        name="ibgvolume"    /> 
    509       <field field_ref="ibgsaltco"        name="ibgsaltco"    /> 
    510       <field field_ref="ibgheatco"        name="ibgheatco"    /> 
    511       <field field_ref="ibgheatfx"        name="ibgheatfx"    /> 
    512  
    513       <!-- global forcings  --> 
    514       <field field_ref="ibgfrcvoltop"     name="ibgfrcvoltop" /> 
    515       <field field_ref="ibgfrcvolbot"     name="ibgfrcvolbot" /> 
    516       <field field_ref="ibgfrctemtop"     name="ibgfrctemtop" /> 
    517       <field field_ref="ibgfrctembot"     name="ibgfrctembot" /> 
    518       <field field_ref="ibgfrcsal"        name="ibgfrcsal"    /> 
    519       <field field_ref="ibgfrchfxtop"     name="ibgfrchfxtop" /> 
    520       <field field_ref="ibgfrchfxbot"     name="ibgfrchfxbot" /> 
    521     </field_group> 
    522  
    523  
     517      <field field_ref="ipbgvol_tot"      name="ipbgvol_tot"  /> 
     518      <field field_ref="ilbgvol_tot"      name="ilbgvol_tot"  /> 
     519    </field_group> 
     520     
     521    <field_group id="ICE_budget"        grid_ref="grid_T_2D" > 
     522      <!-- general --> 
     523      <field field_ref="icemask"          name="simsk"      /> 
     524      <field field_ref="iceconc"          name="siconc"     /> 
     525      <field field_ref="icetemp"          name="sitemp"     /> 
     526      <field field_ref="snwtemp"          name="sntemp"     /> 
     527      <field field_ref="icettop"          name="sittop"     /> 
     528      <field field_ref="icetbot"          name="sitbot"     /> 
     529      <!-- heat fluxes --> 
     530      <field field_ref="qt_oce_ai"        name="qt_oce_ai"  /> 
     531      <field field_ref="qt_atm_oi"        name="qt_atm_oi"  /> 
     532      <field field_ref="qtr_ice_top"      name="qtr_ice_top"/> 
     533      <field field_ref="qtr_ice_bot"      name="qtr_ice_bot"/> 
     534      <field field_ref="qt_ice"           name="qt_ice"     /> 
     535      <field field_ref="qsr_ice"          name="qsr_ice"    /> 
     536      <field field_ref="qns_ice"          name="qns_ice"    /> 
     537      <field field_ref="qemp_ice"         name="qemp_ice"   /> 
     538      <field field_ref="hfxsub"           name="hfxsub"     /> 
     539      <field field_ref="hfxspr"           name="hfxspr"     /> 
     540      <field field_ref="hfxcndtop"        name="hfxcndtop"  /> 
     541      <field field_ref="hfxcndbot"        name="hfxcndbot"  /> 
     542      <field field_ref="hfxsensib"        name="hfxsensib"  /> 
     543      <field field_ref="hfxmelt"          name="hfxmelt"    /> 
     544      <field field_ref="hfxldmelt"        name="hfxldmelt"  /> 
     545      <field field_ref="hfxldgrow"        name="hfxldgrow"  /> 
     546      <!-- salt fluxes --> 
     547      <field field_ref="sfxice"           name="sfxice"     /> 
     548      <!-- mass fluxes --> 
     549      <field field_ref="vfxice"           name="vfxice"     /> 
     550      <field field_ref="vfxsnw"           name="vfxsnw"     /> 
     551      <field field_ref="vfxpnd"           name="vfxpnd"     /> 
     552      <field field_ref="vfxsub"           name="vfxsub"     /> 
     553      <field field_ref="vfxsub_err"       name="vfxsub_err" /> 
     554      <field field_ref="vfxsnw_sub"       name="vfxsnw_sub" /> 
     555      <field field_ref="vfxsnw_pre"       name="vfxsnw_pre" /> 
     556    </field_group> 
     557 
     558     
    524559    <!--============================--> 
    525560    <!-- SIMIP sea ice field groups --> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/field_def_nemo-innerttrc.xml

    r14072 r15574  
    4343    <field id="C14Age"       long_name="Radiocarbon age"                          unit="yr"        grid_ref="grid_T_3D"  /> 
    4444    <field id="RAge"         long_name="Reservoir Age"                            unit="yr"       /> 
    45     <field id="qtr_C14"      long_name="Air-sea flux of C14"                      unit="1/m2/s"   /> 
    46     <field id="qint_C14"     long_name="Cumulative air-sea flux of C14"           unit="1/m2"     /> 
     45    <field id="qtr_c14"      long_name="Air-sea flux of C14"                      unit="1/m2/s"   /> 
     46    <field id="qint_c14"     long_name="Cumulative air-sea flux of C14"           unit="1/m2"     /> 
    4747    <field id="AtmCO2"       long_name="Global atmospheric CO2"                   unit="ppm"      /> 
    4848    <field id="AtmC14"       long_name="Global atmospheric DeltaC14"              unit="permil"   /> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/field_def_nemo-oce.xml

    r14224 r15574  
    131131 
    132132    <!-- additions to diawri.F90 --> 
    133     <field id="socegrad"    long_name="module of salinity gradient"              unit="psu/m"   grid_ref="grid_T_3D"/> 
    134     <field id="socegrad2"   long_name="square of module of salinity gradient"    unit="psu2/m2" grid_ref="grid_T_3D"/> 
     133    <field id="sssgrad"     long_name="module of surface salinity gradient"              unit="1e-3/m"   grid_ref="grid_T_2D_inner"/> 
     134    <field id="sssgrad2"    long_name="square of module of surface salinity gradient"    unit="1e-6/m2"  grid_ref="grid_T_2D_inner"/> 
    135135    <field id="ke"          long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" /> 
    136     <field id="ke_int"      long_name="vertical integration of kinetic energy"   unit="m3/s2"   /> 
    137     <field id="relvor"      long_name="relative vorticity"                       unit="s-1"    grid_ref="grid_T_3D"/> 
    138     <field id="absvor"      long_name="absolute vorticity"                       unit="s-1"    grid_ref="grid_T_3D"/> 
    139     <field id="potvor"      long_name="potential vorticity"                      unit="s-1"    grid_ref="grid_T_3D"/> 
     136    <field id="ke_int"      long_name="vertical integration of kinetic energy"   unit="m3/s2"  grid_ref="grid_T_2D_inner" /> 
    140137    <field id="salt2c"      long_name="Salt content vertically integrated"       unit="1e-3*kg/m2" /> 
    141138 
     
    149146    <field id="sstmax"       long_name="max of sea surface temperature"   field_ref="sst"   operation="maximum"                                                 /> 
    150147    <field id="sstmin"       long_name="min of sea surface temperature"   field_ref="sst"   operation="minimum"                                                 /> 
    151     <field id="sstgrad"      long_name="module of sst gradient"                                                                                 unit="degC/m"   /> 
    152     <field id="sstgrad2"     long_name="square of module of sst gradient"                                                                       unit="degC2/m2" /> 
    153     <field id="sbt"          long_name="sea bottom temperature"                                                                                 unit="degC"     /> 
    154     <field id="tosmint"      long_name="vertical integral of temperature times density"     standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature"  unit="(kg m2) degree_C" /> 
     148    <field id="sstgrad"      long_name="module of sst gradient"                                                                                 unit="degC/m"   grid_ref="grid_T_2D_inner" /> 
     149    <field id="sstgrad2"     long_name="square of module of sst gradient"                                                                       unit="degC2/m2" grid_ref="grid_T_2D_inner" /> 
     150    <field id="sbt"          long_name="sea bottom temperature"                                                                                 unit="degC"     grid_ref="grid_T_2D_inner" /> 
     151    <field id="tosmint"      long_name="vertical integral of temperature times density"     standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature"  unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner" /> 
    155152    <field id="sst_wl"       long_name="Delta SST of warm layer"                                                                                unit="degC"     /> 
    156153    <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                                 unit="degC"     /> 
     
    161158    <field id="sssmax"       long_name="max of sea surface salinity"      field_ref="sss"   operation="maximum"                                                 /> 
    162159    <field id="sssmin"       long_name="min of sea surface salinity"      field_ref="sss"   operation="minimum"                                                 /> 
    163     <field id="sbs"          long_name="sea bottom salinity"                                                                                    unit="0.001"    /> 
    164     <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" /> 
    165  
    166     <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"    /> 
     160    <field id="sbs"          long_name="sea bottom salinity"                                                                                    unit="0.001"   grid_ref="grid_T_2D_inner" /> 
     161    <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" /> 
     162 
     163    <field id="taubot"       long_name="bottom stress module"   unit="N/m2"   grid_ref="grid_T_2D_inner" /> 
    167164 
    168165    <!-- Case EOS = TEOS-10 : output potential temperature --> 
     
    181178    <field id="mldr10_1max"  long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)"   field_ref="mldr10_1"   operation="maximum"                                                                          /> 
    182179    <field id="mldr10_1min"  long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)"   field_ref="mldr10_1"   operation="minimum"                                                                          /> 
    183     <field id="heatc"        long_name="Heat content vertically integrated"                 standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"      /> 
    184     <field id="saltc"        long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3*kg/m2" /> 
    185     <field id="saltc2"       long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3*kg/m2" > saltc * saltc </field > 
     180    <field id="heatc"        long_name="Heat content vertically integrated"                 standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"      grid_ref="grid_T_2D_inner" /> 
     181    <field id="saltc"        long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3.m"    grid_ref="grid_T_2D_inner" /> 
     182    <field id="saltc2"       long_name="square of Salt content vertically integrated"                                                                                                         unit="1e-6.m2"   grid_ref="grid_T_2D_inner" /> 
    186183 
    187184    <!-- EOS --> 
     
    191188 
    192189    <!-- Energy - horizontal divergence --> 
    193     <field id="sKE"          long_name="surface kinetic energy"  standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_2D" /> 
     190    <field id="sKE"          long_name="surface kinetic energy"  standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_2D_inner" /> 
    194191    <field id="hdiv"         long_name="horizontal divergence"                                                          unit="s-1"    grid_ref="grid_T_3D" /> 
    195192 
     
    240237    <field id="mf_Sp"       long_name="plume_salinity"         standard_name="plume_salinity"        unit="1e-3"   grid_ref="grid_T_3D" /> 
    241238    <field id="mf_mf"       long_name="mass flux"              standard_name="mf_mass_flux"          unit="m"      grid_ref="grid_T_3D" /> 
     239 
     240    <!-- fluxes from damping --> 
     241    <field id="sflx_dmp_cea"  long_name="salt flux due to damping"  standard_name="salt_flux_due_to_damping"          unit="g/m2/s"   /> 
     242    <field id="hflx_dmp_cea"  long_name="heat flux due to damping"  standard_name="heat_flux_due_to_damping"          unit="W/m2"     /> 
     243 
     244    <!-- * variable related to ice shelf forcing * --> 
     245 
     246    <!-- * fwf * --> 
     247    <field id="fwfisf_cav"      long_name="Ice shelf fresh water flux ( from isf to oce )"                    unit="kg/m2/s"  /> 
     248    <field id="fwfisf_par"      long_name="Ice shelf fresh water flux ( from isf to oce )"                    unit="kg/m2/s"  /> 
     249    <field id="fwfisf3d_cav"    long_name="3d Ice shelf fresh water flux ( from isf to oce )"                 unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     250    <field id="fwfisf3d_par"    long_name="3d Ice shelf fresh water flux ( from isf to oce )"                 unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     251 
     252    <!-- * heat fluxes * --> 
     253    <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     /> 
     254    <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     /> 
     255    <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     /> 
     256    <field id="qlatisf_par"     long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     /> 
     257    <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     /> 
     258    <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     /> 
     259    <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     260    <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     261    <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     262    <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     263    <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     grid_ref="grid_T_3D" /> 
     264    <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     grid_ref="grid_T_3D" /> 
     265    <field id="qconisf"         long_name="Conductive heat flux through the ice shelf ( from isf to oce )"    unit="W/m2"     /> 
     266 
     267    <!-- top boundary layer properties --> 
     268    <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     /> 
     269    <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     /> 
     270    <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     271    <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     272    <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      /> 
     273    <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      /> 
     274    <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     /> 
     275    <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" /> 
     276    <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
     277    <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      /> 
     278    <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      /> 
     279    <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      /> 
    242280 
    243281  </field_group> <!-- grid_T --> 
     
    276314    <field id="dh"                  long_name="Pycnocline thickness"                     unit=" m"      /> 
    277315    <field id="ibld"                long_name="index of boundary layer depth"            unit="#"       /> 
    278     <field id="imld"                long_name="index of mixed layer depth"            unit="#"       /> 
    279     <field id="zhbl"                long_name="boundary layer depth -grid"                     unit="m"       /> 
    280     <field id="zhml"                long_name="mixed layer depth - grid"                        unit="m"       /> 
     316    <field id="imld"                long_name="index of mixed layer depth"               unit="#"       /> 
     317    <field id="jp_ext"              long_name="flag =1 if pycnocline well resolved"      unit="#"       /> 
     318    <field id="j_ddh"               long_name="index of mixed layer depth"               unit="#"       /> 
     319    <field id="zshear"              long_name="shear production of TKE "                 unit="m^3/s^3" /> 
     320    <field id="zhbl"                long_name="boundary layer depth -grid"               unit="m"       /> 
     321    <field id="zhml"                long_name="mixed layer depth - grid"                 unit="m"       /> 
    281322    <field id="zdh"                 long_name="Pycnocline  depth - grid"                 unit=" m"      /> 
    282323    <field id="zustke"              long_name="magnitude of stokes drift  at T-points"   unit="m/s"     /> 
    283     <field id="us_x"        long_name="i component of active Stokes drift"                      unit="m/s"     /> 
    284     <field id="us_y"        long_name="j component of active Stokes drift"                      unit="m/s"     /> 
     324    <field id="us_x"                long_name="i component of active Stokes drift"       unit="m/s"     /> 
     325    <field id="us_y"                long_name="j component of active Stokes drift"       unit="m/s"     /> 
    285326    <field id="dstokes"             long_name="stokes drift  depth scale"                unit="m"       /> 
    286327    <field id="zwth0"               long_name="surface non-local temperature flux"       unit="deg m/s" /> 
    287328    <field id="zws0"                long_name="surface non-local salinity flux"          unit="psu m/s" /> 
     329    <field id="zwb0"                long_name="surface non-local buoyancy flux"          unit="m^2/s^3" /> 
    288330    <field id="zwstrc"              long_name="convective velocity scale"                unit="m/s"     /> 
    289331    <field id="zustar"              long_name="friction velocity"                        unit="m/s"     /> 
     
    296338 
    297339    <!-- interior BL OSMOSIS diagnostics --> 
    298     <field id="zwthav"              long_name="av turb flux of T in ml"                  unit="deg m/s" /> 
     340    <field id="zwbav"               long_name="av turb flux of buoyancy in ml"           unit="m^2/s^3" /> 
    299341    <field id="zt_ml"               long_name="av T in ml"                               unit="deg"     /> 
    300342    <field id="zhol"                long_name="Hoenekker number"                         unit="#"       /> 
     
    303345    <field id="zwb_ent"            long_name="entrainment turb flux of buoyancy"         unit="m^2/s^-3" /> 
    304346 
    305     <field id="zdt_bl"             long_name="temperature jump at base of BL"                 unit="deg"      /> 
    306     <field id="zds_bl"             long_name="salinity jump at base of BL"                 unit="10^-3"      /> 
    307     <field id="zdb_bl"             long_name="buoyancy jump at base of BL"                 unit="m/s^2"      /> 
    308     <field id="zdu_bl"             long_name="u jump at base of BL"                       unit="m/s"      /> 
    309     <field id="zdv_bl"             long_name="v jump at base of BL"                       unit="m/s"      /> 
    310  
     347    <field id="zdt_bl"             long_name="temperature jump at base of BL"            unit="deg"      /> 
     348    <field id="zds_bl"             long_name="salinity jump at base of BL"               unit="10^-3"    /> 
     349    <field id="zdb_bl"             long_name="buoyancy jump at base of BL"               unit="m/s^2"    /> 
     350    <field id="zdu_bl"             long_name="u jump at base of BL"                      unit="m/s"      /> 
     351    <field id="zdv_bl"             long_name="v jump at base of BL"                      unit="m/s"      /> 
     352    <field id="zdt_ml"             long_name="temperature jump at base of ML"            unit="deg"      /> 
     353    <field id="zds_ml"             long_name="salinity jump at base of ML"               unit="10^-3"    /> 
     354    <field id="zdb_ml"             long_name="buoyancy jump at base of ML"               unit="m/s^2"    /> 
     355    <field id="pb_coup"            long_name="bottom coupling velocity"                  unit="m/s"      /> 
    311356    <!-- extra OSMOSIS diagnostics for debugging --> 
    312357    <field id="zsc_uw_1_0"       long_name="zsc u-momentum flux on T after Stokes"                       unit="m^2/s^2" /> 
     
    315360    <field id="zsc_uw_2_f"       long_name="2nd zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
    316361    <field id="zsc_vw_2_f"       long_name="2nd zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" /> 
    317     <field id="zuw_bse"       long_name="base u-flux T-points"                          unit="m^2/s^2" /> 
    318     <field id="zvw_bse"       long_name="base v-flux T-points"                          unit="m^2/s^2" /> 
    319362 
    320363    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.--> 
     
    375418      <field id="emp_oce"      long_name="Evap minus Precip over ocean"         standard_name="evap_minus_precip_over_sea_water"                                     unit="kg/m2/s"   /> 
    376419      <field id="emp_ice"      long_name="Evap minus Precip over ice"           standard_name="evap_minus_precip_over_sea_ice"                                       unit="kg/m2/s"   /> 
    377       <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="1e-3/m2/s" /> 
     420      <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="g/m2/s"    /> 
    378421      <field id="fmmflx"       long_name="Water flux due to freezing/melting"                                                                                        unit="kg/m2/s"   /> 
    379422      <field id="snowpre"      long_name="Snow precipitation"                   standard_name="snowfall_flux"                                                        unit="kg/m2/s"   /> 
     
    394437      <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 
    395438      <field id="ssh_ib"       long_name="Inverse barometer sea surface height"  standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency"   unit="m"        /> 
    396  
    397       <!-- * variable related to ice shelf forcing * --> 
    398       <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     /> 
    399       <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     /> 
    400       <field id="fwfisf_cav"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
    401       <field id="fwfisf_par"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
    402       <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
    403       <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
    404       <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    405       <field id="qlatisf_par"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    406       <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
    407       <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
    408       <field id="fwfisf3d_cav"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    409       <field id="fwfisf3d_par"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    410       <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    411       <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    412       <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    413       <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    414       <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
    415       <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
    416       <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     /> 
    417       <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" /> 
    418       <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
    419       <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
    420       <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      /> 
    421       <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      /> 
    422       <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
    423       <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      /> 
    424       <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      /> 
    425       <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      /> 
    426       <field id="qconisf"         long_name="Conductive heat flux through the ice shelf"    unit="W/m2"     /> 
    427439 
    428440      <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
     
    475487      <field id="hflx_rain_cea" long_name="heat flux due to rainfall"                                standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water"        unit="W/m2"     /> 
    476488      <field id="hflx_evap_cea" long_name="heat flux due to evaporation"                             standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water"   unit="W/m2"     /> 
     489      <field id="hflx_subl_cea" long_name="heat flux due to sublimation (from atm. forcings)"        standard_name="temperature_flux_due_to_sublimation_expressed_as_heat_flux_out_of_sea_ice"     unit="W/m2"     /> 
    477490      <field id="hflx_prec_cea" long_name="heat flux due to all precip"                              standard_name="temperature_flux_due_to_all_precip_expressed_as_heat_flux_into_sea_water"      unit="W/m2"     /> 
    478491      <field id="hflx_snow_cea" long_name="heat flux due to snow falling"                            standard_name="heat_flux_onto_ocean_and_ice_due_to_snow_thermodynamics"                       unit="W/m2"     /> 
     
    481494      <field id="hflx_ice_cea"  long_name="heat flux due to ice thermodynamics"                      standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics"                        unit="W/m2"     /> 
    482495      <field id="hflx_rnf_cea"  long_name="heat flux due to runoffs"                                 standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water"          unit="W/m2"     /> 
     496      <field id="sflx_rnf_cea"  long_name="salt flux due to runoffs"                                 standard_name="salt_flux_due_to_runoffs"                                                      unit="g/m2/s"   /> 
    483497      <field id="hflx_cal_cea"  long_name="heat flux due to calving"                                 standard_name="heat_flux_into_sea_water_due_to_calving"                                       unit="W/m2"     /> 
    484498      <field id="hflx_icb_cea"  long_name="heat flux due to iceberg"                                 standard_name="heat_flux_into_sea_water_due_to_icebergs"                                      unit="W/m2"     /> 
     
    490504      <field id="ticemel_cea"   long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"  standard_name="tendency_of_sea_ice_amount_due_to_surface_melting"                             unit="kg/m2/s"  /> 
    491505 
     506      <!-- fluxes from relaxation and freshwater budget --> 
     507      <field id="sflx_ssr_cea"  long_name="salt flux due to restoring"    standard_name="salt_flux_due_to_restoring"    unit="g/m2/s"   /> 
     508      <field id="hflx_ssr_cea"  long_name="heat flux due to restoring"    standard_name="heat_flux_due_to_restoring"    unit="W/m2"     /> 
     509      <field id="vflx_ssr_cea"  long_name="volume flux due to restoring"  standard_name="volume_flux_due_to_restoring"  unit="kg/m2/s"  /> 
     510      <field id="hflx_fwb_cea"  long_name="heat flux due to fwb"          standard_name="heat_flux_due_to_fwb"          unit="W/m2"     /> 
     511      <field id="vflx_fwb_cea"  long_name="volume flux due to fwb"        standard_name="volume_flux_due_to_fwb"        unit="kg/m2/s"  /> 
     512       
    492513      <!-- ice field (nn_ice=1)  --> 
    493514      <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            /> 
     
    576597 
    577598    <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             /> 
    578     <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                            /> 
     599    <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"  grid_ref="grid_U_2D_inner" /> 
    579600    <field id="ubar"         long_name="ocean barotropic current along i-axis"                                                              unit="m/s"                             /> 
    580601    <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"                 standard_name="ocean_volume_x_transport"    unit="m3/s"       grid_ref="grid_U_3D" /> 
     
    615636 
    616637    <!-- variables available with diaar5 --> 
    617     <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s"  grid_ref="grid_U_3D" /> 
    618     <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport"    unit="kg/s"                     /> 
    619     <field id="u_heattr"      long_name="ocean eulerian heat transport along i-axis"                      standard_name="ocean_heat_x_transport"                         unit="W"                        /> 
    620     <field id="u_salttr"      long_name="ocean eulerian salt transport along i-axis"                      standard_name="ocean_salt_x_transport"                         unit="1e-3*kg/s"                 /> 
     638    <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s"   grid_ref="grid_U_3D" /> 
     639    <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport"    unit="kg/s"   grid_ref="grid_U_2D_inner" /> 
     640    <field id="u_heattr"      long_name="ocean eulerian heat transport along i-axis"                      standard_name="ocean_heat_x_transport"                         unit="W"      grid_ref="grid_U_2D_inner" /> 
     641    <field id="u_salttr"      long_name="ocean eulerian salt transport along i-axis"                      standard_name="ocean_salt_x_transport"                         unit="1e-3*kg/s" grid_ref="grid_U_2D_inner" /> 
    621642    <field id="uadv_heattr"   long_name="ocean advective heat transport along i-axis"                     standard_name="advectice_ocean_heat_x_transport"               unit="W"                         /> 
    622643    <field id="uadv_salttr"   long_name="ocean advective salt transport along i-axis"                     standard_name="advectice_ocean_salt_x_transport"               unit="1e-3*kg/s"                 /> 
     
    636657    <field id="voce_e3v"     long_name="ocean current along j-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_V_3D"  > voce * e3v </field> 
    637658    <field id="ssv"          long_name="ocean surface current along j-axis"                                                                 unit="m/s"                             /> 
    638     <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"                            /> 
     659    <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"  grid_ref="grid_V_2D_inner" /> 
    639660    <field id="vbar"         long_name="ocean barotropic current along j-axis"                                                              unit="m/s"                             /> 
    640661    <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"                 standard_name="ocean_volume_y_transport"    unit="m3/s"       grid_ref="grid_V_3D" /> 
     
    677698    <!-- variables available with diaar5 --> 
    678699    <field id="v_masstr"      long_name="ocean eulerian mass transport along j-axis"    standard_name="ocean_mass_y_transport"                          unit="kg/s" grid_ref="grid_V_3D" /> 
    679     <field id="v_heattr"      long_name="ocean eulerian heat transport along j-axis"    standard_name="ocean_heat_y_transport"                          unit="W"                        /> 
    680     <field id="v_salttr"      long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_y_transport"                          unit="1e-3*kg/s"                /> 
     700    <field id="v_heattr"      long_name="ocean eulerian heat transport along j-axis"    standard_name="ocean_heat_y_transport"                          unit="W"          grid_ref="grid_V_2D_inner" /> 
     701    <field id="v_salttr"      long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_y_transport"                          unit="1e-3*kg/s"  grid_ref="grid_V_2D_inner" /> 
    681702    <field id="vadv_heattr"   long_name="ocean advective heat transport along j-axis"   standard_name="advectice_ocean_heat_y_transport"                unit="W"                         /> 
    682703    <field id="vadv_salttr"   long_name="ocean advective salt transport along j-axis"   standard_name="advectice_ocean_salt_y_transport"                unit="1e-3*kg/s"                 /> 
     
    737758    <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                unit="s-2" /> 
    738759 
     760    <!-- dissipation diagnostics (note: ediss_k is only available with tke scheme) -->    
     761    <field id="avt_k"        long_name="vertical eddy diffusivity from closure schemes" standard_name="ocean_vertical_eddy_diffusivity"       unit="m2/s" /> 
     762    <field id="avm_k"        long_name="vertical eddy viscosity from closure schemes"   standard_name="ocean_vertical_eddy_viscosity"         unit="m2/s" /> 
     763    <field id="ediss_k"      long_name="Kolmogorov energy dissipation (tke scheme)"     standard_name="Kolmogorov_energy_dissipation"         unit="W/kg" /> 
     764    <field id="eshear_k"     long_name="energy source from vertical shear"              standard_name="energy_source_from_shear"              unit="W/kg" /> 
     765    <field id="estrat_k"     long_name="energy sink from stratification"                standard_name="energy_sink_from_stratification"       unit="W/kg" /> 
     766     
    739767  </field_group> 
    740768 
     
    744772    <field id="e3f_0"        long_name="F-cell thickness"                      standard_name="cell_thickness"         unit="m"   grid_ref="grid_F_3D" /> 
    745773    <field id="hf"           long_name="water column height at F point"        standard_name="water_column_height_F"  unit="m"                     /> 
    746     <field id="ssKEf"        long_name="surface kinetic energy at F point"     standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2" /> 
    747     <field id="ssrelvor"     long_name="surface relative vorticity"            standard_name="relative_vorticity"     unit="1/s"                   /> 
    748     <field id="ssplavor"     long_name="surface planetary vorticity"           standard_name="planetary_vorticity"    unit="1/s"                   /> 
    749     <field id="ssrelpotvor"  long_name="surface relative potential vorticity"  standard_name="relpot_vorticity"       unit="1/m.s"                 /> 
    750     <field id="ssabspotvor"  long_name="surface absolute potential vorticity"  standard_name="abspot_vorticity"       unit="1/m.s"                 /> 
    751     <field id="ssEns"        long_name="surface enstrophy"                     standard_name="enstrophy"              unit="1/m2.s2"               /> 
     774    <field id="ssKEf"        long_name="surface kinetic energy at F point"     standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2" grid_ref="grid_F_2D_inner" /> 
     775    <field id="ssrelvor"     long_name="surface relative vorticity"            standard_name="relative_vorticity"     unit="1/s"       grid_ref="grid_F_2D_inner" /> 
     776    <field id="ssplavor"     long_name="surface planetary vorticity"           standard_name="planetary_vorticity"    unit="1/s"       /> 
     777    <field id="ssrelpotvor"  long_name="surface relative potential vorticity"  standard_name="relpot_vorticity"       unit="1/m.s"     grid_ref="grid_F_2D_inner" /> 
     778    <field id="ssabspotvor"  long_name="surface absolute potential vorticity"  standard_name="abspot_vorticity"       unit="1/m.s"     grid_ref="grid_F_2D_inner" /> 
     779    <field id="ssEns"        long_name="surface enstrophy"                     standard_name="enstrophy"              unit="1/m2.s2"   grid_ref="grid_F_2D_inner" /> 
    752780  </field_group> 
    753781 
     
    10921120  </field_group> 
    10931121 
     1122  <!-- shared variables available with TOP interface --> 
     1123  <field_group id="top_shared" grid_ref="grid_T_3D"> 
     1124    <field id="xeps"           long_name="Broadband light attenuation"                     unit="-"                            /> 
     1125    <field id="Heup"           long_name="Euphotic layer depth"                            unit="m"     grid_ref="grid_T_2D"   /> 
     1126  </field_group> 
     1127 
    10941128 
    10951129  <!-- 
     
    12091243 
    12101244  <!-- 25h diagnostic output --> 
    1211   <field_group id="25h_grid_T" grid_ref="grid_T_3D" operation="instant"> 
     1245  <field_group id="25h_grid_T" grid_ref="grid_T_3D_inner" operation="instant"> 
    12121246    <field id="temper25h"         name="potential temperature 25h mean"    unit="degC" /> 
    12131247    <field id="tempis25h"         name="insitu temperature 25h mean"    unit="degC" /> 
    12141248    <field id="salin25h"          name="salinity 25h mean"                 unit="psu"  /> 
    1215     <field id="ssh25h"            name="sea surface height 25h mean"  grid_ref="grid_T_2D"      unit="m"    /> 
    1216   </field_group> 
    1217  
    1218   <field_group id="25h_grid_U" grid_ref="grid_U_3D" operation="instant" > 
     1249    <field id="ssh25h"            name="sea surface height 25h mean"  grid_ref="grid_T_2D_inner"      unit="m"    /> 
     1250  </field_group> 
     1251 
     1252  <field_group id="25h_grid_U" grid_ref="grid_U_3D_inner" operation="instant" > 
    12191253    <field id="vozocrtx25h"         name="i current 25h mean"    unit="m/s"   /> 
    12201254  </field_group> 
    12211255 
    1222   <field_group id="25h_grid_V" grid_ref="grid_V_3D" operation="instant"> 
     1256  <field_group id="25h_grid_V" grid_ref="grid_V_3D_inner" operation="instant"> 
    12231257    <field id="vomecrty25h"         name="j current 25h mean"    unit="m/s"    /> 
    12241258  </field_group> 
    12251259 
    1226   <field_group id="25h_grid_W" grid_ref="grid_W_3D" operation="instant"> 
     1260  <field_group id="25h_grid_W" grid_ref="grid_W_3D_inner" operation="instant"> 
    12271261    <field id="vovecrtz25h"         name="k current 25h mean"                 unit="m/s"      /> 
    12281262    <field id="avt25h"              name="vertical diffusivity25h mean"       unit="m2/s" /> 
     
    12531287  </field_group> 
    12541288 
     1289  <!--============================--> 
     1290  <!--  CONSERVATION diagnostics  --> 
     1291  <!--============================--> 
     1292  <!-- BE CAREFUL: this group (OCE_budget) cannot be called in file_def.xml as such (unless nn_fsbc=1) 
     1293                   If doing so, the last output (in time) of the netcdf file  
     1294         would be corrupted (NaN values). However calling each of these 
     1295         variables directly in the file_def.xml works. It is probably  
     1296         because there is a mix up of sbc variables with other variables 
     1297    --> 
     1298  <field_group id="OCE_budget"        grid_ref="grid_T_2D" > 
     1299    <field field_ref="sst"                 name="tos"          /> 
     1300    <field field_ref="sss"                 name="sos"          /> 
     1301    <field field_ref="ssh"                 name="zos"          /> 
     1302    <!-- mass flux --> 
     1303    <field field_ref="empmr"               name="empmr"        /> 
     1304    <field field_ref="runoffs"             name="runoffs"      /> 
     1305    <field field_ref="emp_ice"             name="emp_ice"      /> 
     1306    <field field_ref="emp_oce"             name="emp_oce"      /> 
     1307    <field field_ref="iceshelf_cea"        name="iceshelf"     /> 
     1308    <field field_ref="iceberg_cea"         name="iceberg"      /> 
     1309    <field field_ref="calving_cea"         name="calving"      /> 
     1310    <!-- <field field_ref="berg_floating_melt"  name="calving" /> --> 
     1311    <field field_ref="precip"              name="precip"       /> 
     1312    <field field_ref="snowpre"             name="snowpre"      /> 
     1313    <field field_ref="rain"                name="rain"         /> 
     1314    <field field_ref="evap_ao_cea"         name="evap_ao"      /> 
     1315    <field field_ref="subl_ai_cea"         name="subl_ai"      /> 
     1316    <field field_ref="snow_ai_cea"         name="snow_ai"      /> 
     1317    <field field_ref="snow_ao_cea"         name="snow_ao"      /> 
     1318    <!-- heat flux --> 
     1319    <field field_ref="qsr"                 name="qsr"          /> 
     1320    <field field_ref="qns"                 name="qns"          /> 
     1321    <field field_ref="qt_oce"              name="qt_oce"       /> 
     1322    <field field_ref="qemp_oce"            name="qemp_oce"     /> 
     1323    <field field_ref="hflx_rain_cea"       name="hflx_rain"    /> 
     1324    <field field_ref="hflx_evap_cea"       name="hflx_evap"    /> 
     1325    <field field_ref="hflx_snow_cea"       name="hflx_snow"    /> 
     1326    <field field_ref="hflx_snow_ao_cea"    name="hflx_snow_ao" /> 
     1327    <field field_ref="hflx_snow_ai_cea"    name="hflx_snow_ai" /> 
     1328    <field field_ref="hflx_rnf_cea"        name="hflx_rnf"     /> 
     1329    <field field_ref="hflx_icb_cea"        name="hflx_icb"     /> 
     1330    <field field_ref="hflx_isf_cea"        name="hflx_isf"     /> 
     1331    <!-- salt flux (includes ssr) --> 
     1332    <field field_ref="saltflx"             name="saltflx"      /> 
     1333    <field field_ref="sflx_rnf_cea"        name="sflx_rnf"     /> 
     1334    <!-- relaxation and damping --> 
     1335    <field field_ref="hflx_ssr_cea"        name="hflx_ssr"     /> 
     1336    <field field_ref="vflx_ssr_cea"        name="vflx_ssr"     /> 
     1337    <field field_ref="sflx_ssr_cea"        name="sflx_ssr"     /> 
     1338    <field field_ref="hflx_dmp_cea"        name="hflx_dmp"     /> 
     1339    <field field_ref="sflx_dmp_cea"        name="sflx_dmp"     /> 
     1340    <field field_ref="hflx_fwb_cea"        name="hflx_fwb"     /> 
     1341    <field field_ref="vflx_fwb_cea"        name="vflx_fwb"     /> 
     1342  </field_group> 
     1343 
     1344  <field_group id="OCE_globalbudget"  grid_ref="grid_scalar" > 
     1345    <field field_ref="voltot"              name="scvoltot"     /> 
     1346    <field field_ref="saltot"              name="scsaltot"     /> 
     1347    <field field_ref="temptot"             name="sctemtot"     /> 
     1348  </field_group> 
     1349 
     1350   
    12551351</field_definition> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/field_def_nemo-pisces.xml

    r14209 r15574  
    136136    <field id="SedpH"           long_name="PH"                                       unit="1"          /> 
    137137    <field id="SedCO3por"       long_name="Bicarbonates"                             unit="mol/m3" /> 
     138    <field id="Sedligand"       long_name="Ligands"                                  unit="mol/m3" /> 
     139    <field id="SaturCO3"        long_name="CO3 Saturation"                           unit="-" />        
    138140  </field_group> 
    139141 
    140142  <!-- SEDIMENT additional variables on T sediment grid  --> 
    141143  <field_group id="Diag_S" grid_ref="grid_T_2D"> 
    142     <field id="FlxSi"       long_name="Si sediment flux"                        unit="mol/cm2/s"     /> 
    143     <field id="FlxO2"       long_name="O2 sediment flux"                        unit="mol/cm2/s"     /> 
    144     <field id="FlxDIC"      long_name="DIC sediment flux"                       unit="mol/cm2/s"     /> 
    145     <field id="FlxNO3"      long_name="NO3 sediment flux"                       unit="mol/cm2/s"     /> 
    146     <field id="FlxPO4"      long_name="PO4 sediment flux"                       unit="mol/cm2/s"     /> 
    147     <field id="FlxAlkalini" long_name="Alkalinity sediment flux"                unit="mol/cm2/s"     /> 
    148     <field id="FlxNH4"      long_name="Ammonium sediment flux"                  unit="mol/cm2/s"     /> 
    149     <field id="FlxH2S"      long_name="H2S sediment flux"                       unit="mol/cm2/s"     /> 
    150     <field id="FlxSO4"      long_name="SO4 sediment flux"                       unit="mol/cm2/s"     /> 
    151     <field id="FlxFe2"      long_name="Fe2+ sediment flux"                      unit="mol/cm2/s"     /> 
    152     <field id="Flxtot"      long_name="Sediment net burial rate"                unit="cm/s"     /> 
    153     <field id="dzdep"       long_name="Sedimentation rate"                      unit="cm/s"     /> 
    154     <field id="sflxclay"    long_name="Clay sedimentation rate"                 unit="g/m2/s"     /> 
    155     <field id="sflxcal"     long_name="Calcite sedimentation rate"              unit="mol/m2/s"     /> 
    156     <field id="sflxbsi"     long_name="BSi Sedimentation rate"                  unit="mol/m2/s"     /> 
    157     <field id="sflxpoc"     long_name="POC Sedimentation rate"                  unit="mol/m2/s"     /> 
    158     <field id="sflxfeo"     long_name="Fe(OH)3 Sedimentation rate"              unit="mol/m2/s"     /> 
     144      <field id="FlxSi"       long_name="Si sediment flux"                        unit="mol/cm2/s"     /> 
     145       <field id="FlxO2"       long_name="O2 sediment flux"                        unit="mol/cm2/s"     /> 
     146       <field id="FlxDIC"      long_name="DIC sediment flux"                       unit="mol/cm2/s"     /> 
     147       <field id="FlxNO3"      long_name="NO3 sediment flux"                       unit="mol/cm2/s"     /> 
     148       <field id="FlxPO4"      long_name="PO4 sediment flux"                       unit="mol/cm2/s"     /> 
     149       <field id="FlxAlkalini" long_name="Alkalinity sediment flux"                unit="mol/cm2/s"     /> 
     150       <field id="FlxNH4"      long_name="Ammonium sediment flux"                  unit="mol/cm2/s"     /> 
     151       <field id="FlxH2S"      long_name="H2S sediment flux"                       unit="mol/cm2/s"     /> 
     152       <field id="FlxSO4"      long_name="SO4 sediment flux"                       unit="mol/cm2/s"     /> 
     153       <field id="FlxFe2"      long_name="Fe2+ sediment flux"                      unit="mol/cm2/s"     /> 
     154       <field id="FlxLig"      long_name="Ligand sediment flux"                    unit="mol/cm2/s"     /> 
     155       <field id="dzdep"       long_name="Sedimentation rate"                      unit="cm/s"     /> 
     156       <field id="sflxclay"    long_name="Clay sedimentation rate"                 unit="g/cm2/s"     /> 
     157       <field id="sflxbsi"     long_name="BSi sedimentation rate"                  unit="g/cm2/s"     /> 
     158       <field id="sflxpoc"     long_name="POC sedimentation rate"                  unit="g/cm2/s"     /> 
     159       <field id="sflxcal"     long_name="Calcite sedimentation rate"              unit="g/cm2/s"     /> 
     160       <field id="sflxfeo"     long_name="Fe(OH)3 Sedimentation rate"              unit="mol/cm2/s"     /> 
     161       <field id="FlxClay"     long_name="Clay burial rate"                        unit="g/cm2/s"     /> 
     162       <field id="FlxCaCO3"    long_name="Calcite burial rate"                     unit="g/cm2/s"     /> 
     163       <field id="FlxBSi"      long_name="BSi burial rate"                         unit="g/cm2/s"     /> 
     164       <field id="FlxPOC"      long_name="POC burial rate"                         unit="g/cm2/s"     /> 
     165       <field id="FlxFeO"      long_name="Fe(OH)3 burial rate"                     unit="g/cm2/s"     /> 
     166       <field id="FlxFeS"      long_name="FeS burial rate"                         unit="g/cm2/s"     /> 
     167       <field id="FlxPOR"      long_name="POR burial rate"                         unit="g/cm2/s"     /> 
     168       <field id="FlxPOS"      long_name="POS burial rate"                         unit="g/cm2/s"     /> 
     169       <field id="Flxtot"      long_name="total burial flux"                       unit="g/cm2/s"     /> 
     170       <field id="Rstepros"    long_name="Number of iterations"                    unit="-"     /> 
    159171  </field_group> 
    160172 
     
    205217    <field id="SIZEP"       long_name="Mean relative size of picophyto."        unit="-"          grid_ref="grid_T_3D" /> 
    206218    <field id="SIZED"       long_name="Mean relative size of diatoms"           unit="-"          grid_ref="grid_T_3D" /> 
     219    <field id="RASSD"       long_name="Size of the protein machinery (Diat.)"   unit="-"          grid_ref="grid_T_3D" /> 
     220    <field id="RASSN"       long_name="Size of the protein machinery (Nano.)"   unit="-"          grid_ref="grid_T_3D" /> 
     221    <field id="RASSP"       long_name="Size of the protein machinery (Pico.)"   unit="-"          grid_ref="grid_T_3D" /> 
    207222    <field id="Fe3"         long_name="Iron III concentration"                  unit="nmol/m3"    grid_ref="grid_T_3D" /> 
    208223    <field id="FeL1"        long_name="Complexed Iron concentration with L1"    unit="nmol/m3"    grid_ref="grid_T_3D" /> 
     
    241256    <field id="BACT"        long_name="Bacterial Biomass"                       unit="mmol/m3"  grid_ref="grid_T_3D" /> 
    242257    <field id="FEBACT"      long_name="Bacterial uptake of Fe"                  unit="molFe/m3/s"  grid_ref="grid_T_3D" /> 
     258    <field id="FEPREC"      long_name="Precipitation of Fe"                     unit="molFe/m3/s"  grid_ref="grid_T_3D" /> 
    243259    <field id="LPRODR"      long_name="OM remineralisation ligand production rate" unit="nmol-L/m3/s"  grid_ref="grid_T_3D" /> 
    244260    <field id="LPRODP"      long_name="phytoplankton ligand production rate"    unit="nmol-L/m3/s"  grid_ref="grid_T_3D" /> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/grid_def_nemo.xml

    r14072 r15574  
    1212    <domain domain_ref="grid_T" /> 
    1313  </grid> 
     14  <grid id="grid_T_2D_inner" > 
     15    <domain domain_ref="grid_T_inner" /> 
     16  </grid> 
    1417  <!--  --> 
    1518  <grid id="grid_T_ncatice" > 
     
    1720    <axis axis_ref="ncatice" /> 
    1821  </grid> 
     22  <grid id="grid_T_ncatice_inner" > 
     23    <domain domain_ref="grid_T_inner" /> 
     24    <axis axis_ref="ncatice" /> 
     25  </grid> 
    1926  <!--  --> 
    2027  <grid id="grid_T_3D" > 
     
    2229    <axis axis_ref="deptht" /> 
    2330  </grid> 
     31  <grid id="grid_T_3D_inner" > 
     32    <domain domain_ref="grid_T_inner" /> 
     33    <axis axis_ref="deptht" /> 
     34  </grid> 
    2435  <!--  --> 
    2536  <grid id="grid_T_3DS" > 
     
    2738    <axis axis_ref="profsed" /> 
    2839  </grid> 
     40  <grid id="grid_T_3DS_inner" > 
     41    <domain domain_ref="grid_T_inner" /> 
     42    <axis axis_ref="profsed" /> 
     43  </grid> 
    2944  <!--  --> 
    3045  <grid id="grid_U_2D" > 
    3146    <domain domain_ref="grid_U" /> 
    3247  </grid> 
     48  <grid id="grid_U_2D_inner" > 
     49    <domain domain_ref="grid_U_inner" /> 
     50  </grid> 
    3351  <!--  --> 
    3452  <grid id="grid_U_3D" > 
     
    3654    <axis axis_ref="depthu" /> 
    3755  </grid> 
     56  <grid id="grid_U_3D_inner" > 
     57    <domain domain_ref="grid_U_inner" /> 
     58    <axis axis_ref="depthu" /> 
     59  </grid> 
    3860  <!--  --> 
    3961  <grid id="grid_V_2D" > 
    4062    <domain domain_ref="grid_V" /> 
    4163  </grid> 
     64  <grid id="grid_V_2D_inner" > 
     65    <domain domain_ref="grid_V_inner" /> 
     66  </grid> 
    4267  <!--  --> 
    4368  <grid id="grid_V_3D" > 
    4469    <domain domain_ref="grid_V" /> 
     70    <axis axis_ref="depthv" /> 
     71  </grid> 
     72  <grid id="grid_V_3D_inner" > 
     73    <domain domain_ref="grid_V_inner" /> 
    4574    <axis axis_ref="depthv" /> 
    4675  </grid> 
     
    4877  <grid id="grid_W_2D" > 
    4978    <domain domain_ref="grid_W" /> 
     79  </grid> 
     80  <grid id="grid_W_2D_inner" > 
     81    <domain domain_ref="grid_W_inner" /> 
    5082  </grid> 
    5183  <!--  --> 
     
    5486    <axis axis_ref="depthw" /> 
    5587  </grid> 
     88  <grid id="grid_W_3D_inner" > 
     89    <domain domain_ref="grid_W_inner" /> 
     90    <axis axis_ref="depthw" /> 
     91  </grid> 
    5692  <!--  --> 
    5793  <grid id="grid_F_2D" > 
    5894    <domain domain_ref="grid_F" /> 
     95  </grid> 
     96  <grid id="grid_F_2D_inner" > 
     97    <domain domain_ref="grid_F_inner" /> 
    5998  </grid> 
    6099  <!--  --> 
     
    63102    <axis axis_ref="depthf" /> 
    64103  </grid> 
     104  <grid id="grid_F_3D_inner" > 
     105    <domain domain_ref="grid_F_inner" /> 
     106    <axis axis_ref="depthf" /> 
     107  </grid> 
    65108  <!--  --> 
    66109  <grid id="grid_1point" > 
     
    77120  </grid> 
    78121  <!--  --> 
    79  
    80122 
    81123  <grid id="grid_znl_T_2D"> 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/namelist_ice_ref

    r14247 r15574  
    7373!------------------------------------------------------------------------------ 
    7474          ! -- ice_rdgrft_strength -- ! 
    75    ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                       
     75   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A) 
    7676      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2] 
    7777      rn_crhg       =  20.0           !     ice strength conc. parameter (-) 
     78   ln_str_R75       = .false.         !  ice strength param.: Rothrock_75 => P = fn of potential energy 
     79      rn_pe_rdg     =  17.0           !     coef accouting for frictional dissipation 
     80   ln_str_CST       = .false.         !  ice strength param.: Constant 
     81      rn_str        =   0.0           !     ice strength value 
     82   ln_str_smooth    = .true.          !  spatial smoothing of the ice strength 
    7883                   ! -- ice_rdgrft -- ! 
     84   ln_distf_lin     = .true.          !  redistribution function of ridged ice: linear (Hibler 1980) 
     85   ln_distf_exp     = .false.         !  redistribution function of ridged ice: exponential => not coded yet 
     86      rn_murdg      =   3.0           !     e-folding scale of ridged ice (m**.5) 
    7987   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging 
    8088              ! -- ice_rdgrft_prep -- ! 
     
    139147   ln_cndflx        = .false.         !  Use conduction flux as surface boundary conditions (i.e. for Jules coupling) 
    140148      ln_cndemulate = .false.         !     emulate conduction flux (if not provided in the inputs) 
    141    nn_qtrice        =   1             !  Solar flux transmitted thru the surface scattering layer: 
     149   nn_qtrice        =   0             !  Solar flux transmitted thru the surface scattering layer: 
    142150                                      !     = 0  Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow)  
    143151                                      !     = 1  Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 
     
    268276   rn_alb_imlt      =   0.50          !  bare puddled ice albedo :  0.49 -- 0.58 
    269277   rn_alb_dpnd      =   0.27          !  ponded ice albedo       :  0.10 -- 0.30  
     278   rn_alb_hpiv      =   1.00          !  pivotal ice thickness in m (above which albedo is constant) 
    270279/ 
    271280!------------------------------------------------------------------------------ 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/namelist_pisces_ref

    r12845 r15574  
    1414!----------------------------------------------------------------------- 
    1515  ln_p2z    = .false.        !  LOBSTER model used 
    16   ln_p4z    = .true.         !  PISCES model used 
    17   ln_p5z    = .false.        !  PISCES QUOTA model used 
     16  ln_p4z    = .true.        !  PISCES model used 
     17  ln_p5z    = .false.         !  PISCES QUOTA model used 
    1818  ln_ligand = .false.        !  Enable  organic ligands 
    1919  ln_sediment = .false.      !  Enable sediment module 
     
    3434!              !  file name   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    3535!              !              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    36    sn_patm     = 'presatm'    ,     -1            , 'patm'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    37    sn_atmco2   = 'presatmco2' ,     -1            , 'xco2'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     36   sn_patm     = 'presatm.orca2'    ,     24.           , 'presatm'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     37   sn_atmco2   = 'presatmco2' ,     -1.           , 'xco2'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    3838   cn_dir      = './'     !  root directory for the location of the dynamical files 
    3939! 
     
    4444&nampisbio     !   biological parameters 
    4545!----------------------------------------------------------------------- 
    46    nrdttrc    =  1        ! time step frequency for biology 
    47    wsbio      =  2.       ! POC sinking speed 
    48    xkmort     =  2.E-7    ! half saturation constant for mortality 
    49    ferat3     =  10.E-6   ! Fe/C in zooplankton  
    50    wsbio2     =  50.      ! Big particles sinking speed 
    51    wsbio2max  =  50.      ! Big particles maximum sinking speed 
    52    wsbio2scale =  5000.    ! Big particles length scale of sinking 
     46   nrdttrc     =  1       ! time step frequency for biology 
     47   wsbio       =  2.      ! POC sinking speed 
     48   xkmort      =  1.E-7   ! half saturation constant for mortality 
     49   feratz      =  10.E-6  ! Fe/C in zooplankton  
     50   feratm      =  15.E-6  ! Fe/C in mesozooplankton 
     51   wsbio2      =  50.     ! Big particles sinking speed 
     52   wsbio2max   =  50.     ! Big particles maximum sinking speed 
     53   wsbio2scale =  5000.   ! Big particles length scale of sinking 
    5354!                         !  ln_ligand enabled 
    54    ldocp      =  1.E-4    ! Phyto ligand production per unit doc  
    55    ldocz      =  1.E-4    ! Zoo ligand production per unit doc  
    56    lthet      =  1.0      ! Proportional loss of ligands due to Fe uptake  
     55   ldocp       =  1.E-4   ! Phyto ligand production per unit doc  
     56   ldocz       =  1.E-4   ! Zoo ligand production per unit doc  
     57   lthet       =  1.0     ! Proportional loss of ligands due to Fe uptake  
    5758!                         !  ln_p5z enabled 
    58    no3rat3    =  0.182    ! N/C ratio in zooplankton 
    59    po4rat3    =  0.0094  ! P/C ratio in zooplankton 
     59   no3rat3     =  0.151   ! N/C ratio in zooplankton 
     60   po4rat3     =  0.00943 ! P/C ratio in zooplankton 
    6061/ 
    6162!----------------------------------------------------------------------- 
     
    6465   concnno3   =  1.e-6    ! Nitrate half saturation of nanophytoplankton 
    6566   concdno3   =  3.E-6    ! Nitrate half saturation for diatoms 
    66    concnnh4   =  1.E-7    ! NH4 half saturation for phyto 
    67    concdnh4   =  3.E-7    ! NH4 half saturation for diatoms 
    68    concnfer   =  1.E-9    ! Iron half saturation for phyto 
    69    concdfer   =  3.E-9    ! Iron half saturation for diatoms 
    70    concbfe    =  1.E-11   ! Iron half-saturation for DOC remin. 
    71    concbnh4   =  2.E-8    ! NH4 half saturation for DOC remin. 
    72    concbno3   =  2.E-7    ! Nitrate half saturation for DOC remin. 
     67   concnnh4   =  1.E-6    ! NH4 half saturation for phyto 
     68   concdnh4   =  3.E-6    ! NH4 half saturation for diatoms 
     69   concnfer   =  1.7E-9   ! Iron half saturation for phyto 
     70   concdfer   =  5.E-9    ! Iron half saturation for diatoms 
     71   concbfe    =  3.E-11   ! Iron half-saturation for DOC remin. 
     72   concbnh4   =  3.E-7    ! NH4 half saturation for DOC remin. 
     73   concbno3   =  3.E-7    ! Nitrate half saturation for DOC remin. 
    7374   xsizedia   =  1.E-6    ! Minimum size criteria for diatoms 
    7475   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
    7576   xsizern    =  3.0      ! Size ratio for nanophytoplankton 
    76    xsizerd    =  3.0      ! Size ratio for diatoms 
    77    xksi1      =  2.E-6    ! half saturation constant for Si uptake 
     77   xsizerd    =  4.0      ! Size ratio for diatoms 
     78   xksi1      =  8.E-6    ! half saturation constant for Si uptake 
    7879   xksi2      =  20E-6    ! half saturation constant for Si/C 
    7980   xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
    80    qnfelim    =  7.E-6    ! Optimal quota of phyto 
    81    qdfelim    =  7.E-6    ! Optimal quota of diatoms 
     81   qnfelim    =  10.E-6   ! Optimal quota of phyto 
     82   qdfelim    =  10.E-6   ! Optimal quota of diatoms 
     83   caco3r     =  0.2      ! mean rain ratio 
     84   oxymin     =  1.E-6    ! Half-saturation constant for anoxia 
     85/ 
     86!----------------------------------------------------------------------- 
     87&namp5zlim     !   parameters for nutrient limitations PISCES QUOTA    - ln_p5z 
     88!----------------------------------------------------------------------- 
     89   concnno3   =  2e-6     ! Nitrate half saturation of nanophytoplankton 
     90   concpno3   =  7e-7     ! Nitrate half saturation of picophytoplankton 
     91   concdno3   =  3E-6     ! Phosphate half saturation for diatoms 
     92   concnnh4   =  2E-6     ! NH4 half saturation for phyto 
     93   concpnh4   =  7E-7     ! NH4 half saturation for picophytoplankton 
     94   concdnh4   =  3E-6     ! NH4 half saturation for diatoms 
     95   concnpo4   =  2E-6     ! PO4 half saturation for phyto 
     96   concppo4   =  7E-7     ! PO4 half saturation for picophytoplankton 
     97   concdpo4   =  3E-6     ! PO4 half saturation for diatoms 
     98   concnfer   =  3E-9     ! Iron half saturation for phyto 
     99   concpfer   =  1E-9     ! Iron half saturation for picophytoplankton 
     100   concdfer   =  4.5E-9   ! Iron half saturation for diatoms 
     101   concbfe    =  3E-11    ! Half-saturation for Fe limitation of Bacteria 
     102   concbnh4   =  4.E-7    ! NH4 half saturation for phyto 
     103   concbno3   =  4.E-7    ! Phosphate half saturation for diatoms 
     104   concbpo4   =  4.E-7    ! Phosphate half saturation for bacteria 
     105   xsizedia   =  1.E-6    ! Minimum size criteria for diatoms 
     106   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
     107   xsizepic   =  5.E-7    ! Minimum size criteria for picophyto 
     108   xsizern    =  3.0      ! Size ratio for nanophytoplankton 
     109   xsizerp    =  2.0      ! Size ratio for picophytoplankton 
     110   xsizerd    =  4.0      ! Size ratio for diatoms 
     111   xksi1      =  8.E-6    ! half saturation constant for Si uptake 
     112   xksi2      =  20E-6    ! half saturation constant for Si/C 
     113   xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
    82114   caco3r     =  0.3      ! mean rain ratio 
    83115   oxymin     =  1.E-6    ! Half-saturation constant for anoxia 
    84116/ 
    85117!----------------------------------------------------------------------- 
    86 &namp5zlim     !   parameters for nutrient limitations PISCES QUOTA    - ln_p5z 
    87 !----------------------------------------------------------------------- 
    88    concnno3   =  3e-6     ! Nitrate half saturation of nanophytoplankton 
    89    concpno3   =  1e-6 
    90    concdno3   =  4E-6     ! Phosphate half saturation for diatoms 
    91    concnnh4   =  1.5E-6   ! NH4 half saturation for phyto 
    92    concpnh4   =  4E-7 
    93    concdnh4   =  2E-6     ! NH4 half saturation for diatoms 
    94    concnpo4   =  3E-6     ! PO4 half saturation for phyto 
    95    concppo4   =  1.5E-6 
    96    concdpo4   =  4E-6     ! PO4 half saturation for diatoms 
    97    concnfer   =  3E-9   ! Iron half saturation for phyto 
    98    concpfer   =  1.5E-9 
    99    concdfer   =  4E-9   ! Iron half saturation for diatoms 
    100    concbfe    =  1.E-11   ! Half-saturation for Fe limitation of Bacteria 
    101    concbnh4   =  1.E-7    ! NH4 half saturation for phyto 
    102    concbno3   =  5.E-7    ! Phosphate half saturation for diatoms 
    103    concbpo4   =  1E-7     ! Phosphate half saturation for bacteria 
    104    xsizedia   =  1.E-6    ! Minimum size criteria for diatoms 
    105    xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
    106    xsizepic   =  1.E-6 
    107    xsizern    =  1.0      ! Size ratio for nanophytoplankton 
    108    xsizerp    =  1.0 
    109    xsizerd    =  4.0      ! Size ratio for diatoms 
    110    xksi1      =  2.E-6    ! half saturation constant for Si uptake 
    111    xksi2      =  20E-6  ! half saturation constant for Si/C 
    112    xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
    113    caco3r     =  0.35     ! mean rain ratio 
    114    oxymin     =  1.E-6    ! Half-saturation constant for anoxia 
    115 / 
    116 !----------------------------------------------------------------------- 
    117118&namp5zquota   !   parameters for nutrient limitations PISCES quota    - ln_p5z 
    118119!----------------------------------------------------------------------- 
    119    qfnopt     =  7.E-6    ! Optimal Fe quota of nanophyto 
    120    qfpopt     =  7.E-6    ! Optimal Fe quota of picophyto 
    121    qfdopt     =  7.E-6    ! Optimal quota of diatoms 
    122    qnnmin     =  0.29     ! Minimal N quota for nano 
    123    qnnmax     =  1.39     ! Maximal N quota for nano 
    124    qpnmin     =  0.28     ! Minimal P quota for nano 
    125    qpnmax     =  1.06     ! Maximal P quota for nano 
    126    qnpmin     =  0.42     ! Minimal N quota for pico 
     120   qfnopt     =  12.E-6   ! Optimal Fe quota of nanophyto 
     121   qfpopt     =  12.E-6   ! Optimal Fe quota of picophyto 
     122   qfdopt     =  12.E-6   ! Optimal quota of diatoms 
     123   qnnmin     =  0.61     ! Minimal N quota for nano 
     124   qnnmax     =  1.25     ! Maximal N quota for nano 
     125   qpnmin     =  0.24     ! Minimal P quota for nano 
     126   qpnmax     =  1.35     ! Maximal P quota for nano 
     127   qnpmin     =  1.02     ! Minimal N quota for pico 
    127128   qnpmax     =  1.39     ! Maximal N quota for pico 
    128    qppmin     =  0.25     ! Minimal P quota for pico 
    129    qppmax     =  0.7      ! Maximal P quota for pico 
    130    qndmin     =  0.25     ! Minimal N quota for diatoms 
    131    qndmax     =  1.39     ! Maximal N quota for diatoms 
    132    qpdmin     =  0.29     ! Minimal P quota for diatoms 
    133    qpdmax     =  1.32     ! Maximal P quota for diatoms 
    134    qfnmax     =  40E-6    ! Maximal Fe quota for nano 
    135    qfpmax     =  40E-6    ! Maximal Fe quota for pico 
    136    qfdmax     =  40E-6    ! Maximal Fe quota for diatoms 
     129   qppmin     =  0.19     ! Minimal P quota for pico 
     130   qppmax     =  1.15     ! Maximal P quota for pico 
     131   qndmin     =  0.51     ! Minimal N quota for diatoms 
     132   qndmax     =  1.25     ! Maximal N quota for diatoms 
     133   qpdmin     =  0.24     ! Minimal P quota for diatoms 
     134   qpdmax     =  1.525    ! Maximal P quota for diatoms 
     135   qfnmax     =  60E-6    ! Maximal Fe quota for nano 
     136   qfpmax     =  60E-6    ! Maximal Fe quota for pico 
     137   qfdmax     =  60E-6    ! Maximal Fe quota for diatoms 
    137138/ 
    138139!----------------------------------------------------------------------- 
     
    141142!              !  file name       ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    142143!              !                  !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    143    sn_par      = 'par.orca'       ,     24            , 'fr_par'  ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     144   sn_par      = 'par.orca'       ,     24.           , 'fr_par'  ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    144145   cn_dir      = './'      !  root directory for the location of the dynamical files 
    145146   ln_varpar   =  .true.   ! boolean for PAR variable 
    146147   parlux      =  0.43      ! Fraction of shortwave as PAR 
     148   ln_p4z_dcyc =  .false.  ! Diurnal cycle in PISCES 
    147149/  
    148150!----------------------------------------------------------------------- 
     
    157159   chlcnm     =  0.033    ! Maximum Chl/C in nanophytoplankton 
    158160   chlcdm     =  0.05     ! Maximum Chl/C in diatoms 
    159    chlcmin    =  0.004    ! Minimum Chl/c in phytoplankton 
    160    fecnm      =  40E-6    ! Maximum Fe/C in nanophytoplankton 
    161    fecdm      =  40E-6    ! Maximum Fe/C in diatoms 
    162    grosip     =  0.159    ! mean Si/C ratio 
     161   chlcmin    =  0.003    ! Minimum Chl/c in phytoplankton 
     162   fecnm      =  60E-6    ! Maximum Fe/C in nanophytoplankton 
     163   fecdm      =  60E-6    ! Maximum Fe/C in diatoms 
     164   grosip     =  0.13     ! mean Si/C ratio 
    163165/ 
    164166!----------------------------------------------------------------------- 
    165167&namp5zprod    !   parameters for phytoplankton growth for PISCES quota- ln_p5z 
    166168!----------------------------------------------------------------------- 
    167    pislopen   =  3.       ! P-I slope 
    168    pislopep   =  3.       ! P-I slope for picophytoplankton 
    169    pisloped   =  3.       ! P-I slope  for diatoms 
     169   pislopen   =  5        ! P-I slope of nanophytoplankton 
     170   pislopep   =  5        ! P-I slope for picophytoplankton 
     171   pisloped   =  5        ! P-I slope  for diatoms 
    170172   excretn    =  0.05     ! excretion ratio of phytoplankton 
    171173   excretp    =  0.05     ! excretion ratio of picophytoplankton 
     
    173175   xadap      =  0.       ! Adaptation factor to low light 
    174176   bresp      =  0.02     ! Basal respiration rate 
    175    thetannm   =  0.25     ! Maximum Chl/N in nanophytoplankton 
    176    thetanpm   =  0.25     ! Maximum Chl/N in picophytoplankton 
    177    thetandm   =  0.3      ! Maximum Chl/N in diatoms 
    178    chlcmin    =  0.004    ! Minimum Chl/c in phytoplankton 
    179    grosip     =  0.131    ! mean Si/C ratio 
     177   thetannm   =  0.3      ! Maximum Chl/N in nanophytoplankton 
     178   thetanpm   =  0.3      ! Maximum Chl/N in picophytoplankton 
     179   thetandm   =  0.4      ! Maximum Chl/N in diatoms 
     180   chlcmin    =  0.003    ! Minimum Chl/c in phytoplankton 
     181   grosip     =  0.12     ! mean Si/C ratio 
    180182/ 
    181183!----------------------------------------------------------------------- 
    182184&namp4zmort    !   parameters for phytoplankton sinks for PISCES std   - ln_p4z 
    183185!----------------------------------------------------------------------- 
    184    wchl       =  0.01     ! quadratic mortality of phytoplankton 
    185    wchld      =  0.01     ! maximum quadratic mortality of diatoms 
    186    wchldm     =  0.03     ! maximum quadratic mortality of diatoms 
    187    mprat      =  0.01     ! phytoplankton mortality rate 
    188    mprat2     =  0.01     ! Diatoms mortality rate 
     186   wchln      =  0.01     ! quadratic mortality of phytoplankton 
     187   wchld      =  0.03     ! maximum quadratic mortality of diatoms 
     188   mpratn     =  0.01     ! phytoplankton mortality rate 
     189   mpratd     =  0.01     ! Diatoms mortality rate 
    189190/ 
    190191!----------------------------------------------------------------------- 
     
    193194   wchln      =  0.01     ! quadratic mortality of nanophytoplankton 
    194195   wchlp      =  0.01     ! quadratic mortality of picophytoplankton 
    195    wchld      =  0.01     ! maximum quadratic mortality of diatoms 
    196    wchldm     =  0.02     ! maximum quadratic mortality of diatoms 
     196   wchld      =  0.03     ! maximum quadratic mortality of diatoms 
    197197   mpratn     =  0.01     ! nanophytoplankton mortality rate 
    198198   mpratp     =  0.01     ! picophytoplankton mortality rate 
    199    mprat2     =  0.01     ! Diatoms mortality rate 
     199   mpratd     =  0.01     ! Diatoms mortality rate 
    200200/ 
    201201!----------------------------------------------------------------------- 
    202202&namp4zmes     !   parameters for mesozooplankton for PISCES std       - ln_p4z 
    203203!----------------------------------------------------------------------- 
    204    part2      =  0.75     ! part of calcite not dissolved in mesozoo guts 
    205    grazrat2   =  0.75     ! maximal mesozoo grazing rate 
    206    resrat2    =  0.005    ! exsudation rate of mesozooplankton 
    207    mzrat2     =  0.03     ! mesozooplankton mortality rate 
    208    xpref2d    =  1.       ! mesozoo preference for diatoms 
    209    xpref2n    =  0.3      ! mesozoo preference for nanophyto. 
    210    xpref2z    =  1.       ! mesozoo preference for microzoo. 
    211    xpref2c    =  0.3      ! mesozoo preference for poc 
    212    xthresh2zoo = 1E-8     ! zoo feeding threshold for mesozooplankton  
    213    xthresh2dia = 1E-8     ! diatoms feeding threshold for mesozooplankton  
    214    xthresh2phy = 1E-8     ! nanophyto feeding threshold for mesozooplankton  
    215    xthresh2poc = 1E-8     ! poc feeding threshold for mesozooplankton  
    216    xthresh2   =  3E-7     ! Food threshold for grazing 
    217    xkgraz2    =  20.E-6   ! half saturation constant for meso grazing 
    218    epsher2    =  0.35     ! Efficicency of Mesozoo growth 
    219    epsher2min =  0.35     ! Minimum efficiency of mesozoo growth 
    220    sigma2     =  0.6      ! Fraction of mesozoo excretion as DOM 
    221    unass2     =  0.3      ! non assimilated fraction of P by mesozoo 
    222    grazflux   =  3.e3     ! flux-feeding rate 
     204   part2       =  0.75     ! part of calcite not dissolved in mesozoo guts 
     205   grazrat2    =  0.5      ! maximal mesozoo grazing rate 
     206   resrat2     =  0.005    ! exsudation rate of mesozooplankton 
     207   mzrat2      =  0.01     ! mesozooplankton mortality rate 
     208   xpref2d     =  1.       ! mesozoo preference for diatoms 
     209   xpref2n     =  0.3      ! mesozoo preference for nanophyto. 
     210   xpref2z     =  1.       ! mesozoo preference for microzoo. 
     211   xpref2c     =  0.3      ! mesozoo preference for poc 
     212   xthresh2zoo =  1E-8     ! zoo feeding threshold for mesozooplankton  
     213   xthresh2dia =  1E-8     ! diatoms feeding threshold for mesozooplankton  
     214   xthresh2phy =  1E-8     ! nanophyto feeding threshold for mesozooplankton  
     215   xthresh2poc =  1E-8     ! poc feeding threshold for mesozooplankton  
     216   xthresh2    =  3E-7     ! Food threshold for grazing 
     217   xkgraz2     =  20.E-6   ! half saturation constant for meso grazing 
     218   epsher2     =  0.4      ! Efficicency of Mesozoo growth 
     219   epsher2min  =  0.4      ! Minimum efficiency of mesozoo growth 
     220   sigma2      =  0.6      ! Fraction of mesozoo excretion as DOM 
     221   unass2      =  0.3      ! non assimilated fraction of P by mesozoo 
     222   grazflux    =  3.e3     ! flux-feeding rate 
     223   xsigma2     =  0.5      ! Predation window size 
     224   xsigma2del  =  1.0      ! Predation window size scaling 
     225   ln_dvm_meso =  .false.  ! Activates DVM for mesozooplankton 
     226   xfracmig    =  0.3      ! Fraction of mesozooplankton performing DVM 
    223227/ 
    224228!----------------------------------------------------------------------- 
    225229&namp5zmes     !   parameters for mesozooplankton 
    226230!----------------------------------------------------------------------- 
    227    part2      =  0.75     ! part of calcite not dissolved in mesozoo guts 
    228    grazrat2   =  0.85     ! maximal mesozoo grazing rate 
    229    bmetexc2   =  .true.   ! Metabolic use of excess carbon  
    230    resrat2    =  0.005    ! exsudation rate of mesozooplankton 
    231    mzrat2     =  0.02     ! mesozooplankton mortality rate 
    232    xpref2d    =  1.       ! zoo preference for phyto 
    233    xpref2p    =  1.       ! zoo preference for POC 
    234    xpref2z    =  1.       ! zoo preference for zoo 
    235    xpref2m    =  0.2      ! meso preference for zoo 
    236    xpref2c    =  0.3      ! zoo preference for poc 
    237    xthresh2zoo = 1E-8     ! zoo feeding threshold for mesozooplankton 
    238    xthresh2dia = 1E-8     ! diatoms feeding threshold for mesozooplankton 
    239    xthresh2phy = 1E-8     ! nanophyto feeding threshold for mesozooplankton 
    240    xthresh2mes = 1E-8     ! meso feeding threshold for mesozooplankton 
    241    xthresh2poc = 1E-8     ! poc feeding threshold for mesozooplankton 
     231   part2       =  0.75     ! part of calcite not dissolved in mesozoo guts 
     232   grazrat2    =  0.5      ! maximal mesozoo grazing rate 
     233   bmetexc2    =  .true.   ! Metabolic use of excess carbon  
     234   resrat2     =  0.005    ! exsudation rate of mesozooplankton 
     235   mzrat2      =  0.01     ! mesozooplankton mortality rate 
     236   xpref2d     =  1.       ! meso preference for diatoms 
     237   xpref2n     =  0.3      ! meso preference for nano 
     238   xpref2z     =  1.       ! meso preference for zoo 
     239   xpref2m     =  0.       ! meso preference for zoo 
     240   xpref2c     =  0.3      ! meso preference for poc 
     241   xthresh2zoo =  1E-8     ! zoo feeding threshold for mesozooplankton 
     242   xthresh2dia =  1E-8     ! diatoms feeding threshold for mesozooplankton 
     243   xthresh2phy =  1E-8     ! nanophyto feeding threshold for mesozooplankton 
     244   xthresh2mes =  1E-8     ! meso feeding threshold for mesozooplankton 
     245   xthresh2poc =  1E-8     ! poc feeding threshold for mesozooplankton 
    242246   xthresh2    =  3E-7     ! Food threshold for grazing 
    243247   xkgraz2     =  20.E-6   ! half sturation constant for meso grazing 
    244248   epsher2     =  0.5      ! Efficicency of Mesozoo growth 
    245    epsher2min  =  0.2     ! Minimum efficiency of mesozoo growth 
    246    ssigma2     =  0.5     ! Fraction excreted as semi-labile DOM 
    247    srespir2    =  0.2     ! Active respiration 
    248    unass2c     =  0.3     ! non assimilated fraction of P by mesozoo 
    249    unass2n     =  0.3     ! non assimilated fraction of N by mesozoo 
    250    unass2p     =  0.3     ! non assimilated fraction of P by mesozoo 
    251    grazflux   =  3.e3     ! flux-feeding rate 
     249   epsher2min  =  0.5      ! Minimum efficiency of mesozoo growth 
     250   ssigma2     =  0.5      ! Fraction excreted as semi-labile DOM 
     251   srespir2    =  0.2      ! Active respiration 
     252   unass2c     =  0.3      ! non assimilated fraction of C by mesozoo 
     253   unass2n     =  0.3      ! non assimilated fraction of N by mesozoo 
     254   unass2p     =  0.3      ! non assimilated fraction of P by mesozoo 
     255   xsigma2     =  0.5      ! Predation window size 
     256   xsigma2del  =  1.0      ! Predation window size scaling 
     257   grazflux    =  3.e3     ! flux-feeding rate 
     258   ln_dvm_meso =  .false.  ! Activates DVM for mesozooplankton 
     259   xfracmig    =  0.25     ! Fraction of mesozooplankton performing DVM 
    252260/ 
    253261!----------------------------------------------------------------------- 
    254262&namp4zzoo     !   parameters for microzooplankton for PISCES std      - ln_p4z 
    255263!----------------------------------------------------------------------- 
    256    part       =  0.5      ! part of calcite not dissolved in microzoo guts 
    257    grazrat    =  3.0      ! maximal zoo grazing rate 
    258    resrat     =  0.03     ! exsudation rate of zooplankton 
    259    mzrat      =  0.004    ! zooplankton mortality rate 
    260    xprefc     =  0.1      ! Microzoo preference for POM 
     264   part       =  0.75     ! part of calcite not dissolved in microzoo guts 
     265   grazrat    =  2.0      ! maximal zoo grazing rate 
     266   resrat     =  0.02     ! Linear mortality rate of zooplankton 
     267   mzrat      =  0.005    ! zooplankton mortality rate 
     268   xprefc     =  0.15      ! Microzoo preference for POM 
    261269   xprefn     =  1.       ! Microzoo preference for Nanophyto 
    262    xprefd     =  0.6      ! Microzoo preference for Diatoms 
     270   xprefd     =  0.8      ! Microzoo preference for Diatoms 
    263271   xthreshdia =  1.E-8    ! Diatoms feeding threshold for microzooplankton  
    264272   xthreshphy =  1.E-8    ! Nanophyto feeding threshold for microzooplankton  
     
    266274   xthresh    =  3.E-7    ! Food threshold for feeding 
    267275   xkgraz     =  20.E-6   ! half sturation constant for grazing 
    268    epsher     =  0.3      ! Efficiency of microzoo growth 
    269    epshermin  =  0.3      ! Minimum efficiency of microzoo growth 
     276   epsher     =  0.4     ! Efficiency of microzoo growth 
     277   epshermin  =  0.4     ! Minimum efficiency of microzoo growth 
    270278   sigma1     =  0.6      ! Fraction of microzoo excretion as DOM 
    271279   unass      =  0.3      ! non assimilated fraction of phyto by zoo 
     280   xsigma     =  0.5      ! Predation window size 
     281   xsigmadel  =  1.0      ! Predation window size scaling 
    272282/ 
    273283!----------------------------------------------------------------------- 
    274284&namp5zzoo     !   parameters for microzooplankton 
    275285!----------------------------------------------------------------------- 
    276    part       =  0.5      ! part of calcite not dissolved in microzoo gutsa 
    277    grazrat    =  2.75     ! maximal zoo grazing rate 
     286   part       =  0.75     ! part of calcite not dissolved in microzoo gutsa 
     287   grazrat    =  2.0      ! maximal zoo grazing rate 
    278288   bmetexc    =  .true.   ! Metabolic use of excess carbon 
    279    resrat     =  0.03     ! exsudation rate of zooplankton 
     289   resrat     =  0.02     ! exsudation rate of zooplankton 
    280290   mzrat      =  0.005    ! zooplankton mortality rate 
    281    xprefc     =  0.1      ! Microzoo preference for POM 
    282    xprefn     =  1.       ! Microzoo preference for Nanophyto 
    283    xprefp     =  1.6      ! Microzoo preference for picophyto 
     291   xprefc     =  0.15     ! Microzoo preference for POM 
     292   xprefn     =  1.0      ! Microzoo preference for Nanophyto 
     293   xprefp     =  1.0      ! Microzoo preference for picophyto 
    284294   xprefd     =  1.0      ! Microzoo preference for Diatoms 
    285    xprefz     =  0.3      ! Microzoo preference for microzooplankton 
     295   xprefz     =  0.       ! Microzoo preference for microzooplankton 
    286296   xthreshdia =  1.E-8    ! Diatoms feeding threshold for microzooplankton 
    287297   xthreshphy =  1.E-8    ! Nanophyto feeding threshold for microzooplankton 
    288    xthreshpic =  1.E-8 
    289    xthreshzoo =  1.E-8    ! Nanophyto feeding threshold for microzooplankton 
     298   xthreshpic =  1.E-8    ! Picophyto feeding threshold for microzooplankton 
     299   xthreshzoo =  1.E-8    ! Microzoo feeding threshold for microzooplankton 
    290300   xthreshpoc =  1.E-8    ! POC feeding threshold for microzooplankton 
    291301   xthresh    =  3.E-7    ! Food threshold for feeding 
    292    xkgraz     =  20.E-6   ! half sturation constant for grazing 
     302   xkgraz     =  20.E-6   ! half saturation constant for grazing 
    293303   epsher     =  0.5      ! Efficiency of microzoo growth 
    294    epshermin  =  0.2      ! Minimum efficiency of microzoo growth 
     304   epshermin  =  0.5      ! Minimum efficiency of microzoo growth 
    295305   ssigma     =  0.5      ! Fraction excreted as semi-labile DOM 
    296306   srespir    =  0.2      ! Active respiration 
    297307   unassc     =  0.3      ! non assimilated fraction of C by zoo 
    298    unassn     =  0.3      ! non assimilated fraction of C by zoo 
    299    unassp     =  0.3      ! non assimilated fraction of C by zoo 
     308   unassn     =  0.3      ! non assimilated fraction of N by zoo 
     309   unassp     =  0.3      ! non assimilated fraction of P by zoo 
     310   xsigma     =  0.5      ! Predation window size 
     311   xsigmadel  =  1.0      ! Predation window size scaling 
    300312/ 
    301313!----------------------------------------------------------------------- 
    302314&nampisfer     !   parameters for iron chemistry 
    303315!----------------------------------------------------------------------- 
    304    ln_ligvar =  .false.    ! variable ligand concentration 
    305    xlam1     =  0.005     ! scavenging rate of Iron 
    306    xlamdust  =  150.0     ! Scavenging rate of dust 
    307    ligand    =  0.7E-9    ! Ligands concentration  
     316   ln_ligvar =  .false.   ! variable ligand concentration 
     317   xlam1     =  0.02      ! scavenging rate of Iron by biogenic particles 
     318   xlamdust  =  150.0     ! Scavenging rate of Iron by dust 
     319   ligand    =  1E-9     ! Ligands concentration  
    308320   kfep      =  0.01      ! Nanoparticle formation rate constant 
     321   scaveff   =  1.0       ! Fraction of scavenged Fe that goes to POFe 
    309322/ 
    310323!-----------------------------------------------------------------------   
    311324&nampisrem     !   parameters for remineralization 
    312325!----------------------------------------------------------------------- 
    313    xremik    =  0.3       ! remineralization rate of DOC 
    314326   nitrif    =  0.05      ! NH4 nitrification rate 
    315327   xsirem    =  0.003     ! remineralization rate of Si 
    316328   xsiremlab =  0.03      ! fast remineralization rate of Si 
    317329   xsilab    =  0.5       ! Fraction of labile biogenic silica 
    318    feratb    =  10.E-6    ! Fe/C quota in bacteria 
    319    xkferb    =  3E-10     ! Half-saturation constant for bacteria Fe/C 
     330   feratb    =  60.E-6    ! Fe/C quota in bacteria 
     331   xkferb    =  4E-10     ! Half-saturation constant for bacteria Fe/C 
    320332!                         ! ln_p5z 
    321    xremikc   =  0.25      ! remineralization rate of DOC 
    322    xremikn   =  0.35      ! remineralization rate of DON 
    323    xremikp   =  0.4       ! remineralization rate of DOP 
    324 !   feratb    =  20E-6     ! Bacterial Fe/C ratio 
    325 !   xkferb    =  3E-10     ! Half-saturation constant for bact. Fe/C 
     333   xremikc   =  0.4       ! remineralization rate of DOC 
     334   xremikn   =  0.4       ! remineralization rate of DON 
     335   xremikp   =  0.5       ! remineralization rate of DOP 
    326336/ 
    327337!----------------------------------------------------------------------- 
    328338&nampispoc     !   parameters for organic particles 
    329339!----------------------------------------------------------------------- 
    330    xremip    =  0.035     ! remineralisation rate of PON 
     340   xremip    =  0.035     ! remineralisation rate of POC 
    331341   jcpoc     =  15        ! Number of lability classes 
    332342   rshape    =  1.0       ! Shape of the gamma function 
    333343!                         ! ln_p5z 
    334    xremipc   =  0.02      ! remineralisation rate of POC 
    335    xremipn   =  0.025     ! remineralisation rate of PON 
    336    xremipp   =  0.03      ! remineralisation rate of POP 
     344   xremipc   =  0.028     ! remineralisation rate of POC 
     345   xremipn   =  0.03      ! remineralisation rate of PON 
     346   xremipp   =  0.035     ! remineralisation rate of POP 
    337347/ 
    338348!----------------------------------------------------------------------- 
    339349&nampiscal     !   parameters for Calcite chemistry 
    340350!----------------------------------------------------------------------- 
    341    kdca       =  6.       ! calcite dissolution rate constant (1/time) 
    342    nca        =  1.       ! order of dissolution reaction (dimensionless) 
     351   kdca       =  100.     ! calcite dissolution rate constant (1/time) 
     352   nca        =  4.7      ! order of dissolution reaction (dimensionless) 
    343353/ 
    344354!----------------------------------------------------------------------- 
     
    358368   distcoast   =  5.e3     ! Distance off the coast for Iron from sediments 
    359369   mfrac       =  0.035    ! Fe mineral fraction of dust 
    360    wdust       =  2.0      ! Dust sinking speed  
     370   wdust       =  2.0      ! Dust sinking speed 
    361371   icefeinput  =  15.e-9   ! Iron concentration in sea ice 
    362    hratio      =  1.e+7    ! Fe to 3He ratio assumed for vent iron supply  
     372   hratio      =  1.e+7    ! Fe to 3He ratio assumed for vent iron supply 
    363373!                          ! ln_ligand 
    364    lgw_rath    =  0.5      ! Weak ligand ratio from sed hydro sources  
    365 / 
    366 !----------------------------------------------------------------------- 
    367 &nampissed     !   parameters for sediments mobilization 
    368 !----------------------------------------------------------------------- 
    369    nitrfix     =  1.e-7    ! Nitrogen fixation rate 
    370    diazolight  =  50.      ! Diazotrophs sensitivity to light (W/m2) 
     374   lgw_rath    =  0.2      ! Weak ligand ratio from sed hydro sources 
     375/ 
     376!----------------------------------------------------------------------- 
     377&nampislig     !   Namelist parameters for ligands, nampislig 
     378!----------------------------------------------------------------------- 
     379   rlgw        =  300.     ! Lifetime (years) of weak ligands 
     380   rlig        =  1.E-4    ! Remin ligand production per unit C 
     381   prlgw       =  3.E-4    ! Photolysis of weak ligand 
     382   rlgs        =  1.       ! Lifetime (years) of strong ligands 
     383   xklig       =  1.E-9    ! 1/2 saturation constant of photolysis 
     384/ 
     385!----------------------------------------------------------------------- 
     386&nampissed     !   Namelist parameters for sediment mobilisation 
     387!----------------------------------------------------------------------- 
     388   nitrfix     =  2.e-7    ! Nitrogen fixation rate 
     389   diazolight  =  30.      ! Diazotrophs sensitivity to light (W/m2) 
    371390   concfediaz  =  1.e-10   ! Diazotrophs half-saturation Cste for Iron 
    372 / 
    373 !----------------------------------------------------------------------- 
    374 &nampislig     !   Namelist parameters for ligands, nampislig 
    375 !----------------------------------------------------------------------- 
    376    rlgw        =  100.     ! Lifetime (years) of weak ligands 
    377    rlig        =  1.E-4    ! Remin ligand production per unit C 
    378    prlgw       =  1.E-4    ! Photolysis of weak ligand 
    379    rlgs        =  1.       ! Lifetime (years) of strong ligands 
    380391/ 
    381392!----------------------------------------------------------------------- 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/namelist_ref

    r14255 r15574  
    2626!!   namtsd       data: temperature & salinity                          (default: OFF) 
    2727!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T) 
    28 !!   namc1d       1D configuration options                              ("key_c1d") 
    29 !!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
    30 !!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
     28!!   namc1d       1D configuration options                              (ln_c1d =T) 
     29!!   namc1d_dyndmp 1D newtonian damping applied on currents             (ln_c1d =T) 
     30!!   namc1d_uvd   1D data (currents)                                    (ln_c1d =T) 
    3131!!====================================================================== 
    3232! 
     
    7676   ! 
    7777   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
     78   ln_c1d      = .false.   !  Single column domain (1x1pt)              (T => fill namc1d) 
    7879   ! 
    7980   ln_meshmask = .true.   !  =T create a mesh file 
     
    99100!----------------------------------------------------------------------- 
    100101   ln_tile = .false.     !  Use tiling (T) or not (F) 
    101    nn_ltile_i = 10       !  Length of tiles in i 
     102   nn_ltile_i = 99999    !  Length of tiles in i 
    102103   nn_ltile_j = 10       !  Length of tiles in j 
    103104/ 
     
    161162/ 
    162163!----------------------------------------------------------------------- 
    163 &namc1d        !   1D configuration options                             ("key_c1d" default: PAPA station) 
    164 !----------------------------------------------------------------------- 
    165    rn_lat1d    =      50   !  Column latitude 
    166    rn_lon1d    =    -145   !  Column longitude 
    167    ln_c1d_locpt =  .true.   ! Localization of 1D config in a grid (T) or independant point (F) 
    168 / 
    169 !----------------------------------------------------------------------- 
    170 &namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF) 
     164&namc1d        !   1D configuration options                             (ln_c1d =T default: PAPA station) 
     165!----------------------------------------------------------------------- 
     166   rn_lat1d    =    50.1   !  Column latitude 
     167   rn_lon1d    =  -144.9   !  Column longitude 
     168/ 
     169!----------------------------------------------------------------------- 
     170&namc1d_dyndmp !   U & V newtonian damping                              (ln_c1d =T default: OFF) 
    171171!----------------------------------------------------------------------- 
    172172   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    173173/ 
    174174!----------------------------------------------------------------------- 
    175 &namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF) 
     175&namc1d_uvd    !   data: U & V currents                                 (ln_c1d =T default: OFF) 
    176176!----------------------------------------------------------------------- 
    177177   !                       !  =T read U-V fields for: 
     
    531531      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
    532532         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
    533          !                       ! spe = fwfisf is read from a forcing field 
     533         !                       ! spe = fwfisf is read from a forcing field ( melt > 0; freezing < 0 ) 
    534534         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
    535535         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     
    556556      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
    557557         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
    558          !                      ! spe   = fwfisf is read from a forcing field 
     558         !                      ! spe   = fwfisf is read from a forcing field ( melt > 0; freezing < 0 ) 
    559559         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
    560560         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     561         ! 
     562         !* bg03 case 
     563         rn_isfpar_bg03_gt0 = 1.0e-4 ! gamma coeficient used in bg03 paper [m/s] 
     564         ! 
     565         !*** File definition *** 
    561566         ! 
    562567         !* all cases 
     
    566571         sn_isfpar_zmax = 'isfmlt_par',       0.       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    567572         sn_isfpar_zmin = 'isfmlt_par',       0.       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     573         ! 
    568574         !* 'spe' and 'oasis' case 
    569575         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     576         ! 
    570577         !* 'bg03' case 
     578         !* Leff is in [km] 
    571579         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    572580      ! 
     
    644652   rn_test_box             = 108.0,  116.0, -66.0, -58.0 
    645653   ln_use_calving          = .false. ! Use calving data even when nn_test_icebergs > 0 
    646    rn_speed_limit          = 0.      ! CFL speed limit for a berg 
     654   rn_speed_limit          = 0.      ! CFL speed limit for a berg (safe value is 0.4, see #2581) 
    647655   ! 
    648656   ln_M2016                = .false. ! use Merino et al. (2016) modification (use of 3d ocean data instead of only sea surface data) 
     
    686694   ln_vert_remap   = .false. !  use vertical remapping 
    687695   ln_spc_dyn      = .true.  !  use 0 as special value for dynamics 
    688    ln_chk_bathy    = .false. !  =T  check the parent bathymetry 
     696   ln_chk_bathy    = .true. !  =T  check the parent bathymetry 
    689697   rn_sponge_tra   = 0.002   !  coefficient for tracer   sponge layer [] 
    690698   rn_sponge_dyn   = 0.002   !  coefficient for dynamics sponge layer [] 
     
    12341242   !                       !           = 2 roughness uses rn_hsri and is weighted by 1-fr_i 
    12351243   !                       !           = 3 roughness uses rn_hsri and is weighted by 1-MIN(1,4*fr_i) 
     1244   nn_mxlice     =     1   !  mixing under sea ice 
     1245                           !     = 0 No scaling under sea-ice 
     1246                           !     = 1 scaling with constant Ice-ocean roughness (rn_hsri) 
     1247                           !     = 2 scaling with mean sea-ice thickness 
     1248                           !     = 3 scaling with max sea-ice thickness 
    12361249   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
    12371250   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
     
    12651278                               !  = 2:use surface value of SD fit to slope at rn_osm_hblfrac*hbl below surface 
    12661279   ln_zdfosm_ice_shelter = .true.  ! reduce surface SD and depth scale under ice 
    1267    ln_osm_mle = .false.        !  Use integrated FK-OSM model 
     1280   ln_osm_mle = .true.         !  Use integrated FK-OSM model 
    12681281/ 
    12691282!----------------------------------------------------------------------- 
     
    12731286   nn_osm_mle          = 0         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
    12741287   rn_osm_mle_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_osm_mle=0) 
    1275    rn_osm_mle_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_osm_mle=0) 
     1288   rn_osm_mle_time     = 43200.    ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_osm_mle=0) 
    12761289   rn_osm_mle_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
    1277    rn_osm_mle_rho_c =    0.01      ! delta rho criterion used to calculate MLD for FK 
    1278    rn_osm_mle_thresh  = 0.0005     ! delta b criterion used for FK MLE criterion 
    1279    rn_osm_mle_tau     = 172800.    ! time scale for FK-OSM (seconds)  (case rn_osm_mle=0) 
    1280    ln_osm_hmle_limit   = .false.   ! limit hmle to rn_osm_hmle_limit*hbl 
    1281    rn_osm_hmle_limit   = 1.2 
     1290   rn_osm_mle_rho_c    = 0.03      ! delta rho criterion used to calculate MLD for FK 
     1291   rn_osm_mle_thresh   = 0.0001    ! delta b criterion used for FK MLE criterion 
     1292   rn_osm_mle_tau      = 172800.   ! time scale for FK-OSM (seconds)  (case rn_osm_mle=0) 
     1293   ln_osm_hmle_limit   = .true.    ! If true, limit hmle to rn_osm_hmle_limit*hbl 
     1294   rn_osm_hmle_limit   = 1.5 
    12821295   / 
    12831296!----------------------------------------------------------------------- 
     
    14981511   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T 
    14991512   nn_hls      =   1       !  halo width (applies to both rows and columns) 
     1513   nn_comm     =   1       !  comm choice 
    15001514/ 
    15011515!----------------------------------------------------------------------- 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/cfgs/SHARED/namelist_top_ref

    r14032 r15574  
    101101/ 
    102102!----------------------------------------------------------------------- 
     103&namtrc_opt      !  light availability in the water column 
     104!----------------------------------------------------------------------- 
     105!              !  file name       ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     106!              !                  !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     107   sn_par      = 'par.orca'       ,     24            , 'fr_par'  ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     108   cn_dir      = './'        ! root directory for the location of the dynamical files 
     109   ln_varpar   =  .true.     ! Read PAR from file 
     110   parlux      =  0.43       ! Fraction of shortwave as PAR 
     111   light_loc   = 'center'    ! Light location in the water cell ('center', 'integral') 
     112/ 
     113!----------------------------------------------------------------------- 
    103114&namtrc_dmp      !   passive tracer newtonian damping                   (ln_trcdmp=T) 
    104115!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.