- Timestamp:
- 2019-12-05T12:06:36+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/README.rst
r10598 r12065 3 3 *********** 4 4 5 .. todo:: 6 7 8 5 9 .. contents:: 6 10 :local: 7 11 8 12 Output of diagnostics in NEMO is usually done using XIOS. 9 This is an efficient way of writing diagnostics because the time averaging, file writing and even some simple arithmetic or regridding is carried out in parallel to the NEMO model run. 13 This is an efficient way of writing diagnostics because 14 the time averaging, file writing and even some simple arithmetic or regridding is carried out in 15 parallel to the NEMO model run. 10 16 This page gives a basic introduction to using XIOS with NEMO. 11 Much more information is available from the XIOS homepageabove and from the NEMO manual.17 Much more information is available from the :xios:`XIOS homepage<>` above and from the NEMO manual. 12 18 13 Use of XIOS for diagnostics is activated using the pre-compiler key ``key_iomput``. 19 Use of XIOS for diagnostics is activated using the pre-compiler key ``key_iomput``. 14 20 15 21 Extracting and installing XIOS 16 ------------------------------ 22 ============================== 17 23 18 24 1. Install the NetCDF4 library. 19 If you want to use single file output you will need to compile the HDF & NetCDF libraries to allow parallel IO.20 2. Download the version of XIOS that you wish to use. The recommended version is now XIOS 2.5: 21 22 .. code-block:: console 25 If you want to use single file output you will need to compile the HDF & NetCDF libraries to 26 allow parallel IO. 27 2. Download the version of XIOS that you wish to use. 28 The recommended version is now XIOS 2.5: 23 29 24 $ svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/branchs/xios-2.5 xios-2.5 30 .. code-block:: console 25 31 26 and follow the instructions in `XIOS documentation <http://forge.ipsl.jussieu.fr/ioserver/wiki/documentation>`_ to compile it. 27 If you find problems at this stage, support can be found by subscribing to the `XIOS mailing list <http://forge.ipsl.jussieu.fr/mailman/listinfo.cgi/xios-users>`_ and sending a mail message to it. 32 $ svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/branchs/xios-2.5 33 34 and follow the instructions in :xios:`XIOS documentation <wiki/documentation>` to compile it. 35 If you find problems at this stage, support can be found by subscribing to 36 the :xios:`XIOS mailing list <../mailman/listinfo.cgi/xios-users>` and sending a mail message to it. 28 37 29 38 XIOS Configuration files 30 39 ------------------------ 31 40 32 XIOS is controlled using xml input files that should be copied to your model run directory before running the model. 33 Examples of these files can be found in the reference configurations (``cfgs``). The XIOS executable expects to find a file called ``iodef.xml`` in the model run directory. 34 In NEMO we have made the decision to use include statements in the ``iodef.xml`` file to include ``field_def_nemo-oce.xml`` (for physics), ``field_def_nemo-ice.xml`` (for ice), ``field_def_nemo-pisces.xml`` (for biogeochemistry) and ``domain_def.xml`` from the /cfgs/SHARED directory. 35 Most users will not need to modify ``domain_def.xml`` or ``field_def_nemo-???.xml`` unless they want to add new diagnostics to the NEMO code. 36 The definition of the output files is organized into separate ``file_definition.xml`` files which are included in the ``iodef.xml`` file. 41 XIOS is controlled using XML input files that should be copied to 42 your model run directory before running the model. 43 Examples of these files can be found in the reference configurations (:file:`./cfgs`). 44 The XIOS executable expects to find a file called :file:`iodef.xml` in the model run directory. 45 In NEMO we have made the decision to use include statements in the :file:`iodef.xml` file to include: 46 47 - :file:`field_def_nemo-oce.xml` (for physics), 48 - :file:`field_def_nemo-ice.xml` (for ice), 49 - :file:`field_def_nemo-pisces.xml` (for biogeochemistry) and 50 - :file:`domain_def.xml` from the :file:`./cfgs/SHARED` directory. 51 52 Most users will not need to modify :file:`domain_def.xml` or :file:`field_def_nemo-???.xml` unless 53 they want to add new diagnostics to the NEMO code. 54 The definition of the output files is organized into separate :file:`file_definition.xml` files which 55 are included in the :file:`iodef.xml` file. 37 56 38 57 Modes 39 ----- 58 ===== 40 59 41 60 Detached Mode … … 44 63 In detached mode the XIOS executable is executed on separate cores from the NEMO model. 45 64 This is the recommended method for using XIOS for realistic model runs. 46 To use this mode set ``using_server`` to ``true`` at the bottom of the ``iodef.xml`` file:65 To use this mode set ``using_server`` to ``true`` at the bottom of the :file:`iodef.xml` file: 47 66 48 67 .. code-block:: xml 49 68 50 69 <variable id="using_server" type="boolean">true</variable> 51 70 52 Make sure there is a copy (or link to) your XIOS executable in the working directory and in your job submission script allocate processors to XIOS. 71 Make sure there is a copy (or link to) your XIOS executable in the working directory and 72 in your job submission script allocate processors to XIOS. 53 73 54 74 Attached Mode … … 56 76 57 77 In attached mode XIOS runs on each of the cores used by NEMO. 58 This method is less efficient than the detached mode but can be more convenient for testing or with small configurations. 59 To activate this mode simply set ``using_server`` to false in the ``iodef.xml`` file 78 This method is less efficient than the detached mode but can be more convenient for testing or 79 with small configurations. 80 To activate this mode simply set ``using_server`` to false in the :file:`iodef.xml` file 60 81 61 82 .. code-block:: xml 62 83 63 84 <variable id="using_server" type="boolean">false</variable> 64 85 65 86 and don't allocate any cores to XIOS. 66 Note that due to the different domain decompositions between XIOS and NEMO if the total number of cores is larger than the number of grid points in the j direction then the model run will fail. 87 88 .. note:: 89 90 Due to the different domain decompositions between XIOS and NEMO, 91 if the total number of cores is larger than the number of grid points in the ``j`` direction then 92 the model run will fail. 67 93 68 94 Adding new diagnostics 69 ---------------------- 95 ====================== 70 96 71 97 If you want to add a NEMO diagnostic to the NEMO code you will need to do the following: 72 98 73 99 1. Add any necessary code to calculate you new diagnostic in NEMO 74 2. Send the field to XIOS using ``CALL iom_put( 'field_id', variable )`` where ``field_id`` is a unique id for your new diagnostics and variable is the fortran variable containing the data. 75 This should be called at every model timestep regardless of how often you want to output the field. No time averaging should be done in the model code. 76 3. If it is computationally expensive to calculate your new diagnostic you should also use "iom_use" to determine if it is requested in the current model run. For example, 77 78 .. code-block:: fortran 100 2. Send the field to XIOS using ``CALL iom_put( 'field_id', variable )`` where 101 ``field_id`` is a unique id for your new diagnostics and 102 variable is the fortran variable containing the data. 103 This should be called at every model timestep regardless of how often you want to output the field. 104 No time averaging should be done in the model code. 105 3. If it is computationally expensive to calculate your new diagnostic 106 you should also use "iom_use" to determine if it is requested in the current model run. 107 For example, 79 108 80 IF iom_use('field_id') THEN 81 !Some expensive computation 82 !... 83 !... 84 iom_put('field_id', variable) 85 ENDIF 109 .. code-block:: fortran 86 110 87 4. Add a variable definition to the ``field_def_nemo-???.xml`` file. 88 5. Add the variable to the ``iodef.xml`` or ``file_definition.xml`` file. 111 IF iom_use('field_id') THEN 112 !Some expensive computation 113 !... 114 !... 115 iom_put('field_id', variable) 116 ENDIF 117 118 4. Add a variable definition to the :file:`field_def_nemo-???.xml` file. 119 5. Add the variable to the :file:`iodef.xml` or :file:`file_definition.xml` file. -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/domain_def_nemo.xml
r9930 r12065 10 10 </domain> 11 11 12 <domain id="1point" domain_ref="grid_T" >13 <zoom_domain ibegin="1" jbegin="1" ni="1" nj="1"/>14 </domain>15 12 <!-- Eq section --> 16 13 <domain id="EqT" domain_ref="grid_T" > <zoom_domain id="EqT"/> </domain> -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/field_def_nemo-ice.xml
r10578 r12065 32 32 33 33 <!-- general fields --> 34 <field id="icemass" long_name="Sea-ice mass per area"standard_name="sea_ice_amount" unit="kg/m2"/>34 <field id="icemass" long_name="Sea-ice mass per gridcell area" standard_name="sea_ice_amount" unit="kg/m2"/> 35 35 <field id="snwmass" long_name="Snow mass per area" standard_name="liquid_water_content_of_surface_snow" unit="kg/m2"/> 36 <field id="iceconc" long_name="Sea-ice area fraction" standard_name="sea_ice_area_fraction" unit="" /> 37 <field id="icevolu" long_name="Sea-ice volume per area" standard_name="sea_ice_thickness" unit="m" /> 38 <field id="icethic" long_name="Sea-ice thickness per area" standard_name="sea_ice_thickness" unit="m" /> 39 <field id="snwthic" long_name="Snow thickness per area" standard_name="snow_thickness" unit="m" /> 40 <field id="icebrv" long_name="brine volume" unit="%" /> 41 <field id="iceage" long_name="ice age" unit="days"/> 42 <field id="icehnew" long_name="frazil ice collection thickness" unit="m" /> 43 <field id="snwvolu" long_name="snow volume" unit="m" /> 44 <field id="icefrb" long_name="Sea-ice freeboard" standard_name="sea_ice_freeboard" unit="m" /> 45 <field id="icealb" long_name="Sea-ice or snow albedo" standard_name="sea_ice_albedo" unit="" /> 46 <field id="tau_icebfr" long_name="ice friction on ocean bottom for landfast ice" unit="N/2" /> 36 <field id="iceconc" long_name="Sea-ice area fraction" standard_name="sea_ice_area_fraction" unit="" /> 37 <field id="icevolu" long_name="Sea-ice volume per area" standard_name="sea_ice_thickness" unit="m" /> 38 <field id="icethic" long_name="Sea-ice thickness" standard_name="sea_ice_thickness" unit="m" /> 39 <field id="snwthic" long_name="Snow thickness" standard_name="snow_thickness" unit="m" /> 40 <field id="icebrv" long_name="brine volume" unit="%" /> 41 <field id="iceage" long_name="ice age" unit="days" detect_missing_value="true"/> 42 <field id="icehnew" long_name="frazil ice collection thickness" unit="m" /> 43 <field id="snwvolu" long_name="snow volume" unit="m" /> 44 <field id="icefrb" long_name="Sea-ice freeboard" standard_name="sea_ice_freeboard" unit="m" /> 45 <field id="icealb" long_name="Sea-ice or snow albedo" standard_name="sea_ice_albedo" unit="" detect_missing_value="true" /> 47 46 48 47 <!-- melt ponds --> 49 <field id="iceapnd" long_name="melt pond fraction" standard_name="sea_ice_meltpond_fraction" unit="%" /> 48 <field id="iceapnd" long_name="melt pond concentration" standard_name="sea_ice_meltpond_concentration" unit="" /> 49 <field id="icehpnd" long_name="melt pond depth" standard_name="sea_ice_meltpond_depth" unit="m" /> 50 50 <field id="icevpnd" long_name="melt pond volume" standard_name="sea_ice_meltpond_volume" unit="m" /> 51 51 52 52 <!-- heat --> 53 <field id="icetemp" long_name="Mean ice temperature" unit="degC" />54 <field id="snwtemp" long_name="Mean snow temperature" unit="degC" />55 <field id="icettop" long_name="temperature at the ice surface" unit="degC" />56 <field id="icetbot" long_name="temperature at the ice bottom" unit="degC" />57 <field id="icetsni" long_name="temperature at the snow-ice interface" unit="degC" />53 <field id="icetemp" long_name="Mean ice temperature" unit="degC" detect_missing_value="true" /> 54 <field id="snwtemp" long_name="Mean snow temperature" unit="degC" detect_missing_value="true" /> 55 <field id="icettop" long_name="temperature at the ice surface" unit="degC" detect_missing_value="true" /> 56 <field id="icetbot" long_name="temperature at the ice bottom" unit="degC" detect_missing_value="true" /> 57 <field id="icetsni" long_name="temperature at the snow-ice interface" unit="degC" detect_missing_value="true" /> 58 58 <field id="icehc" long_name="ice heat content" unit="J/m2" /> 59 59 <field id="snwhc" long_name="snow heat content" unit="J/m2" /> 60 60 61 61 <!-- salt --> 62 <field id="icesalt" long_name="Sea ice salinity" unit="g/kg" />62 <field id="icesalt" long_name="Sea ice salinity" unit="g/kg" detect_missing_value="true" /> 63 63 <field id="icesalm" long_name="Mass of salt in sea ice per area" standard_name="sea_ice_salt_mass" unit="kg/m2" /> 64 64 … … 71 71 <field id="utau_oi" long_name="X-component of ocean stress on sea ice" standard_name="sea_ice_base_upward_x_stress" unit="N/m2" /> 72 72 <field id="vtau_oi" long_name="Y-component of ocean stress on sea ice" standard_name="sea_ice_base_upward_y_stress" unit="N/m2" /> 73 <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" /> 74 <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" /> 73 75 <field id="isig1" long_name="1st principal stress component for EVP rhg" unit="" /> 74 76 <field id="isig2" long_name="2nd principal stress component for EVP rhg" unit="" /> … … 92 94 93 95 <!-- trends --> 94 <field id="afxthd" long_name="sea-ice area fraction change from thermodynamics" 95 <field id="afxdyn" long_name="sea-ice area fraction change from dynamics" 96 <field id="afxtot" long_name="area tendency (total)" 96 <field id="afxthd" long_name="sea-ice area fraction change from thermodynamics" standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" /> 97 <field id="afxdyn" long_name="sea-ice area fraction change from dynamics" standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" /> 98 <field id="afxtot" long_name="area tendency (total)" unit="s-1" /> 97 99 98 100 <!-- momentum (advection) --> … … 160 162 <field id="hfxdhc" long_name="Heat content variation in snow and ice (neg = ice cooling)" unit="W/m2" /> 161 163 162 <!-- diagnostics of the negative values resulting from the advection scheme -->163 <field id="iceneg_pres" long_name="Fraction of time steps with negative sea ice volume" unit="" />164 <field id="iceneg_volu" long_name="Negative sea ice volume per area arising from advection" unit="m" />165 <field id="iceneg_hfx" long_name="Negative sea ice heat content (eq. heat flux) arising from advection" unit="W/m2" />166 167 164 <!-- sbcssm variables --> 168 165 <field id="sst_m" unit="degC" /> … … 173 170 <field id="e3t_m" unit="m" /> 174 171 <field id="frq_m" unit="-" /> 175 176 <!-- categories -->177 <field id="iceconc_cat" long_name="Sea-ice concentration per category" unit="" grid_ref="grid_T_3D_ncatice" />178 <field id="icethic_cat" long_name="Sea-ice thickness per category" unit="m" grid_ref="grid_T_3D_ncatice" />179 <field id="snwthic_cat" long_name="Snow thickness per category" unit="m" grid_ref="grid_T_3D_ncatice" />180 <field id="icesalt_cat" long_name="Sea-Ice Bulk salinity per category" unit="g/kg" grid_ref="grid_T_3D_ncatice" />181 <field id="icetemp_cat" long_name="Ice temperature per category" unit="degC" grid_ref="grid_T_3D_ncatice" />182 <field id="snwtemp_cat" long_name="Snow temperature per category" unit="degC" grid_ref="grid_T_3D_ncatice" />183 <field id="icettop_cat" long_name="Ice/snow surface temperature per category" unit="degC" grid_ref="grid_T_3D_ncatice" />184 <field id="iceapnd_cat" long_name="Ice melt pond concentration per category" unit="%" grid_ref="grid_T_3D_ncatice" />185 <field id="icehpnd_cat" long_name="Ice melt pond thickness per category" unit="m" grid_ref="grid_T_3D_ncatice" />186 <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category" unit="m" grid_ref="grid_T_3D_ncatice" />187 <field id="icemask_cat" long_name="Fraction of time step with sea ice (per category)" unit="" grid_ref="grid_T_3D_ncatice" />188 <field id="iceage_cat" long_name="Ice age per category" unit="days" grid_ref="grid_T_3D_ncatice" />189 <field id="icebrv_cat" long_name="Brine volume per category" unit="%" grid_ref="grid_T_3D_ncatice" />190 172 191 173 <!-- ================= --> … … 254 236 <field id="xmtrptot" long_name="X-component of sea-ice mass transport" standard_name="sea_ice_x_transport" unit="kg/s" > xmtrpice + xmtrpsnw </field> 255 237 <field id="ymtrptot" long_name="Y-component of sea-ice mass transport" standard_name="sea_ice_y_transport" unit="kg/s" > ymtrpice + ymtrpsnw </field> 256 257 <!-- categories -->258 <field id="iceconc_cat_cmip" long_name="Sea-ice area fractions in thickness categories" standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="" grid_ref="grid_T_3D_ncatice" > iceconc_cat * icemask_cat + $missval * (1.-icemask_cat) </field>259 <field id="icethic_cat_cmip" long_name="Sea-ice thickness in thickness categories" standard_name="sea_ice_thickness_over_categories" detect_missing_value="true" unit="m" grid_ref="grid_T_3D_ncatice" > icethic_cat * icemask_cat + $missval * (1.-icemask_cat) </field>260 <field id="snwthic_cat_cmip" long_name="Snow thickness in thickness categories" standard_name="snow_thickness_over_categories" detect_missing_value="true" unit="m" grid_ref="grid_T_3D_ncatice" > snwthic_cat * icemask_cat + $missval * (1.-icemask_cat) </field>261 <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="%" grid_ref="grid_T_3D_ncatice" > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field>262 238 263 239 </field_group> <!-- SBC_2D --> 264 240 241 <!-- categories --> 242 <field_group id="SBC_3D" grid_ref="grid_T_ncatice" > 243 244 <!-- standard ice fields --> 245 <field id="iceconc_cat" long_name="Sea-ice concentration per category" unit="" /> 246 <field id="icethic_cat" long_name="Sea-ice thickness per category" unit="m" detect_missing_value="true" /> 247 <field id="snwthic_cat" long_name="Snow thickness per category" unit="m" detect_missing_value="true" /> 248 <field id="icesalt_cat" long_name="Sea-Ice Bulk salinity per category" unit="g/kg" detect_missing_value="true" /> 249 <field id="icetemp_cat" long_name="Ice temperature per category" unit="degC" detect_missing_value="true" /> 250 <field id="snwtemp_cat" long_name="Snow temperature per category" unit="degC" detect_missing_value="true" /> 251 <field id="icettop_cat" long_name="Ice/snow surface temperature per category" unit="degC" detect_missing_value="true" /> 252 <field id="iceapnd_cat" long_name="Ice melt pond concentration per category" unit="" /> 253 <field id="icehpnd_cat" long_name="Ice melt pond thickness per category" unit="m" detect_missing_value="true" /> 254 <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category" unit="" /> 255 <field id="icemask_cat" long_name="Fraction of time step with sea ice (per category)" unit="" /> 256 <field id="iceage_cat" long_name="Ice age per category" unit="days" detect_missing_value="true" /> 257 <field id="icebrv_cat" long_name="Brine volume per category" unit="%" detect_missing_value="true" /> 258 <field id="icealb_cat" long_name="Sea-ice or snow albedo" unit="" detect_missing_value="true" /> 259 260 <!-- Add-ons for SIMIP --> 261 <field id="iceconc_cat_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 * icemask_cat + $missval * (1.-icemask_cat) </field> 262 <field id="icethic_cat_cmip" long_name="Sea-ice thickness in thickness categories" standard_name="sea_ice_thickness_over_categories" detect_missing_value="true" unit="m" > icethic_cat * icemask_cat + $missval * (1.-icemask_cat) </field> 263 <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> 264 <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> 265 266 </field_group> <!-- SBC_3D --> 267 265 268 <!-- scalar variables --> 266 <field_group id="SBC_0D" grid_ref="grid_1point" > 269 <field_group id="SBC_scalar" grid_ref="grid_scalar" > 270 <field id="NH_iceextt" long_name="Sea ice extent North" standard_name="sea_ice_extent_n" unit="1e6_km2" /> 271 <field id="SH_iceextt" long_name="Sea ice extent South" standard_name="sea_ice_extent_s" unit="1e6_km2" /> 272 <field id="NH_icevolu" long_name="Sea ice volume North" standard_name="sea_ice_volume_n" unit="1e3_km3" /> 273 <field id="SH_icevolu" long_name="Sea ice volume South" standard_name="sea_ice_volume_s" unit="1e3_km3" /> 274 <field id="NH_icearea" long_name="Sea ice area North" standard_name="sea_ice_area_n" unit="1e6_km2" /> 275 <field id="SH_icearea" long_name="Sea ice area South" standard_name="sea_ice_area_s" unit="1e6_km2" /> 276 267 277 <!-- available with ln_icediaout --> 268 278 <field id="ibgfrcvoltop" long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)" unit="km3" /> … … 286 296 <field id="sbgheat_tot" long_name="global mean snow heat content" unit="1e20J" /> 287 297 288 <field id="NH_iceextt" long_name="Sea ice extent North" standard_name="sea_ice_extent_n" unit="1e6_km2" /> 289 <field id="SH_iceextt" long_name="Sea ice extent South" standard_name="sea_ice_extent_s" unit="1e6_km2" /> 290 <field id="NH_icevolu" long_name="Sea ice volume North" standard_name="sea_ice_volume_n" unit="1e3_km3" /> 291 <field id="SH_icevolu" long_name="Sea ice volume South" standard_name="sea_ice_volume_s" unit="1e3_km3" /> 292 <field id="NH_icearea" long_name="Sea ice area North" standard_name="sea_ice_area_n" unit="1e6_km2" /> 293 <field id="SH_icearea" long_name="Sea ice area South" standard_name="sea_ice_area_s" unit="1e6_km2" /> 294 298 <!-- available later --> 299 <!-- 295 300 <field id="strait_mifl" long_name="Sea ice mass flux through straits" standard_name="sea_ice_mass_transport_across_line" unit="kg/s" grid_ref="grid_4strait" /> 296 301 <field id="strait_arfl" long_name="Sea ice area flux through straits" standard_name="sea_ice_area_transport_across_line" unit="m2/s" grid_ref="grid_4strait" /> 297 302 <field id="strait_msfl" long_name="Sea ice snow flux through straits" standard_name="snow_mass_transport_across_line" unit="kg/s" grid_ref="grid_4strait" /> 298 </field_group> <!-- SBC_0D --> 303 --> 304 </field_group> 299 305 300 306 <!-- … … 304 310 305 311 <field_group id="myvarICE" grid_ref="grid_T_2D" > 306 <field field_ref="icethic" name="sithic" /> 307 <field field_ref="icevolu" name="sivolu" /> 308 <field field_ref="iceconc" name="siconc" /> 312 <!-- ice mask --> 313 <field field_ref="icemask" name="simsk" /> 314 <field field_ref="icemask05" name="simsk05" /> 315 <field field_ref="icemask15" name="simsk15" /> 316 317 <!-- general --> 318 <field field_ref="snwvolu" name="snvolu" /> 319 <field field_ref="snwthic" name="snthic" /> 320 <field field_ref="icethic" name="sithic" /> 321 <!-- 322 <field field_ref="icethic" name="sithic_max" operation="maximum" /> 323 <field field_ref="icethic" name="sithic_min" operation="minimum" /> 324 --> 325 <field field_ref="fasticepres" name="fasticepres" /> 326 <field field_ref="icevolu" name="sivolu" /> 327 <field field_ref="iceconc" name="siconc" /> 328 <field field_ref="icesalt" name="sisali" /> 329 <field field_ref="iceapnd" name="siapnd" /> 330 <field field_ref="icehpnd" name="sihpnd" /> 331 <field field_ref="icevpnd" name="sivpnd" /> 332 <field field_ref="iceage" name="siage" /> 333 <field field_ref="sst_m" name="sst_m" /> 334 <field field_ref="sss_m" name="sss_m" /> 335 336 <!-- heat --> 337 <field field_ref="icetemp" name="sitemp" /> 338 <field field_ref="snwtemp" name="sntemp" /> 339 <field field_ref="icettop" name="sittop" /> 340 <field field_ref="icetbot" name="sitbot" /> 341 <field field_ref="icetsni" name="sitsni" /> 342 343 <!-- momentum --> 344 <field field_ref="uice" name="sivelu" /> 345 <field field_ref="vice" name="sivelv" /> 346 <field field_ref="icevel" name="sivelo" /> 347 <field field_ref="utau_ai" name="utau_ai" /> 348 <field field_ref="vtau_ai" name="vtau_ai" /> 349 <field field_ref="utau_oi" name="utau_oi" /> 350 <field field_ref="vtau_oi" name="vtau_oi" /> 351 352 <!-- rheology --> 353 <field field_ref="icediv" name="sidive" /> 354 <field field_ref="iceshe" name="sishea" /> 355 <field field_ref="icestr" name="sistre" /> 356 <field field_ref="normstr" name="normstr" /> 357 <field field_ref="sheastr" name="sheastr" /> 358 <field field_ref="isig1" name="isig1" /> 359 <field field_ref="isig2" name="isig2" /> 360 <field field_ref="isig3" name="isig3" /> 361 362 <!-- heat fluxes --> 363 <field field_ref="qt_oce_ai" name="qt_oce_ai" /> 364 <field field_ref="qt_atm_oi" name="qt_atm_oi" /> 365 <field field_ref="qtr_ice_top" name="qtr_ice_top"/> 366 <field field_ref="qtr_ice_bot" name="qtr_ice_bot"/> 367 <field field_ref="qt_ice" name="qt_ice" /> 368 <field field_ref="qsr_ice" name="qsr_ice" /> 369 <field field_ref="qns_ice" name="qns_ice" /> 370 <field field_ref="qemp_ice" name="qemp_ice" /> 371 <field field_ref="albedo" name="albedo" /> 372 <field field_ref="icealb" name="icealb" /> 373 374 <field field_ref="hfxcndtop" name="hfxcndtop" /> 375 <field field_ref="hfxcndbot" name="hfxcndbot" /> 376 <field field_ref="hfxsensib" name="hfxsensib" /> 377 378 <field field_ref="icehc" name="icehc" /> 379 <field field_ref="snwhc" name="snwhc" /> 380 <field field_ref="hfxbog" name="hfxbog" /> 381 <field field_ref="hfxbom" name="hfxbom" /> 382 <field field_ref="hfxsum" name="hfxsum" /> 383 <field field_ref="hfxopw" name="hfxopw" /> 384 <field field_ref="hfxdif" name="hfxdif" /> 385 <field field_ref="hfxsnw" name="hfxsnw" /> 386 <field field_ref="hfxerr" name="hfxerr" /> 387 <field field_ref="hfxthd" name="hfxthd" /> 388 <field field_ref="hfxdyn" name="hfxdyn" /> 389 <field field_ref="hfxres" name="hfxres" /> 390 <field field_ref="hfxsub" name="hfxsub" /> 391 <field field_ref="hfxspr" name="hfxspr" /> 392 393 <!-- 394 <field field_ref="hfxdhc" name="hfxbudget" operation="average" freq_op="1d" > @qt_atm_oi - @qt_oce_ai -@hfxdhc </field> 395 --> 396 397 <!-- salt fluxes --> 398 <field field_ref="sfxice" name="sfxice" /> 399 400 <!-- mass fluxes --> 401 <field field_ref="vfxice" name="vfxice" /> 402 <field field_ref="vfxsnw" name="vfxsnw" /> 403 404 <field field_ref="vfxbom" name="vfxbom" /> 405 <field field_ref="vfxbog" name="vfxbog" /> 406 <field field_ref="vfxsum" name="vfxsum" /> 407 <field field_ref="vfxopw" name="vfxopw" /> 408 <field field_ref="vfxlam" name="vfxlam" /> 409 <field field_ref="vfxsni" name="vfxsni" /> 410 <field field_ref="vfxdyn" name="vfxdyn" /> 411 <field field_ref="vfxres" name="vfxres" /> 412 <field field_ref="vfxpnd" name="vfxpnd" /> 413 <field field_ref="vfxsub" name="vfxsub" /> 414 <field field_ref="vfxsub_err" name="vfxsub_err" /> 415 <field field_ref="vfxthin" name="vfxthin" /> 416 417 <!-- diag error for negative ice volume after advection --> 418 <field field_ref="iceneg_pres" name="sineg_pres" /> 419 <field field_ref="iceneg_volu" name="sineg_volu" /> 420 <field field_ref="iceneg_hfx" name="sineg_hfx" /> 421 </field_group> 422 423 424 <field_group id="myvarICE_cat" grid_ref="grid_T_ncatice" > 425 426 <!-- categories --> 427 <field field_ref="icemask_cat" name="simskcat"/> 428 <field field_ref="snwthic_cat" name="snthicat"/> 429 <field field_ref="iceconc_cat" name="siconcat"/> 430 <field field_ref="icethic_cat" name="sithicat"/> 431 <field field_ref="icesalt_cat" name="sisalcat"/> 432 <field field_ref="icetemp_cat" name="sitemcat"/> 433 <field field_ref="snwtemp_cat" name="sntemcat"/> 434 <field field_ref="icettop_cat" name="sitopcat"/> 435 309 436 </field_group> 310 437 311 <field_group id="ICE_globalbudget" grid_ref="grid_1point" > 438 <field_group id="ICE_globalbudget" grid_ref="grid_scalar" > 439 <!-- global contents --> 312 440 <field field_ref="ibgvol_tot" name="ibgvol_tot" /> 313 441 <field field_ref="sbgvol_tot" name="sbgvol_tot" /> … … 317 445 <field field_ref="sbgheat_tot" name="sbgheat_tot" /> 318 446 319 <field field_ref="ibgvolume" name="ibgvolume" /> 447 <!-- global drifts (conservation checks) --> 448 <field field_ref="ibgvolume" name="ibgvolume" /> 320 449 <field field_ref="ibgsaltco" name="ibgsaltco" /> 321 450 <field field_ref="ibgheatco" name="ibgheatco" /> 322 451 <field field_ref="ibgheatfx" name="ibgheatfx" /> 323 452 453 <!-- global forcings --> 324 454 <field field_ref="ibgfrcvoltop" name="ibgfrcvoltop" /> 325 455 <field field_ref="ibgfrcvolbot" name="ibgfrcvolbot" /> … … 338 468 339 469 <!-- SIMIP monthly scalar variables --> 340 <field_group id="SImon_scalars" grid_ref="grid_ 1point" >470 <field_group id="SImon_scalars" grid_ref="grid_scalar" > 341 471 <!-- Integrated quantities --> 342 472 <field field_ref="NH_iceextt" name="siextentn" /> -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/field_def_nemo-oce.xml
r11764 r12065 108 108 <field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" /> 109 109 <field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D" /> 110 111 <!-- variables available with ln_diacfl=.true. --> 112 <field id="cfl_cu" long_name="u-courant number" unit="#" /> 113 <field id="cfl_cv" long_name="v-courant number" unit="#" /> 114 <field id="cfl_cw" long_name="w-courant number" unit="#" /> 115 110 116 </field_group> <!-- grid_T --> 111 117 … … 358 364 </field_group> 359 365 360 <!-- scalar variables -->361 <field_group id="SBC_0D" grid_ref="grid_1point" >362 </field_group>363 366 364 367 </field_group> <!-- SBC --> … … 515 518 <field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/> 516 519 517 <field_group id="scalar" grid_ref="grid_ T_2D" >520 <field_group id="scalar" grid_ref="grid_scalar" > 518 521 <!-- Need to have a salinity reference climatological file : sali_ref_clim_monthly --> 519 522 <field id="voltot" long_name="global total volume" standard_name="sea_water_volume" unit="m3" /> … … 542 545 </field_group> 543 546 544 <!-- variables available with key_float-->547 <!-- variables available with ln_floats --> 545 548 546 549 <field_group id="floatvar" grid_ref="grid_T_nfloat" operation="instant" > -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/field_def_nemo-pisces.xml
r10416 r12065 290 290 </field_group> 291 291 292 <field_group id="tracer_scalar" grid_ref="grid_ T_2D" >292 <field_group id="tracer_scalar" grid_ref="grid_scalar" > 293 293 <!-- PISCES scalar --> 294 294 <field id="pno3tot" long_name="Global mean nitrate concentration" unit="mol/m3" /> -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/grid_def_nemo.xml
r10226 r12065 7 7 </grid> 8 8 <!-- --> 9 <grid id="grid_T_ 3D_ncatice" >9 <grid id="grid_T_ncatice" > 10 10 <domain id="grid_T" /> 11 11 <axis id="ncatice" /> … … 49 49 </grid> 50 50 <!-- --> 51 <grid id="grid_1point" >52 <domain domain_ref="1point"/>53 </grid>54 <!-- -->55 51 <grid id="grid_T_nfloat" > 56 52 <domain id="grid_T" /> 57 53 <axis id="nfloat" /> 58 54 </grid> 55 <!-- scalars --> 56 <grid id="grid_scalar" > 57 <scalar/> 58 </grid> 59 59 60 60 </grid_definition> -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/namelist_ice_ref
r10611 r12065 56 56 rn_ishlat = 2. ! lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 57 57 ln_landfast_L16 = .false. ! landfast: parameterization from Lemieux 2016 58 ln_landfast_home = .false. ! landfast: parameterization from "home made"59 58 rn_depfra = 0.125 ! fraction of ocean depth that ice must reach to initiate landfast 60 ! recommended range: [0.1 ; 0.25] - L16=0.125 - home=0.15 61 rn_icebfr = 15. ! ln_landfast_L16: maximum bottom stress per unit volume [N/m3] 62 ! ln_landfast_home: maximum bottom stress per unit area of contact [N/m2] 63 ! recommended range: ?? L16=15 - home=10 59 ! recommended range: [0.1 ; 0.25] 60 rn_icebfr = 15. ! maximum bottom stress per unit volume [N/m3] 64 61 rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction [s-1] 65 rn_tensile = 0.2 ! ln_landfast_L16: isotropic tensile strength62 rn_tensile = 0.2 ! isotropic tensile strength [0-0.5??] 66 63 / 67 64 !------------------------------------------------------------------------------ … … 104 101 &namdyn_adv ! Ice advection 105 102 !------------------------------------------------------------------------------ 106 ln_adv_Pra = . false. ! Advection scheme (Prather)107 ln_adv_UMx = . true. ! Advection scheme (Ultimate-Macho)103 ln_adv_Pra = .true. ! Advection scheme (Prather) 104 ln_adv_UMx = .false. ! Advection scheme (Ultimate-Macho) 108 105 nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order) 109 106 / … … 154 151 &namthd_do ! Ice growth in open water 155 152 !------------------------------------------------------------------------------ 156 rn_hinew = 0.1 ! thickness for new ice formation in open water (m), must be larger than rn_h newice153 rn_hinew = 0.1 ! thickness for new ice formation in open water (m), must be larger than rn_himin 157 154 ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind) 158 155 rn_maxfraz = 1.0 ! maximum fraction of frazil ice collecting at the ice base … … 178 175 &namthd_pnd ! Melt ponds 179 176 !------------------------------------------------------------------------------ 180 ln_pnd_H12 = .false. ! activate evolutive melt ponds (from Holland et al 2012) 181 ln_pnd_CST = .false. ! activate constant melt ponds 182 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC 183 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC 184 ln_pnd_alb = .false. ! melt ponds affect albedo or not 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 185 183 / 186 184 !------------------------------------------------------------------------------ … … 189 187 ln_iceini = .true. ! activate ice initialization (T) or not (F) 190 188 ln_iceini_file = .false. ! netcdf file provided for initialization (T) or not (F) 191 rn_thres_sst = 2.0 ! max delta temp. above Tfreeze with initial ice = (sst - tfreeze) 192 rn_hts_ini_n = 0.3 ! initial real snow thickness (m), North 189 rn_thres_sst = 2.0 ! max temp. above Tfreeze with initial ice = (sst - tfreeze) 190 rn_hti_ini_n = 3.0 ! initial ice thickness (m), North 191 rn_hti_ini_s = 1.0 ! " " South 192 rn_hts_ini_n = 0.3 ! initial snow thickness (m), North 193 193 rn_hts_ini_s = 0.3 ! " " South 194 rn_hti_ini_n = 3.0 ! initial real ice thickness (m), North195 rn_hti_ini_s = 1.0 ! " " South196 194 rn_ati_ini_n = 0.9 ! initial ice concentration (-), North 197 195 rn_ati_ini_s = 0.9 ! " " South 198 196 rn_smi_ini_n = 6.3 ! initial ice salinity (g/kg), North 199 197 rn_smi_ini_s = 6.3 ! " " South 200 rn_tmi_ini_n = 270. ! initial ice /snw temperature(K), North198 rn_tmi_ini_n = 270. ! initial ice temperature (K), North 201 199 rn_tmi_ini_s = 270. ! " " South 202 200 rn_tsu_ini_n = 270. ! initial surface temperature (K), North 201 rn_tsu_ini_s = 270. ! " " South 202 rn_tms_ini_n = 270. ! initial snw temperature (K), North 203 rn_tms_ini_s = 270. ! " " South 204 rn_apd_ini_n = 0.2 ! initial pond fraction (-), North 205 rn_apd_ini_s = 0.2 ! " " South 206 rn_hpd_ini_n = 0.05 ! initial pond depth (m), North 207 rn_hpd_ini_s = 0.05 ! " " South 208 ! -- for ln_iceini_file = T 203 209 sn_hti = 'Ice_initialization' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' 204 210 sn_hts = 'Ice_initialization' , -12 ,'hts' , .false. , .true., 'yearly' , '' , '', '' 205 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' , '' , '', '' 206 214 sn_tsu = 'Ice_initialization' , -12 ,'tsu' , .false. , .true., 'yearly' , '' , '', '' 207 sn_tmi = 'Ice_initialization' , -12 ,'tmi' , .false. , .true., 'yearly' , '' , '', '' 208 sn_smi = 'Ice_initialization' , -12 ,'smi' , .false. , .true., 'yearly' , '' , '', '' 215 sn_tms = 'NOT USED' , -12 ,'tms' , .false. , .true., 'yearly' , '' , '', '' 216 ! 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' , '' , '', '' 209 219 cn_dir='./' 210 220 / … … 222 232 &namdia ! Diagnostics 223 233 !------------------------------------------------------------------------------ 224 ln_icediachk = .false. ! check online the heat, mass & salt budgets (T) or not (F) 234 ln_icediachk = .false. ! check online heat, mass & salt budgets 235 ! ! rate of ice spuriously gained/lost at each time step => rn_icechk=1 <=> 1.e-6 m/hour 236 rn_icechk_cel = 100. ! check at each gridcell (1.e-4m/h)=> stops the code if violated (and writes a file) 237 rn_icechk_glo = 1. ! check over the entire ice cover (1.e-6m/h)=> only prints warnings 225 238 ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F) 226 239 ln_icectl = .false. ! ice points output for debug (T or F) -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/namelist_pisces_ref
r10721 r12065 34 34 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 35 35 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 36 sn_patm = 'presatm' , -1 37 sn_atmco2 = 'presatmco2' , -1 36 sn_patm = 'presatm' , -1. , 'patm' , .true. , .true. , 'yearly' , '' , '' , '' 37 sn_atmco2 = 'presatmco2' , -1. , 'xco2' , .true. , .true. , 'yearly' , '' , '' , '' 38 38 cn_dir = './' ! root directory for the location of the dynamical files 39 39 ! … … 141 141 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 142 142 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 143 sn_par = 'par.orca' , 24 143 sn_par = 'par.orca' , 24. , 'fr_par' , .true. , .true. , 'yearly' , '' , '' , '' 144 144 cn_dir = './' ! root directory for the location of the dynamical files 145 145 ln_varpar = .true. ! boolean for PAR variable … … 347 347 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 348 348 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 349 sn_dust = 'dust.orca' , -1 350 sn_solub = 'solubility.orca' , -12 351 sn_riverdic = 'river.orca' , 120 352 sn_riverdoc = 'river.orca' , 120 353 sn_riverdin = 'river.orca' , 120 354 sn_riverdon = 'river.orca' , 120 355 sn_riverdip = 'river.orca' , 120 356 sn_riverdop = 'river.orca' , 120 357 sn_riverdsi = 'river.orca' , 120 358 sn_ndepo = 'ndeposition.orca', -12 359 sn_ironsed = 'bathy.orca' , -12 360 sn_hydrofe = 'hydrofe.orca' , -12 349 sn_dust = 'dust.orca' , -1. , 'dust' , .true. , .true. , 'yearly' , '' , '' , '' 350 sn_solub = 'solubility.orca' , -12. , 'solubility1' , .false. , .true. , 'yearly' , '' , '' , '' 351 sn_riverdic = 'river.orca' , 120. , 'riverdic' , .true. , .true. , 'yearly' , '' , '' , '' 352 sn_riverdoc = 'river.orca' , 120. , 'riverdoc' , .true. , .true. , 'yearly' , '' , '' , '' 353 sn_riverdin = 'river.orca' , 120. , 'riverdin' , .true. , .true. , 'yearly' , '' , '' , '' 354 sn_riverdon = 'river.orca' , 120. , 'riverdon' , .true. , .true. , 'yearly' , '' , '' , '' 355 sn_riverdip = 'river.orca' , 120. , 'riverdip' , .true. , .true. , 'yearly' , '' , '' , '' 356 sn_riverdop = 'river.orca' , 120. , 'riverdop' , .true. , .true. , 'yearly' , '' , '' , '' 357 sn_riverdsi = 'river.orca' , 120. , 'riverdsi' , .true. , .true. , 'yearly' , '' , '' , '' 358 sn_ndepo = 'ndeposition.orca', -12. , 'ndep' , .false. , .true. , 'yearly' , '' , '' , '' 359 sn_ironsed = 'bathy.orca' , -12. , 'bathy' , .false. , .true. , 'yearly' , '' , '' , '' 360 sn_hydrofe = 'hydrofe.orca' , -12. , 'epsdb' , .false. , .true. , 'yearly' , '' , '' , '' 361 361 ! 362 362 cn_dir = './' ! root directory for the location of the dynamical files -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/namelist_ref
r12037 r12065 37 37 cn_exp = "ORCA2" ! experience name 38 38 nn_it000 = 1 ! first time step 39 nn_itend = 5 475 ! last time step (std 5475)39 nn_itend = 5840 ! last time step (std 5840) 40 40 nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 41 41 nn_time0 = 0 ! initial time of day in hhmm … … 50 50 cn_ocerst_indir = "." ! directory from which to read input ocean restarts 51 51 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 52 cn_ocerst_outdir = "." 52 cn_ocerst_outdir = "." ! directory in which to write output ocean restarts 53 53 ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model 54 54 nn_istate = 0 ! output the initial state (1) or not (0) 55 55 ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 56 nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1) 56 nn_stock = 0 ! used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1) 57 ! ! = 0 force to write restart files only at the end of the run 58 ! ! = -1 do not do any restart 57 59 nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 58 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) 59 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 60 nn_write = 0 ! used only if key_iomput is not defined: output frequency (modulo referenced to nn_it000) 61 ! ! = 0 force to write output files only at the end of the run 62 ! ! = -1 do not do any output file 63 ln_mskland = .false. ! mask land points in NetCDF outputs 60 64 ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard 61 65 ln_clobber = .true. ! clobber (overwrite) an existing file … … 70 74 rn_isfhmin = 1.00 ! treshold [m] to discriminate grounding ice from floating ice 71 75 ! 72 rn_rdt = 5 760. ! time step for the dynamics and tracer76 rn_rdt = 5400. ! time step for the dynamics and tracer 73 77 rn_atfp = 0.1 ! asselin time filter parameter 74 78 ! … … 107 111 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 108 112 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 109 sn_tem = 'data_1m_potential_temperature_nomask', -1 110 sn_sal = 'data_1m_salinity_nomask' , -1 113 sn_tem = 'data_1m_potential_temperature_nomask', -1. , 'votemper', .true. , .true. , 'yearly' , '' , '' , '' 114 sn_sal = 'data_1m_salinity_nomask' , -1. , 'vosaline', .true. , .true. , 'yearly' , '' , '' , '' 111 115 / 112 116 !----------------------------------------------------------------------- 113 117 &namwad ! Wetting and Drying (WaD) (default: OFF) 114 118 !----------------------------------------------------------------------- 115 ln_wd_il = .false 119 ln_wd_il = .false. ! T/F activation of iterative limiter 116 120 ln_wd_dl = .false. ! T/F activation of directional limiter 117 121 ln_wd_dl_bc = .false. ! T/F Directional limiteer Baroclinic option … … 163 167 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 164 168 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 165 sn_ucur = 'ucurrent' , -1 166 sn_vcur = 'vcurrent' , -1 169 sn_ucur = 'ucurrent' , -1. ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 170 sn_vcur = 'vcurrent' , -1. ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 167 171 / 168 172 … … 189 193 &namsbc ! Surface Boundary Condition manager (default: NO selection) 190 194 !----------------------------------------------------------------------- 191 nn_fsbc = 5! frequency of SBC module call195 nn_fsbc = 2 ! frequency of SBC module call 192 196 ! ! (control sea-ice & iceberg model call) 193 197 ! Type of air-sea fluxes … … 239 243 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 240 244 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 241 sn_utau = 'utau' , 24 242 sn_vtau = 'vtau' , 24 243 sn_qtot = 'qtot' , 24 244 sn_qsr = 'qsr' , 24 245 sn_emp = 'emp' , 24 245 sn_utau = 'utau' , 24. , 'utau' , .false. , .false., 'yearly' , '' , '' , '' 246 sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false., 'yearly' , '' , '' , '' 247 sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false., 'yearly' , '' , '' , '' 248 sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false., 'yearly' , '' , '' , '' 249 sn_emp = 'emp' , 24. , 'emp' , .false. , .false., 'yearly' , '' , '' , '' 246 250 / 247 251 !----------------------------------------------------------------------- … … 268 272 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 269 273 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 270 sn_wndi = 'u_10.15JUNE2009_fill' , 6 271 sn_wndj = 'v_10.15JUNE2009_fill' , 6 272 sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24 273 sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24 274 sn_tair = 't_10.15JUNE2009_fill' , 6 275 sn_humi = 'q_10.15JUNE2009_fill' , 6 276 sn_prec = 'ncar_precip.15JUNE2009_fill', -1 277 sn_snow = 'ncar_precip.15JUNE2009_fill', -1 278 sn_slp = 'slp.15JUNE2009_fill' , 6 274 sn_wndi = 'u_10.15JUNE2009_fill' , 6. , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , '' 275 sn_wndj = 'v_10.15JUNE2009_fill' , 6. , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , '' 276 sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24. , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 277 sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24. , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 278 sn_tair = 't_10.15JUNE2009_fill' , 6. , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 279 sn_humi = 'q_10.15JUNE2009_fill' , 6. , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 280 sn_prec = 'ncar_precip.15JUNE2009_fill', -1. , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 281 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 282 sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 279 283 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 280 284 / … … 286 290 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 287 291 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 288 289 292 !_____________!__________________________!____________!_____________!______________________!________! 290 293 ! ! description ! multiple ! vector ! vector ! vector ! … … 323 326 sn_rcv_wper = 'none' , 'no' , '' , '' , '' 324 327 sn_rcv_wnum = 'none' , 'no' , '' , '' , '' 325 sn_rcv_w strf= 'none' , 'no' , '' , '' , ''328 sn_rcv_wfreq = 'none' , 'no' , '' , '' , '' 326 329 sn_rcv_wdrag = 'none' , 'no' , '' , '' , '' 327 330 sn_rcv_ts_ice = 'none' , 'no' , '' , '' , '' … … 343 346 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 344 347 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 345 sn_usp = 'sas_grid_U' , 120 346 sn_vsp = 'sas_grid_V' , 120 347 sn_tem = 'sas_grid_T' , 120 348 sn_sal = 'sas_grid_T' , 120 349 sn_ssh = 'sas_grid_T' , 120 350 sn_e3t = 'sas_grid_T' , 120 351 sn_frq = 'sas_grid_T' , 120 348 sn_usp = 'sas_grid_U' , 120. , 'uos' , .true. , .true. , 'yearly' , '' , '' , '' 349 sn_vsp = 'sas_grid_V' , 120. , 'vos' , .true. , .true. , 'yearly' , '' , '' , '' 350 sn_tem = 'sas_grid_T' , 120. , 'sosstsst', .true. , .true. , 'yearly' , '' , '' , '' 351 sn_sal = 'sas_grid_T' , 120. , 'sosaline', .true. , .true. , 'yearly' , '' , '' , '' 352 sn_ssh = 'sas_grid_T' , 120. , 'sossheig', .true. , .true. , 'yearly' , '' , '' , '' 353 sn_e3t = 'sas_grid_T' , 120. , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , '' 354 sn_frq = 'sas_grid_T' , 120. , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , '' 352 355 / 353 356 !----------------------------------------------------------------------- … … 377 380 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 378 381 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 379 sn_chl ='chlorophyll' , -1 382 sn_chl ='chlorophyll' , -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , '' 380 383 / 381 384 !----------------------------------------------------------------------- … … 394 397 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 395 398 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 396 sn_sst = 'sst_data' , 24 397 sn_sss = 'sss_data' , -1 399 sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false., 'yearly' , '' , '' , '' 400 sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' , '' , '' , '' 398 401 / 399 402 !----------------------------------------------------------------------- … … 416 419 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 417 420 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 418 sn_rnf = 'runoff_core_monthly' , -1 419 sn_cnf = 'runoff_core_monthly' , 0 420 sn_s_rnf = 'runoffs' , 24 421 sn_t_rnf = 'runoffs' , 24 422 sn_dep_rnf = 'runoffs' , 0 421 sn_rnf = 'runoff_core_monthly' , -1. , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , '' 422 sn_cnf = 'runoff_core_monthly' , 0. , 'socoefr0', .false. , .true. , 'yearly' , '' , '' , '' 423 sn_s_rnf = 'runoffs' , 24. , 'rosaline', .true. , .true. , 'yearly' , '' , '' , '' 424 sn_t_rnf = 'runoffs' , 24. , 'rotemper', .true. , .true. , 'yearly' , '' , '' , '' 425 sn_dep_rnf = 'runoffs' , 0. , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , '' 423 426 / 424 427 !----------------------------------------------------------------------- … … 433 436 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 434 437 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 435 sn_apr = 'patm' , -1 438 sn_apr = 'patm' , -1. ,'somslpre' , .true. , .true. , 'yearly' , '' , '' , '' 436 439 / 437 440 !----------------------------------------------------------------------- … … 460 463 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 461 464 !* nn_isf = 4 case 462 sn_fwfisf = 'rnfisf' , -12 465 sn_fwfisf = 'rnfisf' , -12. ,'sowflisf' , .false. , .true. , 'yearly' , '' , '' , '' 463 466 !* nn_isf = 3 case 464 sn_rnfisf = 'rnfisf' , -12 467 sn_rnfisf = 'rnfisf' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 465 468 !* nn_isf = 2 and 3 cases 466 sn_depmax_isf ='rnfisf' , -12,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , ''467 sn_depmin_isf ='rnfisf' , -12,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , ''469 sn_depmax_isf ='rnfisf' , -12. ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 470 sn_depmin_isf ='rnfisf' , -12. ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' 468 471 !* nn_isf = 2 case 469 sn_Leff_isf = 'rnfisf' , -12 472 sn_Leff_isf = 'rnfisf' , -12. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 470 473 / 471 474 !----------------------------------------------------------------------- … … 483 486 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 484 487 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 485 sn_cdg = 'sdw_ecwaves_orca2' , 6 486 sn_usd = 'sdw_ecwaves_orca2' , 6 487 sn_vsd = 'sdw_ecwaves_orca2' , 6 488 sn_hsw = 'sdw_ecwaves_orca2' , 6 489 sn_wmp = 'sdw_ecwaves_orca2' , 6 490 sn_wfr = 'sdw_ecwaves_orca2' , 6 491 sn_wnum = 'sdw_ecwaves_orca2' , 6 492 sn_tauwoc = 'sdw_ecwaves_orca2' , 6 493 sn_tauwx = 'sdw_ecwaves_orca2' , 6 494 sn_tauwy = 'sdw_ecwaves_orca2' , 6 488 sn_cdg = 'sdw_ecwaves_orca2' , 6. , 'drag_coeff' , .true. , .true. , 'yearly' , '' , '' , '' 489 sn_usd = 'sdw_ecwaves_orca2' , 6. , 'u_sd2d' , .true. , .true. , 'yearly' , '' , '' , '' 490 sn_vsd = 'sdw_ecwaves_orca2' , 6. , 'v_sd2d' , .true. , .true. , 'yearly' , '' , '' , '' 491 sn_hsw = 'sdw_ecwaves_orca2' , 6. , 'hs' , .true. , .true. , 'yearly' , '' , '' , '' 492 sn_wmp = 'sdw_ecwaves_orca2' , 6. , 'wmp' , .true. , .true. , 'yearly' , '' , '' , '' 493 sn_wfr = 'sdw_ecwaves_orca2' , 6. , 'wfr' , .true. , .true. , 'yearly' , '' , '' , '' 494 sn_wnum = 'sdw_ecwaves_orca2' , 6. , 'wave_num' , .true. , .true. , 'yearly' , '' , '' , '' 495 sn_tauwoc = 'sdw_ecwaves_orca2' , 6. , 'wave_stress', .true. , .true. , 'yearly' , '' , '' , '' 496 sn_tauwx = 'sdw_ecwaves_orca2' , 6. , 'wave_stress', .true. , .true. , 'yearly' , '' , '' , '' 497 sn_tauwy = 'sdw_ecwaves_orca2' , 6. , 'wave_stress', .true. , .true. , 'yearly' , '' , '' , '' 495 498 / 496 499 !----------------------------------------------------------------------- … … 500 503 ! 501 504 ! ! diagnostics: 502 ln_bergdia = .true. ! Calculate budgets503 nn_verbose_level = 0 ! Turn on more verbose output if level > 0504 nn_verbose_write = 15 ! Timesteps between verbose messages505 nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage505 ln_bergdia = .true. ! Calculate budgets 506 nn_verbose_level = 0 ! Turn on more verbose output if level > 0 507 nn_verbose_write = 15 ! Timesteps between verbose messages 508 nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage 506 509 ! 507 510 ! ! iceberg setting: 508 ! ! Initial mass required for an iceberg of each class511 ! ! Initial mass required for an iceberg of each class 509 512 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 510 ! ! Proportion of calving mass to apportion to each class513 ! ! Proportion of calving mass to apportion to each class 511 514 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 512 ! ! Ratio between effective and real iceberg mass (non-dim)513 ! ! i.e. number of icebergs represented at a point514 rn_mass_scaling = 2000 , 200, 50, 20, 10, 5, 2, 1, 1, 1515 ! thickness of newly calved bergs (m)515 ! ! Ratio between effective and real iceberg mass (non-dim) 516 ! ! i.e. number of icebergs represented at a point 517 rn_mass_scaling = 2000., 200., 50., 20., 10., 5., 2., 1., 1., 1. 518 ! thickness of newly calved bergs (m) 516 519 rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 517 520 ! 518 rn_rho_bergs = 850. ! Density of icebergs519 rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs520 ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics521 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits522 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)523 ln_passive_mode = .false. ! iceberg - ocean decoupling524 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no)525 ! ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)521 rn_rho_bergs = 850. ! Density of icebergs 522 rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs 523 ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics 524 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits 525 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 526 ln_passive_mode = .false. ! iceberg - ocean decoupling 527 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no) 528 ! ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 526 529 rn_test_box = 108.0, 116.0, -66.0, -58.0 527 ln_use_calving = .false. ! Use calving data even when nn_test_icebergs > 0528 rn_speed_limit = 0. ! CFL speed limit for a berg530 ln_use_calving = .false. ! Use calving data even when nn_test_icebergs > 0 531 rn_speed_limit = 0. ! CFL speed limit for a berg 529 532 530 533 cn_dir = './' ! root directory for the calving data location … … 532 535 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 533 536 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 534 sn_icb = 'calving' , -1 537 sn_icb = 'calving' , -1. ,'calvingmask', .true. , .true. , 'yearly' , '' , '' , '' 535 538 / 536 539 … … 602 605 nn_ice_dta = 0 ! = 0, bdy data are equal to the initial state 603 606 ! ! = 1, bdy data are read in 'bdydata .nc' files 604 rn_ice_tem = 270. ! si3 only: arbitrary temperature of incoming sea ice605 rn_ice_sal = 10. ! si3 only: -- salinity --606 rn_ice_age = 30. ! si3 only: -- age --607 607 ! 608 608 ln_tra_dmp =.false. ! open boudaries conditions for tracers 609 609 ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities 610 610 rn_time_dmp = 1. ! Damping time scale in days 611 rn_time_dmp_out = 1. 611 rn_time_dmp_out = 1. ! Outflow damping time scale 612 612 nn_rimwidth = 10 ! width of the relaxation zone 613 613 ln_vol = .false. ! total volume correction (see nn_volctl parameter) 614 614 nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero 615 nb_jpk_bdy = -1 ! number of levels in the bdy data (set < 0 if consistent with planned run)616 615 / 617 616 !----------------------------------------------------------------------- 618 617 &nambdy_dta ! open boundaries - external data (see nam_bdy) 619 618 !----------------------------------------------------------------------- 620 ln_full_vel = .false. ! ??? 621 619 ln_zinterp = .false. ! T if a vertical interpolation is required. Variables gdep[tuv] and e3[tuv] must exist in the file 620 ! ! automatically defined to T if the number of vertical levels in bdy dta /= jpk 621 ln_full_vel = .false. ! T if [uv]3d are "full" velocities and not only its baroclinic components 622 ! ! in this case, baroclinic and barotropic velocities will be recomputed -> [uv]2d not needed 623 ! 622 624 cn_dir = 'bdydta/' ! root directory for the BDY data location 623 625 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 624 626 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 625 627 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 626 bn_ssh = 'amm12_bdyT_u2d' , 24 627 bn_u2d = 'amm12_bdyU_u2d' , 24 628 bn_v2d = 'amm12_bdyV_u2d' , 24 629 bn_u3d = 'amm12_bdyU_u3d' , 24 630 bn_v3d = 'amm12_bdyV_u3d' , 24 631 bn_tem = 'amm12_bdyT_tra' , 24 632 bn_sal = 'amm12_bdyT_tra' , 24 628 bn_ssh = 'amm12_bdyT_u2d' , 24. , 'sossheig', .true. , .false., 'daily' , '' , '' , '' 629 bn_u2d = 'amm12_bdyU_u2d' , 24. , 'vobtcrtx', .true. , .false., 'daily' , '' , '' , '' 630 bn_v2d = 'amm12_bdyV_u2d' , 24. , 'vobtcrty', .true. , .false., 'daily' , '' , '' , '' 631 bn_u3d = 'amm12_bdyU_u3d' , 24. , 'vozocrtx', .true. , .false., 'daily' , '' , '' , '' 632 bn_v3d = 'amm12_bdyV_u3d' , 24. , 'vomecrty', .true. , .false., 'daily' , '' , '' , '' 633 bn_tem = 'amm12_bdyT_tra' , 24. , 'votemper', .true. , .false., 'daily' , '' , '' , '' 634 bn_sal = 'amm12_bdyT_tra' , 24. , 'vosaline', .true. , .false., 'daily' , '' , '' , '' 633 635 !* for si3 634 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra', .true. , .false., 'daily' , '' , '' , '' 635 ! bn_h_i = 'amm12_bdyT_ice' , 24 , 'iicethic', .true. , .false., 'daily' , '' , '' , '' 636 ! bn_h_s = 'amm12_bdyT_ice' , 24 , 'isnowthi', .true. , .false., 'daily' , '' , '' , '' 636 bn_a_i = 'amm12_bdyT_ice' , 24. , 'siconc' , .true. , .false., 'daily' , '' , '' , '' 637 bn_h_i = 'amm12_bdyT_ice' , 24. , 'sithic' , .true. , .false., 'daily' , '' , '' , '' 638 bn_h_s = 'amm12_bdyT_ice' , 24. , 'snthic' , .true. , .false., 'daily' , '' , '' , '' 639 bn_t_i = 'NOT USED' , 24. , 'sitemp' , .true. , .false., 'daily' , '' , '' , '' 640 bn_t_s = 'NOT USED' , 24. , 'sntemp' , .true. , .false., 'daily' , '' , '' , '' 641 bn_tsu = 'NOT USED' , 24. , 'sittop' , .true. , .false., 'daily' , '' , '' , '' 642 bn_s_i = 'NOT USED' , 24. , 'sisalt' , .true. , .false., 'daily' , '' , '' , '' 643 ! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd) 644 bn_aip = 'NOT USED' , 24. , 'siapnd' , .true. , .false., 'daily' , '' , '' , '' 645 bn_hip = 'NOT USED' , 24. , 'sihpnd' , .true. , .false., 'daily' , '' , '' , '' 646 ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 647 rn_ice_tem = 270. ! arbitrary temperature of incoming sea ice 648 rn_ice_sal = 10. ! -- salinity -- 649 rn_ice_age = 30. ! -- age -- 650 rn_ice_apnd = 0.2 ! -- pond fraction = a_ip/a_i -- 651 rn_ice_hpnd = 0.05 ! -- pond depth -- 637 652 / 638 653 !----------------------------------------------------------------------- … … 874 889 ln_dynvor_eeT = .false. ! energy conserving scheme (een using e3t) 875 890 ln_dynvor_een = .false. ! energy & enstrophy scheme 876 nn_een_e3f = 1! =0 e3f = mi(mj(e3t))/4891 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 877 892 ! ! =1 e3f = mi(mj(e3t))/mi(mj( tmask)) 878 893 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) ==>>> PLEASE DO NOT ACTIVATE … … 947 962 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 948 963 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 949 sn_tem = 'dyna_grid_T' , 120 950 sn_sal = 'dyna_grid_T' , 120 951 sn_mld = 'dyna_grid_T' , 120 952 sn_emp = 'dyna_grid_T' , 120 953 sn_fmf = 'dyna_grid_T' , 120 954 sn_ice = 'dyna_grid_T' , 120 955 sn_qsr = 'dyna_grid_T' , 120 956 sn_wnd = 'dyna_grid_T' , 120 957 sn_uwd = 'dyna_grid_U' , 120 958 sn_vwd = 'dyna_grid_V' , 120 959 sn_wwd = 'dyna_grid_W' , 120 960 sn_avt = 'dyna_grid_W' , 120 961 sn_ubl = 'dyna_grid_U' , 120 962 sn_vbl = 'dyna_grid_V' , 120 964 sn_tem = 'dyna_grid_T' , 120. , 'votemper' , .true. , .true. , 'yearly' , '' , '' , '' 965 sn_sal = 'dyna_grid_T' , 120. , 'vosaline' , .true. , .true. , 'yearly' , '' , '' , '' 966 sn_mld = 'dyna_grid_T' , 120. , 'somixhgt' , .true. , .true. , 'yearly' , '' , '' , '' 967 sn_emp = 'dyna_grid_T' , 120. , 'sowaflup' , .true. , .true. , 'yearly' , '' , '' , '' 968 sn_fmf = 'dyna_grid_T' , 120. , 'iowaflup' , .true. , .true. , 'yearly' , '' , '' , '' 969 sn_ice = 'dyna_grid_T' , 120. , 'soicecov' , .true. , .true. , 'yearly' , '' , '' , '' 970 sn_qsr = 'dyna_grid_T' , 120. , 'soshfldo' , .true. , .true. , 'yearly' , '' , '' , '' 971 sn_wnd = 'dyna_grid_T' , 120. , 'sowindsp' , .true. , .true. , 'yearly' , '' , '' , '' 972 sn_uwd = 'dyna_grid_U' , 120. , 'uocetr_eff', .true. , .true. , 'yearly' , '' , '' , '' 973 sn_vwd = 'dyna_grid_V' , 120. , 'vocetr_eff', .true. , .true. , 'yearly' , '' , '' , '' 974 sn_wwd = 'dyna_grid_W' , 120. , 'wocetr_eff', .true. , .true. , 'yearly' , '' , '' , '' 975 sn_avt = 'dyna_grid_W' , 120. , 'voddmavs' , .true. , .true. , 'yearly' , '' , '' , '' 976 sn_ubl = 'dyna_grid_U' , 120. , 'sobblcox' , .true. , .true. , 'yearly' , '' , '' , '' 977 sn_vbl = 'dyna_grid_V' , 120. , 'sobblcoy' , .true. , .true. , 'yearly' , '' , '' , '' 963 978 / 964 979 … … 1105 1120 !! namdiu Cool skin and warm layer models (default: OFF) 1106 1121 !! namdiu Cool skin and warm layer models (default: OFF) 1107 !! namflo float parameters ( "key_float")1108 !! nam_diaharm Harmonic analysis of tidal constituents ( "key_diaharm")1109 !! nam dct transports through some sections ("key_diadct")1122 !! namflo float parameters (default: OFF) 1123 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF) 1124 !! nam_diadct transports through some sections (default: OFF) 1110 1125 !! nam_diatmb Top Middle Bottom Output (default: OFF) 1111 1126 !! nam_dia25h 25h Mean Output (default: OFF) … … 1142 1157 &namhsb ! Heat and salt budgets (default: OFF) 1143 1158 !----------------------------------------------------------------------- 1144 ln_diahsb = .false. ! checkthe heat and salt budgets (T) or not (F)1159 ln_diahsb = .false. ! output the heat and salt budgets (T) or not (F) 1145 1160 / 1146 1161 !----------------------------------------------------------------------- … … 1151 1166 / 1152 1167 !----------------------------------------------------------------------- 1153 &namflo ! float parameters ("key_float") 1154 !----------------------------------------------------------------------- 1155 jpnfl = 1 ! total number of floats during the run 1156 jpnnewflo = 0 ! number of floats for the restart 1157 ln_rstflo = .false. ! float restart (T) or not (F) 1158 nn_writefl = 75 ! frequency of writing in float output file 1159 nn_stockfl = 5475 ! frequency of creation of the float restart file 1160 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 1161 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 1162 ! ! or computed with Blanke' scheme (F) 1163 ln_ariane = .true. ! Input with Ariane tool convention(T) 1164 ln_flo_ascii = .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 1165 / 1166 !----------------------------------------------------------------------- 1167 &nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm") 1168 !----------------------------------------------------------------------- 1169 nit000_han = 1 ! First time step used for harmonic analysis 1170 nitend_han = 75 ! Last time step used for harmonic analysis 1171 nstep_han = 15 ! Time step frequency for harmonic analysis 1172 tname(1) = 'M2' ! Name of tidal constituents 1173 tname(2) = 'K1' 1174 / 1175 !----------------------------------------------------------------------- 1176 &namdct ! transports through some sections ("key_diadct") 1177 !----------------------------------------------------------------------- 1178 nn_dct = 15 ! time step frequency for transports computing 1179 nn_dctwri = 15 ! time step frequency for transports writing 1180 nn_secdebug = 112 ! 0 : no section to debug 1181 ! ! -1 : debug all section 1182 ! ! 0 < n : debug section number n 1168 &namflo ! float parameters (default: OFF) 1169 !----------------------------------------------------------------------- 1170 ln_floats = .false. ! activate floats or not 1171 jpnfl = 1 ! total number of floats during the run 1172 jpnnewflo = 0 ! number of floats for the restart 1173 ln_rstflo = .false. ! float restart (T) or not (F) 1174 nn_writefl = 75 ! frequency of writing in float output file 1175 nn_stockfl = 5475 ! frequency of creation of the float restart file 1176 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 1177 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 1178 ! ! or computed with Blanke' scheme (F) 1179 ln_ariane = .true. ! Input with Ariane tool convention(T) 1180 ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 1181 / 1182 !----------------------------------------------------------------------- 1183 &nam_diaharm ! Harmonic analysis of tidal constituents (default: OFF) 1184 !----------------------------------------------------------------------- 1185 ln_diaharm = .false. ! Choose tidal harmonic output or not 1186 nit000_han = 1 ! First time step used for harmonic analysis 1187 nitend_han = 75 ! Last time step used for harmonic analysis 1188 nstep_han = 15 ! Time step frequency for harmonic analysis 1189 tname(1) = 'M2' ! Name of tidal constituents 1190 tname(2) = 'K1' ! --- 1191 / 1192 !----------------------------------------------------------------------- 1193 &nam_diadct ! transports through some sections (default: OFF) 1194 !----------------------------------------------------------------------- 1195 ln_diadct = .false. ! Calculate transport thru sections or not 1196 nn_dct = 15 ! time step frequency for transports computing 1197 nn_dctwri = 15 ! time step frequency for transports writing 1198 nn_secdebug = 112 ! 0 : no section to debug 1199 ! ! -1 : debug all section 1200 ! ! 0 < n : debug section number n 1183 1201 / 1184 1202 !----------------------------------------------------------------------- … … 1231 1249 ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there 1232 1250 ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs 1251 ln_bound_reject = .false. ! Logical to remove obs near boundaries in LAMs. 1233 1252 ln_sla_fp_indegs = .true. ! Logical for SLA: T=> averaging footprint is in degrees, F=> in metres 1234 1253 ln_sst_fp_indegs = .true. ! Logical for SST: T=> averaging footprint is in degrees, F=> in metres … … 1297 1316 &nammpp ! Massively Parallel Processing ("key_mpp_mpi") 1298 1317 !----------------------------------------------------------------------- 1299 cn_mpi_send = 'I' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, 1300 ! ! buffer blocking send or immediate non-blocking sends, resp. 1301 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 1318 ln_listonly = .false. ! do nothing else than listing the best domain decompositions (with land domains suppression) 1319 ! ! if T: the largest number of cores tested is defined by max(mppsize, jpni*jpnj) 1302 1320 ln_nnogather = .true. ! activate code to avoid mpi_allgather use at the northfold 1303 jpni = 0 ! jpni number of processors following i (set automatically if < 1)1304 jpnj = 0 ! jpnj number of processors following j (set automatically if < 1)1321 jpni = 0 ! number of processors following i (set automatically if < 1), see also ln_listonly = T 1322 jpnj = 0 ! number of processors following j (set automatically if < 1), see also ln_listonly = T 1305 1323 / 1306 1324 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/namelist_top_ref
r10375 r12065 59 59 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 60 60 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 61 sn_trcdta(1) = 'data_TRC_nomask' , -12 61 sn_trcdta(1) = 'data_TRC_nomask' , -12. , 'TRC' , .false. , .true. , 'yearly' , '' , '' , '' 62 62 ! 63 63 cn_dir = './' ! root directory for the location of the data files
Note: See TracChangeset
for help on using the changeset viewer.