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 13710 for NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2020-11-02T10:56:42+01:00 (4 years ago)
Author:
emanuelaclementi
Message:

branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves: merge with trunk@13708, see #2155 and #2339

Location:
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@13559        sette 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/cfgs/SHARED/field_def_nemo-ice.xml

    r12377 r13710  
    4949          <field id="icehpnd"      long_name="melt pond depth"                                         standard_name="sea_ice_meltpond_depth"                    unit="m" />  
    5050          <field id="icevpnd"      long_name="melt pond volume"                                        standard_name="sea_ice_meltpond_volume"                   unit="m" />  
     51          <field id="icehlid"      long_name="melt pond lid depth"                                     standard_name="sea_ice_meltpondlid_depth"                 unit="m" />  
     52          <field id="icevlid"      long_name="melt pond lid volume"                                    standard_name="sea_ice_meltpondlid_volume"                unit="m" />  
    5153      
    5254     <!-- heat --> 
     
    7375          <field id="utau_bi"      long_name="X-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_x_stress"              unit="N/m2" /> 
    7476          <field id="vtau_bi"      long_name="Y-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_y_stress"              unit="N/m2" /> 
    75           <field id="isig1"        long_name="1st principal stress component for EVP rhg"                                                                        unit=""     /> 
    76           <field id="isig2"        long_name="2nd principal stress component for EVP rhg"                                                                        unit=""     /> 
    77           <field id="isig3"        long_name="convergence measure for EVP rheology (must be around 1)"                                                           unit=""     /> 
     77          <field id="sig1_pnorm"   long_name="P-normalized 1st principal stress component"                                                                       unit=""     /> 
     78          <field id="sig2_pnorm"   long_name="P-normalized 2nd principal stress component"                                                                       unit=""     /> 
    7879          <field id="normstr"      long_name="Average normal stress in sea ice"                        standard_name="average_normal_stress"                     unit="N/m"  /> 
    7980          <field id="sheastr"      long_name="Maximum shear stress in sea ice"                         standard_name="maximum_shear_stress"                      unit="N/m"  /> 
     
    8182          <field id="icediv"       long_name="Divergence of the sea-ice velocity field"                standard_name="divergence_of_sea_ice_velocity"            unit="s-1"  /> 
    8283          <field id="iceshe"       long_name="Maximum shear of sea-ice velocity field"                 standard_name="maximum_shear_of_sea_ice_velocity"         unit="s-1"  /> 
    83       
     84          <field id="beta_evp"     long_name="Relaxation parameter of ice rheology (beta)"             standard_name="relaxation_parameter_of_ice_rheology"      unit=""  />    
     85  
    8486     <!-- surface heat fluxes --> 
    8587          <field id="qt_ice"       long_name="total heat flux at ice surface"                          standard_name="surface_downward_heat_flux_in_air"         unit="W/m2" /> 
     
    162164 
    163165          <!-- diags --> 
    164           <field id="hfxdhc"       long_name="Heat content variation in snow and ice (neg = ice cooling)"   unit="W/m2" /> 
     166          <field id="hfxdhc"        long_name="Heat content variation in snow and ice (neg = ice cooling)"   unit="W/m2" /> 
     167          <!-- available if ln_icediachk=T --> 
     168          <field id="icedrift_mass" long_name="Ice mass drift (conservation check)"   unit="kg/m2/s" /> 
     169          <field id="icedrift_salt" long_name="Ice salt drift (conservation check)"   unit="kg/m2/s" /> 
     170          <field id="icedrift_heat" long_name="Ice heat drift (conservation check)"   unit="W/m2"    /> 
    165171 
    166172     <!-- sbcssm variables --> 
     
    172178          <field id="e3t_m"    unit="m"    /> 
    173179          <field id="frq_m"    unit="-"    /> 
     180 
     181          <!-- rheology convergence tests --> 
     182          <field id="uice_cvg"   long_name="sea ice velocity convergence"      standard_name="sea_ice_velocity_convergence"      unit="m/s" /> 
    174183 
    175184     <!-- ================= --> 
     
    211220          <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" /> 
    212221          <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" /> 
     222          <field id="dmilam"       long_name="sea-ice mass change through lateral melting"             standard_name="tendency_of_sea_ice_amount_due_to_lateral_melting"                       unit="kg/m2/s" /> 
    213223          <field id="dmsspr"       long_name="snow mass change through snow fall"                      standard_name="snowfall_flux"                                                           unit="kg/m2/s" /> 
    214224          <field id="dmsmel"       long_name="snow mass change through melt"                           standard_name="surface_snow_melt_flux"                                                  unit="kg/m2/s" /> 
     
    289299          <field id="iceapnd_cat"  long_name="Ice melt pond concentration per category"          unit=""        />  
    290300          <field id="icehpnd_cat"  long_name="Ice melt pond thickness per category"              unit="m"       detect_missing_value="true" />  
     301          <field id="icehlid_cat"  long_name="Ice melt pond lid thickness per category"          unit="m"       detect_missing_value="true" />  
    291302          <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category"               unit=""        />  
     303          <field id="iceaepnd_cat" long_name="Ice melt pond effective fraction per category"     unit=""        />  
    292304          <field id="icemask_cat"  long_name="Fraction of time step with sea ice (per category)" unit=""        /> 
    293305          <field id="iceage_cat"   long_name="Ice age per category"                              unit="days"    detect_missing_value="true" /> 
     
    300312          <field id="snwthic_cat_cmip"     long_name="Snow thickness in thickness categories"          standard_name="snow_thickness_over_categories"        detect_missing_value="true" unit="m"  > snwthic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
    301313          <field id="iceconc_cat_pct_cmip" long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="%"  > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field> 
     314 
     315          <!-- heat diffusion convergence tests --> 
     316          <field id="tice_cvgerr" long_name="sea ice temperature convergence error"      standard_name="sea_ice_temperature_convergence_err" unit="K" /> 
     317          <field id="tice_cvgstp" long_name="sea ice temperature convergence iterations" standard_name="sea_ice_temperature_convergence_stp" unit=""  /> 
    302318 
    303319   </field_group> <!-- SBC_3D --> 
     
    387403     <field field_ref="normstr"          name="normstr" /> 
    388404     <field field_ref="sheastr"          name="sheastr" /> 
    389      <field field_ref="isig1"            name="isig1"   /> 
    390      <field field_ref="isig2"            name="isig2"   /> 
    391      <field field_ref="isig3"            name="isig3"   /> 
     405     <field field_ref="sig1_pnorm"       name="sig1_pnorm"/> 
     406     <field field_ref="sig2_pnorm"       name="sig2_pnorm"/> 
    392407      
    393408     <!-- heat fluxes --> 
     
    446461          <field field_ref="vfxthin"          name="vfxthin" /> 
    447462    
    448      <!-- diag error for negative ice volume after advection --> 
    449      <field field_ref="iceneg_pres"      name="sineg_pres" /> 
    450      <field field_ref="iceneg_volu"      name="sineg_volu" /> 
    451      <field field_ref="iceneg_hfx"       name="sineg_hfx"  /> 
    452463        </field_group> 
    453464 
     
    560571          <field field_ref="dmisum"           name="sidmassmelttop"   /> 
    561572          <field field_ref="dmibom"           name="sidmassmeltbot"   /> 
     573          <field field_ref="dmilam"           name="sidmassmeltlat"   /> 
    562574          <field field_ref="dmsspr"           name="sndmasssnf"       /> 
    563575          <field field_ref="dmsmel"           name="sndmassmelt"      /> 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/cfgs/SHARED/field_def_nemo-oce.xml

    r12377 r13710  
    1 <?xml version="1.0"?>  
     1<?xml version="1.0"?> 
    22    <!-- $id$ --> 
    33 
     
    1616                         Configuration of multiple-linear-regression analysis (diamlr) 
    1717       ===================================================================================================== 
    18         
     18 
    1919       This field group configures diamlr for tidal harmonic analysis of field 
    2020       ssh: in addition to a regressor for fitting the mean value (diamlr_r101), 
     
    7373 
    7474    </field_group> 
    75      
    76     <!--  
     75 
     76    <!-- 
    7777============================================================================================================ 
    7878=                                  definition of all existing variables                                    = 
     
    101101    </field_group> 
    102102 
    103     <!--  
     103    <!-- 
    104104============================================================================================================ 
    105105                                  Physical ocean model variables 
     
    108108 
    109109      <!-- T grid --> 
    110        
     110 
    111111      <field_group id="grid_T" grid_ref="grid_T_2D" > 
    112112        <field id="e3t"          long_name="T-cell thickness"                    standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_3D" /> 
    113113        <field id="e3ts"         long_name="T-cell thickness"   field_ref="e3t"  standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_SFC"/> 
    114114        <field id="e3t_0"        long_name="Initial T-cell thickness"            standard_name="ref_cell_thickness"    unit="m"   grid_ref="grid_T_3D" /> 
    115         <field id="e3tb"         long_name="bottom T-cell thickness"             standard_name="bottom_cell_thickness" unit="m"   grid_ref="grid_T_2D"/>  
     115        <field id="e3tb"         long_name="bottom T-cell thickness"             standard_name="bottom_cell_thickness" unit="m"   grid_ref="grid_T_2D"/> 
    116116        <field id="e3t_300"      field_ref="e3t"                grid_ref="grid_T_zoom_300"       detect_missing_value="true" /> 
    117117        <field id="e3t_vsum300"  field_ref="e3t_300"            grid_ref="grid_T_vsum"   detect_missing_value="true" /> 
    118118   <field id="masscello"    long_name="Sea Water Mass per unit area"   standard_name="sea_water_mass_per_unit_area"   unit="kg/m2"   grid_ref="grid_T_3D"/> 
    119         <field id="volcello"     long_name="Ocean Volume"                   standard_name="ocean_volume"   unit="m3"       grid_ref="grid_T_3D"/>  
     119        <field id="volcello"     long_name="Ocean Volume"                   standard_name="ocean_volume"   unit="m3"       grid_ref="grid_T_3D"/> 
    120120        <field id="toce"         long_name="temperature"                         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
    121121        <field id="toce_e3t"     long_name="temperature (thickness weighted)"                                                      unit="degC"     grid_ref="grid_T_3D" > toce * e3t </field > 
     
    129129        <!-- AGRIF sponge --> 
    130130        <field id="agrif_spt"         long_name=" AGRIF t-sponge coefficient"   unit=" " /> 
     131    
     132   <!-- additions to diawri.F90 --> 
     133        <field id="socegrad"    long_name="module of salinity gradient"              unit="psu/m"   grid_ref="grid_T_3D"/> 
     134        <field id="socegrad2"   long_name="square of module of salinity gradient"    unit="psu2/m2" grid_ref="grid_T_3D"/> 
     135        <field id="ke"          long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" /> 
     136        <field id="ke_int"      long_name="vertical integration of kinetic energy"   unit="m3/s2"   /> 
     137        <field id="relvor"      long_name="relative vorticity"                       unit="s-1"    grid_ref="grid_T_3D"/> 
     138        <field id="absvor"      long_name="absolute vorticity"                       unit="s-1"    grid_ref="grid_T_3D"/> 
     139        <field id="potvor"      long_name="potential vorticity"                      unit="s-1"    grid_ref="grid_T_3D"/> 
     140        <field id="salt2c"      long_name="Salt content vertically integrated"       unit="1e-3*kg/m2" /> 
    131141 
    132142        <!-- t-eddy viscosity coefficients (ldfdyn) --> 
     
    146156        <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                                 unit="degC"     /> 
    147157   <field id="temp_3m"      long_name="temperature at 3m"                                                                                      unit="degC"     /> 
    148          
     158 
    149159        <field id="sss"          long_name="sea surface salinity"                               standard_name="sea_surface_salinity"                unit="1e-3"     /> 
    150160        <field id="sss2"         long_name="square of sea surface salinity"                                                                         unit="1e-6"      > sss * sss </field > 
     
    152162        <field id="sssmin"       long_name="min of sea surface salinity"      field_ref="sss"   operation="minimum"                                                 /> 
    153163        <field id="sbs"          long_name="sea bottom salinity"                                                                                    unit="0.001"    /> 
    154         <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" />  
    155  
    156         <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     />  
     164        <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" /> 
     165 
     166        <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     /> 
    157167 
    158168         <!-- Case EOS = TEOS-10 : output potential temperature --> 
     
    177187        <field id="alpha"        long_name="thermal expansion"                                                         unit="degC-1" grid_ref="grid_T_3D" /> 
    178188        <field id="beta"         long_name="haline contraction"                                                        unit="1e3"    grid_ref="grid_T_3D" /> 
    179         <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                                           unit="s-1"    grid_ref="grid_T_3D" /> 
    180189        <field id="rhop"         long_name="potential density (sigma0)"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D" /> 
    181190 
    182191        <!-- Energy - horizontal divergence --> 
    183         <field id="eken"         long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" /> 
    184192        <field id="hdiv"         long_name="horizontal divergence"                                                          unit="s-1"    grid_ref="grid_T_3D" /> 
    185193 
     
    295303        <field id="us_y"        long_name="j component of Stokes drift"                      unit="m/s"     /> 
    296304      </field_group> 
    297        
     305 
    298306      <!-- SBC --> 
    299307      <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> 
     
    311319          <field id="precip"       long_name="Total precipitation"                  standard_name="precipitation_flux"                                                   unit="kg/m2/s"   /> 
    312320          <field id="wclosea"      long_name="closed sea empmr correction"          standard_name="closea_empmr"                                                         unit="kg/m2/s"   /> 
    313       
     321 
    314322          <field id="qt"           long_name="Net Downward Heat Flux"                standard_name="surface_downward_heat_flux_in_sea_water"                              unit="W/m2"                           /> 
    315323          <field id="qns"          long_name="non solar Downward Heat Flux"                                                                                               unit="W/m2"                           /> 
     
    321329          <field id="taum"         long_name="wind stress module"                    standard_name="magnitude_of_surface_downward_stress"                                 unit="N/m2"                           /> 
    322330          <field id="wspd"         long_name="wind speed module"                     standard_name="wind_speed"                                                           unit="m/s"                            /> 
    323            
     331 
    324332          <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 
    325333          <field id="ssh_ib"       long_name="Inverse barometer sea surface height"  standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency"   unit="m"        /> 
     
    369377          <field id="taum_oce"     long_name="wind stress module over open ocean"           standard_name="magnitude_of_surface_downward_stress"               unit="N/m2"  /> 
    370378 
     379          <!-- variables computed by the bulk parameterization algorithms (ln_blk) --> 
     380          <field id="Cd_oce"      long_name="Drag coefficient over open ocean"              standard_name="drag_coefficient_water"                unit=""  /> 
     381          <field id="Ce_oce"      long_name="Evaporaion coefficient over open ocean"        standard_name="evap_coefficient_water"                unit=""  /> 
     382          <field id="Ch_oce"      long_name="Sensible heat coefficient over open ocean"     standard_name="sensible_heat_coefficient_water"       unit=""  /> 
     383          <field id="theta_zt"    long_name="Potential air temperature at z=zt"             standard_name="potential_air_temperature_at_zt"       unit="degC" /> 
     384          <field id="q_zt"        long_name="Specific air humidity at z=zt"                 standard_name="specific_air_humidity_at_zt"           unit="kg/kg" /> 
     385          <field id="theta_zu"    long_name="Potential air temperature at z=zu"             standard_name="potential_air_temperature_at_zu"       unit="degC" /> 
     386          <field id="q_zu"        long_name="Specific air humidity at z=zu"                 standard_name="specific_air_humidity_at_zu"           unit="kg/kg" /> 
     387          <field id="ssq"         long_name="Saturation specific humidity of air at z=0"    standard_name="surface_air_saturation_spec_humidity"  unit="kg/kg" /> 
     388          <field id="wspd_blk"    long_name="Bulk wind speed at z=zu"                       standard_name="bulk_wind_speed_at_zu"                 unit="m/s"   /> 
     389          <!-- ln_blk + key_si3 --> 
     390          <field id="Cd_ice"      long_name="Drag coefficient over ice"                     standard_name="drag_coefficient_ice"                 unit=""  /> 
     391          <field id="Ce_ice"      long_name="Evaporaion coefficient over ice"               standard_name="evap_coefficient_ice"                 unit=""  /> 
     392          <field id="Ch_ice"      long_name="Sensible heat coefficient over ice"            standard_name="sensible_heat_coefficient_ice"        unit=""  /> 
     393 
    371394          <!-- available key_oasis3 --> 
    372395          <field id="snow_ao_cea"  long_name="Snow over ice-free ocean (cell average)"   standard_name="snowfall_flux"                             unit="kg/m2/s"  /> 
     
    405428          <!-- ice field (nn_ice=1)  --> 
    406429          <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            /> 
    407            
     430 
    408431          <!-- dilution --> 
    409432          <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                                                                                              unit="kg*degC/m2/s" /> 
    410           <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" />         
     433          <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" /> 
    411434          <field id="rnf_x_sst"    long_name="Runoff term on SST"                                                                                                              unit="kg*degC/m2/s" /> 
    412435          <field id="rnf_x_sss"    long_name="Runoff term on SSS"                                                                                                              unit="kg*1e-3/m2/s" /> 
    413       
     436 
    414437     <!-- sbcssm variables --> 
    415438          <field id="sst_m"    unit="degC" /> 
     
    422445 
    423446   </field_group> 
    424     
     447 
    425448 
    426449      </field_group> <!-- SBC --> 
    427        
     450 
    428451      <!-- ABL --> 
    429452      <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 
     
    440463          <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
    441464          <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
    442           <field id="coeft"      long_name="ABL nudging coefficient"       standard_name="coeft"          unit=""         /> 
     465          <field id="u_geo"      long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     466          <field id="v_geo"      long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
    443467          <field id="tke_abl"    long_name="ABL turbulent kinetic energy"  standard_name="abl_tke"        unit="m2/s2"    /> 
    444468          <field id="avm_abl"    long_name="ABL turbulent viscosity"       standard_name="abl_avm"        unit="m2/s"     /> 
    445469          <field id="avt_abl"    long_name="ABL turbulent diffusivity"     standard_name="abl_avt"        unit="m2/s"     /> 
    446           <field id="mxl_abl"    long_name="ABL mixing length"             standard_name="abl_mxl"        unit="m"        /> 
     470          <field id="mxlm_abl"   long_name="ABL master mixing length"      standard_name="abl_mxlm"       unit="m"        /> 
     471          <field id="mxld_abl"   long_name="ABL dissipative mixing length" standard_name="abl_mxld"       unit="m"        /> 
    447472   </field_group> 
    448473 
     
    456481          <field id="uz1_dta"    long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
    457482          <field id="vz1_dta"    long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
    458           <field id="uvz1_dta"   long_name="DTA wind speed module"         standard_name="dta_wind_speed" unit="m/s"       > sqrt( uz1_dta^2 + vz1_dta^2 ) </field>  
     483          <field id="uvz1_dta"   long_name="DTA wind speed module"         standard_name="dta_wind_speed" unit="m/s"       > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 
    459484          <field id="tz1_dta"    long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
    460485          <field id="qz1_dta"    long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
     
    462487          <field id="uz1_geo"    long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
    463488          <field id="vz1_geo"    long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
    464           <field id="uvz1_geo"   long_name="GEO wind speed module"         standard_name="geo_wind_speed" unit="m/s"       > sqrt( uz1_geo^2 + vz1_geo^2 ) </field>  
     489          <field id="uvz1_geo"   long_name="GEO wind speed module"         standard_name="geo_wind_speed" unit="m/s"       > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 
    465490   </field_group> 
    466491 
    467492      </field_group> <!-- ABL --> 
    468493 
    469        
     494 
    470495      <!-- U grid --> 
    471        
     496 
    472497      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    473498        <field id="e2u"           long_name="U-cell width in meridional direction"                   standard_name="cell_width"                  unit="m"                               /> 
     
    478503      <field id="uoce_e3u"      long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
    479504      <field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical"  field_ref="uoce_e3u"    unit="m3/s"       grid_ref="grid_U_vsum"/> 
    480         <field id="uocetr_vsum"   long_name="ocean transport along i-axis  summed on the vertical"         field_ref="e2u"       unit="m3/s"> this * uoce_e3u_vsum  </field>  
     505        <field id="uocetr_vsum"   long_name="ocean transport along i-axis  summed on the vertical"         field_ref="e2u"       unit="m3/s"> this * uoce_e3u_vsum  </field> 
    481506 
    482507        <field id="uocetr_vsum_op"    long_name="ocean current along i-axis * e3u * e2u summed on the vertical"  read_access="true"  freq_op="1mo"    field_ref="e2u"       unit="m3/s"> @uocetr_vsum </field> 
    483         <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo"  unit="m3/s" />  
    484         <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 
     508        <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo"  unit="m3/s" /> 
     509        <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rho0 </field> 
    485510 
    486511 
     
    534559        <field id="udiff_salttr"  long_name="ocean diffusion salt transport along i-axis"                     standard_name="ocean_salt_x_transport_due_to_diffusion"        unit="1e-3*kg/s"                 /> 
    535560      </field_group> 
    536        
     561 
    537562      <!-- V grid --> 
    538        
     563 
    539564      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
    540565        <field id="e1v"          long_name="V-cell width in longitudinal direction"                 standard_name="cell_width"                  unit="m"                              /> 
     
    593618        <field id="vdiff_salttr"  long_name="ocean diffusion salt transport along j-axis"   standard_name="ocean_salt_y_transport_due_to_diffusion"         unit="1e-3*kg/s"                 /> 
    594619      </field_group> 
    595        
     620 
    596621      <!-- W grid --> 
    597        
     622 
    598623      <field_group id="grid_W" grid_ref="grid_W_3D"> 
    599624        <field id="e3w"          long_name="W-cell thickness"                              standard_name="cell_thickness"                         unit="m"    /> 
    600625   <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  /> 
    601    <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field>   
     626   <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field> 
    602627        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                                                   unit="m3/s" /> 
    603628 
     
    609634 
    610635   <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    611         <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field>      
     636        <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field> 
    612637        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"         standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    613638        <field id="avm"          long_name="vertical eddy viscosity"                        standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
    614         <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field>  
     639        <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field> 
    615640 
    616641        <!-- avs: /= avt with ln_zdfddm=T --> 
    617642        <field id="avs"          long_name="salt vertical eddy diffusivity"                 standard_name="ocean_vertical_salt_diffusivity"       unit="m2/s" /> 
    618         <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field>  
     643        <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field> 
    619644   <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    620645 
    621646        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    622647        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    623         <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field>  
     648        <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field> 
    624649   <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    625650 
     
    634659        <field id="wstokes"      long_name="Stokes Drift vertical velocity"                 standard_name="upward_StokesDrift_velocity"   unit="m/s" /> 
    635660 
    636         <!-- variables available with diaar5 -->    
     661        <!-- variables available with diaar5 --> 
    637662        <field id="w_masstr"     long_name="vertical mass transport"                        standard_name="upward_ocean_mass_transport"             unit="kg/s"   /> 
    638663        <field id="w_masstr2"    long_name="square of vertical mass transport"              standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" /> 
    639664 
     665        <!-- EOS --> 
     666        <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                unit="s-2" /> 
     667 
    640668      </field_group> 
    641        
     669 
    642670      <!-- F grid --> 
    643671      <!-- AGRIF sponge --> 
     
    683711         <field id="uocetr_vsum_section"  long_name="Total 2D transport in i-direction"               field_ref="uoce_e3u_ave_vsum"    grid_ref="grid_U_scalar"  detect_missing_value="true"> this * e2u </field> 
    684712         <field id="uocetr_strait"        long_name="Total transport across lines in i-direction"     field_ref="uocetr_vsum_section"  grid_ref="grid_U_4strait" /> 
    685          <field id="u_masstr_strait"      long_name="Sea water transport across line in i-direction"  field_ref="uocetr_strait"        grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rau0 </field> 
     713         <field id="u_masstr_strait"      long_name="Sea water transport across line in i-direction"  field_ref="uocetr_strait"        grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rho0 </field> 
    686714 
    687715         <field id="voce_e3v_ave"         long_name="Monthly average of v*e3v"                        field_ref="voce_e3v"                    freq_op="1mo"   freq_offset="_reset_" > @voce_e3v </field> 
     
    689717         <field id="vocetr_vsum_section"  long_name="Total 2D transport of in j-direction"            field_ref="voce_e3v_ave_vsum"    grid_ref="grid_V_scalar"  detect_missing_value="true"> this * e1v </field> 
    690718         <field id="vocetr_strait"        long_name="Total transport across lines in j-direction"     field_ref="vocetr_vsum_section"  grid_ref="grid_V_4strait"  /> 
    691          <field id="v_masstr_strait"      long_name="Sea water transport across line in j-direction"  field_ref="vocetr_strait"        grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rau0 </field> 
     719         <field id="v_masstr_strait"      long_name="Sea water transport across line in j-direction"  field_ref="vocetr_strait"        grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rho0 </field> 
    692720 
    693721         <field id="masstr_strait"        long_name="Sea water transport across line"                                                  grid_ref="grid_4strait"  > u_masstr_strait + v_masstr_strait </field> 
    694722      </field_group> 
    695723 
    696        
    697724      <!-- variables available with ln_floats --> 
    698725 
     
    709736      <!-- variables available with iceberg trajectories --> 
    710737 
    711       <field_group id="icbvar" domain_ref="grid_T"  >  
     738      <field_group id="icbvar" domain_ref="grid_T"  > 
    712739        <field id="berg_melt"          long_name="icb melt rate of icebergs"                       unit="kg/m2/s"                    /> 
    713740        <field id="berg_melt_hcflx"    long_name="icb heat flux to ocean due to melting heat content"   unit="J/m2/s"                /> 
     
    727754      </field_group> 
    728755 
    729       <!-- Poleward transport : ptr -->      
    730       <field_group id="diaptr" >   
     756      <!-- Poleward transport : ptr --> 
     757      <field_group id="diaptr" > 
    731758        <field id="zomsf"         long_name="Overturning Stream-Function : All basins"                     unit="Sv"         grid_ref="grid_znl_W_3D" /> 
    732759        <field id="zotem"         long_name="Zonal Mean Temperature : All basins"                          unit="degree_C"   grid_ref="grid_znl_T_3D" /> 
     
    736763        <field id="sopstove"      long_name="Overturning Salt Transport: All basins"                       unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    737764        <field id="sophtbtr"      long_name="Barotropic Heat Transport: All basins"                        unit="PW"         grid_ref="grid_znl_T_2D" /> 
    738         <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" />  
     765        <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    739766        <field id="sophtadv"      long_name="Advective Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" /> 
    740767        <field id="sopstadv"      long_name="Advective Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     
    742769        <field id="sopstldf"      long_name="Diffusive Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    743770        <field id="sophtvtr"      long_name="Heat Transport : All basins"                                  unit="PW"         grid_ref="grid_znl_T_2D" /> 
    744         <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" />   
     771        <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
    745772        <field id="sophteiv"      long_name="Heat Transport from mesoscale eddy advection: All basins"     unit="PW"         grid_ref="grid_znl_T_2D" /> 
    746773        <field id="sopsteiv"      long_name="Salt Transport from mesoscale eddy advection : All basins"    unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     
    758785 
    759786 
    760     <!--  
     787    <!-- 
    761788============================================================================================================ 
    762789                  Physical ocean model trend diagnostics : temperature, KE, PE, momentum 
     
    899926     <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                    unit="W/s^3"                        /> 
    900927     <field id="ketrd_tau"     long_name="ke-trend: wind stress "                           unit="W/s^3"   grid_ref="grid_T_2D" /> 
    901      <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        />    
    902      <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        />    
    903      <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        />   
     928     <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        /> 
     929     <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        /> 
     930     <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        /> 
    904931     <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"      unit="W/s^3"                        /> 
    905      <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        />    
     932     <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        /> 
    906933 
    907934     <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) --> 
    908      <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" />    
    909      <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" />    
     935     <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" /> 
     936     <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" /> 
    910937     <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 
    911      <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" />    
     938     <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" /> 
    912939 
    913940     <!-- variables available with ln_PE_trd --> 
     
    926953     <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"         unit="W/m^3"                        /> 
    927954     <field id="petrd_atf"     long_name="pe-trend: asselin time filter"        unit="W/m^3"                        /> 
    928      <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            />    
    929      <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       />    
    930      <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          />    
     955     <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            /> 
     956     <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       /> 
     957     <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          /> 
    931958   </field_group> 
    932959 
     
    945972     <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                    unit="m/s^2"                        /> 
    946973     <field id="utrd_tau"       long_name="i-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_U_2D" /> 
    947      <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        />    
    948      <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        />    
    949      <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        />    
    950      <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        />    
     974     <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        /> 
     975     <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        /> 
     976     <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        /> 
     977     <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        /> 
    951978   </field_group> 
    952979 
     
    965992     <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                    unit="m/s^2"                        /> 
    966993     <field id="vtrd_tau"       long_name="j-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_V_2D" /> 
    967      <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        />    
    968      <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        />    
    969      <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        />    
    970      <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        />    
     994     <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        /> 
     995     <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        /> 
     996     <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        /> 
     997     <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        /> 
    971998   </field_group> 
    972999 
    9731000 
    974     <!--  
     1001    <!-- 
    9751002============================================================================================================ 
    9761003                                        Definitions for iodef_demo.xml 
     
    9901017      <field field_ref="strd_zdfp_li"    name="osaltdiff" /> 
    9911018    </field_group> 
    992      
     1019 
    9931020    <field_group id="mooring" > 
    9941021      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"      /> 
     
    9991026      <field field_ref="avt"          name="difvho"   long_name="ocean_vertical_heat_diffusivity"      /> 
    10001027      <field field_ref="avm"          name="difvmo"   long_name="ocean_vertical_momentum_diffusivity"  /> 
    1001        
     1028 
    10021029      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       /> 
    10031030      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             /> 
     
    10471074      <field field_ref="BLT"          name="blt"      long_name="Barrier Layer Thickness"                       /> 
    10481075    </field_group> 
    1049      
     1076 
    10501077    <field_group id="groupU" > 
    10511078      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    10521079      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    10531080    </field_group> 
    1054      
     1081 
    10551082    <field_group id="groupV" > 
    10561083      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    10571084      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    10581085    </field_group> 
    1059      
     1086 
    10601087    <field_group id="groupW" > 
    10611088      <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  /> 
     
    11001127    </field_group> 
    11011128 
    1102     <!--  
     1129    <!-- 
    11031130============================================================================================================ 
    11041131    --> 
    1105     <!-- output variables for my configuration (example) -->  
    1106      
     1132    <!-- output variables for my configuration (example) --> 
     1133 
    11071134    <field_group id="myvarOCE" > 
    1108       <!-- grid T -->  
     1135      <!-- grid T --> 
    11091136      <field field_ref="e3t"          name="e3t"      long_name="vertical scale factor"           /> 
    11101137      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"         /> 
    11111138      <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"            /> 
    11121139      <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"  /> 
    1113        
    1114       <!-- grid U -->  
     1140 
     1141      <!-- grid U --> 
    11151142      <field field_ref="e3u"          name="e3u"     long_name="vertical scale factor"            /> 
    11161143      <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"           /> 
    1117        
    1118       <!-- grid V -->  
     1144 
     1145      <!-- grid V --> 
    11191146      <field field_ref="e3v"          name="e3v"     long_name="vertical scale factor"            /> 
    1120       <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           />      
    1121     </field_group>     
     1147      <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           /> 
     1148    </field_group> 
    11221149 
    11231150   </field_definition> 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/cfgs/SHARED/namelist_ice_ref

    r12377 r13710  
    4343   ln_cat_usr       = .false.         !  ice categories are defined by rn_catbnd below (m) 
    4444      rn_catbnd     =   0.,0.45,1.1,2.1,3.7,6.0   
    45    rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
     45   rn_himin         =   0.1           !  minimum ice thickness (m) allowed 
     46   rn_himax         =  99.0           !  maximum ice thickness (m) allowed 
    4647/ 
    4748!------------------------------------------------------------------------------ 
     
    5657   rn_ishlat        =   2.            !  lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 
    5758   ln_landfast_L16  = .false.         !  landfast: parameterization from Lemieux 2016 
    58       rn_depfra     =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast 
     59      rn_lf_depfra  =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast 
    5960                                      !          recommended range: [0.1 ; 0.25] 
    60       rn_icebfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
    61       rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
    62       rn_tensile    =   0.05          !        isotropic tensile strength [0-0.5??] 
     61      rn_lf_bfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
     62      rn_lf_relax   =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
     63      rn_lf_tensile =   0.05          !        isotropic tensile strength [0-0.5??] 
    6364/ 
    6465!------------------------------------------------------------------------------ 
     
    9192!------------------------------------------------------------------------------ 
    9293   ln_rhg_EVP       = .true.          !  EVP rheology 
    93       ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
     94      ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
    9495      rn_creepl     =   2.0e-9        !     creep limit [1/s] 
    9596      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve           
    96       nn_nevp       = 120             !     number of EVP subcycles                              
     97      nn_nevp       = 100             !     number of EVP subcycles                              
    9798      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
    98                                       !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     99                                      !        advised value: 1/3 (nn_nevp=100) or 1/9 (nn_nevp=300) 
     100   nn_rhg_chkcvg    =   0             !  check convergence of rheology 
     101                                      !     = 0  no check 
     102                                      !     = 1  check at the main time step (output xml: uice_cvg) 
     103                                      !     = 2  check at both main and rheology time steps (additional output: ice_cvg.nc) 
     104                                      !          this option 2 asks a lot of communications between cpu 
    99105/ 
    100106!------------------------------------------------------------------------------ 
    101107&namdyn_adv     !   Ice advection 
    102108!------------------------------------------------------------------------------ 
    103    ln_adv_Pra       = .true.         !  Advection scheme (Prather) 
    104    ln_adv_UMx       = .false.          !  Advection scheme (Ultimate-Macho) 
     109   ln_adv_Pra       = .true.          !  Advection scheme (Prather) 
     110   ln_adv_UMx       = .false.         !  Advection scheme (Ultimate-Macho) 
    105111      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
    106112/ 
     
    109115!------------------------------------------------------------------------------ 
    110116   rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-) 
    111    rn_blow_s        =   0.66          !  mesure of snow blowing into the leads 
     117   nn_snwfra        =   2             !  calculate the fraction of ice covered by snow (for zdf and albedo) 
     118                                      !     = 0  fraction = 1 (if snow) or 0 (if no snow) 
     119                                      !     = 1  fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation] 
     120                                      !     = 2  fraction = hsnw / (hsnw+0.02)    [CICE formulation] 
     121   rn_snwblow       =   0.66          !  mesure of snow blowing into the leads 
    112122                                      !     = 1 => no snow blowing, < 1 => some snow blowing 
    113123   nn_flxdist       =  -1             !  Redistribute heat flux over ice categories 
     
    118128   ln_cndflx        = .false.         !  Use conduction flux as surface boundary conditions (i.e. for Jules coupling) 
    119129      ln_cndemulate = .false.         !     emulate conduction flux (if not provided in the inputs) 
     130   nn_qtrice        =   1             !  Solar flux transmitted thru the surface scattering layer: 
     131                                      !     = 0  Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow)  
     132                                      !     = 1  Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 
    120133/ 
    121134!------------------------------------------------------------------------------ 
     
    126139   ln_icedO         = .true.          !  activate ice growth in open-water (T) or not (F) 
    127140   ln_icedS         = .true.          !  activate brine drainage (T) or not (F) 
     141   ! 
     142   ln_leadhfx       = .true.          !  heat in the leads is used to melt sea-ice before warming the ocean 
    128143/ 
    129144!------------------------------------------------------------------------------ 
     
    135150   rn_cnd_s         =   0.31          !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 
    136151                                      !     Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 
    137    rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice [1/m] 
     152   rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice                     [1/m] 
     153   rn_kappa_s       =  10.0           !  nn_qtrice = 0: radiation attenuation coefficient in snow         [1/m] 
     154   rn_kappa_smlt    =   7.0           !  nn_qtrice = 1: radiation attenuation coefficient in melting snow [1/m] 
     155   rn_kappa_sdry    =  10.0           !                 radiation attenuation coefficient in dry snow     [1/m] 
     156   ln_zdf_chkcvg    = .false.         !  check convergence of heat diffusion scheme (outputs: tice_cvgerr, tice_cvgstp) 
    138157/ 
    139158!------------------------------------------------------------------------------ 
     
    175194&namthd_pnd     !   Melt ponds 
    176195!------------------------------------------------------------------------------ 
    177    ln_pnd           = .false.         !  activate melt ponds or not 
    178      ln_pnd_H12     = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
    179      ln_pnd_CST     = .false.         !  activate constant  melt ponds 
    180        rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC 
    181        rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC 
    182      ln_pnd_alb     = .false.         !  melt ponds affect albedo or not 
     196   ln_pnd            = .true.         !  activate melt ponds or not 
     197      ln_pnd_LEV     = .true.         !  level ice melt ponds (from Flocco et al 2007,2010 & Holland et al 2012) 
     198         rn_apnd_min =   0.15         !     minimum ice fraction that contributes to melt pond. range: 0.0 -- 0.15 ?? 
     199         rn_apnd_max =   0.85         !     maximum ice fraction that contributes to melt pond. range: 0.7 -- 0.85 ?? 
     200      ln_pnd_CST     = .false.        !  constant  melt ponds 
     201         rn_apnd     =   0.2          !     prescribed pond fraction, at Tsu=0 degC 
     202         rn_hpnd     =   0.05         !     prescribed pond depth,    at Tsu=0 degC 
     203      ln_pnd_lids    = .true.         !  frozen lids on top of the ponds (only for ln_pnd_LEV) 
     204      ln_pnd_alb     = .true.         !  effect of melt ponds on ice albedo 
    183205/ 
    184206!------------------------------------------------------------------------------ 
     
    186208!------------------------------------------------------------------------------ 
    187209   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
    188    ln_iceini_file   = .false.         !  netcdf file provided for initialization (T) or not (F) 
     210   nn_iceini_file   =   0             !     0 = Initialise sea ice based on SSTs 
     211                                      !     1 = Initialise sea ice from single category netcdf file 
     212                                      !     2 = Initialise sea ice from multi category restart file 
    189213   rn_thres_sst     =   2.0           !  max temp. above Tfreeze with initial ice = (sst - tfreeze) 
    190214   rn_hti_ini_n     =   3.0           !  initial ice thickness       (m), North 
     
    206230   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North 
    207231   rn_hpd_ini_s     =   0.05          !        "            "             South 
    208    ! -- for ln_iceini_file = T 
    209    sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    210    sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    211    sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    212    sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    213    sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    214    sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    215    sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     232   rn_hld_ini_n     =   0.0           !  initial pond lid depth      (m), North 
     233   rn_hld_ini_s     =   0.0           !        "            "             South 
     234   ! -- for nn_iceini_file = 1 
     235   sn_hti = 'Ice_initialization'    , -12. ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     236   sn_hts = 'Ice_initialization'    , -12. ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     237   sn_ati = 'Ice_initialization'    , -12. ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     238   sn_smi = 'Ice_initialization'    , -12. ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     239   sn_tmi = 'Ice_initialization'    , -12. ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     240   sn_tsu = 'Ice_initialization'    , -12. ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     241   sn_tms = 'NOT USED'              , -12. ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    216242   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 
    217    sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    218    sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     243   sn_apd = 'NOT USED'              , -12. ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     244   sn_hpd = 'NOT USED'              , -12. ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     245   sn_hld = 'NOT USED'              , -12. ,'hld'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    219246   cn_dir='./' 
    220247/ 
     
    238265   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F) 
    239266   ln_icectl        = .false.         !  ice points output for debug (T or F) 
    240    iiceprt          =  10             !  i-index for debug 
    241    jiceprt          =  10             !  j-index for debug 
    242 / 
     267      iiceprt       =  10             !     i-index for debug 
     268      jiceprt       =  10             !     j-index for debug 
     269/ 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/cfgs/SHARED/namelist_pisces_ref

    r12377 r13710  
    352352! 
    353353   cn_dir      = './'      !  root directory for the location of the dynamical files 
    354    ln_ironsed  =  .true.   ! boolean for Fe input from sediments 
    355    ln_ironice  =  .true.   ! boolean for Fe input from sea ice 
    356    ln_hydrofe  =  .true.   ! boolean for from hydrothermal vents 
     354   ln_ironsed  =  .false.   ! boolean for Fe input from sediments 
     355   ln_ironice  =  .false.   ! boolean for Fe input from sea ice 
     356   ln_hydrofe  =  .false.   ! boolean for from hydrothermal vents 
    357357   sedfeinput  =  2.e-9    ! Coastal release of Iron 
    358358   distcoast   =  5.e3     ! Distance off the coast for Iron from sediments 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/cfgs/SHARED/namelist_ref

    r12991 r13710  
    259259      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
    260260      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
    261       !                     !  - module of the mean stress" data 
     261      ln_crt_fbk = .false.     !  Add surface current feedback to the wind stress (Renault et al. 2020, doi: 10.1029/2019MS001715) 
     262         rn_stau_a = -2.9e-3   !     Alpha from eq. 10: Stau = Alpha * Wnd + Beta 
     263         rn_stau_b =  8.0e-3   !     Beta  
    262264      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
    263265      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
    264       rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity 
    265       !                     !  used to calculate the wind stress 
    266       !                     ! (0. => absolute or 1. => relative winds) 
    267266      ln_skin_cs = .false.  !  use the cool-skin parameterization 
    268267      ln_skin_wl = .false.  !  use the warm-layer parameterization 
     
    271270      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
    272271      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     272      ln_tpot     = .true.  !!GS: compute potential temperature or not 
    273273   ! 
    274274   cn_dir      = './'      !  root directory for the bulk data location 
     
    282282   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    283283   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    284    sn_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
    285    sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    286284   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    287285   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    288286   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     287   sn_uoatm    = 'NOT USED'                   ,    6.        , 'UOATM'   ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , 'Uoceatm', '' 
     288   sn_voatm    = 'NOT USED'                   ,    6.        , 'VOATM'   ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , 'Voceatm', '' 
     289   sn_cc       = 'NOT USED'                   ,   24.        , 'CC'      ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     290   sn_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
     291   sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    289292/ 
    290293!----------------------------------------------------------------------- 
     
    299302   cn_ablrst_outdir = "."             !  directory to write output abl restarts 
    300303 
     304   ln_rstart_abl  = .false. 
    301305   ln_hpgls_frc   = .false. 
    302306   ln_geos_winds  = .false. 
    303    nn_dyn_restore = 2         ! restoring option for dynamical ABL variables: = 0 no restoring 
     307   ln_smth_pblh   = .false. 
     308   nn_dyn_restore = 0         ! restoring option for dynamical ABL variables: = 0 no restoring 
    304309                              !                                               = 1 equatorial restoring 
    305310                              !                                               = 2 global restoring 
    306    rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
    307    rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
    308    rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
    309    rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
     311   rn_ldyn_min   =  4.5       ! dynamics nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     312   rn_ldyn_max   =  1.5       ! dynamics nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
     313   rn_ltra_min   =  4.5       ! tracers  nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     314   rn_ltra_max   =  1.5       ! tracers  nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
    310315   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale 
    311316                              !                = 1 length-scale based on the distance to the PBL height 
    312317                              !                = 2 Bougeault & Lacarrere 89 length-scale 
    313    rn_Cm         = 0.0667     ! 0.126 in MesoNH 
    314    rn_Ct         = 0.1667     ! 0.143 in MesoNH 
    315    rn_Ce         = 0.4        ! 0.4   in MesoNH 
    316    rn_Ceps       = 0.7        ! 0.85  in MesoNH 
    317    rn_Rod        = 0.15       ! c0 in RMCA17 mixing length formulation (not yet implemented) 
    318    rn_Ric        = 0.139      !  Critical Richardson number (to compute PBL height and diffusivities) 
     318                              ! CBR00  ! CCH02  ! MesoNH ! 
     319   rn_Cm          = 0.0667    ! 0.0667 ! 0.1260 ! 0.1260 ! 
     320   rn_Ct          = 0.1667    ! 0.1667 ! 0.1430 ! 0.1430 ! 
     321   rn_Ce          = 0.40      ! 0.40   ! 0.34   ! 0.40   ! 
     322   rn_Ceps        = 0.700     ! 0.700  ! 0.845  ! 0.850  ! 
     323   rn_Ric         = 0.139     ! 0.139  ! 0.143  !   ?    ! Critical Richardson number (to compute PBL height and diffusivities) 
     324   rn_Rod         = 0.15      ! c0 in RMCA17 mixing length formulation (not yet implemented) 
    319325/ 
    320326!----------------------------------------------------------------------- 
    321327&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    322328!----------------------------------------------------------------------- 
    323    nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data 
    324    ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models 
    325    !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    326    nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
     329   nn_cplmodel       =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data 
     330   ln_usecplmask     = .false. !  use a coupling mask file to merge data received from several models 
     331   !                           !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
     332   ln_scale_ice_flux = .false. !  use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 
     333   nn_cats_cpl       =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
    327334   !_____________!__________________________!____________!_____________!______________________!________! 
    328335   !             !        description       !  multiple  !    vector   !       vector         ! vector ! 
     
    445452!----------------------------------------------------------------------- 
    446453   ln_rnf_mouth = .false.   !  specific treatment at rivers mouths 
    447       rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
    448       rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
    449    rn_rfact    =   1.e0    !  multiplicative factor for runoff 
     454      rn_hrnf     =  15.e0     !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
     455      rn_avt_rnf  =   1.e-3    !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
     456   rn_rfact     =   1.e0    !  multiplicative factor for runoff 
    450457   ln_rnf_depth = .false.   !  read in depth information for runoff 
    451    ln_rnf_tem  = .false.   !  read in temperature information for runoff 
    452    ln_rnf_sal  = .false.   !  read in salinity information for runoff 
    453    ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file 
    454       rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
    455       rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
    456       nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0) 
    457  
    458    cn_dir      = './'      !  root directory for the runoff data location 
     458   ln_rnf_tem   = .false.   !  read in temperature information for runoff 
     459   ln_rnf_sal   = .false.   !  read in salinity information for runoff 
     460   ln_rnf_icb   = .false.   !  read iceberg flux 
     461   ln_rnf_depth_ini = .false.  ! compute depth at initialisation from runoff file 
     462      rn_rnf_max  = 5.735e-4   !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
     463      rn_dep_max  = 150.       !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
     464      nn_rnf_depth_file = 0    !  create (=1) a runoff depth file or not (=0) 
     465 
     466   cn_dir       = './'      !  root directory for the runoff data location 
    459467   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    460468   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    461469   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    462470   sn_rnf      = 'runoff_core_monthly'   ,        -1.        , 'sorunoff',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    463    sn_cnf      = 'runoff_core_monthly'   ,         0.        , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
     471   sn_cnf      = 'runoff_core_monthly'   ,       -12.        , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
    464472   sn_s_rnf    = 'runoffs'               ,        24.        , 'rosaline',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    465473   sn_t_rnf    = 'runoffs'               ,        24.        , 'rotemper',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    466    sn_dep_rnf  = 'runoffs'               ,         0.        , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
     474   sn_i_rnf    = 'NOT USED'              ,        24.        , 'xxxxxxxx',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     475   sn_dep_rnf  = 'runoffs'               ,       -12.        , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
    467476/ 
    468477!----------------------------------------------------------------------- 
     
    532541         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    533542         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    534          sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    535          sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     543         sn_isfpar_zmax = 'isfmlt_par',       0.       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     544         sn_isfpar_zmin = 'isfmlt_par',       0.       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    536545         !* 'spe' and 'oasis' case 
    537          sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     546         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    538547         !* 'bg03' case 
    539          sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     548         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    540549      ! 
    541550      ! ---------------- ice sheet coupling ------------------------------- 
     
    572581   sn_wmp      =  'sdw_ecwaves_orca2'    ,        6.         , 'wmp'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    573582   sn_wnum     =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_num'   ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     583   sn_tauoc    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    574584/ 
    575585!----------------------------------------------------------------------- 
     
    641651&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    642652!----------------------------------------------------------------------- 
    643    ln_agrif_2way = .true.  !  activate two way nesting 
    644    ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    645    rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    646    rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
    647    rn_trelax_tra = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
    648    rn_trelax_dyn = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
    649    ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
     653   ln_agrif_2way   = .true.  !  activate two way nesting 
     654   ln_init_chfrpar = .false. !  initialize child grids from parent 
     655   ln_spc_dyn      = .true.  !  use 0 as special value for dynamics 
     656   rn_sponge_tra   = 0.002   !  coefficient for tracer   sponge layer [] 
     657   rn_sponge_dyn   = 0.002   !  coefficient for dynamics sponge layer [] 
     658   rn_trelax_tra   = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
     659   rn_trelax_dyn   = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
     660   ln_chk_bathy    = .false. !  =T  check the parent bathymetry 
    650661/ 
    651662!----------------------------------------------------------------------- 
     
    729740   bn_aip      = 'NOT USED'              ,         24.       , 'siapnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    730741   bn_hip      = 'NOT USED'              ,         24.       , 'sihpnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     742   bn_hil      = 'NOT USED'              ,         24.       , 'sihlid'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    731743   ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 
    732744   rn_ice_tem  = 270.         !  arbitrary temperature               of incoming sea ice 
     
    735747   rn_ice_apnd = 0.2          !       --   pond fraction = a_ip/a_i            -- 
    736748   rn_ice_hpnd = 0.05         !       --   pond depth                          -- 
     749   rn_ice_hlid = 0.0          !       --   pond lid depth                      -- 
    737750/ 
    738751!----------------------------------------------------------------------- 
     
    747760!!                                                                    !! 
    748761!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    749 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    750 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     762!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     763!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    751764!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    752765!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
     
    756769&namdrg        !   top/bottom drag coefficient                          (default: NO selection) 
    757770!----------------------------------------------------------------------- 
    758    ln_OFF      = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
     771   ln_drg_OFF  = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
    759772   ln_lin      = .false.   !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top) 
    760773   ln_non_lin  = .false.   !  non-linear  drag: Cd = Cd0 |U| 
     
    762775   ! 
    763776   ln_drgimp   = .true.    !  implicit top/bottom friction flag 
    764 / 
    765 !----------------------------------------------------------------------- 
    766 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
     777      ln_drgice_imp = .true. ! implicit ice-ocean drag 
     778/ 
     779!----------------------------------------------------------------------- 
     780&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
    767781!----------------------------------------------------------------------- 
    768782   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
     
    775789/ 
    776790!----------------------------------------------------------------------- 
    777 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     791&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    778792!----------------------------------------------------------------------- 
    779793   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
     
    828842                                 ! 
    829843   !                     ! S-EOS coefficients (ln_seos=T): 
    830    !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     844   !                             !  rd(T,S,Z)*rho0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    831845   rn_a0       =  1.6550e-1      !  thermal expension coefficient 
    832846   rn_b0       =  7.6554e-1      !  saline  expension coefficient 
     
    10411055   ln_dynrnf       =  .false.    !  runoffs option enabled (T) or not (F) 
    10421056   ln_dynrnf_depth =  .false.    !  runoffs is spread in vertical (T) or not (F) 
    1043 !   fwbcorr        = 3.786e-06   !  annual global mean of empmr for ssh correction 
     1057   fwbcorr         =    0.0      !  annual global mean of empmr for ssh correction 
    10441058 
    10451059   cn_dir      = './'      !  root directory for the ocean data location 
     
    11321146   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it) 
    11331147   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
    1134    nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
     1148   nn_mxl      =   3       !  mixing length: = 0 bounded by the distance to surface and bottom 
    11351149   !                       !                 = 1 bounded by the local vertical scale factor 
    11361150   !                       !                 = 2 first vertical derivative of mixing length bounded by 1 
    11371151   !                       !                 = 3 as =2 with distinct dissipative an mixing length scale 
    11381152   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
     1153      nn_mxlice    = 2        ! type of scaling under sea-ice 
     1154                              !    = 0 no scaling under sea-ice 
     1155                              !    = 1 scaling with constant sea-ice thickness 
     1156                              !    = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 
     1157                              !    = 3 scaling with maximum sea-ice thickness 
     1158      rn_mxlice   = 10.       ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 
    11391159   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
    11401160   ln_mxhsw    = .false.   !  surface mixing length scale = F(wave height) 
    1141    ln_drg      = .false.   !  top/bottom friction added as boundary condition of TKE 
    11421161   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    11431162      rn_lc       =   0.15    !  coef. associated to Langmuir cells 
     
    11501169                              !        = 0  constant 10 m length scale 
    11511170                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    1152       rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 
     1171   nn_eice     =   1       !  attenutaion of langmuir & surface wave breaking under ice 
     1172   !                       !           = 0 no impact of ice cover on langmuir & surface wave breaking 
     1173   !                       !           = 1 weigthed by 1-TANH(10*fr_i) 
     1174   !                       !           = 2 weighted by 1-fr_i 
     1175   !                       !           = 3 weighted by 1-MIN(1,4*fr_i)    
    11531176   nn_bc_surf   =     1    !  surface condition (0/1=Dir/Neum) ! Only applicable for wave coupling (ln_cplwave=1) 
    11541177   nn_bc_bot    =     1    !  bottom condition (0/1=Dir/Neum) ! Only applicable for wave coupling (ln_cplwave=1) 
     
    11651188   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
    11661189   rn_hsro       =  0.02   !  Minimum surface roughness 
     1190   rn_hsri       =  0.03   !  Ice-ocean roughness 
    11671191   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1) 
    11681192   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3) 
    1169    !                             ! =3 requires ln_wave=T 
     1193   !                       !     = 3 requires ln_wave=T 
     1194   nn_z0_ice     =   1     !  attenutaion of surface wave breaking under ice 
     1195   !                       !           = 0 no impact of ice cover 
     1196   !                       !           = 1 roughness uses rn_hsri and is weigthed by 1-TANH(10*fr_i) 
     1197   !                       !           = 2 roughness uses rn_hsri and is weighted by 1-fr_i 
     1198   !                       !           = 3 roughness uses rn_hsri and is weighted by 1-MIN(1,4*fr_i) 
    11701199   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
    11711200   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
     
    12031232   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    12041233   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    1205    sn_mpb      = 'NOT USED'              , -12               , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , '' 
    1206    sn_mpp      = 'NOT USED'              , -12               , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , '' 
    1207    sn_mpc      = 'NOT USED'              , -12               , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , '' 
    1208    sn_dsb      = 'NOT USED'              , -12               , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
    1209    sn_dsc      = 'NOT USED'              , -12               , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1234   sn_mpb      = 'NOT USED'              , -12.              , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1235   sn_mpp      = 'NOT USED'              , -12.              , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1236   sn_mpc      = 'NOT USED'              , -12.              , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1237   sn_dsb      = 'NOT USED'              , -12.              , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1238   sn_dsc      = 'NOT USED'              , -12.              , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
    12101239/ 
    12111240!!====================================================================== 
     
    12131242!!                                                                    !! 
    12141243!!   namtrd       dynamics and/or tracer trends                         (default: OFF) 
    1215 !!   namptr       Poleward Transport Diagnostics                        (default: OFF) 
    12161244!!   namhsb       Heat and salt budgets                                 (default: OFF) 
    12171245!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
     
    13961424   jpni        =   0       !  number of processors following i (set automatically if < 1), see also ln_listonly = T 
    13971425   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T 
     1426   nn_hls      =   1       !  halo width (applies to both rows and columns) 
    13981427/ 
    13991428!----------------------------------------------------------------------- 
    14001429&namctl        !   Control prints                                       (default: OFF) 
    14011430!----------------------------------------------------------------------- 
    1402    sn_cfctl%l_glochk = .FALSE.    ! Range sanity checks are local (F) or global (T). Set T for debugging only 
    1403    sn_cfctl%l_allon  = .FALSE.    ! IF T activate all options. If F deactivate all unless l_config is T 
    1404      sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
    1405        sn_cfctl%l_runstat = .TRUE.  ! switches and which areas produce reports with the proc integer settings. 
    1406        sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
    1407        sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
    1408        sn_cfctl%l_layout  = .FALSE. ! 
    1409        sn_cfctl%l_prtctl  = .FALSE. ! 
    1410        sn_cfctl%l_prttrc  = .FALSE. ! 
    1411        sn_cfctl%l_oasout  = .FALSE. ! 
    1412        sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
    1413        sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
    1414        sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
    1415        sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info 
    1416    nn_print    =    0      !  level of print (0 no extra print) 
    1417    nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
    1418    nn_ictle    =    0      !  end   i indice of control sum        multi processor runs 
    1419    nn_jctls    =    0      !  start j indice of control               over a subdomain) 
    1420    nn_jctle    =    0      !  end   j indice of control 
    1421    nn_isplt    =    1      !  number of processors in i-direction 
    1422    nn_jsplt    =    1      !  number of processors in j-direction 
    1423    ln_timing   = .false.   !  timing by routine write out in timing.output file 
    1424    ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
     1431   sn_cfctl%l_runstat = .TRUE.    ! switches and which areas produce reports with the proc integer settings. 
     1432   sn_cfctl%l_trcstat = .FALSE.   ! The default settings for the proc integers should ensure 
     1433   sn_cfctl%l_oceout  = .FALSE.   ! that  all areas report. 
     1434   sn_cfctl%l_layout  = .FALSE.   ! 
     1435   sn_cfctl%l_prtctl  = .FALSE.   ! 
     1436   sn_cfctl%l_prttrc  = .FALSE.   ! 
     1437   sn_cfctl%l_oasout  = .FALSE.   ! 
     1438   sn_cfctl%procmin   = 0         ! Minimum area number for reporting [default:0] 
     1439   sn_cfctl%procmax   = 1000000   ! Maximum area number for reporting [default:1000000] 
     1440   sn_cfctl%procincr  = 1         ! Increment for optional subsetting of areas [default:1] 
     1441   sn_cfctl%ptimincr  = 1         ! Timestep increment for writing time step progress info 
     1442   nn_ictls    =    0             !  start i indice of control sum (use to compare mono versus 
     1443   nn_ictle    =    0             !  end   i indice of control sum        multi processor runs 
     1444   nn_jctls    =    0             !  start j indice of control               over a subdomain) 
     1445   nn_jctle    =    0             !  end   j indice of control 
     1446   nn_isplt    =    1             !  number of processors in i-direction 
     1447   nn_jsplt    =    1             !  number of processors in j-direction 
     1448   ln_timing   = .false.          !  timing by routine write out in timing.output file 
     1449   ln_diacfl   = .false.          !  CFL diagnostics write out in cfl_diagnostics.ascii 
    14251450/ 
    14261451!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.