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 8881 – NEMO

Changeset 8881


Ignore:
Timestamp:
2017-12-01T18:42:21+01:00 (6 years ago)
Author:
flavoni
Message:

dev_CNRS_2017 branch: merged dev_r7881_ENHANCE09_RK3 with trunk r8864

Location:
branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED
Files:
1 deleted
5 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED/field_def_nemo-lim.xml

    r7767 r8881  
    44    <!--  
    55============================================================================================================ 
    6 =                                  definition of all existing variables                                    = 
     6=                                  definition of all existing ESIM variables                               = 
    77=                                            DO NOT CHANGE                                                 = 
    88============================================================================================================ 
     
    1212      <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc --> 
    1313 
    14          <!-- LIM2(only) fields --> 
    15          <field id="qsr_ai_cea"   long_name="Air-Ice downward solar heat flux (cell average)"              standard_name="surface_downwelling_shortwave_flux_in_air"          unit="W/m2"         /> 
    16          <field id="qns_ai_cea"   long_name="Air-Ice downward non-solar heat flux (cell average)"                                                                             unit="W/m2"         /> 
    17          <field id="qla_ai_cea"   long_name="Air-Ice downward Latent heat flux (cell average)"             standard_name="surface_downward_latent_heat_flux"                  unit="W/m2"         />          
    18          <field id="qsr_io_cea"   long_name="Ice-Oce downward solar heat flux (cell average)"              standard_name="net_downward_shortwave_flux_at_sea_water_surface"   unit="W/m2"         /> 
    19          <field id="qns_io_cea"   long_name="Ice-Oce downward non-solar heat flux (cell average)"                                                                             unit="W/m2"         /> 
    20          <field id="iceprod_cea"  long_name="Ice production (cell average)"                                                                                                   unit="m/s"          /> 
    21          <field id="iiceconc"     long_name="Ice concentration"                                            standard_name="sea_ice_area_fraction"                              unit="1"            /> 
    22          <field id="ice_pres"     long_name="Ice presence"                                                                                                                    unit=""             /> 
    23          <field id="ist_cea"      long_name="Ice surface temperature (cell average)"                       standard_name="surface_temperature"                                unit="degC"         /> 
    24          <field id="ist_ipa"      long_name="Ice surface temperature (ice presence average)"               standard_name="surface_temperature"                                unit="degC"         /> 
    25          <field id="u_imasstr"    long_name="Sea-ice mass transport along i-axis"                          standard_name="sea_ice_x_transport"                                unit="kg/s"         /> 
    26          <field id="v_imasstr"    long_name="Sea-ice mass transport along j-axis"                          standard_name="sea_ice_y_transport"                                unit="kg/s"         /> 
    27  
    28          <!-- LIM3 fields -->   
     14    <!-- ================ --> 
     15         <!-- SIMIP ice fields --> 
     16         <!-- ================ --> 
     17    <field id="icepres"      long_name="Fraction of time steps with sea ice"                     standard_name="sea_ice_time_fraction"                      unit=""             /> 
     18         <field id="icethick"     long_name="Sea-ice thickness per area"                              standard_name="sea_ice_thickness"                          unit="m"   /> 
     19         <field id="icethic"      long_name="Sea-ice thickness"                                       standard_name="sea_ice_thickness"                          unit="m"   /> 
     20         <field id="snothic"      long_name="Snow thickness"                                          standard_name="surface_snow_thickness"                     unit="m"            /> 
     21         <field id="icevolu"      long_name="Sea-ice volume per area"                                 standard_name="sea_ice_thickness"                          unit="m"            /> 
     22         <field id="icefb"        long_name="Sea-ice freeboard"                                       standard_name="sea_ice_freeboard"                          unit="m"            /> 
     23         <field id="iceconc"      long_name="Sea-ice area fraction"                                   standard_name="sea_ice_area_fraction"                      unit=""             /> 
     24         <field id="iceconc_pct"  long_name="Sea-ice area fraction in percent"                        standard_name="sea_ice_area_fraction_in_percent"           unit="%"             > iceconc * 100. </field > 
     25         <field id="icemass"      long_name="Sea-ice mass per area"                                   standard_name="sea_ice_amount"                             unit="kg/m2"        /> 
     26         <field id="snomass"      long_name="Snow mass per area"                                      standard_name="liquid_water_content_of_surface_snow"       unit="kg/m2"        /> 
     27         <field id="iceage"       long_name="Age of sea ice"                                          standard_name="age_of_sea_ice"                             unit="s"            /> 
     28         <field id="icealb"       long_name="Sea-ice or snow albedo"                                  standard_name="sea_ice_albedo"                             unit=""             /> 
     29         <field id="icesmass"     long_name="Mass of salt in sea ice per area"                        standard_name="sea_ice_salt_mass"                          unit="kg/m2"        /> 
     30         <field id="icesal"       long_name="Sea ice salinity"                                        standard_name="sea_ice_salinity"                           unit="g/kg"         /> 
     31 
     32    <!-- heat --> 
     33         <field id="icestK"       long_name="Surface temperature of sea ice"                          standard_name="sea_ice_surface_temperature"                unit="K"            />  
     34         <field id="icesntK"      long_name="Temperature at snow-ice interface"                       standard_name="sea_ice_snow_interface_temperature"         unit="K"            /> 
     35         <field id="icebotK"      long_name="Temperature at ice-ocean interface"                      standard_name="sea_ice_bottom_temperature"                 unit="K"            /> 
     36         <field id="icehcneg"     long_name="Sea-ice heat content per unit area"                      standard_name="integral_of_sea_ice_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" > (-icehc) </field> 
     37         <field id="isnhcneg"     long_name="Snow-heat content per unit area"                         standard_name="thermal_energy_content_of_surface_snow"                              unit="J/m2" /> 
     38 
     39    <!-- fluxes (mass, heat, salt) --> 
     40         <field id="wfxsum"       long_name="Freshwater flux from sea-ice surface"                    standard_name="freshwater_flux_from_ice_surface"             unit="kg/m2/s"      /> 
     41         <field id="wfxtot"       long_name="Freshwater flux from sea ice"                            standard_name="freshwater_flux_from_ice"                     unit="kg/m2/s" />  
     42         <field id="hfxsenso"     long_name="Net sensible heat flux under sea ice"                    standard_name="ice_ocean_heat_flux"                          unit="W/m2" /> 
     43         <field id="hfxconsu"     long_name="Net conductive heat flux in ice at the surface"          standard_name="conductive_heat_flux_at_sea_ice_surface"      unit="W/m2" />  
     44         <field id="hfxconbo"     long_name="Net conductive heat fluxes in ice at the bottom"         standard_name="conductive_heat_flux_at_sea_ice_bottom"       unit="W/m2" /> 
     45         <field id="sfx_mv"       long_name="Salt flux from sea ice"                                  standard_name="salt_flux_from_ice"                           unit="kg/m2/s" /> 
     46 
     47    <!-- trends --> 
     48         <field id="afxthd"       long_name="sea-ice area fraction change from thermodynamics"        standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics"  unit="s-1"   /> 
     49         <field id="afxdyn"       long_name="sea-ice area fraction change from dynamics"              standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics"  unit="s-1"   /> 
     50 
     51         <field id="dmithd"       long_name="sea-ice mass change from thermodynamics"                 standard_name="tendency_of_sea_ice_amount_due_to_thermodynamics"   unit="kg/m2/s" /> 
     52         <field id="dmidyn"       long_name="sea-ice mass change from dynamics"                       standard_name="tendency_of_sea_ice_amount_due_to_dynamics"         unit="kg/m2/s" /> 
     53         <field id="dmiopw"       long_name="sea-ice mass change through growth in supercooled open water (aka frazil)"   standard_name="tendency_of_sea_ice_amount_due_to_freezing_in_open_water"   unit="kg/m2/s" /> 
     54         <field id="dmibog"       long_name="sea-ice mass change through basal growth"                standard_name="tendency_of_sea_ice_amount_due_to_congelation_ice_accumulation"  unit="kg/m2/s" /> 
     55         <field id="dmisni"       long_name="sea-ice mass change through snow-to-ice conversion"      standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion"               unit="kg/m2/s" /> 
     56         <field id="dmtsub"       long_name="snow and sea ice mass change through sublimation"        standard_name="tendency_of_snow_and_ice_amount_due_to_sublimation"              unit="kg/m2/s" /> 
     57         <field id="dmssub"       long_name="sea-ice mass change through evaporation and sublimation" standard_name="water_evaporation_flux"                                     unit="kg/m2/s" /> 
     58         <field id="dmisub"       long_name="snow mass change through evaporation or sublimation"     standard_name="surface_snow_sublimation_flux"           unit="kg/m2/s" /> 
     59         <field id="dmisum"       long_name="sea-ice mass change through surface melting"             standard_name="tendency_of_sea_ice_amount_due_to_surface_melting"               unit="kg/m2/s" /> 
     60         <field id="dmibom"       long_name="sea-ice mass change through bottom melting"              standard_name="tendency_of_sea_ice_amount_due_to_basal_melting"                 unit="kg/m2/s" /> 
     61 
     62         <field id="dmsspr"       long_name="snow mass change through snow fall"                      standard_name="snowfall_flux"                                                   unit="kg/m2/s" /> 
     63         <field id="dmsmel"       long_name="snow mass change through melt"                           standard_name="surface_snow_melt_flux"                                          unit="kg/m2/s" /> 
     64         <field id="dmsdyn"       long_name="snow mass change through advection by sea-ice dynamics"  standard_name="tendency_of_snow_mass_due_to_sea_ice_dynamics"               unit="kg/m2/s" /> 
     65         <field id="dmsssi"       long_name="snow mass change through snow-to-ice conversion"         standard_name="tendency_of_snow_mass_due_to_snow_to_ice_conversion"             unit="kg/m2/s" /> 
     66 
     67    <!-- momentum (rheology) --> 
     68         <field id="uice_mv"      long_name="X-component of sea ice velocity"                         standard_name="sea_ice_x_velocity"                         unit="m/s" />  
     69         <field id="vice_mv"      long_name="Y-component of sea ice velocity"                         standard_name="sea_ice_y_velocity"                         unit="m/s" />       
     70         <field id="icevel_mv"    long_name="Sea-ice speed"                                           standard_name="sea_ice_speed"                              unit="m/s" />     
     71         <field id="utau_ice"     long_name="X-component of atmospheric stress on sea ice"            standard_name="surface_downward_x_stress"                    unit="N/m2"    /> 
     72         <field id="vtau_ice"     long_name="Y-component of atmospheric stress on sea ice"            standard_name="surface_downward_y_stress"                    unit="N/m2"    /> 
     73         <field id="utau_oi"      long_name="X-component of ocean stress on sea ice"                  standard_name="sea_ice_base_upward_x_stress"                 unit="N/m2"    /> 
     74         <field id="vtau_oi"      long_name="Y-component of ocean stress on sea ice"                  standard_name="sea_ice_base_upward_y_stress"                 unit="N/m2"    /> 
     75         <field id="dssh_dx"      long_name="Sea-surface tilt term in force balance (x-component)"    standard_name="sea_surface_tilt_force_on_sea_ice_x"       unit="N/m2"    /> 
     76         <field id="dssh_dy"      long_name="Sea-surface tilt term in force balance (y-component)"    standard_name="sea_surface_tilt_force_on_sea_ice_y"       unit="N/m2"    /> 
     77         <field id="corstrx"      long_name="Coriolis force term in force balance (x-component)"      standard_name="coriolis_force_on_sea_ice_x"               unit="N/m2"    /> 
     78         <field id="corstry"      long_name="Coriolis force term in force balance (y-component)"      standard_name="coriolis_force_on_sea_ice_y"               unit="N/m2"    /> 
     79         <field id="intstrx"      long_name="Internal stress term in force balance (x-component)"     standard_name="internal_stress_in_sea_ice_x"              unit="N/m2"    /> 
     80         <field id="intstry"      long_name="Internal stress term in force balance (y-component)"     standard_name="internal_stress_in_sea_ice_y"              unit="N/m2"    /> 
     81         <field id="icestr"       long_name="Compressive sea ice strength"                            standard_name="compressive_strength_of_sea_ice"              unit="N/m"     /> 
     82         <field id="normstr"      long_name="Average normal stress in sea ice"                        standard_name="average_normal_stress"                        unit="N/m"     /> 
     83         <field id="sheastr"      long_name="Maximum shear stress in sea ice"                         standard_name="maximum_shear_stress"                         unit="N/m"     /> 
     84         <field id="idive"        long_name="Divergence of the sea-ice velocity field"                standard_name="divergence_of_sea_ice_velocity"               unit="s-1"     /> 
     85         <field id="ishear"       long_name="Maximum shear of sea-ice velocity field"                 standard_name="maximum_shear_of_sea_ice_velocity"            unit="s-1"     /> 
     86 
     87    <!-- momentum (advection) --> 
     88    <field id="xmtrpice"     long_name="X-component of ice mass transport"                       standard_name="ice_x_transport"                              unit="kg/s" /> 
     89         <field id="ymtrpice"     long_name="Y-component of ice mass transport"                       standard_name="ice_y_transport"                              unit="kg/s" /> 
     90         <field id="xmtrpsnw"     long_name="X-component of snw mass transport"                       standard_name="snw_x_transport"                              unit="kg/s" /> 
     91         <field id="ymtrpsnw"     long_name="Y-component of snw mass transport"                       standard_name="snw_y_transport"                              unit="kg/s" /> 
     92         <field id="xatrp"        long_name="X-component of ice area transport"                       standard_name="area_x_transport"                             unit="m2/s" /> 
     93         <field id="yatrp"        long_name="Y-component of ice area transport"                       standard_name="area_y_transport"                             unit="m2/s" /> 
     94         <field id="xmtrptot"     long_name="X-component of sea-ice mass transport"                   standard_name="sea_ice_x_transport"                          unit="kg/s" > xmtrpice + xmtrpsnw </field> 
     95         <field id="ymtrptot"     long_name="Y-component of sea-ice mass transport"                   standard_name="sea_ice_y_transport"                          unit="kg/s" > ymtrpice + ymtrpsnw </field> 
     96 
     97    <!-- categories --> 
     98         <field id="iceconc_cat_mv"     long_name="Sea-ice area fractions in thickness categories"    standard_name="sea_ice_area_fraction_over_categories"      unit=""   grid_ref="grid_T_3D_ncatice" /> 
     99         <field id="iceconc_cat_pct_mv" long_name="Sea-ice area fractions in thickness categories"    standard_name="sea_ice_area_fraction_over_categories"      unit="%"  grid_ref="grid_T_3D_ncatice"  > iceconc_cat_mv * 100. </field > 
     100         <field id="icethic_cat_mv"     long_name="Sea-ice thickness in thickness categories"         standard_name="sea_ice_thickness_over_categories"          unit="m"  grid_ref="grid_T_3D_ncatice" /> 
     101         <field id="snowthic_cat_mv"    long_name="Snow thickness in thickness categories"            standard_name="snow_thickness_over_categories"             unit="m"  grid_ref="grid_T_3D_ncatice" /> 
     102 
     103         <!-- =================== --> 
     104         <!-- standard ice fields -->   
     105         <!-- =================== --> 
    29106         <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            /> 
    30   
    31          <field id="snowthic_cea" long_name="Snow thickness (cell average)"                                standard_name="surface_snow_thickness"                             unit="m"            /> 
    32          <field id="icethic_cea"  long_name="Ice thickness (cell average)"                                 standard_name="sea_ice_thickness"                                  unit="m"            /> 
    33  
    34          <field id="uice_ipa"     long_name="Ice velocity along i-axis at I-point (ice presence average)"  standard_name="sea_ice_x_velocity"                                 unit="m/s"          />       
    35          <field id="vice_ipa"     long_name="Ice velocity along j-axis at I-point (ice presence average)"  standard_name="sea_ice_y_velocity"                                 unit="m/s"          />                
    36          <field id="utau_ice"     long_name="Wind stress along i-axis over the ice at i-point"             standard_name="surface_downward_x_stress"                          unit="N/m2"         /> 
    37          <field id="vtau_ice"     long_name="Wind stress along j-axis over the ice at i-point"             standard_name="surface_downward_y_stress"                          unit="N/m2"         /> 
    38           
    39          <field id="iceconc"      long_name="ice concentration"                                            standard_name="sea_ice_area_fraction"                              unit="%"            /> 
    40           <field id="isst"         long_name="sea surface temperature"                                      standard_name="sea_surface_temperature"                            unit="degC"         /> 
    41          <field id="isss"         long_name="sea surface salinity"                                         standard_name="sea_surface_salinity"                               unit="1e-3"         />  
     107 
     108    <field id="miceage"      long_name="Age of sea ice"                                               unit="s"        /> 
     109         <field id="micesalt"     long_name="Sea ice salinity"                                             unit="g/kg"     /> 
     110         <field id="ibrinv"       long_name="brine volume"                                                 unit="%"        /> 
     111         <field id="icecolf"      long_name="frazil ice collection thickness"                              unit="m"        /> 
     112         <field id="snowvol"      long_name="snow volume"                                                  unit="m"        /> 
     113         <field id="tau_icebfr"   long_name="ice friction on ocean bottom for landfast ice"                unit="N/2"      /> 
     114 
     115    <!-- melt ponds --> 
     116    <field id="iceamp"       long_name="melt pond fraction"                                           standard_name="sea_ice_meltpond_fraction"                          unit="%"            />  
     117         <field id="icevmp"       long_name="melt pond volume"                                             standard_name="sea_ice_meltpond_volume"                            unit="m"            />  
     118     
     119    <!-- momentum (rheology) --> 
     120         <field id="uice"         long_name="X-component of sea ice velocity"                              standard_name="sea_ice_x_velocity"                                 unit="m/s"          />  
     121         <field id="vice"         long_name="Y-component of sea ice velocity"                              standard_name="sea_ice_y_velocity"                                 unit="m/s"          />       
     122         <field id="icevel"       long_name="Sea-ice speed"                                                standard_name="sea_ice_speed"                                      unit="m/s"          /> 
     123         <field id="isig1"        long_name="1st principal stress component for EVP rhg"                   unit="unitless" /> 
     124         <field id="isig2"        long_name="2nd principal stress component for EVP rhg"                   unit="unitless" /> 
     125         <field id="isig3"        long_name="convergence measure for EVP rheology (must be =1)"            unit="unitless" /> 
     126 
     127    <!-- heat --> 
     128         <field id="micet"        long_name="Mean ice temperature"                                         unit="degC"     /> 
     129         <field id="icehc"        long_name="ice total heat content"                                       unit="J/m2"     />  
     130         <field id="isnowhc"      long_name="snow total heat content"                                      unit="J/m2"     /> 
     131         <field id="icest"        long_name="ice surface temperature"                                      unit="degC"     /> 
     132     
     133    <!-- surface heat fluxes --> 
    42134         <field id="qt_oce"       long_name="total flux at ocean surface"                                  standard_name="surface_downward_heat_flux_in_sea_water"            unit="W/m2"         /> 
    43135         <field id="qsr_oce"      long_name="solar heat flux at ocean surface"                             standard_name="net_downward_shortwave_flux_at_sea_water_surface"   unit="W/m2"         /> 
     
    48140         <field id="qtr_ice"      long_name="solar heat flux transmitted through ice: sum over categories"                                                                    unit="W/m2"         /> 
    49141         <field id="qemp_ice"     long_name="Downward Heat Flux from E-P over ice"                                                                                            unit="W/m2"         /> 
    50          <field id="micesalt"     long_name="Mean ice salinity"                                                                                                               unit="1e-3"         /> 
    51          <field id="miceage"      long_name="Mean ice age"                                                                                                                    unit="years"        /> 
    52          <field id="alb_ice"      long_name="Mean albedo over sea ice"                                                                                                        unit=""             /> 
    53142         <field id="albedo"       long_name="Mean albedo over sea ice and ocean"                                                                                              unit=""             /> 
    54143 
    55          <field id="iceage_cat"   long_name="Ice age for categories"                                       unit="days"   grid_ref="grid_T_3D_ncatice" /> 
    56          <field id="iceconc_cat"  long_name="Ice concentration for categories"                             unit="%"      grid_ref="grid_T_3D_ncatice" /> 
    57          <field id="icethic_cat"  long_name="Ice thickness for categories"                                 unit="m"      grid_ref="grid_T_3D_ncatice" /> 
    58          <field id="snowthic_cat" long_name="Snow thicknessi for categories"                               unit="m"      grid_ref="grid_T_3D_ncatice" /> 
    59          <field id="salinity_cat" long_name="Sea-Ice Bulk salinity for categories"                         unit="g/kg"   grid_ref="grid_T_3D_ncatice" /> 
    60          <field id="brinevol_cat" long_name="Brine volume for categories"                                  unit="%"      grid_ref="grid_T_3D_ncatice" /> 
    61          <field id="icetemp_cat"  long_name="Ice temperature for categories"                               unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
    62          <field id="snwtemp_cat"  long_name="Snow temperature for categories"                              unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
    63  
    64          <field id="micet"        long_name="Mean ice temperature"                                         unit="degC"     /> 
    65          <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"    />  
    66          <field id="isnowhc"      long_name="snow total heat content"                                      unit="10^9J"    /> 
    67          <field id="icest"        long_name="ice surface temperature"                                      unit="degC"     /> 
    68          <field id="ibrinv"       long_name="brine volume"                                                 unit="%"        /> 
    69          <field id="icecolf"      long_name="frazil ice collection thickness"                              unit="m"        /> 
    70          <field id="icestr"       long_name="ice strength"                                                 unit="N/m"      /> 
    71          <field id="icevel"       long_name="ice velocity"                                                 unit="m/s"      /> 
    72          <field id="idive"        long_name="divergence"                                                   unit="1e-8s-1"  /> 
    73          <field id="ishear"       long_name="shear"                                                        unit="1e-8s-1"  /> 
    74          <field id="icevolu"      long_name="ice volume"                                                   unit="m"        /> 
    75          <field id="snowvol"      long_name="snow volume"                                                  unit="m"        /> 
    76          <field id="tau_icebfr"   long_name="ice friction on ocean bottom for landfast ice"                unit="N/2"      /> 
    77  
    78          <field id="icetrp"       long_name="ice volume transport"                                         unit="m/day"          /> 
    79          <field id="snwtrp"       long_name="snw volume transport"                                         unit="m/day"          /> 
    80          <field id="saltrp"       long_name="salt content transport"                                       unit="1e-3*kg/m2/day" /> 
     144    <!-- trends --> 
     145    <field id="afxtot"       long_name="area tendency (total)"                                        unit="s-1"   /> 
     146 
     147    <!-- momentum (advection) --> 
     148         <field id="icetrp"       long_name="ice mass transport"                                           unit="kg/m2/s"          /> 
     149         <field id="snwtrp"       long_name="snw mass transport"                                           unit="kg/m2/s"          /> 
     150         <field id="saltrp"       long_name="salt     transport"                                           unit="1e-3*kg/m2/s" /> 
    81151         <field id="deitrp"       long_name="advected ice enthalpy"                                        unit="W/m2"           /> 
    82152         <field id="destrp"       long_name="advected snw enthalpy"                                        unit="W/m2"           /> 
    83153 
    84          <field id="sfxbri"       long_name="brine salt flux"                                              unit="1e-3*kg/m2/day" /> 
    85          <field id="sfxdyn"       long_name="salt flux from ridging rafting"                               unit="1e-3*kg/m2/day" /> 
    86          <field id="sfxres"       long_name="salt flux from lipupdate (resultant)"                         unit="1e-3*kg/m2/day" /> 
    87          <field id="sfxbog"       long_name="salt flux from bot growth"                                    unit="1e-3*kg/m2/day" /> 
    88          <field id="sfxbom"       long_name="salt flux from bot melt"                                      unit="1e-3*kg/m2/day" /> 
    89          <field id="sfxsum"       long_name="salt flux from surf melt"                                     unit="1e-3*kg/m2/day" /> 
    90          <field id="sfxlam"       long_name="salt flux from lateral melt"                                  unit="1e-3*kg/m2/day" /> 
    91          <field id="sfxsni"       long_name="salt flux from snow-ice formation"                            unit="1e-3*kg/m2/day" /> 
    92          <field id="sfxopw"       long_name="salt flux from open water ice formation"                      unit="1e-3*kg/m2/day" /> 
    93          <field id="sfxsub"       long_name="salt flux from sublimation"                                   unit="1e-3*kg/m2/day" /> 
    94          <field id="sfx"          long_name="salt flux total"                                              unit="1e-3*kg/m2/day" /> 
    95  
    96          <field id="vfxbog"       long_name="daily bottom thermo ice prod."                                unit="m/day"   /> 
    97          <field id="vfxdyn"       long_name="daily  dynamic ice prod."                                     unit="m/day"   /> 
    98          <field id="vfxopw"       long_name="daily lateral thermo ice prod."                               unit="m/day"   /> 
    99          <field id="vfxsni"       long_name="daily snowice ice prod."                                      unit="m/day"   /> 
    100          <field id="vfxsum"       long_name="surface melt"                                                 unit="m/day"   /> 
    101          <field id="vfxlam"       long_name="lateral melt"                                                 unit="m/day"   /> 
    102          <field id="vfxbom"       long_name="bottom melt"                                                  unit="m/day"   /> 
    103          <field id="vfxres"       long_name="daily resultant ice prod./melting from limupdate"             unit="m/day"   /> 
    104          <field id="vfxice"       long_name="ice melt/growth"                                              unit="m/day"   /> 
    105          <field id="vfxsnw"       long_name="snw melt/growth"                                              unit="m/day"   /> 
    106          <field id="vfxsub"       long_name="snw sublimation"                                              unit="m/day"   /> 
    107          <field id="vfxsub_err"   long_name="excess of snw sublimation sent to ocean"                      unit="m/day"   /> 
    108          <field id="vfxspr"       long_name="snw precipitation on ice"                                     unit="m/day"   /> 
    109          <field id="vfxthin"      long_name="daily thermo ice prod. for thin ice(20cm) + open water"      unit="m/day"   /> 
    110  
    111          <field id="afxtot"       long_name="area tendency (total)"                                        unit="day-1"   /> 
    112          <field id="afxdyn"       long_name="area tendency (dynamics)"                                     unit="day-1"   /> 
    113          <field id="afxthd"       long_name="area tendency (thermo)"                                       unit="day-1"   /> 
    114  
     154    <!-- salt fluxes --> 
     155         <field id="sfxbri"       long_name="salt flux from brines"                                        unit="1e-3*kg/m2/s" /> 
     156         <field id="sfxdyn"       long_name="salt flux from ridging rafting"                               unit="1e-3*kg/m2/s" /> 
     157         <field id="sfxres"       long_name="salt flux from lipupdate (resultant)"                         unit="1e-3*kg/m2/s" /> 
     158         <field id="sfxbog"       long_name="salt flux from bot growth"                                    unit="1e-3*kg/m2/s" /> 
     159         <field id="sfxbom"       long_name="salt flux from bot melt"                                      unit="1e-3*kg/m2/s" /> 
     160         <field id="sfxsum"       long_name="salt flux from surf melt"                                     unit="1e-3*kg/m2/s" /> 
     161         <field id="sfxlam"       long_name="salt flux from lateral melt"                                  unit="1e-3*kg/m2/s" /> 
     162         <field id="sfxsni"       long_name="salt flux from snow-ice formation"                            unit="1e-3*kg/m2/s" /> 
     163         <field id="sfxopw"       long_name="salt flux from open water ice formation"                      unit="1e-3*kg/m2/s" /> 
     164         <field id="sfxsub"       long_name="salt flux from sublimation"                                   unit="1e-3*kg/m2/s" /> 
     165         <field id="sfx"          long_name="Salt flux from sea ice"                                       unit="1e-3*kg/m2/s" /> 
     166 
     167    <!-- mass fluxes --> 
     168         <field id="vfxbog"       long_name="bottom thermo ice prod."                                      unit="kg/m2/s"   /> 
     169         <field id="vfxdyn"       long_name="dynamic ice prod."                                            unit="kg/m2/s"   /> 
     170         <field id="vfxopw"       long_name="lateral thermo ice prod."                                     unit="kg/m2/s"   /> 
     171         <field id="vfxsni"       long_name="snowice ice prod."                                            unit="kg/m2/s"   /> 
     172         <field id="vfxsum"       long_name="surface melt"                                                 unit="kg/m2/s"   /> 
     173         <field id="vfxlam"       long_name="lateral melt"                                                 unit="kg/m2/s"   /> 
     174         <field id="vfxbom"       long_name="bottom melt"                                                  unit="kg/m2/s"   /> 
     175         <field id="vfxres"       long_name="resultant ice prod./melting"                                  unit="kg/m2/s"   /> 
     176         <field id="vfxice"       long_name="ice melt/growth"                                              unit="kg/m2/s"   /> 
     177         <field id="vfxsnw"       long_name="snw melt/growth"                                              unit="kg/m2/s"   /> 
     178         <field id="vfxsub"       long_name="snw sublimation"                                              unit="kg/m2/s"   /> 
     179         <field id="vfxsub_err"   long_name="excess of snw sublimation sent to ocean"                      unit="kg/m2/s"   /> 
     180         <field id="vfxspr"       long_name="snw precipitation on ice"                                     unit="kg/m2/s"   /> 
     181         <field id="vfxthin"      long_name="thermo ice prod. for thin ice(20cm) + open water"             unit="kg/m2/s"   /> 
     182         <field id="vfxpnd"       long_name="melt pond water flux to ocean"                                unit="kg/m2/s"   /> 
     183 
     184    <!-- heat fluxes --> 
    115185         <field id="hfxsum"       long_name="heat fluxes causing surface ice melt"                         unit="W/m2"  /> 
    116186         <field id="hfxbom"       long_name="heat fluxes causing bottom ice melt"                          unit="W/m2"  /> 
     
    124194         <field id="hfxin"        long_name="total heat fluxes at the ice/ocean surface"                   unit="W/m2"  /> 
    125195 
    126          <!-- heat flux associated with mass exchange --> 
     196         <!-- heat fluxes associated with mass exchange --> 
    127197         <field id="hfxthd"       long_name="heat fluxes from ice-ocean mass exchange during thermo"       unit="W/m2"  /> 
    128198         <field id="hfxdyn"       long_name="heat fluxes from ice-ocean mass exchange during dynamic"      unit="W/m2"  /> 
     
    144214         <field id="frq_m"    unit="-"    /> 
    145215 
     216    <!-- categories --> 
     217    <field id="iceconc_cat"  long_name="Sea-ice concentration in thickness categories"                unit=""       grid_ref="grid_T_3D_ncatice" /> 
     218         <field id="icethic_cat"  long_name="Sea-ice thickness in thickness categories"                    unit="m"      grid_ref="grid_T_3D_ncatice" /> 
     219         <field id="snowthic_cat" long_name="Snow thickness in thickness categories"                       unit="m"      grid_ref="grid_T_3D_ncatice" /> 
     220         <field id="iceage_cat"   long_name="Ice age for categories"                                       unit="days"   grid_ref="grid_T_3D_ncatice" /> 
     221         <field id="salinity_cat" long_name="Sea-Ice Bulk salinity for categories"                         unit="g/kg"   grid_ref="grid_T_3D_ncatice" /> 
     222         <field id="brinevol_cat" long_name="Brine volume for categories"                                  unit="%"      grid_ref="grid_T_3D_ncatice" /> 
     223         <field id="icetemp_cat"  long_name="Ice temperature for categories"                               unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
     224         <field id="snwtemp_cat"  long_name="Snow temperature for categories"                              unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
     225         <field id="iceamp_cat"   long_name="Ice melt pond concentration for categories"                   unit="%"      grid_ref="grid_T_3D_ncatice" />  
     226         <field id="icevmp_cat"   long_name="Ice melt pond volume for categories"                          unit="m"      grid_ref="grid_T_3D_ncatice" />  
     227         <field id="icehmp_cat"   long_name="Ice melt pond thickness for categories"                       unit="m"      grid_ref="grid_T_3D_ncatice" />  
     228         <field id="iceafp_cat"   long_name="Ice melt pond fraction for categories"                        unit="m"      grid_ref="grid_T_3D_ncatice" />  
     229     
    146230      </field_group> 
    147231 
     
    149233      <field_group id="SBC_scalar"  grid_ref="grid_T_2D" > 
    150234         <!-- available with ln_limdiaout --> 
    151          <field id="ibgfrcvoltop"    long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)"   unit="km3"       /> 
    152          <field id="ibgfrcvolbot"    long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3"       /> 
    153          <field id="ibgfrctemtop"    long_name="global mean heat on top of ice/snw/ocean-atm "                                             unit="1e20J"     /> 
    154          <field id="ibgfrctembot"    long_name="global mean heat below ice (on top of ocean) "                                             unit="1e20J"     /> 
    155          <field id="ibgfrcsal"       long_name="global mean ice/snow forcing (salt equivalent ocean volume)"                               unit="pss*km3"   /> 
    156          <field id="ibgfrchfxtop"    long_name="global mean heat flux on top of ice/snw/ocean-atm "                                        unit="W/m2"      /> 
    157          <field id="ibgfrchfxbot"    long_name="global mean heat flux below ice (on top of ocean) "                                        unit="W/m2"      /> 
     235         <field id="ibgfrcvoltop"    long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)"   unit="km3"        /> 
     236         <field id="ibgfrcvolbot"    long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3"        /> 
     237         <field id="ibgfrctemtop"    long_name="global mean heat on top of ice/snw/ocean-atm "                                             unit="1e20J"      /> 
     238         <field id="ibgfrctembot"    long_name="global mean heat below ice (on top of ocean) "                                             unit="1e20J"      /> 
     239         <field id="ibgfrcsal"       long_name="global mean ice/snow forcing (salt equivalent ocean volume)"                               unit="pss*km3"    /> 
     240         <field id="ibgfrchfxtop"    long_name="global mean heat flux on top of ice/snw/ocean-atm "                                        unit="W/m2"       /> 
     241         <field id="ibgfrchfxbot"    long_name="global mean heat flux below ice (on top of ocean) "                                        unit="W/m2"       /> 
    158242  
    159          <field id="ibgvolume"       long_name="drift in ice/snow volume (equivalent ocean volume)"            unit="km3"        /> 
    160          <field id="ibgsaltco"       long_name="drift in ice salt content (equivalent ocean volume)"           unit="pss*km3"    /> 
    161          <field id="ibgheatco"       long_name="drift in ice/snow heat content"                                unit="1e20J"      /> 
    162          <field id="ibgheatfx"       long_name="drift in ice/snow heat flux"                                   unit="W/m2"       /> 
    163  
    164          <field id="ibgvol_tot"      long_name="global mean ice volume"                                        unit="km3"        /> 
    165          <field id="sbgvol_tot"      long_name="global mean snow volume"                                       unit="km3"        /> 
    166          <field id="ibgarea_tot"     long_name="global mean ice area"                                          unit="km2"        /> 
    167          <field id="ibgsalt_tot"     long_name="global mean ice salt content"                                  unit="1e-3*km3"   /> 
    168          <field id="ibgheat_tot"     long_name="global mean ice heat content"                                  unit="1e20J"      /> 
    169          <field id="sbgheat_tot"     long_name="global mean snow heat content"                                 unit="1e20J"      /> 
     243         <field id="ibgvolume"       long_name="drift in ice/snow volume (equivalent ocean volume)"                                        unit="km3"        /> 
     244         <field id="ibgsaltco"       long_name="drift in ice salt content (equivalent ocean volume)"                                       unit="pss*km3"    /> 
     245         <field id="ibgheatco"       long_name="drift in ice/snow heat content"                                                            unit="1e20J"      /> 
     246         <field id="ibgheatfx"       long_name="drift in ice/snow heat flux"                                                               unit="W/m2"       /> 
     247 
     248         <field id="ibgvol_tot"      long_name="global mean ice volume"                                                                    unit="km3"        /> 
     249         <field id="sbgvol_tot"      long_name="global mean snow volume"                                                                   unit="km3"        /> 
     250         <field id="ibgarea_tot"     long_name="global mean ice area"                                                                      unit="km2"        /> 
     251         <field id="ibgsalt_tot"     long_name="global mean ice salt content"                                                              unit="1e-3*km3"   /> 
     252         <field id="ibgheat_tot"     long_name="global mean ice heat content"                                                              unit="1e20J"      /> 
     253         <field id="sbgheat_tot"     long_name="global mean snow heat content"                                                             unit="1e20J"      /> 
     254 
     255         <field id="NH_iceextt"      long_name="Sea ice extent North"                            standard_name="sea_ice_extent_n"          unit="1e6_km2"    /> 
     256         <field id="SH_iceextt"      long_name="Sea ice extent South"                            standard_name="sea_ice_extent_s"          unit="1e6_km2"    /> 
     257         <field id="NH_icevolu"      long_name="Sea ice volume North"                            standard_name="sea_ice_volume_n"          unit="1e3_km3"    /> 
     258         <field id="SH_icevolu"      long_name="Sea ice volume South"                            standard_name="sea_ice_volume_s"          unit="1e3_km3"    /> 
     259         <field id="NH_icearea"      long_name="Sea ice area North"                              standard_name="sea_ice_area_n"            unit="1e6_km2"    /> 
     260         <field id="SH_icearea"      long_name="Sea ice area South"                              standard_name="sea_ice_area_s"            unit="1e6_km2"    /> 
     261 
     262         <field id="strait_mifl"     long_name="Sea ice mass flux through straits"      standard_name="sea_ice_mass_transport_across_line" unit="kg/s"  grid_ref="grid_4strait"   /> 
     263         <field id="strait_arfl"     long_name="Sea ice area flux through straits"      standard_name="sea_ice_area_transport_across_line" unit="m2/s"  grid_ref="grid_4strait"   />   
     264         <field id="strait_msfl"     long_name="Sea ice snow flux through straits"      standard_name="snow_mass_transport_across_line"    unit="kg/s"  grid_ref="grid_4strait"   /> 
     265 
    170266      </field_group> 
    171267   
     
    176272       
    177273      <field_group id="myvarICE" > 
    178          <field field_ref="icethic_cea"      name="sithic"     long_name="sea_ice_thickness"        /> 
     274         <field field_ref="icethick"         name="sithic"  /> 
    179275         <field field_ref="icevolu"          name="sivolu"  /> 
    180276         <field field_ref="iceconc"          name="siconc"  /> 
    181277      </field_group>     
    182        
     278 
     279      <field_group id="ICE_globalbudget"> 
     280 
     281      <field field_ref="ibgvol_tot"     grid_ref="grid_1point"  name="ibgvol_tot"   /> 
     282   <field field_ref="sbgvol_tot"     grid_ref="grid_1point"  name="sbgvol_tot"   /> 
     283      <field field_ref="ibgarea_tot"    grid_ref="grid_1point"  name="ibgarea_tot"  /> 
     284      <field field_ref="ibgsalt_tot"    grid_ref="grid_1point"  name="ibgsalt_tot"  /> 
     285      <field field_ref="ibgheat_tot"    grid_ref="grid_1point"  name="ibgheat_tot"  /> 
     286      <field field_ref="sbgheat_tot"    grid_ref="grid_1point"  name="sbgheat_tot"  /> 
     287 
     288      <field field_ref="ibgvolume"      grid_ref="grid_1point"  name="ibgvolume"    /> 
     289      <field field_ref="ibgsaltco"      grid_ref="grid_1point"  name="ibgsaltco"    /> 
     290      <field field_ref="ibgheatco"      grid_ref="grid_1point"  name="ibgheatco"    /> 
     291        <field field_ref="ibgheatfx"      grid_ref="grid_1point"  name="ibgheatfx"    /> 
     292 
     293      <field field_ref="ibgfrcvoltop"   grid_ref="grid_1point"  name="ibgfrcvoltop" /> 
     294      <field field_ref="ibgfrcvolbot"   grid_ref="grid_1point"  name="ibgfrcvolbot" /> 
     295      <field field_ref="ibgfrctemtop"   grid_ref="grid_1point"  name="ibgfrctemtop" /> 
     296      <field field_ref="ibgfrctembot"   grid_ref="grid_1point"  name="ibgfrctembot" /> 
     297      <field field_ref="ibgfrcsal"      grid_ref="grid_1point"  name="ibgfrcsal"    /> 
     298        <field field_ref="ibgfrchfxtop"   grid_ref="grid_1point"  name="ibgfrchfxtop" /> 
     299        <field field_ref="ibgfrchfxbot"   grid_ref="grid_1point"  name="ibgfrchfxbot" /> 
     300 
     301      </field_group> 
     302 
     303 
     304      <!--============================--> 
     305      <!-- SIMIP sea ice field groups --> 
     306      <!-- Notz et al 2016            --> 
     307      <!--============================--> 
     308 
     309      <!-- SIMIP monthly scalar variables --> 
     310      <field_group id="SImon_scalars"> 
     311        <!-- Integrated quantities --> 
     312        <field field_ref="NH_iceextt"     grid_ref="grid_1point"  name="siextentn"    /> 
     313        <field field_ref="SH_iceextt"     grid_ref="grid_1point"  name="siextents"    /> 
     314        <field field_ref="NH_icevolu"     grid_ref="grid_1point"  name="sivoln"       /> 
     315        <field field_ref="SH_icevolu"     grid_ref="grid_1point"  name="sivols"       /> 
     316        <field field_ref="NH_icearea"     grid_ref="grid_1point"  name="siarean"      /> 
     317        <field field_ref="SH_icearea"     grid_ref="grid_1point"  name="siareas"      /> 
     318      </field_group> 
     319 
     320      <!-- SIMIP daily fields --> <field_group id="SIday_fields">  
     321          <field field_ref="icepres"          name="sitimefrac"   />  
     322          <field field_ref="iceconc_pct"      name="siconc"       />  
     323          <field field_ref="icethic"          name="sithick"      default_value="1.0e20" detect_missing_value="true" operation="average" />  
     324          <field field_ref="snothic"          name="sisnthick"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     325          <field field_ref="icestK"           name="sitemptop"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     326          <field field_ref="uice_mv"          name="siu"          default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     327          <field field_ref="vice_mv"          name="siv"          default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     328          <field field_ref="icevel_mv"        name="sispeed"      default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     329      </field_group> 
     330 
     331      <!-- SIMIP monthly fields --> 
     332      <field_group id="SImon_fields"> 
     333          <!-- Sea-ice state variables --> 
     334          <field field_ref="icepres"          name="sitimefrac"   /> 
     335          <field field_ref="iceconc_pct"      name="siconc"       /> 
     336          <field field_ref="icemass"          name="simass"       /> 
     337          <field field_ref="icethic"          name="sithick"      default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     338          <field field_ref="icevolu"          name="sivol"        /> 
     339          <field field_ref="snomass"          name="sisnmass"     default_value="1.0e20" detect_missing_value="true" operation="average" />  
     340          <field field_ref="snothic"          name="sisnthick"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     341    
     342          <!-- category-dependent fields --> 
     343          <field field_ref="iceconc_cat_pct_mv"  name="siitdconc"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     344          <field field_ref="icethic_cat_mv"      name="siitdthick"   default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     345          <field field_ref="snowthic_cat_mv"     name="siitdsnthick" default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     346 
     347          <!-- additional state variables--> 
     348          <field field_ref="icestK"           name="sitemptop"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     349          <field field_ref="icesntK"          name="sitempsnic"   default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     350          <field field_ref="icebotK"          name="sitempbot"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     351          <field field_ref="iceage"           name="siage"        default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     352          <field field_ref="icesmass"         name="sisaltmass"   /> 
     353          <field field_ref="icesal"           name="sisali"       default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     354          <field field_ref="icefb"            name="sifb"         default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     355          <field field_ref="icehcneg"         name="sihc"         /> 
     356          <field field_ref="isnhcneg"         name="sisnhc"       default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     357 
     358          <!-- freshwater fluxes --> 
     359          <field field_ref="wfxtot"           name="siflfwbot"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     360          <field field_ref="wfxsum"           name="siflfwdrain"  default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     361 
     362          <!-- area balance --> 
     363          <field field_ref="afxthd"           name="sidconcth"    /> 
     364          <field field_ref="afxdyn"           name="sidconcdyn"   /> 
     365 
     366          <!-- mass balance --> 
     367          <field field_ref="dmithd"           name="sidmassth"        /> 
     368          <field field_ref="dmidyn"           name="sidmassdyn"       /> 
     369          <field field_ref="dmiopw"           name="sidmassgrowthwat" /> 
     370          <field field_ref="dmibog"           name="sidmassgrowthbot" /> 
     371          <field field_ref="dmisni"           name="sidmasssi"        /> 
     372          <field field_ref="dmisub"           name="sidmassevapsubl"  /> 
     373          <field field_ref="dmisum"           name="sidmassmelttop"   /> 
     374          <field field_ref="dmibom"           name="sidmassmeltbot"   /> 
     375          <field field_ref="dmsspr"           name="sndmasssnf"       /> 
     376          <field field_ref="dmsmel"           name="sndmassmelt"      /> 
     377          <field field_ref="dmssub"           name="sndmasssubl"      /> 
     378          <field field_ref="dmsdyn"           name="sndmassdyn"       /> 
     379          <field field_ref="dmsssi"           name="sndmasssi"        /> 
     380 
     381          <!-- heat balance --> 
     382          <field field_ref="hfxsenso"         name="siflsensupbot"    default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     383          <field field_ref="hfxconsu"         name="siflcondtop"      default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     384          <field field_ref="hfxconbo"         name="siflcondbot"      default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     385 
     386          <!-- salt balance --> 
     387          <field field_ref="sfx_mv"           name="siflsaltbot"      default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     388 
     389          <!-- sea-ice dynamics --> 
     390          <field field_ref="uice_mv"          name="siu"              default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     391          <field field_ref="vice_mv"          name="siv"              default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     392          <field field_ref="icevel_mv"        name="sispeed"          default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     393 
     394          <field field_ref="utau_ice"         name="sistrxdtop"       default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     395          <field field_ref="vtau_ice"         name="sistrydtop"       default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     396 
     397          <field field_ref="utau_oi"          name="sistrxubot"       default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     398          <field field_ref="vtau_oi"          name="sistryubot"       default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     399 
     400          <field field_ref="icestr"           name="sicompstren"      default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     401 
     402          <field field_ref="dssh_dx"          name="siforcetiltx"     default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     403          <field field_ref="dssh_dy"          name="siforcetilty"     default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     404 
     405          <field field_ref="corstrx"          name="siforcecoriolx"   default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     406          <field field_ref="corstry"          name="siforcecorioly"   default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     407 
     408          <field field_ref="intstrx"          name="siforceintstrx"   default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     409          <field field_ref="intstry"          name="siforceintstry"   default_value="1.0e20" detect_missing_value="true" operation="average" /> 
     410 
     411          <field field_ref="xmtrptot"         name="sidmasstranx"     /> 
     412          <field field_ref="ymtrptot"         name="sidmasstrany"     /> 
     413 
     414          <field field_ref="normstr"          name="sistresave"       default_value="1.0e20" detect_missing_value="true" operation="instant" /> 
     415          <field field_ref="sheastr"          name="sistremax"        default_value="1.0e20" detect_missing_value="true" operation="instant" /> 
     416 
     417          <field field_ref="idive"            name="sidivvel"         default_value="1.0e20" detect_missing_value="true" operation="instant" /> 
     418          <field field_ref="ishear"           name="sishevel"         default_value="1.0e20" detect_missing_value="true" operation="instant" /> 
     419 
     420          <!-- sea ice fluxes across straits: 2D arrays to be post-processed  --> 
     421          <field field_ref="xmtrpice"         name="simassacrossline_x" /> 
     422          <field field_ref="ymtrpice"         name="simassacrossline_y" /> 
     423          <field field_ref="xmtrpsnw"         name="snmassacrossline_x" /> 
     424          <field field_ref="ymtrpsnw"         name="snmassacrossline_y" /> 
     425          <field field_ref="xatrp"            name="siareaacrossline_x" /> 
     426          <field field_ref="yatrp"            name="siareaacrossline_y" /> 
     427 
     428      </field_group> 
     429 
    183430   </field_definition> 
  • branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED/field_def_nemo-opa.xml

    r8698 r8881  
    299299         <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
    300300 
    301          <!-- u-eddy coefficients (ldftra) --> 
     301         <!-- u-eddy diffusivity coefficients (available if ln_traldf_NONE=F) --> 
    302302         <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
    303303         <field id="ahtu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s or m4/s"      grid_ref="grid_U_3D"/> 
     304         <!-- u-eiv diffusivity coefficients (available if ln_ldfeiv=F) --> 
    304305         <field id="aeiu_2d"      long_name=" surface u-EIV coefficient"                unit="m2/s" /> 
    305306         <field id="aeiu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s"              grid_ref="grid_U_3D"/> 
    306307 
    307          <!-- variables available with MLE --> 
     308         <!-- variables available with MLE (ln_mle=T) --> 
    308309         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" /> 
    309310 
    310          <!-- uoce_eiv: available EIV --> 
     311         <!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 
    311312         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"   standard_name="bolus_sea_water_x_velocity"   unit="m/s"   grid_ref="grid_U_3D" /> 
    312313 
    313          <!-- uoce_eiv: available with key_trabbl --> 
     314         <!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 
    314315         <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"    unit="m/s"  /> 
     316         <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 --> 
    315317         <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"   unit="m3/s" /> 
    316318 
     
    347349         <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
    348350 
    349          <!-- v-eddy coefficients (ldftra, ldfdyn) --> 
     351         <!-- v-eddy diffusivity coefficients (available if ln_traldf_NONE=F) --> 
    350352         <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
    351353         <field id="ahtv_3d"      long_name=" 3D v-eddy diffusivity coefficient"          unit="m2/s or (m4/s)^1/2"           grid_ref="grid_V_3D"/> 
     354         <!-- v-eiv diffusivity coefficients (available if ln_ldfeiv=F) --> 
    352355         <field id="aeiv_2d"      long_name=" surface v-EIV coefficient"                  unit="m2/s" /> 
    353356         <field id="aeiv_3d"      long_name=" 3D v-EIV coefficient"                       unit="m2/s"                         grid_ref="grid_V_3D" /> 
    354357 
    355          <!-- variables available with MLE --> 
     358         <!-- variables available with MLE (ln_mle=T) --> 
    356359         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"   unit="m3/s"   grid_ref="grid_V_3D" /> 
    357360 
    358          <!-- voce_eiv: available with EIV --> 
     361         <!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 
    359362         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"   standard_name="bolus_sea_water_y_velocity"   unit="m/s"   grid_ref="grid_V_3D" /> 
    360363 
    361          <!-- voce_eiv: available with key_trabbl --> 
     364         <!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 
    362365         <field id="voce_bbl"     long_name="BBL ocean current along j-axis"    unit="m/s"  /> 
     366         <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 --> 
    363367         <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"   unit="m3/s" /> 
    364368 
     
    387391         <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                               unit="m3/s" /> 
    388392 
    389          <!-- woce_eiv: available with EIV --> 
     393         <!-- woce_eiv: available with EIV  (ln_ldfeiv=T and ln_ldfeiv_dia=T)  --> 
    390394         <field id="woce_eiv"     long_name="EIV ocean vertical velocity"   standard_name="bolus_upward_sea_water_velocity"   unit="m/s" /> 
    391395 
     
    394398         <field id="avm"          long_name="vertical eddy viscosity"     standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
    395399 
    396          <!-- avs: available with key_zdfddm --> 
     400         <!-- avs: /= avt with ln_zdfddm=T --> 
    397401         <field id="avs"          long_name="salt vertical eddy diffusivity"   standard_name="ocean_vertical_salt_diffusivity"   unit="m2/s" /> 
    398402         <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     
    402406         <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    403407 
    404          <!-- avt_tide: available with key_zdftmx --> 
    405          <field id="av_tide"      long_name="tidal vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_tides"   unit="m2/s" /> 
    406  
    407          <!-- variables available with key_zdftmx_new --> 
    408          <field id="av_ratio"     long_name="S over T diffusivity ratio"            standard_name="salinity_over_temperature_diffusivity_ratio"                     unit="1"    /> 
    409          <field id="av_wave"      long_name="wave-induced vertical diffusivity"     standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves"         unit="m2/s" /> 
    410          <field id="bflx_tmx"     long_name="wave-induced buoyancy flux"            standard_name="buoyancy_flux_due_to_internal_waves"                             unit="W/kg" /> 
    411          <field id="pcmap_tmx"    long_name="power consumed by wave-driven mixing"  standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing"   unit="W/m2"      grid_ref="grid_W_2D" /> 
    412          <field id="emix_tmx"     long_name="power density available for mixing"    standard_name="power_available_for_mixing_from_breaking_internal_waves"         unit="W/kg" /> 
     408         <!-- avt_tide: available with ln_zdfiwm=T --> 
     409         <field id="av_ratio"     long_name="S over T diffusivity ratio"                     standard_name="salinity_over_temperature_diffusivity_ratio"                     unit="1"    /> 
     410         <field id="av_wave"      long_name="internal wave-induced vertical diffusivity"     standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves"         unit="m2/s" /> 
     411         <field id="bflx_iwm"     long_name="internal wave-induced buoyancy flux"            standard_name="buoyancy_flux_due_to_internal_waves"                             unit="W/kg" /> 
     412         <field id="pcmap_iwm"    long_name="power consumed by wave-driven mixing"           standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing"   unit="W/m2"      grid_ref="grid_W_2D" /> 
     413         <field id="emix_iwm"     long_name="power density available for mixing"             standard_name="power_available_for_mixing_from_breaking_internal_waves"         unit="W/kg" /> 
    413414 
    414415         <!-- variables available with diaar5 -->    
     
    704705       <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        />    
    705706 
     707      <!-- variables available when explicit lateral mixing is used (ln_dynldf_NONE=F) --> 
     708      <field id="dispkexyfo"    long_name="KE-trend: lateral  mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction"                   unit="W/m^2" grid_ref="grid_T_2D" />    
     709      <field id="dispkevfo"     long_name="KE-trend: vertical mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction"             unit="W/m^2" grid_ref="grid_T_2D" />    
     710      <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 
     711      <field id="eketrd_eiv"    long_name="EKE-trend due to parameterized eddy advection"   standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection"   unit="W/m^2" grid_ref="grid_T_2D" />    
     712 
    706713       <!-- variables available with ln_PE_trd --> 
    707714       <field id="petrd_xad"     long_name="pe-trend: i-advection"                unit="W/m^3"                        /> 
    708715       <field id="petrd_yad"     long_name="pe-trend: j-advection"                unit="W/m^3"                        /> 
    709716       <field id="petrd_zad"     long_name="pe-trend: k-advection"                unit="W/m^3"                        /> 
    710        <field id="petrd_sad"     long_name="pe-trend: surface adv. (linssh true)"      unit="W/m^3"   grid_ref="grid_T_2D" /> 
     717       <field id="petrd_sad"     long_name="pe-trend: surface adv. (linssh true)" unit="W/m^3"   grid_ref="grid_T_2D" /> 
    711718       <field id="petrd_ldf"     long_name="pe-trend: lateral  diffusion"         unit="W/m^3"                        /> 
    712719       <field id="petrd_zdf"     long_name="pe-trend: vertical diffusion"         unit="W/m^3"                        /> 
  • branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r8599 r8881  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! LIM3 namelist:   
    3 !!              1 - Generic parameters                 (namicerun) 
    4 !!              2 - Diagnostics                        (namicediag) 
    5 !!              3 - Ice initialization                 (namiceini) 
    6 !!              4 - Ice discretization                 (namiceitd) 
    7 !!              5 - Ice dynamics and transport         (namicedyn) 
    8 !!              6 - Ice diffusion                      (namicehdf) 
    9 !!              7 - Ice thermodynamics                 (namicethd) 
    10 !!              8 - Ice salinity                       (namicesal) 
    11 !!              9 - Ice mechanical redistribution      (namiceitdme) 
     2!! ESIM namelist:   
     3!!              1 - Generic parameters                 (nampar) 
     4!!              2 - Ice thickness discretization       (namitd) 
     5!!              3 - Ice dynamics                       (namdyn) 
     6!!              4 - Ice ridging/rafting                (namdyn_rdgrft) 
     7!!              5 - Ice rheology                       (namdyn_rhg) 
     8!!              6 - Ice advection                      (namdyn_adv) 
     9!!              7 - Ice surface forcing                (namforcing) 
     10!!              8 - Ice thermodynamics                 (namthd) 
     11!!              9 - Ice heat diffusion                 (namthd_zdf) 
     12!!             10 - Ice lateral melting                (namthd_da) 
     13!!             11 - Ice growth in open water           (namthd_do) 
     14!!             12 - Ice salinity                       (namthd_sal) 
     15!!             13 - Ice melt ponds                     (namthd_pnd) 
     16!!             14 - Ice initialization                 (namini) 
     17!!             15 - Ice/snow albedos                   (namalb) 
     18!!             16 - Ice diagnostics                    (namdia) 
    1219!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1320! 
    1421!------------------------------------------------------------------------------ 
    15 &namicerun     !   Generic parameters 
    16 !------------------------------------------------------------------------------ 
    17    jpl              =    5          !  number of ice  categories 
    18    nlay_i           =    2          !  number of ice  layers 
    19    nlay_s           =    1          !  number of snow layers (only 1 is working) 
    20    rn_amax_n        =   0.997       !  maximum tolerated ice concentration NH 
    21    rn_amax_s        =   0.997       !  maximum tolerated ice concentration SH 
    22    cn_icerst_in     = "restart_ice" !  suffix of ice restart name (input) 
    23    cn_icerst_out    = "restart_ice" !  suffix of ice restart name (output) 
    24    cn_icerst_indir  = "."           !  directory to read   input ice restarts 
    25    cn_icerst_outdir = "."           !  directory to write output ice restarts 
    26    ln_limthd        =  .true.       !  ice thermo   (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    27    ln_limdyn        =  .true.       !  ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    28    nn_limdyn        =   2           !     (ln_limdyn=T) switch for ice dynamics    
    29                                     !      2: total 
    30                                     !      1: advection only (no diffusion, no ridging/rafting) 
    31                                     !      0: advection only (as 1 but with prescribed velocity, bypass rheology) 
    32    rn_uice          =   0.00001     !     (nn_limdyn=0) ice u-velocity 
    33    rn_vice          =  -0.00001     !     (nn_limdyn=0) ice v-velocity 
    34 / 
    35 !------------------------------------------------------------------------------ 
    36 &namicediag    !   Diagnostics 
    37 !------------------------------------------------------------------------------ 
    38    ln_limdiachk   =  .false.        !  check online the heat, mass & salt budgets (T) or not (F) 
    39    ln_limdiahsb   =  .false.        !  output the heat, mass & salt budgets (T) or not (F) 
    40    ln_limctl      =  .false.        !  ice points output for debug (T or F) 
    41    iiceprt        =    10           !  i-index for debug 
    42    jiceprt        =    10           !  j-index for debug 
    43 / 
    44 !------------------------------------------------------------------------------ 
    45 &namiceini     !   Ice initialization 
    46 !------------------------------------------------------------------------------ 
    47                   ! -- limistate -- ! 
    48    ln_limini      = .true.          !  activate ice initialization (T) or not (F) 
    49    ln_limini_file = .false.         !  netcdf file provided for initialization (T) or not (F) 
    50    rn_thres_sst   =  2.0            !  maximum water temperature with initial ice (degC) 
    51    rn_hts_ini_n   =  0.3            !  initial real snow thickness (m), North 
    52    rn_hts_ini_s   =  0.3            !        "            "             South 
    53    rn_hti_ini_n   =  3.0            !  initial real ice thickness  (m), North 
    54    rn_hti_ini_s   =  1.0            !        "            "             South 
    55    rn_ati_ini_n   =  0.9            !  initial ice concentration   (-), North 
    56    rn_ati_ini_s   =  0.9            !        "            "             South 
    57    rn_smi_ini_n   =  6.3            !  initial ice salinity     (g/kg), North 
    58    rn_smi_ini_s   =  6.3            !        "            "             South 
    59    rn_tmi_ini_n   =  270.           !  initial ice/snw temperature (K), North 
    60    rn_tmi_ini_s   =  270.           !        "            "             South 
    61 / 
    62 !------------------------------------------------------------------------------ 
    63 &namiceitd     !   Ice discretization 
    64 !------------------------------------------------------------------------------ 
    65    nn_catbnd      =    2           !  computation of ice category boundaries based on 
    66                                    !      1: tanh function 
    67                                    !      2: h^(-alpha), function of rn_himean 
    68    rn_himean      =    2.0         !  expected domain-average ice thickness (m), nn_catbnd = 2 only 
    69 / 
    70 !------------------------------------------------------------------------------ 
    71 &namicedyn     !   Ice dynamics and transport 
    72 !------------------------------------------------------------------------------ 
    73                   ! -- limtrp & limadv -- ! 
    74    nn_limadv      =    0            !  choose the advection scheme (-1=Prather ; 0=Ultimate-Macho) 
    75    nn_limadv_ord  =    5            !  choose the order of the advection scheme (if nn_limadv=0) 
    76                   ! -- limitd_me -- ! 
    77    nn_icestr      =    0            !  ice strength parameteriztaion                       
    78                                     !     0: Hibler_79     P = pstar*<h>*exp(-c_rhg*A) 
    79                                     !     1: Rothrock_75   P = Cf*coeff*integral(wr.h^2)     
    80    rn_pe_rdg      =   17.0          !     (nn_icestr=1) ridging work divided by pot. energy change in ridging 
    81    rn_pstar       =    2.0e+04      !     (nn_icestr=0) ice strength thickness parameter (N/m2)  
    82    rn_crhg        =   20.0          !     (nn_icestr=0) ice strength conc. parameter (-) 
    83    ln_icestr_bvf  =    .false.      !     ice strength function brine volume (T) or not (F) 
    84                                     ! 
    85             ! -- limdyn & limrhg -- ! 
    86    rn_cio         =    5.0e-03      !  ice-ocean drag coefficient (-) 
    87    rn_creepl      =    1.0e-12      !  creep limit (s-1) 
    88    rn_ecc         =    2.0          !  eccentricity of the elliptical yield curve           
    89    nn_nevp        =  120            !  number of EVP subcycles                              
    90    rn_relast      =    0.333        !  ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
    91                                     !     advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
    92    ln_landfast    =  .false.        !  landfast ice parameterization (T or F)                            
    93    rn_gamma       =    0.15         !     (ln_landfast=T)  fraction of ocean depth that ice must reach to initiate landfast 
    94                                     !                      recommended range: [0.1 ; 0.25] 
    95    rn_icebfr      =    10.          !     (ln_landfast=T)  maximum bottom stress per unit area of contact (N/m2)                  
    96                                     !                      a very large value ensures ice velocity=0 even with a small contact area 
    97                                     !                      recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2) 
    98    rn_lfrelax     =    1.e-5        !     (ln_landfast=T)  relaxation time scale to reach static friction (s-1)                  
    99 / 
    100 !------------------------------------------------------------------------------ 
    101 &namicehdf     !   Ice horizontal diffusion 
    102 !------------------------------------------------------------------------------ 
    103                      ! -- limhdf -- ! 
    104    nn_ahi0        =    -1           !  horizontal diffusivity computation 
    105                                     !    -1: no diffusion (bypass limhdf) 
    106                                     !     0: use rn_ahi0_ref 
    107                                     !     1: use rn_ahi0_ref x mean grid cell length / ( 2deg mean grid cell length ) 
    108                                     !     2: use rn_ahi0_ref x grid cell length      / ( 2deg mean grid cell length ) 
    109    rn_ahi0_ref    = 350.0           !  horizontal sea ice diffusivity (m2/s)  
    110                                     !     if nn_ahi0 > 0, rn_ahi0_ref is the reference value at a nominal 2 deg resolution 
    111 / 
    112 !------------------------------------------------------------------------------ 
    113 &namicethd     !   Ice thermodynamics 
    114 !------------------------------------------------------------------------------ 
    115                  ! -- limthd_dif -- ! 
    116    rn_kappa_i     = 1.0             !  radiation attenuation coefficient in sea ice (m-1) 
    117    nn_conv_dif    = 50              !  maximal number of iterations for heat diffusion computation 
    118    rn_terr_dif    = 1.0e-04         !  maximum temperature after heat diffusion (degC) 
    119    nn_ice_thcon   = 1               !  sea ice thermal conductivity 
    120                                     !     0: k = k0 + beta.S/T            (Untersteiner, 1964) 
    121                                     !     1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 
    122    ln_it_qnsice   = .true.          !  iterate the surface non-solar flux with surface temperature (T) or not (F) 
    123    nn_monocat     = 0               !  virtual ITD mono-category parameterizations (1, jpl = 1 only) or not (0) 
    124                                     !     2: simple piling instead of ridging    --- temporary option 
    125                                     !     3: activate G(he) only                 --- temporary option 
    126                                     !     4: activate extra lateral melting only --- temporary option 
    127    rn_cdsn     = 0.31              !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 
    128                                    !  Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 
    129                   ! -- limthd_dh -- ! 
    130    ln_limdH       = .true.          !  activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    131    rn_betas       = 0.66            !  exponent in lead-ice repratition of snow precipitation 
    132                                     !     betas = 1 -> equipartition, betas < 1 -> more on leads 
    133                   ! -- limthd_da -- ! 
    134    ln_limdA       = .true.          !  activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    135    rn_beta        = 1.0             !     (ln_latmelt=T) coef. beta for lateral melting param. Recommended range=[0.8-1.2] 
    136                                     !      => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2) 
    137                                     !         0.3 = best fit for western Fram Strait and Antarctica 
    138                                     !         1.4 = best fit for eastern Fram Strait       
    139    rn_dmin        = 8.              !     (ln_latmelt=T) minimum floe diameter for lateral melting param. Recommended range=[6-10] 
    140                                     !      => 6  vs 8m = +40% melting at the peak (A~0.5) 
    141                                     !         10 vs 8m = -20% melting 
    142                  ! -- limthd_lac -- ! 
    143    ln_limdO       = .true.          !  activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    144    rn_hnewice     = 0.1             !  thickness for new ice formation in open water (m) 
    145    ln_frazil      = .false.         !  Frazil ice parameterization (ice collection as a function of wind) 
    146    rn_maxfrazb    = 1.0             !     (ln_frazil=T) maximum fraction of frazil ice collecting at the ice base 
    147    rn_vfrazb      = 0.417           !     (ln_frazil=T) thresold drift speed for frazil ice collecting at the ice bottom (m/s) 
    148    rn_Cfrazb      = 5.0             !     (ln_frazil=T) squeezing coefficient for frazil ice collecting at the ice bottom 
    149                   ! -- limitd_th -- ! 
    150    rn_himin       = 0.1             !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 
    151 / 
    152 !------------------------------------------------------------------------------ 
    153 &namicesal     !   Ice salinity 
    154 !------------------------------------------------------------------------------ 
    155                  ! -- limthd_sal -- ! 
    156    ln_limdS       = .true.          !  activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    157    nn_icesal      =  2              !  ice salinity option 
    158                                     !     1: constant ice salinity (S=rn_icesal) 
    159                                     !     2: varying salinity parameterization S(z,t) 
    160                                     !     3: prescribed salinity profile S(z), Schwarzacher, 1959 
    161    rn_icesal      =  4.             !    (nn_icesal=1) ice salinity (g/kg) 
    162    rn_sal_gd      =  5.             !  restoring ice salinity, gravity drainage (g/kg) 
    163    rn_time_gd     =  1.73e+6        !  restoring time scale, gravity drainage  (s) 
    164    rn_sal_fl      =  2.             !  restoring ice salinity, flushing (g/kg) 
    165    rn_time_fl     =  8.64e+5        !  restoring time scale, flushing (s) 
    166    rn_simax       = 20.             !  maximum tolerated ice salinity (g/kg) 
    167    rn_simin       =  0.1            !  minimum tolerated ice salinity (g/kg) 
    168 / 
    169 !------------------------------------------------------------------------------ 
    170 &namiceitdme   !   Ice mechanical redistribution (ridging and rafting) 
    171 !------------------------------------------------------------------------------ 
    172                   ! -- limitd_me -- ! 
    173    rn_cs          =   0.5           !  fraction of shearing energy contributing to ridging 
    174    nn_partfun     =   1             !  type of ridging participation function 
    175                                     !     0: linear      (Thorndike et al, 1975) 
    176                                     !     1: exponential (Lipscomb, 2007) 
    177    rn_gstar       =   0.15          !     (nn_partfun = 0) fractional area of thin ice being ridged  
    178    rn_astar       =   0.03          !     (nn_partfun = 1) exponential measure of ridging ice fraction [set to 0.05 if hstar=100] 
    179    ln_ridging     =   .true.        !  ridging activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    180    rn_hstar       =  25.0           !     (ln_ridging = T) determines the maximum thickness of ridged ice (m) (Hibler, 1980) 
    181    rn_por_rdg     =   0.3           !     (ln_ridging = T) porosity of newly ridged ice (Lepparanta et al., 1995) 
    182    rn_fsnowrdg    =   0.5           !     (ln_ridging = T) snow volume fraction that survives in ridging 
    183    ln_rafting     =   .true.        !  rafting activated (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    184    rn_hraft       =   0.75          !     (ln_rafting = T) threshold thickness for rafting (m) 
    185    rn_craft       =   5.0           !     (ln_rafting = T) squeezing coefficient used in the rafting function 
    186    rn_fsnowrft    =   0.5           !     (ln_rafting = T) snow volume fraction that survives in rafting 
    187 / 
     22&nampar         !   Generic parameters 
     23!------------------------------------------------------------------------------ 
     24   jpl              =   5             !  number of ice  categories 
     25   nlay_i           =   2             !  number of ice  layers 
     26   nlay_s           =   1             !  number of snow layers (only 1 is working) 
     27   nn_monocat       =   0             !  virtual ITD mono-category parameterizations (1-4 => jpl = 1 only) or not (0) 
     28                                      !     2: simple piling instead of ridging    --- temporary option 
     29                                      !     3: activate G(he) only                 --- temporary option 
     30                                      !     4: activate extra lateral melting only --- temporary option 
     31   ln_icedyn        = .true.          !  ice dynamics (T) or not (F) 
     32   ln_icethd        = .true.          !  ice thermo   (T) or not (F) 
     33   rn_amax_n        =   0.997         !  maximum tolerated ice concentration NH 
     34   rn_amax_s        =   0.997         !  maximum tolerated ice concentration SH 
     35   cn_icerst_in     = "restart_ice"   !  suffix of ice restart name (input) 
     36   cn_icerst_out    = "restart_ice"   !  suffix of ice restart name (output) 
     37   cn_icerst_indir  = "."             !  directory to read   input ice restarts 
     38   cn_icerst_outdir = "."             !  directory to write output ice restarts 
     39/ 
     40!------------------------------------------------------------------------------ 
     41&namitd         !   Ice discretization 
     42!------------------------------------------------------------------------------ 
     43   ln_cat_hfn       = .true.          !  ice categories are defined by a function following rn_himean**(-0.05) 
     44      rn_himean     =   2.0           !  expected domain-average ice thickness (m) 
     45   ln_cat_usr       = .false.         !  ice categories are defined by rn_catbnd below (m) 
     46      rn_catbnd     =   0.,0.45,1.1,2.1,3.7,6.0   
     47   rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
     48/ 
     49!------------------------------------------------------------------------------ 
     50&namdyn         !   Ice dynamics 
     51!------------------------------------------------------------------------------ 
     52   ln_dynFULL       = .true.          !  dyn.: full ice dynamics               (rheology + advection + ridging/rafting + correction) 
     53   ln_dynRHGADV     = .false.         !  dyn.: no ridge/raft & no corrections  (rheology + advection) 
     54   ln_dynADV        = .false.         !  dyn.: only advection w prescribed vel.(rn_uvice + advection) 
     55      rn_uice       =   0.00001       !        prescribed ice u-velocity 
     56      rn_vice       =   0.            !        prescribed ice v-velocity 
     57   rn_ishlat        =   2.            !  free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 
     58   ln_landfast      = .false.         !  landfast ice parameterization (T or F)                            
     59      rn_gamma      =   0.15          !     fraction of ocean depth that ice must reach to initiate landfast 
     60                                      !        recommended range: [0.1 ; 0.25] 
     61      rn_icebfr     =  10.            !     maximum bottom stress per unit area of contact [N/m2]                  
     62                                      !        a very large value ensures ice velocity=0 even with a small contact area 
     63                                      !        recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2) 
     64      rn_lfrelax    =   1.e-5         !     relaxation time scale to reach static friction [s-1] 
     65/ 
     66!------------------------------------------------------------------------------ 
     67&namdyn_rdgrft  !   Ice ridging/rafting 
     68!------------------------------------------------------------------------------ 
     69          ! -- ice_rdgrft_strength -- ! 
     70   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                       
     71      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2] 
     72      rn_crhg       =   20.0          !     ice strength conc. parameter (-) 
     73                   ! -- ice_rdgrft -- ! 
     74   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging 
     75              ! -- ice_rdgrft_prep -- ! 
     76   ln_partf_lin     = .false.         !  Linear ridging participation function (Thorndike et al, 1975) 
     77      rn_gstar      =   0.15          !     fractional area of thin ice being ridged  
     78   ln_partf_exp     = .true.          !  Exponential ridging participation function (Lipscomb, 2007) 
     79      rn_astar      =   0.03          !     exponential measure of ridging ice fraction [set to 0.05 if hstar=100] 
     80   ln_ridging       = .true.          !  ridging activated (T) or not (F) 
     81      rn_hstar      =  25.0           !     determines the maximum thickness of ridged ice [m] (Hibler, 1980) 
     82      rn_porordg    =   0.3           !     porosity of newly ridged ice (Lepparanta et al., 1995) 
     83      rn_fsnwrdg    =   0.5           !     snow volume fraction that survives in ridging 
     84      rn_fpndrdg    =   1.0           !     pond fraction that survives in ridging (small a priori) 
     85   ln_rafting       = .true.          !  rafting activated (T) or not (F) 
     86      rn_hraft      =   0.75          !     threshold thickness for rafting [m] 
     87      rn_craft      =   5.0           !     squeezing coefficient used in the rafting function 
     88      rn_fsnwrft    =   0.5           !     snow volume fraction that survives in rafting 
     89      rn_fpndrft    =   1.0           !     pond fraction that survives in rafting (0.5 a priori) 
     90/ 
     91!------------------------------------------------------------------------------ 
     92&namdyn_rhg     !   Ice rheology 
     93!------------------------------------------------------------------------------ 
     94   ln_rhg_EVP       = .true.          !  EVP rheology 
     95      ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
     96      rn_creepl     =   2.0e-9        !     creep limit [1/s] 
     97      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve           
     98      nn_nevp       = 120             !     number of EVP subcycles                              
     99      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
     100                                      !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     101/ 
     102!------------------------------------------------------------------------------ 
     103&namdyn_adv     !   Ice advection 
     104!------------------------------------------------------------------------------ 
     105   ln_adv_Pra       = .false.         !  Advection scheme (Prather) 
     106   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho) 
     107      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
     108/ 
     109!------------------------------------------------------------------------------ 
     110&namforcing     !   Ice surface forcing 
     111!------------------------------------------------------------------------------ 
     112   rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-) 
     113   rn_blow_s        =   0.66          !  mesure of snow blowing into the leads 
     114                                      !     = 1 => no snow blowing, < 1 => some snow blowing 
     115   nn_iceflx        =  -1             !  Redistribute heat flux over ice categories 
     116                                      !     =-1  Do nothing (needs N(cat) fluxes) 
     117                                      !          ==> forced mode only 
     118                                      !     = 0  Average N(cat) fluxes then apply the average over the N(cat) ice 
     119                                      !          ==> forced and coupled modes 
     120                                      !     = 1  Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity 
     121                                      !          ==> forced mode only 
     122                                      !     = 2  Redistribute a single flux over categories 
     123                                      !          ==> coupled mode only 
     124   nice_jules       = 0               !  Jules coupling (0=OFF, 1=EMULATED, 2=ACTIVE) 
     125/ 
     126!------------------------------------------------------------------------------ 
     127&namthd         !   Ice thermodynamics 
     128!------------------------------------------------------------------------------ 
     129   ln_icedH         = .true.          !  activate ice thickness change from growing/melting (T) or not (F) 
     130   ln_icedA         = .true.          !  activate lateral melting param. (T) or not (F) 
     131   ln_icedO         = .true.          !  activate ice growth in open-water (T) or not (F) 
     132   ln_icedS         = .true.          !  activate brine drainage (T) or not (F) 
     133/ 
     134!------------------------------------------------------------------------------ 
     135&namthd_zdf     !   Ice heat diffusion 
     136!------------------------------------------------------------------------------ 
     137   ln_zdf_BL99      = .true.          !  Heat diffusion follows Bitz and Lipscomb 1999 
     138   ln_cndi_U64      = .false.         !  sea ice thermal conductivity: k = k0 + beta.S/T            (Untersteiner, 1964) 
     139   ln_cndi_P07      = .true.          !  sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 
     140   rn_cnd_s         =   0.31          !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 
     141                                      !     Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 
     142   rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice [1/m] 
     143/ 
     144!------------------------------------------------------------------------------ 
     145&namthd_da     !   Ice lateral melting 
     146!------------------------------------------------------------------------------ 
     147   rn_beta          =   1.0           !  coef. beta for lateral melting param. Recommended range=[0.8-1.2] 
     148                                      !   => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2) 
     149                                      !      0.3 = best fit for western Fram Strait and Antarctica 
     150                                      !      1.4 = best fit for eastern Fram Strait       
     151   rn_dmin          =   8.            !  minimum floe diameter for lateral melting param. Recommended range=[6-10] 
     152                                      !   => 6  vs 8m = +40% melting at the peak (A~0.5) 
     153                                      !      10 vs 8m = -20% melting 
     154/ 
     155!------------------------------------------------------------------------------ 
     156&namthd_do      !   Ice growth in open water 
     157!------------------------------------------------------------------------------ 
     158   rn_hinew         =   0.1           !  thickness for new ice formation in open water (m), must be larger than rn_hnewice 
     159   ln_frazil        = .false.         !  Frazil ice parameterization (ice collection as a function of wind) 
     160      rn_maxfraz    =   1.0           !     maximum fraction of frazil ice collecting at the ice base 
     161      rn_vfraz      =   0.417         !     thresold drift speed for frazil ice collecting at the ice bottom (m/s) 
     162      rn_Cfraz      =   5.0           !     squeezing coefficient for frazil ice collecting at the ice bottom 
     163/ 
     164!------------------------------------------------------------------------------ 
     165&namthd_sal     !   Ice salinity 
     166!------------------------------------------------------------------------------ 
     167   nn_icesal        =   2             !  ice salinity option 
     168                                      !     1: constant ice salinity (S=rn_icesal) 
     169                                      !     2: varying salinity parameterization S(z,t) 
     170                                      !     3: prescribed salinity profile S(z), Schwarzacher, 1959 
     171   rn_icesal        =   4.            !    (nn_icesal=1) ice salinity (g/kg) 
     172   rn_sal_gd        =   5.            !  restoring ice salinity, gravity drainage (g/kg) 
     173   rn_time_gd       =   1.73e+6       !  restoring time scale, gravity drainage  (s) 
     174   rn_sal_fl        =   2.            !  restoring ice salinity, flushing (g/kg) 
     175   rn_time_fl       =   8.64e+5       !  restoring time scale, flushing (s) 
     176   rn_simax         =  20.            !  maximum tolerated ice salinity (g/kg) 
     177   rn_simin         =   0.1           !  minimum tolerated ice salinity (g/kg) 
     178/ 
     179!------------------------------------------------------------------------------ 
     180&namthd_pnd     !   Melt ponds 
     181!------------------------------------------------------------------------------ 
     182   ln_pnd_H12       = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
     183      ln_pnd_fwb    = .false.         !     melt ponds store freshwater or not 
     184   ln_pnd_CST       = .false.         !  activate constant melt ponds 
     185      rn_apnd       =   0.2           !     prescribed pond fraction, at Tsu=0 
     186      rn_hpnd       =   0.05          !     prescribed pond depth, at Tsu=0 
     187   ln_pnd_alb       = .false.         !  melt ponds affect albedo or not 
     188/ 
     189!------------------------------------------------------------------------------ 
     190&namini         !   Ice initialization 
     191!------------------------------------------------------------------------------ 
     192   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
     193   ln_iceini_file   = .false.         !  netcdf file provided for initialization (T) or not (F) 
     194   rn_thres_sst     =   2.0           !  max delta temp. above Tfreeze with initial ice = (sst - tfreeze) 
     195   rn_hts_ini_n     =   0.3           !  initial real snow thickness (m), North 
     196   rn_hts_ini_s     =   0.3           !        "            "             South 
     197   rn_hti_ini_n     =   3.0           !  initial real ice thickness  (m), North 
     198   rn_hti_ini_s     =   1.0           !        "            "             South 
     199   rn_ati_ini_n     =   0.9           !  initial ice concentration   (-), North 
     200   rn_ati_ini_s     =   0.9           !        "            "             South 
     201   rn_smi_ini_n     =   6.3           !  initial ice salinity     (g/kg), North 
     202   rn_smi_ini_s     =   6.3           !        "            "             South 
     203   rn_tmi_ini_n     = 270.            !  initial ice/snw temperature (K), North 
     204   rn_tmi_ini_s     = 270.            !        "            "             South 
     205 
     206   sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     207   sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     208   sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     209   sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     210   sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     211   sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     212   cn_dir='./' 
     213/ 
     214!------------------------------------------------------------------------------ 
     215&namalb         !   albedo parameters 
     216!------------------------------------------------------------------------------ 
     217   !                                  !                          !  obs range (cloud-sky) 
     218   rn_alb_sdry      =   0.85          !  dry snow albedo         :  0.85 -- 0.87 
     219   rn_alb_smlt      =   0.75          !  melting snow albedo     :  0.72 -- 0.82 
     220   rn_alb_idry      =   0.60          !  dry ice albedo          :  0.54 -- 0.65 
     221   rn_alb_imlt      =   0.50          !  bare puddled ice albedo :  0.49 -- 0.58 
     222   rn_alb_dpnd      =   0.27          !  ponded ice albedo       :  0.10 -- 0.30  
     223/ 
     224!------------------------------------------------------------------------------ 
     225&namdia         !   Diagnostics 
     226!------------------------------------------------------------------------------ 
     227   ln_icediachk     = .false.         !  check online the heat, mass & salt budgets (T) or not (F) 
     228   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F) 
     229   ln_icectl        = .false.         !  ice points output for debug (T or F) 
     230   iiceprt          =  10             !  i-index for debug 
     231   jiceprt          =  10             !  j-index for debug 
     232/ 
  • branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED/namelist_ref

    r8599 r8881  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !!                            namelist_ref 
     2!!                            namelist_ref                            !! 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!! NEMO/OPA  :  1 - run manager      (namrun) 
    5 !! namelists    2 - Domain           (namcfg, namzgr, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
    6 !!              3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_sas) 
    7 !!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    8 !!                                    namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave) 
     5!! namelists    2 - Domain           (namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
     6!!              3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 
     7!!                                    namsbc_sas, namtra_qsr, namsbc_rnf, 
     8!!                                    namberg, namsbc_isf, namsbc_iscpl,  
     9!!                                    namsbc_apr, namsbc_ssr, namsbc_wave) 
    910!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    10 !!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
     11!!              5 - bottom  boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl) 
    1112!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    1213!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    13 !!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx, namzdf_tmx_new) 
    14 !!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    15 !!             10 - miscellaneous    (nammpp, namctl) 
     14!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm) 
     15!!              9 - miscellaneous    (nammpp, namctl) 
     16!!             10 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    1617!!             11 - Obs & Assim      (namobs, nam_asminc) 
    1718!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1819 
    1920!!====================================================================== 
    20 !!                   ***  Run management namelists  *** 
     21!!                   ***  Run management namelists  ***               !! 
    2122!!====================================================================== 
    2223!!   namrun       parameters of the run 
     
    6263!!   namwad       Wetting and drying                                    (default F) 
    6364!!   namtsd       data: temperature & salinity 
    64 !!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs") 
     65!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T) 
    6566!!   namc1d       1D configuration options                              ("key_c1d") 
    6667!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
     
    6970! 
    7071!----------------------------------------------------------------------- 
    71 &namcfg        !   parameters of the configuration 
     72&namcfg        !   parameters of the configuration                     !   (default: user defined GYRE) 
    7273!----------------------------------------------------------------------- 
    7374   ln_read_cfg = .false.   !  (=T) read the domain configuration file 
     
    8283/ 
    8384!----------------------------------------------------------------------- 
    84 &namdom        !   space and time domain (bathymetry, mesh, timestep) 
     85&namdom        !   time and space domain 
    8586!----------------------------------------------------------------------- 
    8687   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
     
    9091   rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice 
    9192   ! 
    92    rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
     93   rn_rdt      = 5760.     !  time step for the dynamics and tracer 
    9394   rn_atfp     =    0.1    !  asselin time filter parameter 
    9495   ! 
    95    ln_crs      = .false.   !  Logical switch for coarsening module 
     96   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    9697/ 
    9798!----------------------------------------------------------------------- 
     
    117118/ 
    118119!----------------------------------------------------------------------- 
    119 &namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs") 
     120&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    120121!----------------------------------------------------------------------- 
    121122   nn_factx    = 3         !  Reduction factor of x-direction 
     
    170171!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    171172!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
    172 !!   namsbc_alb      albedo parameters 
    173173!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    174174!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    183183   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    184184   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    185    ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
     185   ln_blk      = .false.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    186186                     ! Type of coupling (Ocean/Ice/Atmosphere) : 
    187187   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
    188188   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 ) 
    189189   nn_components = 0       !  configuration of the opa-sas OASIS coupling 
    190                            !  =0 no opa-sas OASIS coupling: default single executable configuration 
    191                            !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component 
    192                            !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component  
    193    nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
    194                            !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
    195                            !  = 0  Average per-category fluxes (forced and coupled mode) 
    196                            !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
    197                            !  = 2  Redistribute a single flux over categories (coupled mode only) 
     190                           !  =0 no opa-sas OASIS coupling: default single executable config. 
     191                           !  =1 opa-sas OASIS coupling: multi executable config., OPA component 
     192                           !  =2 opa-sas OASIS coupling: multi executable config., SAS component  
    198193                     ! Sea-ice : 
    199    nn_ice      = 3         !  =0 no ice boundary condition   , 
     194   nn_ice      = 2         !  =0 no ice boundary condition   , 
    200195                           !  =1 use observed ice-cover      , 
    201                            !  =2 to 4 :  ice-model used (LIM2, LIM3 or CICE)                         ("key_lim3", "key_lim2", or "key_cice") 
    202    nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    203                            !  =1 levitating ice with mass and salt exchange but no presure effect 
    204                            !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
     196                           !  =2 or 3 automatically for LIM3 or CICE    ("key_lim3" or "key_cice") 
     197                           !          except in AGRIF zoom where it has to be specified 
     198   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
     199                           !  =F levitating ice (no pressure, mass and salt exchanges) 
    205200                     ! Misc. options of sbc :  
    206201   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
     
    231226   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    232227   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    233  
     228   ! 
    234229   cn_dir      = './'      !  root directory for the location of the flux files 
    235230/ 
    236231!----------------------------------------------------------------------- 
    237 &namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk = T) 
     232&namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk =T) 
    238233!----------------------------------------------------------------------- 
    239234!              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                              ! rotation ! land/sea mask ! 
     
    261256   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    262257   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    263    rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity 
    264                            !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    265    ln_Cd_L12   = .false.   !  Modify the drag ice-atm and oce-atm depending on ice concentration 
    266                            !  This parameterization is from Lupkes et al. (JGR 2012) 
     258   rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
     259   !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
     260   ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
     261   ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    267262/ 
    268263!----------------------------------------------------------------------- 
     
    307302/ 
    308303!----------------------------------------------------------------------- 
    309 &namsbc_sas    !   Stand Alone Surface boundary condition 
     304&namsbc_sas    !   Stand-Alone Surface boundary condition 
    310305!----------------------------------------------------------------------- 
    311306!              !  file name  ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    312307!              !             !  (if <0  months)  !   name    !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    313    l_sasread   = .TRUE.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
    314    sn_usp      = 'sas_grid_U',     120           , 'vozocrtx',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    315    sn_vsp      = 'sas_grid_V',     120           , 'vomecrty',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     308   sn_usp      = 'sas_grid_U',     120           , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     309   sn_vsp      = 'sas_grid_V',     120           , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    316310   sn_tem      = 'sas_grid_T',     120           , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    317311   sn_sal      = 'sas_grid_T',     120           , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     
    320314   sn_frq      = 'sas_grid_T',     120           , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    321315 
    322    ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field 
     316   l_sasread   = .true.    !  =T Read the above fields in a file, =F initialize to 0. in sbcssm.F90 
     317   ln_3d_uve   = .false.   !  specify whether we are supplying a 3D u,v and e3 field 
    323318   ln_read_frq = .false.   !  specify whether we must read frq or not 
    324319   cn_dir      = './'      !  root directory for the location of the bulk files are 
    325320/ 
    326321!----------------------------------------------------------------------- 
    327 &namtra_qsr    !   penetrative solar radiation                          (ln_traqsr=T) 
     322&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T) 
    328323!----------------------------------------------------------------------- 
    329324!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    339334   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction 
    340335   rn_si1      =   23.0    !  2 bands: longest depth of extinction 
    341    ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3 
    342 / 
    343 !----------------------------------------------------------------------- 
    344 &namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf=T) 
     336/ 
     337!----------------------------------------------------------------------- 
     338&namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf =T) 
    345339!----------------------------------------------------------------------- 
    346340!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    418412/ 
    419413!----------------------------------------------------------------------- 
    420 &namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr=T) 
     414&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
    421415!----------------------------------------------------------------------- 
    422416!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    428422   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0) 
    429423   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2) 
    430                            !  or to SSS only (=1) or no damping term (=0) 
     424   !                       !  or to SSS only (=1) or no damping term (=0) 
    431425   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K] 
    432426   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
    433427   ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    434428   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    435 / 
    436 !----------------------------------------------------------------------- 
    437 &namsbc_alb    !   albedo parameters 
    438 !----------------------------------------------------------------------- 
    439    nn_ice_alb   =    1   !  parameterization of ice/snow albedo 
    440                          !     0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo 
    441                          !     1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004), 
    442                          !        giving cloud-sky albedo 
    443    rn_alb_sdry  =  0.85  !  dry snow albedo         : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky) 
    444    rn_alb_smlt  =  0.75  !  melting snow albedo     : 0.65 ( '' )          ; 0.75 ( '' )          ; obs 0.72-0.82 ( '' ) 
    445    rn_alb_idry  =  0.60  !  dry ice albedo          : 0.72 ( '' )          ; 0.60 ( '' )          ; obs 0.54-0.65 ( '' ) 
    446    rn_alb_imlt  =  0.50  !  bare puddled ice albedo : 0.53 ( '' )          ; 0.50 ( '' )          ; obs 0.49-0.58 ( '' ) 
    447429/ 
    448430!----------------------------------------------------------------------- 
     
    497479 
    498480!!====================================================================== 
    499 !!               ***  Lateral boundary condition  *** 
     481!!               ***  Lateral boundary condition  ***                 !! 
    500482!!====================================================================== 
    501483!!   namlbc        lateral momentum boundary condition 
     
    580562   bn_tem      = 'amm12_bdyT_tra',         24        , 'votemper',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    581563   bn_sal      = 'amm12_bdyT_tra',         24        , 'vosaline',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    582 ! for lim2 
    583 !   bn_frld    = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    584 !   bn_hicif   = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    585 !   bn_hsnif   = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    586564! for lim3 
    587565!   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    588 !   bn_ht_i    = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    589 !   bn_ht_s    = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
     566!   bn_h_i     = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
     567!   bn_h_s     = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    590568 
    591569   cn_dir      = 'bdydta/' !  root directory for the location of the bulk files 
     
    601579 
    602580!!====================================================================== 
    603 !!                 ***  Bottom boundary condition  *** 
    604 !!====================================================================== 
    605 !!   nambfr        bottom friction 
    606 !!   nambbc        bottom temperature boundary condition 
    607 !!   nambbl        bottom boundary layer scheme                         ("key_trabbl") 
    608 !!====================================================================== 
    609 ! 
    610 !----------------------------------------------------------------------- 
    611 &nambfr        !   bottom friction                                      (default: linear) 
    612 !----------------------------------------------------------------------- 
    613    nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    614                            !                              = 2 : nonlinear friction 
    615    rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case) 
    616    rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
    617    rn_bfri2_max=    1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
    618    rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    619    rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T 
    620    ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    621    rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
    622    rn_tfri1    =    4.e-4  !  top drag coefficient (linear case) 
    623    rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
    624    rn_tfri2_max=    1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
    625    rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2) 
    626    rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T 
    627    ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file ) 
    628    rn_tfrien   =   50.     !  local multiplying factor of tfr (ln_tfr2d=T) 
    629  
    630    ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
    631    ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
     581!!                ***  Top/Bottom boundary condition  ***             !! 
     582!!====================================================================== 
     583!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
     584!!   namdrg_top    top    friction                                      (ln_isfcav=T) 
     585!!   namdrg_bot    bottom friction                                       
     586!!   nambbc        bottom temperature boundary condition                (default: NO) 
     587!!   nambbl        bottom boundary layer scheme                         (default: NO) 
     588!!====================================================================== 
     589! 
     590!----------------------------------------------------------------------- 
     591&namdrg            !   top/bottom drag coefficient                      (default: NO selection) 
     592!----------------------------------------------------------------------- 
     593   ln_NONE    = .false.    !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
     594   ln_lin     = .false.    !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top) 
     595   ln_non_lin = .false.    !  non-linear  drag: Cd = Cd0 |U| 
     596   ln_loglayer= .false.    !  logarithmic drag: Cd = vkarmn/log(z/z0) |U| 
     597   ! 
     598   ln_drgimp  = .true.     !  implicit top/bottom friction flag 
     599/ 
     600!----------------------------------------------------------------------- 
     601&namdrg_top        !   TOP friction                                     (ln_isfcav=T) 
     602!----------------------------------------------------------------------- 
     603   rn_Cd0     =  1.e-3     !  drag coefficient [-] 
     604   rn_Uc0     =  0.4       !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     605   rn_Cdmax   =  0.1       !  drag value maximum [-] (logarithmic drag) 
     606   rn_ke0     =  2.5e-3    !  background kinetic energy  [m2/s2] (non-linear cases) 
     607   rn_z0      =  3.0e-3    !  roughness [m] (ln_loglayer=T) 
     608   ln_boost   = .false.    !  =T regional boost of Cd0 ; =F constant 
     609      rn_boost=  50.          !  local boost factor  [-] 
     610/ 
     611!----------------------------------------------------------------------- 
     612&namdrg_bot        !   BOTTOM friction                                   
     613!----------------------------------------------------------------------- 
     614   rn_Cd0     =  1.e-3    !  drag coefficient [-] 
     615   rn_Uc0     =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     616   rn_Cdmax   =  0.1      !  drag value maximum [-] (logarithmic drag) 
     617   rn_ke0     =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     618   rn_z0      =  3.e-3    !  roughness [m] (ln_loglayer=T) 
     619   ln_boost   = .false.   !  =T regional boost of Cd0 ; =F constant 
     620      rn_boost=  50.         !  local boost factor  [-] 
    632621/ 
    633622!----------------------------------------------------------------------- 
     
    646635/ 
    647636!----------------------------------------------------------------------- 
    648 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    649 !----------------------------------------------------------------------- 
     637&nambbl        !   bottom boundary layer scheme                         (default: NO) 
     638!----------------------------------------------------------------------- 
     639   ln_trabbl   = .false.   !  Bottom Boundary Layer parameterisation flag 
    650640   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    651641   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     
    655645 
    656646!!====================================================================== 
    657 !!                        Tracer (T & S ) namelists 
     647!!                        Tracer (T & S) namelists 
    658648!!====================================================================== 
    659649!!   nameos           equation of state 
     
    666656! 
    667657!----------------------------------------------------------------------- 
    668 &nameos        !   ocean physical parameters 
    669 !----------------------------------------------------------------------- 
    670    ln_teos10   = .false.         !  = Use TEOS-10 equation of state 
    671    ln_eos80    = .false.         !  = Use EOS80 equation of state 
    672    ln_seos     = .false.         !  = Use simplified equation of state (S-EOS) 
     658&nameos        !   ocean Equation Of Seawater                           (default: NO) 
     659!----------------------------------------------------------------------- 
     660   ln_teos10   = .false.         !  = Use TEOS-10 
     661   ln_eos80    = .false.         !  = Use EOS80 
     662   ln_seos     = .false.         !  = Use S-EOS (simplified Eq.) 
    673663                                 ! 
    674664   !                     ! S-EOS coefficients (ln_seos=T): 
     
    683673/ 
    684674!----------------------------------------------------------------------- 
    685 &namtra_adv    !   advection scheme for tracer                          (default: NO advection) 
    686 !----------------------------------------------------------------------- 
     675&namtra_adv    !   advection scheme for tracer                          (default: NO selection) 
     676!----------------------------------------------------------------------- 
     677   ln_traadv_NONE= .false. !  No tracer advection 
    687678   ln_traadv_cen = .false. !  2nd order centered scheme 
    688679      nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     
    691682      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order  
    692683      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  
    693       nn_fct_zts =  0            !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
    694       !                          !        (number of sub-timestep = nn_fct_zts) 
    695684   ln_traadv_mus = .false. !  MUSCL scheme 
    696685      ln_mus_ups = .false.       !  use upstream scheme near river mouths 
     
    713702/ 
    714703!----------------------------------------------------------------------- 
    715 &namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO diffusion) 
     704&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection) 
    716705!----------------------------------------------------------------------- 
    717706   !                       !  Operator type: 
    718    !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     707   ln_traldf_NONE  =  .false.  !  No explicit diffusion 
    719708   ln_traldf_lap   =  .false.  !    laplacian operator 
    720709   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     
    748737&namtra_ldfeiv !   eddy induced velocity param.                         (default: NO) 
    749738!----------------------------------------------------------------------- 
    750    ln_ldfeiv     =.false. ! use eddy induced velocity parameterization 
    751    ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities 
    752    rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    753    nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
    754    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    755    !                                !   =  0           constant  
    756    !                                !   = 10 F(k)      =ldf_c1d  
    757    !                                !   = 20 F(i,j)    =ldf_c2d  
    758    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    759    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     739   ln_ldfeiv     = .false. ! use eddy induced velocity parameterization 
     740      rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     741      nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     742      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     743      !                             !   =  0           constant  
     744      !                             !   = 10 F(k)      =ldf_c1d  
     745      !                             !   = 20 F(i,j)    =ldf_c2d  
     746      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     747      !                             !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     748      ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities 
    760749/ 
    761750!----------------------------------------------------------------------- 
    762751&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    763752!----------------------------------------------------------------------- 
    764    ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    765    nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column 
    766                            !                   =1 no damping in the mixing layer (kz  criteria) 
    767                            !                   =2 no damping in the mixed  layer (rho crieria) 
    768    cn_resto    ='resto.nc' !  Name of file containing restoration coeff. field (use dmp_tools to create this) 
    769 / 
    770  
    771 !!====================================================================== 
    772 !!                      ***  Dynamics namelists  *** 
    773 !!====================================================================== 
     753   ln_tradmp   =  .true.   !  add a damping term 
     754      nn_zdmp     =    0      !  vertical shape =0    damping throughout the water column 
     755      !                       !                 =1 no damping in the mixing layer (kz  criteria) 
     756      !                       !                 =2 no damping in the mixed  layer (rho crieria) 
     757      cn_resto    ='resto.nc' !  Name of file containing restoration coeff. field (use dmp_tools to create this) 
     758/ 
     759 
     760!!====================================================================== 
     761!!                      ***  Dynamics namelists  ***                  !! 
     762!!====================================================================== 
     763!!   nam_vvl       vertical coordinate options 
    774764!!   namdyn_adv    formulation of the momentum advection 
    775765!!   namdyn_vor    advection scheme 
     
    780770! 
    781771!----------------------------------------------------------------------- 
    782 &namdyn_adv    !   formulation of the momentum advection                (default: vector form) 
    783 !----------------------------------------------------------------------- 
    784    ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    785    nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction 
    786    ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    787    ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    788    ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection 
    789 / 
    790 !----------------------------------------------------------------------- 
    791 &nam_vvl    !   vertical coordinate options                             (default: zstar) 
    792 !----------------------------------------------------------------------- 
    793    ln_vvl_zstar  = .true.           !  zstar vertical coordinate 
    794    ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
     772&nam_vvl       !   vertical coordinate options                             (default: z-star) 
     773!----------------------------------------------------------------------- 
     774   ln_vvl_zstar  = .true.           !  z-star vertical coordinate 
     775   ln_vvl_ztilde = .false.          !  z-tilde vertical coordinate: only high frequency variations 
    795776   ln_vvl_layer  = .false.          !  full layer vertical coordinate 
    796777   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar 
    797778   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator 
    798    rn_ahe3       = 0.0e0            !  thickness diffusion coefficient 
    799    rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days] 
    800    rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days] 
    801    rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation 
     779   rn_ahe3       =  0.0             !  thickness diffusion coefficient 
     780   rn_rst_e3t    = 30.0             !  ztilde to zstar restoration timescale [days] 
     781   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days] 
     782   rn_zdef_max   =  0.9             !  maximum fractional e3t deformation 
    802783   ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
     784/ 
     785!----------------------------------------------------------------------- 
     786&namdyn_adv    !   formulation of the momentum advection                (default: NO selection) 
     787!----------------------------------------------------------------------- 
     788   ln_dynadv_NONE= .false. !  linear dynamics (no momentum advection) 
     789   ln_dynadv_vec = .false. !  vector form - 2nd centered scheme 
     790     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction 
     791   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
     792   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    803793/ 
    804794!----------------------------------------------------------------------- 
     
    813803/ 
    814804!----------------------------------------------------------------------- 
    815 &namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: zps) 
     805&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: NO selection) 
    816806!----------------------------------------------------------------------- 
    817807   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
     
    837827/ 
    838828!----------------------------------------------------------------------- 
    839 &namdyn_ldf    !   lateral diffusion on momentum                        (default: NO) 
     829&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO selection) 
    840830!----------------------------------------------------------------------- 
    841831   !                       !  Type of the operator : 
    842    !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     832   ln_dynldf_NONE=  .false.    !  No operator (i.e. no explicit diffusion) 
    843833   ln_dynldf_lap =  .false.    !    laplacian operator 
    844834   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     
    866856   rn_maxfac     = 1.0         !  multiplier of theorectical upper limit 
    867857/ 
    868  
    869 !!====================================================================== 
    870 !!             Tracers & Dynamics vertical physics namelists 
     858! 
     859!!====================================================================== 
     860!!                     vertical physics namelists                     !! 
    871861!!====================================================================== 
    872862!!    namzdf        vertical physics 
    873 !!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric") 
    874 !!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke") 
    875 !!    namzdf_gls    GLS vertical mixing                                 ("key_zdfgls") 
    876 !!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm") 
    877 !!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx") 
    878 !!====================================================================== 
    879 ! 
    880 !----------------------------------------------------------------------- 
    881 &namzdf        !   vertical physics 
    882 !----------------------------------------------------------------------- 
    883    rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    884    rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    885    nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0) 
    886    nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0) 
    887    ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F) 
    888       nn_evdm     =    0        ! evd apply on tracer (=0) or on tracer and momentum (=1) 
    889       rn_avevd    =  100.       !  evd mixing coefficient [m2/s] 
    890    ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F) 
    891       nn_npc      =    1        ! frequency of application of npc 
    892       nn_npcp     =  365        ! npc control print frequency 
    893    ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping 
    894       nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T 
    895    ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 
    896 / 
    897 !----------------------------------------------------------------------- 
    898 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
     863!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T) 
     864!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T) 
     865!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T) 
     866!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T) 
     867!!====================================================================== 
     868! 
     869!----------------------------------------------------------------------- 
     870&namzdf        !   vertical physics                                     (default: NO selection) 
     871!----------------------------------------------------------------------- 
     872   !                       ! type of vertical closure (required) 
     873   ln_zdfcst   = .false.      !  constant mixing 
     874   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     875   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     876   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     877   ! 
     878   !                       ! convection 
     879   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     880      nn_evdm     =    0         ! apply on tracer (=0) or on tracer and momentum (=1) 
     881      rn_evd      =  100.        ! mixing coefficient [m2/s] 
     882   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     883      nn_npc      =    1         ! frequency of application of npc 
     884      nn_npcp     =  365         ! npc control print frequency 
     885   ! 
     886   ln_zdfddm   = .false.   ! double diffusive mixing 
     887      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     888      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     889   ! 
     890   !                       ! gravity wave-driven vertical mixing 
     891   ln_zdfiwm   = .false.      ! internal wave-induced mixing            (T =>   fill namzdf_iwm) 
     892   ln_zdfswm   = .false.      ! surface  wave-induced mixing            (T => ln_wave=ln_sdw=T ) 
     893   ! 
     894   !                       ! coefficients 
     895   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     896   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     897   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     898   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
     899/ 
     900!----------------------------------------------------------------------- 
     901&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T) 
    899902!----------------------------------------------------------------------- 
    900903   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity 
    901904   rn_alp      =    5.     !  coefficient of the parameterization 
    902905   nn_ric      =    2      !  coefficient of the parameterization 
    903    rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation 
    904    rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m) 
    905    rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m) 
    906    rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
    907    rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
    908    ln_mldw     =  .true.   !  Flag to use or not the mixed layer depth param. 
    909 / 
    910 !----------------------------------------------------------------------- 
    911 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
     906   ln_mldw     =  .false.  !  enhanced mixing in the Ekman layer 
     907      rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation 
     908      rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m) 
     909      rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m) 
     910      rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
     911      rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
     912/ 
     913!----------------------------------------------------------------------- 
     914&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
    912915!----------------------------------------------------------------------- 
    913916   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     
    917920   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2] 
    918921   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it) 
     922   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
    919923   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
    920                            !                 = 1 bounded by the local vertical scale factor 
    921                            !                 = 2 first vertical derivative of mixing length bounded by 1 
    922                            !                 = 3 as =2 with distinct disspipative an mixing length scale 
    923    nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
     924   !                       !                 = 1 bounded by the local vertical scale factor 
     925   !                       !                 = 2 first vertical derivative of mixing length bounded by 1 
     926   !                       !                 = 3 as =2 with distinct disspipative an mixing length scale 
    924927   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
    925928   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
     929   ln_drg      = .false.   !  top/bottom friction added as boundary condition of TKE 
    926930   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    927    rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    928    nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to near intertial waves 
    929                            !        = 0 no penetration 
    930                            !        = 1 add a tke source below the ML 
    931                            !        = 2 add a tke source just at the base of the ML 
    932                            !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
    933    rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
    934    nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
    935                            !        = 0  constant 10 m length scale 
    936                            !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    937 / 
    938 !----------------------------------------------------------------------- 
    939 &namzdf_gls    !   GLS vertical diffusion                               ("key_zdfgls") 
     931      rn_lc       =   0.15    !  coef. associated to Langmuir cells 
     932   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs 
     933                              !        = 0 none ; = 1 add a tke source below the ML 
     934                              !        = 2 add a tke source just at the base of the ML 
     935                              !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
     936      rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
     937      nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
     938                              !        = 0  constant 10 m length scale 
     939                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
     940/ 
     941!----------------------------------------------------------------------- 
     942&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T) 
    940943!----------------------------------------------------------------------- 
    941944   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2] 
     
    947950   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
    948951   rn_hsro       =  0.02   !  Minimum surface roughness 
    949    rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2) 
     952   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1) 
    950953   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3) 
    951954   !                             ! =3 requires ln_wave=T 
     
    956959/ 
    957960!----------------------------------------------------------------------- 
    958 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    959 !----------------------------------------------------------------------- 
    960    rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity) 
    961    rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio 
    962 / 
    963 !----------------------------------------------------------------------- 
    964 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    965 !----------------------------------------------------------------------- 
    966    rn_htmx     = 500.      !  vertical decay scale for turbulence (meters) 
    967    rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1) 
    968    rn_tfe      = 0.333     !  tidal dissipation efficiency 
    969    rn_me       = 0.2       !  mixing efficiency 
    970    ln_tmx_itf  = .true.    !  ITF specific parameterisation 
    971    rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency 
    972 / 
    973 !----------------------------------------------------------------------- 
    974 &namzdf_tmx_new !   internal wave-driven mixing parameterization        ("key_zdftmx_new" & "key_zdfddm") 
     961&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T) 
    975962!----------------------------------------------------------------------- 
    976963   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 
     
    981968!!                  ***  Miscellaneous namelists  *** 
    982969!!====================================================================== 
    983 !!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
     970!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi") 
    984971!!   namctl            Control prints  
    985972!!   namsto            Stochastic parametrization of EOS 
     
    987974! 
    988975!----------------------------------------------------------------------- 
    989 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     976&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    990977!----------------------------------------------------------------------- 
    991978   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send, 
    992                            !  buffer blocking send or immediate non-blocking sends, resp. 
     979   !                       !  buffer blocking send or immediate non-blocking sends, resp. 
    993980   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
    994981   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold 
     
    1008995   nn_isplt    =    1      !  number of processors in i-direction 
    1009996   nn_jsplt    =    1      !  number of processors in j-direction 
    1010    nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    1011    nn_diacfl   =    0      !  Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 
     997   ln_timing   = .false.   !  timing by routine write out in timing.output file 
     998   ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
    1012999/ 
    10131000!----------------------------------------------------------------------- 
  • branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/SHARED/namelist_top_ref

    r8599 r8881  
    6161/ 
    6262!----------------------------------------------------------------------- 
    63 &namtrc_adv      !   advection scheme for passive tracer  
     63&namtrc_adv      !   advection scheme for passive tracer                (default: NO selection) 
    6464!----------------------------------------------------------------------- 
     65   ln_trcadv_NONE=  .false.  !  No passive tracer advection 
    6566   ln_trcadv_cen =  .false.  !  2nd order centered scheme 
    6667      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     
    6970      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
    7071      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
    71       nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
    72       !                             !        (number of sub-timestep = nn_fct_zts) 
    7372   ln_trcadv_mus =  .false.  !  MUSCL scheme 
    7473      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     
    9392   ! 
    9493   rn_fact_lap      =     1.     !     enhanced zonal eddy diffusivity 
    95 / 
    96 !----------------------------------------------------------------------- 
    97 &namtrc_zdf      !   vertical physics 
    98 !----------------------------------------------------------------------- 
    99    ln_trczdf_exp =  .false.  !  split explicit (T) or implicit (F) time stepping 
    100    nn_trczdf_exp =   3       !  number of sub-timestep for ln_trczdfexp=T 
    10194/ 
    10295!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.