Changeset 2061 for CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM
- Timestamp:
- 04/26/13 13:01:38 (10 years ago)
- Location:
- CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM
- Files:
-
- 3 added
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/iodef.xml
r1437 r2061 1 1 <?xml version="1.0"?> 2 2 <simulation> 3 4 <context id="nemo">5 3 4 <context id="nemo" time_origin="1900-01-01 00:00:00" > 5 6 6 <!-- $id$ --> 7 7 … … 12 12 ============================================================================================================ 13 13 --> 14 15 <field_definition level="0" prec="4" operation="ave(X)" enabled=".TRUE."> <!-- time step automaticaly defined --> 16 17 <!-- T grid --> 18 19 <group id="grid_T" axis_ref="none" grid_ref="grid_T"> 20 <field id="toce" description="temperature" unit="degC" axis_ref="deptht" /> 21 <field id="soce" description="salinity" unit="psu" axis_ref="deptht" /> 22 <field id="sst" description="sea surface temperature" unit="degC" /> 23 <field id="sst2" description="square of sea surface temperature" unit="degC2" /> 24 <field id="|sstgrad|" description="module of sst gradient" unit="degC/m" /> 25 <field id="|sstgrad|2" description="square of module of sst gradient" unit="degC2/m2" /> 26 <field id="sss" description="sea surface salinity" unit="psu" /> 27 <field id="sss2" description="square of sea surface salinity" unit="psu2" /> 28 <field id="ssh" description="sea surface height" unit="m" /> 29 <field id="ssh2" description="square of sea surface height" unit="m2" /> 30 <field id="mldkz5" description="mixing layer depth (Turbocline)" unit="m" /> 31 <field id="mldr10_1" description="Mixed Layer Depth 0.01 ref.10m" unit="m" /> 32 <field id="rhop" description="potential density (sigma0)" unit="kg/m3" axis_ref="deptht" /> 33 <!-- next variables available with key_diahth --> 34 <field id="mlddzt" description="Thermocline Depth (max dT/dz)" unit="m" /> 35 <field id="mldr10_3" description="Mixed Layer Depth dr=0.03 (ref.10m)" unit="m" /> 36 <field id="mldr0_1" description="Mixed Layer Depth dr=0.01 (ref.surf)" unit="m" /> 37 <field id="mldr0_3" description="Mixed Layer Depth dr=0.03 (ref.surf)" unit="m" /> 38 <field id="mld|dt|" description="Mixed Layer Depth |dt|=0.2 (ref.10m)" unit="m" /> 39 <field id="topthdep" description="Top of the thermocline dt=-0.2 (ref.10m)" unit="m" /> 40 <field id="pycndep" description="Pycnocline depth dr~dt=-0.2 (ref.10m)" unit="m" /> 41 <field id="BLT" description="Barrier Layer Thickness" unit="m" /> 42 <field id="tinv" description="Max of vertical invertion of temperature" unit="degC" /> 43 <field id="depti" description="Depth of max. vert. inv. of temperature" unit="m" /> 44 <field id="20d" description="Depth of 20C isotherm" unit="m" /> 45 <field id="28d" description="Depth of 28C isotherm" unit="m" /> 46 <field id="hc300" description="Heat content 300 m" unit="W" /> 47 <!-- variables available with key_diaar5 --> 48 <field id="botpres" description="Pressure at sea floor" unit="dbar" /> 49 <field id="cellthc" description="Cell thickness" unit="m" axis_ref="deptht" /> 50 </group> 51 52 <!-- SBC --> 53 54 <group id="SBC" axis_ref="none" grid_ref="grid_T" > <!-- time step automaticaly defined based on nn_fsbc --> 55 56 <field id="emp" description="Net Upward Water Flux" unit="kg/m2/s" /> 57 <field id="emps" description="concentration/dilution water flux" unit="kg/m2/s" /> 58 <field id="snowpre" description="Snow precipitation" unit="kg/m2/s" /> 59 <field id="runoffs" description="River Runoffs" unit="kg/m2/s" /> 60 61 <field id="qns+qsr" description="Net Downward Heat Flux" unit="W/m2" /> 62 <field id="qns" description="non solar Downward Heat Flux" unit="W/m2" /> 63 <field id="qsr" description="Shortwave Radiation" unit="W/m2" /> 64 <field id="qsr3d" description="Shortwave Radiation 3D distribution" axis_ref="deptht" unit="W/m2" /> 65 <field id="qrp" description="Surface Heat Flux: Damping" unit="W/m2" /> 66 <field id="erp" description="Surface Water Flux: Damping" unit="kg/m2/s" /> 67 <field id="taum" description="wind stress module" unit="N/m2" /> 68 <field id="wspd" description="Wind speed module at 10 m" unit="m/s" /> 69 70 <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 71 <field id="qns_oce" description="Non solar Downward Heat Flux over open ocean" unit="W/m2" /> 72 <field id="qlw_oce" description="Longwave Downward Heat Flux over open ocean" unit="W/m2" /> 73 <field id="qsb_oce" description="Sensible Downward Heat Flux over open ocean" unit="W/m2" /> 74 <field id="qla_oce" description="Latent Downward Heat Flux over open ocean" unit="W/m2" /> 75 <field id="taum_oce" description="wind stress module over open ocean" unit="N/m2" /> 76 77 <field id="ice_cover" description="Ice fraction" unit="1" /> 78 79 <field id="ioceflxb" description="Oceanic flux at the ice base" unit="W/m2" /> 80 <field id="qsr_ai_cea" description="Air-Ice downward solar heat flux (cell average)" unit="W/m2" /> 81 <field id="qns_ai_cea" description="Air-Ice downward non-solar heat flux (cell average)" unit="W/m2" /> 82 <field id="qla_ai_cea" description="Air-Ice downward Latent heat flux (cell average)" unit="W/m2" /> 83 84 <field id="qsr_io_cea" description="Ice-Oce downward solar heat flux (cell average)" unit="W/m2" /> 85 <field id="qns_io_cea" description="Ice-Oce downward non-solar heat flux (cell average)" unit="W/m2" /> 86 87 <field id="snowthic_cea" description="Snow thickness (cell average)" unit="m" /> 88 <field id="icethic_cea" description="Ice thickness (cell average)" unit="m" /> 89 <field id="iceprod_cea" description="Ice production (cell average)" unit="m/s" /> 90 <field id="sim" description="Sea Ice Plus Surface Snow Amount" unit="kg/m2" /> 91 92 <field id="ice_pres" description="Ice presence" unit="-" /> 93 <field id="ist_cea" description="Ice surface temperature (cell average)" unit="degC" /> 94 <field id="ist_ipa" description="Ice surface temperature (ice presence average)" unit="degC" /> 95 <field id="uice_ipa" description="Ice velocity along i-axis at I-point (ice presence average)" unit="m/s" /> 96 <field id="vice_ipa" description="Ice velocity along j-axis at I-point (ice presence average)" unit="m/s" /> 97 98 <field id="utau_ice" description="Wind stress along i-axis over the ice at i-point" unit="N/m2" /> 99 <field id="vtau_ice" description="Wind stress along j-axis over the ice at i-point" unit="N/m2" /> 100 101 <field id="u_imasstr" description="Sea-ice mass transport along i-axis" unit="kg/s" /> 102 <field id="v_imasstr" description="Sea-ice mass transport along j-axis" unit="kg/s" /> 103 104 <!-- available key_coupled --> 105 <field id="snow_ao_cea" description="Snow over ice-free ocean (cell average)" unit="kg/m2/s" /> 106 <field id="snow_ai_cea" description="Snow over sea-ice (cell average)" unit="kg/m2/s" /> 107 <field id="subl_ai_cea" description="Sublimation over sea-ice (cell average)" unit="kg/m2/s" /> 108 <field id="icealb_cea" description="Ice albedo (cell average)" unit="1" /> 109 <field id="calving" description="Calving" unit="kg/m2/s" /> 110 <!-- available if key_coupled + conservative method --> 111 <field id="rain" description="Liquid precipitation" unit="kg/m2/s" /> 112 <field id="evap_ao_cea" description="Evaporation over ice-free ocean (cell average)" unit="kg/m2/s" /> 113 <!-- variables available with key_diaar5 --> 114 <field id="isnwmlt_cea" description="Snow over Ice melting (cell average)" unit="kg/m2/s" /> 115 <field id="fsal_virt_cea" description="Virtual salt flux due to ice formation (cell average)" unit="kg/m2/s" /> 116 <field id="fsal_real_cea" description="Real salt flux due to ice formation (cell average)" unit="kg/m2/s" /> 117 <field id="hflx_rain_cea" description="heat flux due to rainfall" unit="W/m2" /> 118 <field id="hflx_evap_cea" description="heat flux due to evaporation" unit="W/m2" /> 119 <field id="hflx_snow_cea" description="heat flux due to snow falling over ice-free ocean" unit="W/m2" /> 120 <field id="hflx_ice_cea" description="heat flux due to ice thermodynamics" unit="W/m2" /> 121 <field id="hflx_rnf_cea" description="heat flux due to runoffs" unit="W/m2" /> 122 <field id="hflx_cal_cea" description="heat flux due to calving" unit="W/m2" /> 123 <field id="bicemel_cea" description="Rate of Melt at Sea Ice Base (cell average)" unit="kg/m2/s" /> 124 <field id="licepro_cea" description="Lateral Sea Ice Growth Rate (cell average)" unit="kg/m2/s" /> 125 <field id="snowmel_cea" description="Snow Melt Rate (cell average)" unit="kg/m2/s" /> 126 <field id="sntoice_cea" description="Snow-Ice Formation Rate (cell average)" unit="kg/m2/s" /> 127 <field id="ticemel_cea" description="Rate of Melt at Upper Surface of Sea Ice (cell average)" unit="kg/m2/s" /> 128 129 </group> 130 131 <!-- U grid --> 132 133 <group id="grid_U" axis_ref="depthu" grid_ref="grid_U"> 134 <field id="utau" description="Wind Stress along i-axis" unit="N/m2" axis_ref="none" /> 135 <field id="uoce" description="ocean current along i-axis" unit="m/s" /> 136 <field id="uoce_eff" description="Effective ocean current along i-axis" unit="m/s" /> 137 <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 138 <field id="uoce_eiv" description="EIV ocean current along i-axis" unit="m/s" /> 139 <!-- uoce_eiv: available with key_trabbl_adv --> 140 <field id="uoce_bbl" description="BBL ocean current along i-axis" unit="m/s" /> 141 <!-- variables available with key_diaar5 --> 142 <field id="u_masstr" description="ocean eulerian mass transport along i-axis" unit="kg/s" /> 143 <field id="u_heattr" description="ocean eulerian heat transport along i-axis" unit="W" axis_ref="none" /> 144 <field id="ueiv_heattr" description="ocean bolus heat transport along i-axis" unit="W" axis_ref="none" /> 145 <field id="udiff_heattr" description="ocean diffusion heat transport along i-axis" unit="W" axis_ref="none" /> 146 </group> 147 148 <!-- V grid --> 149 150 <group id="grid_V" axis_ref="depthv" grid_ref="grid_V"> 151 <field id="vtau" description="Wind Stress along j-axis" unit="N/m2" axis_ref="none" /> 152 <field id="voce" description="ocean current along j-axis" unit="m/s" /> 153 <field id="voce_eff" description="Effective ocean current along j-axis" unit="m/s" /> 154 <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 155 <field id="voce_eiv" description="EIV ocean current along j-axis" unit="m/s" /> 156 <!-- voce_eiv: available with key_trabbl_adv --> 157 <field id="voce_bbl" description="BBL ocean current along j-axis" unit="m/s" /> 158 <!-- variables available with key_diaar5 --> 159 <field id="v_masstr" description="ocean eulerian mass transport along j-axis" unit="kg/s" /> 160 <field id="v_heattr" description="ocean eulerian heat transport along j-axis" unit="W" axis_ref="none" /> 161 <field id="veiv_heattr" description="ocean bolus heat transport along j-axis" unit="W" axis_ref="none" /> 162 <field id="vdiff_heattr" description="ocean diffusion heat transport along j-axis" unit="W" axis_ref="none" /> 163 </group> 164 165 <!-- W grid --> 166 167 <group id="grid_W" axis_ref="depthw" grid_ref="grid_W"> 168 <field id="woce" description="ocean vertical velocity" unit="m/s" /> 169 <field id="woce_eff" description="effective ocean vertical velocity" unit="m/s" /> 170 <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 171 <field id="woce_eiv" description="EIV ocean vertical velocity" unit="m/s" /> 172 <!-- woce_eiv: available with key_trabbl_adv --> 173 <field id="woce_bbl" description="BBL ocean vertical velocity" unit="m/s" /> 174 <field id="avt" description="vertical eddy diffusivity" unit="m2/s" /> 175 <field id="avm" description="vertical eddy viscosity" unit="m2/s" /> 176 <!-- avs: available with key_zdfddm --> 177 <field id="avs" description="salt vertical eddy diffusivity" unit="m2/s" /> 178 <!-- avt_evd and avm_evd: available with ln_zdfevd --> 179 <field id="avt_evd" description="enhanced vertical diffusivity" unit="m2/s" /> 180 <field id="avm_evd" description="enhanced vertical viscosity" unit="m2/s" /> 181 <!-- aht2d and aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d --> 182 <field id="aht2d" description="lateral eddy diffusivity" unit="m2/s" axis_ref="none" /> 183 <field id="aht2d_eiv" description="EIV lateral eddy diffusivity" unit="m2/s" axis_ref="none" /> 184 <!-- avt_tide: available with key_zdftmx --> 185 <field id="av_tide" description="tidal vertical diffusivity" unit="m2/s" /> 186 <!-- variables available with key_diaar5 --> 187 <field id="w_masstr" description="vertical mass trasport" unit="kg/s" /> 188 <field id="w_masstr2" description="square of vertical mass trasport" unit="kg2/s2" /> 189 </group> 190 191 <!-- scalar --> 192 193 <!-- variables available with key_diaar5 --> 194 <group id="scalar" axis_ref="none" grid_ref="scalarpoint" zoom_ref="1point" > 195 <field id="voltot" description="global mean volume" unit="m3" /> 196 <field id="sshtot" description="global mean ssh" unit="m" /> 197 <field id="sshsteric" description="global mean ssh steric" unit="m" /> 198 <field id="sshthster" description="global mean ssh thermosteric" unit="m" /> 199 <field id="masstot" description="global mean mass" unit="kg" /> 200 <field id="temptot" description="global mean temperature" unit="degC" /> 201 <field id="saltot" description="global mean salinity" unit="psu" /> 202 <field id="fram_trans" description="Sea Ice Mass Transport Through Fram Strait" unit="kg/s" /> 203 </group> 204 205 <!-- ptrc on T grid --> 206 207 <group id="ptrc_T" axis_ref="deptht" grid_ref="grid_T"> 208 <field id="DIC" description="Dissolved inorganic Concentration" unit="molC/L" /> 209 <field id="Alkalini" description="Total Alkalinity Concentration" unit="eq/L" /> 210 <field id="O2" description="Oxygen Concentration" unit="molO2/L" /> 211 <field id="CaCO3" description="Calcite Concentration" unit="molC/L" /> 212 <field id="PO4" description="Phosphate Concentration" unit="molC/L" /> 213 <field id="POC" description="Small organic carbon Concentration" unit="molC/L" /> 214 <field id="Si" description="Silicate Concentration" unit="molSi/L" /> 215 <field id="PHY" description="Nanophytoplankton Concentration" unit="molC/L" /> 216 <field id="ZOO" description="Microzooplankton Concentration" unit="molC/L" /> 217 <field id="DOC" description="Dissolved organic Concentration" unit="molC/L" /> 218 <field id="PHY2" description="Diatoms Concentration" unit="molC/L" /> 219 <field id="ZOO2" description="Mesozooplankton Concentration" unit="molC/L" /> 220 <field id="BSi" description="Diatoms Silicate Concentration" unit="molC/L" /> 221 <field id="Fer" description="Dissolved Iron Concentration" unit="molFe/L" /> 222 <field id="BFe" description="Big iron particles Concentration" unit="molFe/L" /> 223 <field id="GOC" description="Big organic carbon Concentration" unit="molC/L" /> 224 <field id="SFe" description="Small iron particles Concentration" unit="molFe/L" /> 225 <field id="DFe" description="Diatoms iron Concentration" unit="molFe/L" /> 226 <field id="DSi" description="Sinking biogenic Silicate Concentration" unit="molC/L" /> 227 <field id="NFe" description="Nano iron Concentration" unit="molC/L" /> 228 <field id="NCHL" description="Nano chlorophyl Concentration" unit="gChl/L" /> 229 <field id="DCHL" description="Diatoms chlorophyl Concentration" unit="gChl/L" /> 230 <field id="NO3" description="Nitrates Concentration" unit="molC/L" /> 231 <field id="NH4" description="Ammonium Concentration" unit="molC/L" /> 232 </group> 233 234 <!-- diad on T grid : : variables available with key_trc_diaadd && key_trc_dia3d --> 235 236 <group id="diad_T" axis_ref="none" grid_ref="grid_T"> 237 <field id="PH" description="PH" unit="-" axis_ref="deptht" /> 238 <field id="CO3" description="Bicarbonates" unit="mol/L" axis_ref="deptht" /> 239 <field id="CO3sat" description="CO3 saturation" unit="mol/L" axis_ref="deptht" /> 240 <field id="PAR" description="Photosynthetically Available Radiation" unit="W/m2" axis_ref="deptht" /> 241 <field id="PPPHY" description="Primary production of nanophyto" unit="molC/m3/s" axis_ref="deptht" /> 242 <field id="PPPHY2" description="Primary production of diatoms" unit="molC/m3/s" axis_ref="deptht" /> 243 <field id="PPNEWN" description="New Primary production of nanophyto" unit="molC/m3/s" axis_ref="deptht" /> 244 <field id="PPNEWD" description="New Primary production of diatoms" unit="molC/m3/s" axis_ref="deptht" /> 245 <field id="PBSi" description="Primary production of Si diatoms" unit="molSi/m3/s" axis_ref="deptht" /> 246 <field id="PFeN" description="Primary production of nano iron" unit="molFe/m3/s" axis_ref="deptht" /> 247 <field id="PFeD" description="Primary production of diatoms iron" unit="molFe/m3/s" axis_ref="deptht" /> 248 <field id="PCAL" description="Calcite production" unit="molC/m3/s" axis_ref="deptht" /> 249 <field id="DCAL" description="Calcite dissolution" unit="molC/m3/s" axis_ref="deptht" /> 250 <field id="GRAZ" description="Grazing by zooplankton" unit="molC/m3/s" axis_ref="deptht" /> 251 <field id="Nfix" description="Nitrogen fixation at surface" unit="molN/m2/s" /> 252 <field id="EPC100" description="Export of carbon particles at 100 m" unit="molC/m2/s" /> 253 <field id="EPFE100" description="Export of biogenic iron at 100 m" unit="molFe/m2/s" /> 254 <field id="EPSI100" description="Export of Silicate at 100 m" unit="molSi/m2/s" /> 255 <field id="EPCAL100" description="Export of Calcite at 100 m" unit="molC/m2/s" /> 256 <field id="Cflx" description="DIC flux" unit="molC/m2/s" /> 257 <field id="Oflx" description="Oxygen flux" unit="molC/m2/s" /> 258 <field id="Kg" description="Gas transfer" unit="molC/m2/s/uatm" /> 259 <field id="Dpco2" description="Delta CO2" unit="uatm" /> 260 <field id="Dpo2" description="Delta O2" unit="uatm" /> 261 <field id="Heup" description="Euphotic layer depth" unit="m" /> 262 <field id="Irondep" description="Iron deposition" unit="molFe/m2/s" /> 263 </group> 264 265 <!-- dbio_T on T grid : variables available with key_diaar5 --> 266 <group id="dbio_T" axis_ref="none" grid_ref="grid_T"> 267 <field id="PHYT" description="Total phytoplankton Concentration" unit="molC/L" axis_ref="deptht" /> 268 <field id="ZOOT" description="Total zooplankton Concentration" unit="molC/L" axis_ref="deptht" /> 269 <field id="CHLT" description="Total chlorophyll Concentration" unit="gChl/L" axis_ref="deptht" /> 270 <field id="POCT" description="Total carbon particles Concentration" unit="molC/L" axis_ref="deptht" /> 271 <field id="PFET" description="Total biogenic Iron Concentration" unit="molFe/L" axis_ref="deptht" /> 272 <field id="TPP" description="Total Primary production of phyto" unit="molC/m3/s" axis_ref="deptht" /> 273 <field id="TPNEW" description="New Primary production of phyto" unit="molC/m3/s" axis_ref="deptht" /> 274 <field id="TPBFE" description="Total biogenic iron production" unit="molFe/m3/s" axis_ref="deptht" /> 275 <field id="EXPC" description="Export of carbon" unit="molC/m2/s" axis_ref="deptht" /> 276 <field id="EXPFE" description="Export of biogenic iron" unit="molFe/m2/s" axis_ref="deptht" /> 277 <field id="EXPSI" description="Export of Silicate" unit="molSi/m2/s" axis_ref="deptht" /> 278 <field id="EXPCAL" description="Export of Calcite" unit="molC/m2/s" axis_ref="deptht" /> 279 <field id="PHSFC" description="PH at surface" unit="-" /> 280 <field id="INTDIC" description="DIC content" unit="kg/m2" /> 281 <field id="O2MIN" description="Oxygen minimum concentration" unit="molC/L" /> 282 <field id="ZO2MIN" description="Depth of oxygen minimum concentration" unit="m" /> 283 <field id="INTNFIX" description="Nitrogen fixation rate : vert. integrated" unit="molN/m2/s" /> 284 <field id="INTPPPHY" description="Vertically integrated primary production by nanophy" unit="molC/m2/s" /> 285 <field id="INTPPPHY2" description="Vertically integrated primary production by diatom" unit="molC/m2/s" /> 286 <field id="INTPP" description="Vertically integrated primary production by phyto" unit="molC/m2/s" /> 287 <field id="INTPNEW" description="Vertically integrated new primary production" unit="molC/m2/s" /> 288 <field id="INTPBFE" description="Vertically integrated of biogenic iron production" unit="molFe/m2/s" /> 289 <field id="INTPBSI" description="Vertically integrated of biogenic Si production" unit="molSi/m2/s" /> 290 <field id="INTPCAL" description="Vertically integrated of calcite production" unit="molSi/m2/s" /> 291 </group> 292 293 </field_definition> 294 14 <field_definition src="./field_def.xml"/> 295 15 <!-- 296 16 ============================================================================================================ 297 = =298 17 = output files definition = 299 = = 300 = http://www2-pcmdi.llnl.gov/svn/repository/cmor/trunk/Tables/CMIP5_Oclim = 301 = http://www2-pcmdi.llnl.gov/svn/repository/cmor/trunk/Tables/CMIP5_Omon = 302 = http://www2-pcmdi.llnl.gov/svn/repository/cmor/trunk/Tables/CMIP5_Oimon = 303 = http://www.clivar.org/organization/wgomd/references/WGOMD_CMIP5_ocean_fields.June2009.pdf = 304 = = 18 = Define your own files = 19 = put the variables you want... = 305 20 ============================================================================================================ 306 21 --> 307 22 308 <file_definition output_level="<OUTPUT_LEVEL>" > 309 <!-- 310 +++++++++++++++++++++++++++++++++++++++++++++++ daily ++++++++++++++++++++++++++++++++++++++++++++++++++ 311 --> 312 <group id="1d" output_freq="86400" enabled="<1D_ENABLE>"> <!-- 1d files --> 313 <!-- 314 .............................................. grid T ................................................. 315 --> 316 <file id="ORCA2_1d_grid_T" description="table 2.2" > <!-- grid T Table 2.2 --> 317 <field ref="sst" name="sosstsst" description="sea_surface_temperature" level="2" /> 318 <field ref="sst2" name="tossq" description="square_of_sea_surface_temperature" level="2" /> 319 <field ref="mldkz5" name="omldamax" description="ocean_mixed_layer_thickness_defined_by_mixing_scheme" level="2" operation="t_max(X)" /> 320 <!-- IPSL variables --> 321 <field ref="sss" name="sosaline" description="sea_surface_salinity" level="3" /> 322 <field ref="ssh" name="sossheig" description="sea_surface_height_above_geoid" level="3" /> 323 <field ref="toce" name="votemper" description="sea_water_potential_temperature" level="3" /> 324 <field ref="soce" name="vosaline" description="sea_water_salinity" level="3" /> 325 </file> 326 <!-- 327 .............................................. grid U ................................................. 328 --> 329 <file id="ORCA2_1d_grid_U" description="ocean U grid variables : table 2.3" > <!-- grid U --> 330 331 <!-- table 2.3 --> 332 <field ref="uoce" name="vozocrtx" description="sea_water_x_velocity" level="3" /> 333 334 </file> 335 336 <!-- 337 .............................................. grid V ................................................. 338 --> 339 <file id="ORCA2_1d_grid_V" description="ocean V grid variables : table 2.3" > <!-- grid V --> 340 341 <!-- table 2.3 --> 342 <field ref="voce" name="vomecrty" description="sea_water_y_velocity" level="3" /> 343 344 </file> 23 <file_definition type="multiple_file" sync_freq="1d" min_digits="4"> 24 25 <file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE."/> <!-- 1h files --> 26 <file_group id="2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files --> 27 <file_group id="3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files --> 28 <file_group id="4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files --> 29 <file_group id="6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files --> 30 31 <file_group id="1d" output_freq="1d" output_level="10" enabled=".TRUE."> <!-- 1d files --> 32 33 <file id="1d_grid_T" name="auto" description="ocean T grid variables" > 34 <field field_ref="sst" name="sosstsst" /> 35 <field field_ref="sss" name="sosaline" /> 36 <field field_ref="ssh" name="sossheig" /> 37 </file> 38 39 <file id="1d_grid_U" name="auto" description="ocean U grid variables" > 40 <field field_ref="suoce" name="vozocrtx" /> 41 </file> 42 43 <file id="1d_grid_V" name="auto" description="ocean V grid variables" > 44 <field field_ref="svoce" name="vomecrty" /> 45 </file> 46 47 </file_group> 48 <file_group id="3d" output_freq="3d" output_level="10" enabled=".TRUE."/> <!-- 3d files --> 49 <file_group id="5d" output_freq="5d" output_level="10" enabled=".TRUE."> <!-- 5d files --> 50 51 <file id="5d_grid_T" name="auto" description="ocean T grid variables" > 52 <field field_ref="toce" name="votemper" /> 53 <field field_ref="soce" name="vosaline" /> 54 <field field_ref="sst" name="sosstsst" /> 55 <field field_ref="sss" name="sosaline" /> 56 <field field_ref="ssh" name="sossheig" /> 57 <field field_ref="empmr" name="sowaflup" /> 58 <field field_ref="qsr" name="soshfldo" /> 59 <field field_ref="saltflx" name="sosfldow" /> 60 <field field_ref="qt" name="sohefldo" /> 61 <field field_ref="mldr10_1" name="somxl010" /> 62 <field field_ref="mldkz5" name="somixhgt" /> 63 <field field_ref="ice_cover" name="soicecov" /> 64 <field field_ref="wspd" name="sowindsp" /> 65 <field field_ref="erp" name="sowafldp" /> 66 <field field_ref="ist_ipa" name="soicetem" /> 67 <field field_ref="icealb_cea" name="soicealb" /> 68 </file> 69 70 <file id="5d_grid_U" name="auto" description="ocean U grid variables" > 71 <field field_ref="uoce" name="vozocrtx" /> 72 <field field_ref="uoce_eiv" name="vozoeivu" /> 73 <field field_ref="utau" name="sozotaux" /> 74 </file> 75 76 <file id="5d_grid_V" name="auto" description="ocean V grid variables" > 77 <field field_ref="voce" name="vomecrty" /> 78 <field field_ref="voce_eiv" name="vomeeivv" /> 79 <field field_ref="vtau" name="sometauy" /> 80 </file> 81 82 <file id="5d_grid_W" name="auto" description="ocean W grid variables" > 83 <field field_ref="woce" name="vovecrtz" /> 84 <field field_ref="avt" name="votkeavt" /> 85 <field field_ref="avs" name="voddmavs" /> 86 <field field_ref="aht2d_eiv" name="soleaeiw" /> 87 </file> 88 89 <file id="5d_icemod" name="auto" description="ice variables" > 90 <field field_ref="ice_pres" /> 91 <field field_ref="snowthic_cea" name="isnowthi" /> 92 <field field_ref="icethic_cea" name="iicethic" /> 93 <field field_ref="iceprod_cea" name="iiceprod" /> 94 <field field_ref="ist_ipa" name="iicetemp" /> 95 <field field_ref="ioceflxb" name="ioceflxb" /> 96 <field field_ref="uice_ipa" name="iicevelu" /> 97 <field field_ref="vice_ipa" name="iicevelv" /> 98 <field field_ref="utau_ice" name="iicestru" /> 99 <field field_ref="vtau_ice" name="iicestrv" /> 100 <field field_ref="qsr_io_cea" name="iicesflx" /> 101 <field field_ref="qns_io_cea" name="iicenflx" /> 102 <field field_ref="snowpre" name="isnowpre" /> 103 </file> 104 105 </file_group> 106 107 <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files --> 108 109 <file id="1m_ptrc_T" name="auto" description="pisces sms variables" > 110 <field field_ref="DIC" /> 111 <field field_ref="Alkalini" /> 112 <field field_ref="O2" /> 113 <field field_ref="PO4" /> 114 <field field_ref="Si" /> 115 <field field_ref="Fer" /> 116 <field field_ref="NCHL" /> 117 <field field_ref="DCHL" /> 118 <field field_ref="NO3" /> 119 </file> 120 121 <file id="1m_diad_T" name="auto" description="additional pisces diagnostics" > 122 <field field_ref="Cflx" /> 123 <field field_ref="Dpco2" /> 124 </file> 125 126 </file_group> 127 <file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files --> 128 <file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files --> 129 <file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files --> 130 <file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files --> 131 132 <file_group id="1y" output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files --> 133 134 <file id="1y_ptrc_T" name="auto" description="pisces sms variables" > 135 <field field_ref="DIC" /> 136 <field field_ref="Alkalini" /> 137 <field field_ref="O2" /> 138 <field field_ref="CaCO3" /> 139 <field field_ref="PO4" /> 140 <field field_ref="POC" /> 141 <field field_ref="Si" /> 142 <field field_ref="PHY" /> 143 <field field_ref="ZOO" /> 144 <field field_ref="DOC" /> 145 <field field_ref="PHY2" /> 146 <field field_ref="ZOO2" /> 147 <field field_ref="DSi" /> 148 <field field_ref="Fer" /> 149 <field field_ref="BFe" /> 150 <field field_ref="GOC" /> 151 <field field_ref="SFe" /> 152 <field field_ref="DFe" /> 153 <field field_ref="GSi" /> 154 <field field_ref="NFe" /> 155 <field field_ref="NCHL" /> 156 <field field_ref="DCHL" /> 157 <field field_ref="NO3" /> 158 <field field_ref="NH4" /> 159 </file> 160 161 <file id="1y_diad_T" name="auto" description="additional pisces diagnostics" > 162 <field field_ref="PH" /> 163 <field field_ref="CO3" /> 164 <field field_ref="CO3sat" /> 165 <field field_ref="PAR" /> 166 <field field_ref="PPPHY" /> 167 <field field_ref="PPPHY2" /> 168 <field field_ref="PPNEWN" /> 169 <field field_ref="PPNEWD" /> 170 <field field_ref="PBSi" /> 171 <field field_ref="PFeN" /> 172 <field field_ref="PFeD" /> 173 <field field_ref="xfracal" /> 174 <field field_ref="PCAL" /> 175 <field field_ref="DCAL" /> 176 <field field_ref="GRAZ1" /> 177 <field field_ref="GRAZ2" /> 178 <field field_ref="EPC100" /> 179 <field field_ref="EPFE100" /> 180 <field field_ref="EPSI100" /> 181 <field field_ref="EPCAL100" /> 182 <field field_ref="Cflx" /> 183 <field field_ref="Oflx" /> 184 <field field_ref="Kg" /> 185 <field field_ref="Dpco2" /> 186 <field field_ref="Dpo2" /> 187 <field field_ref="Heup" /> 188 <field field_ref="Irondep" /> 189 <field field_ref="Ironsed" /> 190 <field field_ref="Ironice" /> 191 <field field_ref="Nfix" /> 192 <field field_ref="MuN" /> 193 <field field_ref="MuD" /> 194 <field field_ref="LNnut" /> 195 <field field_ref="LDnut" /> 196 <field field_ref="LNFe" /> 197 <field field_ref="LDFe" /> 198 <field field_ref="LNlight" /> 199 <field field_ref="LDlight" /> 200 <field field_ref="pdust" /> 201 <field field_ref="Fe2" /> 202 <field field_ref="Fe3" /> 203 <field field_ref="FeL1" /> 204 <field field_ref="FeL2" /> 205 <field field_ref="FeP" /> 206 <field field_ref="TL1" /> 207 <field field_ref="TL2" /> 208 <field field_ref="Sdenit" /> 209 <field field_ref="Totlig" /> 210 </file> 211 212 </file_group> 213 <file_group id="2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files --> 214 <file_group id="5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files --> 215 <file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files --> 216 217 </file_definition> 218 219 <!-- 220 ============================================================================================================ 221 = grid definition = = DO NOT CHANGE = 222 ============================================================================================================ 223 --> 224 225 <axis_definition> 226 <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" /> 227 <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" /> 228 <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" /> 229 <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" /> 230 <axis id="nfloat" long_name="Float number" unit="-" /> 231 </axis_definition> 232 233 <domain_definition src="./domain_def.xml"/> 234 235 <grid_definition> 236 <grid id="grid_T_2D" domain_ref="grid_T"/> 237 <grid id="grid_T_3D" domain_ref="grid_T" axis_ref="deptht"/> 238 <grid id="grid_U_2D" domain_ref="grid_U"/> 239 <grid id="grid_U_3D" domain_ref="grid_U" axis_ref="depthu"/> 240 <grid id="grid_V_2D" domain_ref="grid_V"/> 241 <grid id="grid_V_3D" domain_ref="grid_V" axis_ref="depthv"/> 242 <grid id="grid_W_2D" domain_ref="grid_W"/> 243 <grid id="grid_W_3D" domain_ref="grid_W" axis_ref="depthw"/> 244 </grid_definition> 345 245 346 <!-- 347 .............................................. icemod ................................................. 348 --> 349 <file id="ORCA2_1d_icemod" description=" ice variables" > <!-- icemod --> 350 351 <!-- table 2.2 --> 352 <field ref="uice_ipa" name="iicevelu" description="Ice velocity along i-axis at I-point (ice presence average)" level="2" /> 353 <field ref="vice_ipa" name="iicevelv" description="Ice velocity along j-axis at I-point (ice presence average)" level="2" /> 354 <field ref="ice_cover" name="ileadfra" description="sea_ice_area_fraction" level="2" /> 355 <field ref="icethic_cea" name="iicethic" description="sea_ice_thickness" level="2" /> 356 </file> 357 358 </group> 359 <!-- 360 ++++++++++++++++++++++++++++++++++++++++++++++ monthly +++++++++++++++++++++++++++++++++++++++++++++++++ 361 --> 362 <group id="1m" output_freq="-1" enabled="<1M_ENABLE>"> <!-- real monthly files --> 363 <!-- 364 .............................................. grid T ................................................. 365 --> 366 <file id="ORCA2_1m_grid_T" description="ocean T grid variables table 2.2 2.5 2.6 2.7"> <!-- grid T --> 367 368 <!-- table 2.2 --> 369 <field ref="botpres" name="pbo" description="sea_water_pressure_at_sea_floor" /> 370 <!-- pso : sea_water_pressure_at_sea_water_surface = 0 --> 371 <field ref="ssh" name="sossheig" description="sea_surface_height_above_geoid" /> 372 <field ref="ssh2" name="zossq" description="square_of_sea_surface_height_above_geoid" level="2" /> 373 <!-- masscello : sea_water_mass_per_unit_area = cellthc*rau0 no time changes --> 374 <field ref="cellthc" name="thkcello" description="cell_thickness" /> <!-- no time changes --> 375 <field ref="toce" name="votemper" description="sea_water_potential_temperature" /> 376 <field ref="sst" name="sosstsst" description="sea_surface_temperature" level="1" /> 377 <field ref="sst2" name="tossq" description="square_of_sea_surface_temperature" level="2" /> 378 <field ref="soce" name="vosaline" description="sea_water_salinity" /> 379 <field ref="sss" name="sosaline" description="sea_surface_salinity" level="1" /> 380 <field ref="rhop" name="rhopoto" description="sea_water_potential_density" level="2" /> 381 <!-- no agessc : sea_water_age_since_surface_contact --> 382 <!-- no cfc11 : moles_per_unit_mass_of_cfc11_in_sea_water --> 383 <!-- msftbarot : ocean_barotropic_mass_streamfunction : offline --> 384 <!-- mlotst : ocean_mixed_layer_thickness_defined_by_sigma_t : must be done offline --> 385 <!-- mlotstsq : square_of_ocean_mixed_layer_thickness_defined_by_sigma_t : must be done offline --> 386 <field ref="mldkz5" name='omlmax' description="ocean_mixed_layer_thickness_defined_by_mixing_scheme" level="2" operation="t_max(X)" /> 387 388 <!-- table 2.5 --> 389 <field ref="rain" name="pr" description="rainfall_flux" level="1" /> 390 <field ref="snow_ao_cea" name="prsn" description="snowfall_flux" level="1" /> 391 <field ref="evap_ao_cea" name="evs" description="water_evaporation_flux" level="1" /> 392 <field ref="runoffs" name="friver" description="water_flux_into_sea_water_from_rivers" level="1" /> 393 <field ref="calving" name="ficeberg" description="water_flux_into_sea_water_from_icebergs" level="1" /> 394 <field ref="isnwmlt_cea" name="fsitherm" description="water_flux_into_sea_water_due_to_sea_ice_thermodynamics" level="1" /> 395 <field ref="emp" name="sowaflup" description="water_flux_into_sea_water" level="1" /> 396 <!-- wfonocorr : water_flux_into_sea_water_without_flux_correction : emp - erp --> 397 <field ref="erp" name="wfcorr" description="water_flux_correction" level="1" /> <!-- usually = 0 --> 398 399 <!-- table2.6 --> 400 <!-- vsfpr : virtual_salt_flux_into_sea_water_due_to_rainfall = 0 --> 401 <!-- vsfevap : virtual_salt_flux_into_sea_water_due_to_evaporation = 0 --> 402 <!-- vsfriver : virtual_salt_flux_into_sea_water_from_rivers = 0 --> 403 <field ref="fsal_virt_cea" name="vsfsit" description="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics" level="1" /> 404 <!-- vsf : virtual_salt_flux_into_sea_water = fsal_virtual + fsal_real --> 405 <!-- wfcorr : virtual_salt_flux_correction = 0 --> 406 <field ref="fsal_virt_cea" name="sfdsi" description="downward_sea_ice_basal_salt_flux" level="1" /> 407 <!-- sfriver : salt_flux_into_sea_water_from_rivers = 0 --> 408 409 <!-- table 2.7 --> 410 <!-- hfgeou : upward_geothermal_heat_flux_at_sea_floor : cte, see nambbc and trabbc.F90 --> 411 <field ref="hflx_rain_cea" name="hfrainds" description="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water" level="1" /> 412 <field ref="hflx_evap_cea" name="hfevapds" description="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water" level="1" /> 413 <field ref="hflx_rnf_cea" name="hfrunoffds" description="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water" level="1" /> 414 <field ref="hflx_snow_cea" name="hfsnthermds" description="heat_flux_into_sea_water_due_to_snow_thermodynamics" level="1" /> 415 <field ref="hflx_ice_cea" name="hfsithermds" description="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics" level="1" /> 416 <field ref="hflx_cal_cea" name="hfibthermds" description="heat_flux_into_sea_water_due_to_iceberg_thermodynamics" level="1" /> 417 <!-- rlds : surface_net_downward_longwave_flux : not available --> 418 <!-- hfls : surface_downward_latent_heat_flux : not available --> 419 <!-- hfss : surface_downward_sensible_heat_flux: not available --> 420 <field ref="qns" name="nshfls" description="surface_net_downward_non_solar_flux" level="1" /> 421 <field ref="qsr" name="soshfldo" description="surface_net_downward_shortwave_flux" level="1" /> 422 <field ref="qsr3d" name="rsds" description="downwelling_shortwave_flux_in_sea_water" level="1" /> 423 <field ref="qrp" name="hfcorr" description="heat_flux_correction" level="1" /> 424 425 <!-- IPSL variables --> 426 <field ref="qns+qsr" name="sohefldo" description="Net Downward Heat Flux" /> 427 <field ref="20d" name="so20chgt" description="Depth of 20C isotherm" /> 428 <field ref="28d" name="so28chgt" description="Depth of 28C isotherm" /> 429 <field ref="hc300" name="sohtc300" description="Heat content 300 m" /> 430 <field ref="mldr10_1" name="somxl010" description="Mixed Layer Depth 0.01 ref.10m" /> 431 <field ref="mldr10_3" name="somx3010" description="Mixed Layer Depth dr=0.03 (ref.10m)" /> 432 <field ref="mldkz5" name="somixhgt" description="mixing layer depth (Turbocline)" /> 433 <field ref="emps" name="sowaflcd" description="concentration/dilution water flux" /> 434 <field ref="taum" name="sotaumod" description="Wind stress module" /> 435 <field ref="wspd" name="sowindsp" description="Wind speed module at 10 m" /> 436 <field ref="ice_cover" name="soicecov" description="sea_ice_area_fraction" /> 437 <field ref="BLT" name="blt" description="Barrier Layer Thickness" level="3" /> 438 <field ref="mld|dt|" name="mld_dt02" description="Mixed Layer Depth |dt|=0.2 (ref.10m)" level="3" /> 439 440 </file> 441 <!-- 442 .............................................. grid U ................................................. 443 --> 444 <file id="ORCA2_1m_grid_U" description="ocean U grid variables : table 2.3 2.8" > <!-- grid U --> 445 446 <!-- table 2.3 --> 447 <field ref="uoce" name="vozocrtx" description="sea_water_x_velocity" /> 448 <field ref="u_masstr" name="umo" description="ocean_mass_x_transport" level="1" /> 449 <field ref="u_heattr" name="hfx" description="ocean_heat_x_transport" level="1" /> 450 <field ref="ueiv_heattr" name="hfxba" description="ocean_heat_x_transport_due_to_bolus_advection" level="2" /> 451 <field ref="udiff_heattr" name="hfxdiff" description="ocean_heat_x_transport_due_to_diffusion" level="2" /> 452 453 <!-- table 2.8 --> 454 <field ref="utau" name="sozotaux" description="surface_downward_x_stress" level="1" /> 455 <!-- tauucorr : surface_downward_x_stress_correction = 0 --> 456 457 <!-- IPSL variables --> 458 <field ref="uoce_eiv" name="vozoeivu" description="EIV ocean current along i-axis" /> 459 460 </file> 461 462 <!-- 463 .............................................. grid V ................................................. 464 --> 465 <file id="ORCA2_1m_grid_V" description="ocean V grid variables : table 2.3 2.8" > <!-- grid V --> 466 467 <!-- table 2.3" --> 468 <field ref="voce" name="vomecrty" description="sea_water_y_velocity" /> 469 <field ref="v_masstr" name="vmo" description="ocean_mass_y_transport" level="1" /> 470 <field ref="v_heattr" name="hfy" description="ocean_heat_y_transport" level="1" /> 471 <field ref="veiv_heattr" name="hfyba" description="ocean_heat_y_transport_due_to_bolus_advection" level="2" /> 472 <field ref="vdiff_heattr" name="hfydiff" description="ocean_heat_y_transport_due_to_diffusion" level="2" /> 473 474 <!-- table 2.8 --> 475 <field ref="vtau" name="sometauy" description="surface_downward_y_stress" level="1" /> 476 <!-- tauvcorr : surface_downward_y_stress_correction = 0 --> 477 478 <!-- IPSL variables --> 479 <field ref="voce_eiv" name="vomeeivv" description="EIV ocean current along j-axis" /> 480 481 </file> 482 483 <!-- 484 .............................................. grid W ................................................. 485 --> 486 <file id="ORCA2_1m_grid_W" description="ocean W grid variables table 2.3 2.9 2.10" > <!-- grid W --> 487 488 <!-- table 2.3 --> 489 <field ref="w_masstr" name="wmo" description="upward_ocean_mass_transport" /> 490 <field ref="w_masstr2" name="wmosq" description="square_pf_upward_ocean_mass_transport" /> 491 492 <!-- table 2.9 --> 493 <field ref="avt" name="votkeavt" description="ocean_vertical_heat_diffusivity" level="2" /> 494 <field ref="avs" name="voddmavs" description="ocean_vertical_salt_diffusivity" level="2" /> 495 <!-- difvtrbo : ocean_vertical_tracer_diffusivity_due_to_background : cte with time, see namelist parameters nn_avb and nn_havtb --> 496 <field ref="av_tide" name="difvtrto" description="ocean_vertical_tracer_diffusivity_due_to_tides" level="2" /> 497 <!-- tnpeo : tendency_of_ocean_potential_energy_content : not available --> 498 <!-- tnpeot : tendency_of_ocean_potential_energy_content_due_to_tides : not available --> 499 <!-- tnpeotb : tendency_of_ocean_potential_energy_content_due_to_background : not available --> 500 <field ref="avm" name="difvmo" description="ocean_vertical_momentum_diffusivity" level="2" /> 501 <!-- difvmbo : ocean_vertical_momentum_diffusivity_due_to_background : cte with time, see namelist parameters nn_avb --> 502 <field ref="av_tide" name="difvmto" description="ocean_vertical_momentum_diffusivity_due_to_tides" level="2" /> <!-- same as tracer --> 503 <!-- difvmfdo : ocean_vertical_momentum_diffusivity_due_to_form_drag : ??? --> 504 <!-- dispkevfo : ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction : not available --> 505 506 <!-- table 2.10 --> 507 <!-- if ln_traldf_lap = .true. --> 508 <field ref="aht2d_eiv" name="soleaeiw" description="ocean_tracer_bolus_laplacian_diffusivity" level="2" /> 509 <!-- diftrelo : ocean_tracer_epineutral_laplacian_diffusivity : cte with time, see ln_traldf_iso --> 510 <!-- diftrxylo : ocean_tracer_xy_laplacian_diffusivity : cte with time --> 511 <field ref="aht2d" name="soleahtw" description="ocean_tracer_xy_laplacian_diffusivity" level="2" /> 512 <!-- if ln_traldf_bilap = .true. --> 513 <!-- field ref="diftrbbo" name="aht2d_eiv" description="ocean_tracer_bolus_biharmonic_diffusivity" level="2" /--> 514 <!-- diftrebo : ocean_tracer_epineutral_biharmonic_diffusivity : cte with time, see ln_traldf_iso --> 515 <!-- diftrxybo : ocean_tracer_xy_biharmonic_diffusivity : cte with time --> 516 <!-- tnkebto : tendency_of_ocean_eddy_kinetic_energy_content_due_to_bolus_transport : not available --> 517 <!-- difmxylo : ocean_momentum_xy_laplacian_diffusivity : cte with time, see ln_dynldf_lap --> 518 <!-- difmxybo : ocean_momentum_xy_biharmonic_diffusivity : cte with time, see ln_dynldf_bilap --> 519 <!-- dispkexyfo : ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction : not available --> 520 521 <!-- IPSL variables --> 522 <field ref="woce" name="vovecrtz" description="ocean vertical velocity" /> 523 <field ref="woce_eiv" name="voveeivw" description="EIV ocean vertical velocity" /> 524 525 526 </file> 527 528 <!-- 529 .............................................. icemod ................................................. 530 --> 531 <file id="ORCA2_1m_icemod" description=" ice variables" > <!-- icemod --> 532 <!-- table 2.2 --> 533 <field ref="ice_pres" /> 534 <field ref="ice_cover" name="soicecov" description="sea_ice_area_fraction" /> 535 <field ref="icethic_cea" name="iicethic" description="sea_ice_thickness" /> 536 <field ref="subl_ai_cea" name="evap" description="water_evaporation_flux" /> 537 <field ref="snowthic_cea" name="isnowthi" description="surface_snow_thickness" /> 538 <field ref="sim" name="sim" description="sea_ice_and_surface_snow_amount" /> 539 <!-- snc : surface_snow_area_fraction : same as sic --> 540 <field ref="icealb_cea" name="ialb" description="Bare Sea Ice Albedo" /> 541 <!-- ssi : Sea Ice Salinity : cte = sice = 6, expect in baltic sea = 2 --> 542 <field ref="ist_cea" name="tsice" description="Surface Temperature of Sea Ice" /> 543 <!-- tsnint : Temperature at Interface Between Sea Ice and Snow : not available --> 544 <!-- pr : Surface Rainfall Rate into the Sea Ice Portion of the Grid Cell = 0 --> 545 <field ref="snow_ai_cea" name="prsn" description="Surface Snowfall Rate into the Sea Ice Portion of the Grid Cell" /> 546 <!-- ageice : Age of Sea Ice : not available ??? --> 547 <!-- grFrazil : Frazil Sea Ice Growth (Leads) Rate ??? --> 548 <!-- grCongel : Congelation Sea Ice Growth Rate ??? --> 549 <field ref="licepro_cea" name="grLateral" description="Lateral Sea Ice Growth Rate" /> 550 <field ref="sntoice_cea" name="snoToIce" description="Snow-Ice Formation Rate" /> 551 <field ref="snowmel_cea" name="snomelt" description="Snow Melt Rate" /> 552 <field ref="ticemel_cea" name="tmelt" description="Rate of Melt at Upper Surface of Sea Ice" /> 553 <field ref="bicemel_cea" name="bmelt" description="Rate of Melt at Sea Ice Base" /> 554 <!-- hcice : Sea Ice Total Heat Content ??? --> 555 <field ref="qsr_ai_cea" name="rsdssi" description="surface_downwelling_shortwave_flux_in_air" /> 556 <!-- rsussi : surface_upwelling_shortwave_flux_in_air ??? --> 557 <!-- rldssi : surface_downwelling_longwave_flux_in_air ??? --> 558 <!-- rlussi : surface_upwelling_longwave_flux_in_air ??? --> 559 <!-- hfssi : surface_upward_sensible_heat_flux ??? --> 560 <!-- hflssi : surface_upward_latent_heat_flux ??? --> 561 <field ref="subl_ai_cea" name="sblsi" description="surface_snow_and_ice_sublimation_flux" /> 562 <field ref="u_imasstr" name="transix" description="Eastward Sea Ice Transport" /> 563 <field ref="v_imasstr" name="transiy" description="Northward Sea Ice Transport" /> 564 <field ref="fram_trans" name="transifs" description="Sea Ice Mass Transport Through Fram Strait" /> 565 <field ref="utau_ice" name="strairx" description="Eastward Atmospheric Stress On Sea Ice" /> 566 <field ref="vtau_ice" name="strairy" description="Northward Atmospheric Stress On Sea Ice" /> 567 <!-- strocnx : Eastward Ocean Stress On Sea Ice ??? --> 568 <!-- strocny : Northward Ocean Stress On Sea Ice ??? --> 569 <!-- streng : Compressive Sea Ice Strength ??? --> 570 <!-- divice : Strain Rate Divergence of Sea Ice ??? --> 571 <!-- shrice : Strain Rate Shear of Sea Ice ??? --> 572 <!-- ridgice : Sea Ice Ridging Rate ??? --> 573 <!-- IPSL variables --> 574 <field ref="ist_ipa" name="iicetemp" description="Ice surface temperature (ice presence average)" /> 575 <field ref="uice_ipa" name="iicevelu" description="Ice velocity along i-axis at I-point (ice presence average)" /> 576 <field ref="vice_ipa" name="iicevelv" description="Ice velocity along j-axis at I-point (ice presence average)" /> 577 <field ref="ice_cover" name="ileadfra" description="sea_ice_area_fraction" /> 578 579 </file> 580 581 <!-- 582 .............................................. scalar ................................................. 583 --> 584 <file id="ORCA2_1m_scalar" description="Ocean scalar" > <!-- scalar --> 585 <!-- table 2.2 --> 586 <field ref="masstot" name="masso" description="sea_water_mass" /> 587 <field ref="voltot" name="volo" description="sea_water_volume" /> 588 <field ref="sshtot" name="zosga" description="global_average_sea_level_change" /> 589 <field ref="sshsteric" name="zossga" description="global_average_steric_sea_level_change" /> 590 <field ref="sshthster" name="zostoga" description="global_average_thermosteric_sea_level_change" /> 591 <field ref="temptot" name="thetaoga" description="sea_water_potential_temperature" /> 592 <field ref="saltot" name="soga" description="sea_water_salinity" /> 593 </file> 594 595 <!-- 596 .............................................. ptrc_T ................................................. 597 --> 598 <file id="ORCA2_1m_ptrc_T" description="Marine Biogeochemical 3-D Fields" > <!-- ptrc_T --> 599 <field ref="DIC" /> 600 <field ref="Alkalini" /> 601 <field ref="O2" /> 602 <field ref="CaCO3" /> 603 <field ref="PO4" /> 604 <field ref="POC" /> 605 <field ref="Si" /> 606 <field ref="PHY" /> 607 <field ref="ZOO" /> 608 <field ref="DOC" /> 609 <field ref="PHY2" /> 610 <field ref="ZOO2" /> 611 <field ref="BSi" /> 612 <field ref="Fer" /> 613 <field ref="BFe" /> 614 <field ref="GOC" /> 615 <field ref="SFe" /> 616 <field ref="DFe" /> 617 <field ref="DSi" /> 618 <field ref="NFe" /> 619 <field ref="NCHL" /> 620 <field ref="DCHL" /> 621 <field ref="NO3" /> 622 <field ref="NH4" /> 623 </file> 624 625 <!-- 626 .............................................. diad_T ................................................. 627 --> 628 <file id="ORCA2_1m_diad_T" description="Marine Biogeochemical Fields : Additional diagnostics" > 629 <field ref="PH" /> 630 <field ref="CO3" /> 631 <field ref="CO3sat" /> 632 <field ref="PAR" /> 633 <field ref="PPPHY" /> 634 <field ref="PPPHY2" /> 635 <field ref="PPNEWN" /> 636 <field ref="PPNEWD" /> 637 <field ref="PBSi" /> 638 <field ref="PFeN" /> 639 <field ref="PFeD" /> 640 <field ref="PCAL" /> 641 <field ref="DCAL" /> 642 <field ref="GRAZ" /> 643 <field ref="EPC100" /> 644 <field ref="EPFE100" /> 645 <field ref="EPSI100" /> 646 <field ref="EPCAL100" /> 647 <field ref="Cflx" /> 648 <field ref="Oflx" /> 649 <field ref="Kg" /> 650 <field ref="Dpco2" /> 651 <field ref="Dpo2" /> 652 <field ref="Heup" /> 653 <field ref="Irondep" /> 654 <field ref="Nfix" /> 655 </file> 656 657 <!-- 658 .............................................. dbio_T ................................................. 659 --> 660 <file id="ORCA2_1m_dbio_T" description="Marine Biogeochemical Fields : Additional AR5 variables" enabled="<DBIO_ENABLE>" > 661 <field ref="PHYT" /> 662 <field ref="ZOOT" /> 663 <field ref="CHLT" /> 664 <field ref="POCT" /> 665 <field ref="PFET" /> 666 <field ref="TPP" /> 667 <field ref="TPNEW" /> 668 <field ref="TPBFE" /> 669 <field ref="EXPC" /> 670 <field ref="EXPFE" /> 671 <field ref="EXPSI" /> 672 <field ref="EXPCAL" /> 673 <field ref="INTDIC" /> 674 <field ref="O2MIN" /> 675 <field ref="ZO2MIN" /> 676 <field ref="INTNFIX" /> 677 <field ref="INTPPPHY" /> 678 <field ref="INTPPPHY2" /> 679 <field ref="INTPP" /> 680 <field ref="INTPNEW" /> 681 <field ref="INTPBFE" /> 682 <field ref="INTPBSI" /> 683 <field ref="INTPCAL" /> 684 </file> 685 686 </group> 687 688 <!-- 689 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 690 --> 691 </file_definition> 692 693 <!-- 694 ============================================================================================================ 695 = grid definition = 696 = DO NOT CHANGE = 697 ============================================================================================================ 698 --> 699 700 <axis_definition> 701 <axis id="deptht" description="Vertical T levels" unit="m" positive=".false." /> 702 <axis id="depthu" description="Vertical U levels" unit="m" positive=".false." /> 703 <axis id="depthv" description="Vertical V levels" unit="m" positive=".false." /> 704 <axis id="depthw" description="Vertical W levels" unit="m" positive=".false." /> 705 <axis id="none" description="axe non defini" unit="none" size="1" /> 706 </axis_definition> 707 708 <grid_definition> 709 <grid id="grid_T" description="grid T" /> 710 <grid id="grid_U" description="grid U" /> 711 <grid id="grid_V" description="grid V" /> 712 <grid id="grid_W" description="grid W" /> 713 <grid id="scalarpoint" description="scalar" > 714 <zoom id="1point" ibegin="1" jbegin="1" ni="1" nj="1" /> 715 </grid> 716 </grid_definition> 717 246 </context> 247 248 249 <context id="xios"> 250 251 <variable_definition> 252 253 <!-- 254 We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size) 255 --> 256 <variable id="buffer_size" type="integer">50000000</variable> 257 <variable id="buffer_server_factor_size" type="integer">2</variable> 258 <variable id="info_level" type="integer">0</variable> 259 <variable id="using_server" type="boolean">false</variable> 260 <variable id="using_oasis" type="boolean">false</variable> 261 <variable id="oasis_codes_id" type="string" >oceanx</variable> 262 263 </variable_definition> 264 718 265 </context> 719 266 -
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2
r970 r2061 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 2 !! NEMO/OPA : 1 - run manager (namrun) 3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom )3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom, namtsd) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb) 5 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 !! namsbc_apr, namsbc_ssr, namsbc_alb) 6 7 !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 7 8 !! 5 - bottom boundary (nambfr, nambbc, nambbl) … … 9 10 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 10 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol, nammpp, nammpp_dyndist, namctl) 12 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namptr, namhsb) 13 !! 10 - miscellaneous (namsol, nammpp, nammpp_dyndist, namctl) 14 !! 11 - Obs & Assim (namobs, nam_asminc) 13 15 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE.15 16 16 17 !!====================================================================== 17 18 !! *** Run management namelists *** 18 19 !!====================================================================== 19 !! namrun 20 !!====================================================================== 21 20 !! namrun parameters of the run 21 !!====================================================================== 22 ! 22 23 !----------------------------------------------------------------------- 23 24 &namrun ! parameters of the run 24 25 !----------------------------------------------------------------------- 25 nn_no = 0 ! job number 26 cn_exp = "ORCA2" ! AUTO - experience name 26 nn_no = 0 ! job number (no more used...) 27 cn_exp = "ORCA2" ! AUTO - experience name 27 28 nn_it000 = 1 ! AUTO - first time step 28 nn_itend = 5475 ! AUTO - last time step (std 5475)29 nn_date0 = 010101 ! AUTO - initial calendar date yymmdd (used if nn_rstctl=1)29 nn_itend = 5475 ! AUTO - last time step (std 5475) 30 nn_date0 = 010101 ! AUTO - date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 30 31 nn_leapy = 0 ! AUTO - Leap year calendar (1) or not (0) 32 ln_rstart = .false. ! AUTO - start from rest (F) or from a restart file (T) 33 nn_rstctl = 0 ! AUTO - restart control => activated only if ln_rstart = T 34 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 35 ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 36 ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 37 cn_ocerst_in = "restartopa" ! suffix of ocean restart name (input) 38 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 31 39 nn_istate = 0 ! output the initial state (1) or not (0) 32 40 nn_stock = 5475 ! AUTO - frequency of creation of a restart file (modulo referenced to 1) 33 41 nn_write = 5475 ! AUTO - frequency of write in the output file (modulo referenced to nn_it000) 34 42 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) 35 ln_mskland = .true. 43 ln_mskland = .true. ! mask land points in NetCDF outputs (costly: + ~15%) 36 44 ln_clobber = .false. ! clobber (overwrite) an existing file 37 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (working only with iom_nf90 routines) 38 ln_rstart = .false. ! AUTO - start from rest (F) or from a restart file (T) 39 nn_rstctl = 0 ! AUTO - restart control = 0 nn_it000 is not compared to the restart file value 40 ! = 1 use nn_date0 in namelist (not the value in the restart file) 41 ! = 2 calendar parameters read in the restart file 42 cn_ocerst_in = "restartopa"! suffix of ocean restart name (input) 43 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 44 / 45 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 46 / 47 45 48 !!====================================================================== 46 49 !! *** Domain namelists *** … … 49 52 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 50 53 !! namdom space and time domain (bathymetry, mesh, timestep) 51 !!====================================================================== 52 54 !! namtsd data: temperature & salinity 55 !!====================================================================== 56 ! 53 57 !----------------------------------------------------------------------- 54 58 &namzgr ! vertical coordinate … … 61 65 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 62 66 !----------------------------------------------------------------------- 63 rn_sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) 64 rn_sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 67 ln_s_sh94 = .true. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 68 ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 69 ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 70 ! stretching coefficients for all functions 71 rn_sbot_min = 10.0 ! minimum depth of s-bottom surface (>0) (m) 72 rn_sbot_max = 7000.0 ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 73 rn_hc = 150.0 ! critical depth for transition to stretched coordinates 74 !!!!!!! Envelop bathymetry 75 rn_rmax = 0.3 ! maximum cut-off r-value allowed (0<r_max<1) 76 !!!!!!! SH94 stretching coefficients (ln_s_sh94 = .true.) 65 77 rn_theta = 6.0 ! surface control parameter (0<=theta<=20) 66 rn_thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) 67 rn_rmax = 0.15 ! maximum cut-off r-value allowed (0<r_max<1) 68 ln_s_sigma = .false. ! hybrid s-sigma coordinates 69 rn_bb = 0.8 ! stretching with s-sigma 70 rn_hc = 150.0 ! critical depth with s-sigma 78 rn_bb = 0.8 ! stretching with SH94 s-sigma 79 !!!!!!! SF12 stretching coefficient (ln_s_sf12 = .true.) 80 rn_alpha = 4.4 ! stretching with SF12 s-sigma 81 rn_efold = 0.0 ! efold length scale for transition to stretched coord 82 rn_zs = 1.0 ! depth of surface grid box 83 ! bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 84 rn_zb_a = 0.024 ! bathymetry scaling factor for calculating Zb 85 rn_zb_b = -0.2 ! offset for calculating Zb 86 !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 87 rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) 71 88 / 72 89 !----------------------------------------------------------------------- 73 90 &namdom ! space and time domain (bathymetry, mesh, timestep) 74 91 !----------------------------------------------------------------------- 75 nn_bathy = 1 ! compute (=0) or read(=1) the bathymetry file 76 nn_closea = 1 ! closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain 77 nn_msh = 0 ! AUTO create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 78 rn_e3zps_min= 20. ! the thickness of the partial step is set larger than the minimum 79 rn_e3zps_rat= 0.1 ! of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 92 nn_bathy = 1 ! compute (=0) or read (=1) the bathymetry file 93 nn_closea = 0 ! remove (=0) or keep (=1) closed seas and lakes (ORCA) 94 nn_msh = 0 ! AUTO - create (=1) a mesh file or not (=0) 95 rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) 96 rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of 97 rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 80 98 ! 81 rn_rdt = 5760. ! time step for the dynamics (and tracer if n acc=0) ==> 576082 nn_baro = 64 ! number of barotropic time step (for the split explicit algorithm)("key_dynspg_ts")99 rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) 100 nn_baro = 64 ! number of barotropic time step ("key_dynspg_ts") 83 101 rn_atfp = 0.1 ! asselin time filter parameter 84 102 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k) 85 103 ! =0, not used, rdt = rdttra 86 rn_rdtmin = 28800. ! minimum time step on tracers (used if nacc=1) 87 rn_rdtmax = 28800. ! maximum time step on tracers (used if nacc=1) 88 rn_rdth = 800. ! depth variation of tracer time step (used if nacc=1) 104 rn_rdtmin = 28800. ! minimum time step on tracers (used if nn_acc=1) 105 rn_rdtmax = 28800. ! maximum time step on tracers (used if nn_acc=1) 106 rn_rdth = 800. ! depth variation of tracer time step (used if nn_acc=1) 107 / 108 !----------------------------------------------------------------------- 109 &namtsd ! data : Temperature & Salinity 110 !----------------------------------------------------------------------- 111 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! 112 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 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' , '' , ' ' 115 ! 116 cn_dir = './' ! root directory for the location of the runoff files 117 ln_tsd_init = .true. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 118 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T &S input data (T) or not (F) 89 119 / 90 120 !!====================================================================== 91 121 !! *** Surface Boundary Condition namelists *** 92 122 !!====================================================================== 93 !! namsbc surface boundary condition 94 !! namsbc_ana analytical formulation 95 !! namsbc_flx flux formulation 96 !! namsbc_clio CLIO bulk formulea formulation 97 !! namsbc_core CORE bulk formulea formulation 98 !! namsbc_cpl CouPLed formulation ("key_coupled") 99 !! namtra_qsr penetrative solar radiation 100 !! namsbc_rnf river runoffs 101 !! namsbc_ssr sea surface restoring term (for T and/or S) 102 !! namsbc_alb albedo parameters 103 !!====================================================================== 104 123 !! namsbc surface boundary condition 124 !! namsbc_ana analytical formulation 125 !! namsbc_flx flux formulation 126 !! namsbc_clio CLIO bulk formulae formulation 127 !! namsbc_core CORE bulk formulae formulation 128 !! namsbc_mfs MFS bulk formulae formulation 129 !! namsbc_cpl CouPLed formulation ("key_coupled") 130 !! namtra_qsr penetrative solar radiation 131 !! namsbc_rnf river runoffs 132 !! namsbc_apr Atmospheric Pressure 133 !! namsbc_ssr sea surface restoring term (for T and/or S) 134 !! namsbc_alb albedo parameters 135 !!====================================================================== 136 ! 105 137 !----------------------------------------------------------------------- 106 138 &namsbc ! Surface Boundary Condition (surface module) 107 139 !----------------------------------------------------------------------- 108 nn_fsbc = 5 ! frequency of surface boundary condition computation 109 ! (= the frequency of sea-ice model call) 110 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 111 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 112 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 113 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 114 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl ) 140 nn_fsbc = 5 ! frequency of surface boundary condition computation 141 ! (also = the frequency of sea-ice model call) 142 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 143 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 144 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) 145 ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) 146 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 147 ln_cpl = .true. ! Coupled formulation (T => fill namsbc_cpl ) 148 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 115 149 nn_ice = 2 ! =0 no ice boundary condition , 116 150 ! =1 use observed ice-cover , 117 ! =2 ice-model used ("key_lim3" or "key_lim2) 118 nn_ico_cpl = 0 ! ice-ocean coupling : =0 each nn_fsbc 119 ! =1 stresses recomputed each ocean time step ("key_lim3" only) 120 ! =2 combination of 0 and 1 cases ("key_lim3" only) 121 ln_dm2dc = .false. ! daily mean to diurnal cycle short wave (qsr) 122 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 123 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 124 nn_fwb = 0 ! FreshWater Budget: =0 unchecked 125 ! =1 global mean of e-p-r set to zero at each time step 126 ! =2 annual global mean of e-p-r set to zero 127 ! =3 global emp set to zero and spread out over erp area 151 ! =2 ice-model used ("key_lim3" or "key_lim2) 152 nn_ice_embd = 0 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 153 ! =1 levitating ice with mass and salt exchange but no presure effect 154 ! =2 embedded sea-ice (full salt and mass exchanges and pressure) 155 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 156 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 157 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 158 nn_fwb = 0 ! FreshWater Budget: =0 unchecked 159 ! =1 global mean of e-p-r set to zero at each time step 160 ! =2 annual global mean of e-p-r set to zero 161 ! =3 global emp set to zero and spread out over erp area 162 ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave) 163 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave) 164 ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave) 128 165 / 129 166 !----------------------------------------------------------------------- … … 140 177 &namsbc_flx ! surface boundary condition : flux formulation 141 178 !----------------------------------------------------------------------- 142 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !143 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !144 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , ''145 sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , ''146 sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , ''147 sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , ''148 sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , ''149 ! 179 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 180 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 181 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' 182 sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , '' 183 sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , '' 184 sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , '' 185 sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , '' 186 150 187 cn_dir = './' ! root directory for the location of the flux files 151 / 152 !----------------------------------------------------------------------- 153 &namsbc_clio ! namsbc_clio CLIO bulk formul ea154 !----------------------------------------------------------------------- 155 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !156 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !157 sn_utau = 'taux_1m' , -1 , 'sozotaux' , .true. , .true., 'yearly' , '' , ''158 sn_vtau = 'tauy_1m' , -1 , 'sometauy' , .true. , .true., 'yearly' , '' , ''159 sn_wndm = 'flx' , -1 , 'socliowi' , .true. , .true., 'yearly' , '' , ''160 sn_tair = 'flx' , -1 , 'socliot2' , .true. , .true., 'yearly' , '' , ''161 sn_humi = 'flx' , -1 , 'socliohu' , .true. , .true., 'yearly' , '' , ''162 sn_ccov = 'flx' , -1 , 'socliocl' , .false. , .true., 'yearly' , '' , ''163 sn_prec = 'flx' , -1 , 'socliopl' , .false. , .true., 'yearly' , '' , ''164 ! 188 / 189 !----------------------------------------------------------------------- 190 &namsbc_clio ! namsbc_clio CLIO bulk formulae 191 !----------------------------------------------------------------------- 192 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 193 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 194 sn_utau = 'taux_1m' , -1 , 'sozotaux', .true. , .true. , 'yearly' , '' , '' 195 sn_vtau = 'tauy_1m' , -1 , 'sometauy', .true. , .true. , 'yearly' , '' , '' 196 sn_wndm = 'flx' , -1 , 'socliowi', .true. , .true. , 'yearly' , '' , '' 197 sn_tair = 'flx' , -1 , 'socliot2', .true. , .true. , 'yearly' , '' , '' 198 sn_humi = 'flx' , -1 , 'socliohu', .true. , .true. , 'yearly' , '' , '' 199 sn_ccov = 'flx' , -1 , 'socliocl', .false. , .true. , 'yearly' , '' , '' 200 sn_prec = 'flx' , -1 , 'socliopl', .false. , .true. , 'yearly' , '' , '' 201 165 202 cn_dir = './' ! root directory for the location of the bulk files are 166 203 / 167 204 !----------------------------------------------------------------------- 168 &namsbc_core ! namsbc_core CORE bulk formul ea169 !----------------------------------------------------------------------- 170 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights! rotation !171 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename! pairing !172 sn_wndi = 'u 10_core' , -1 , 'u10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'U1'173 sn_wndj = 'v 10_core' , -1 , 'v10' , .true. , .true. , 'yearly' ,'bicubic_weights_orca2.nc' , 'V1'174 sn_qsr = ' qsw_core' , -1 , 'swdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''175 sn_qlw = ' qlw_core' , -1 , 'lwdn' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''176 sn_tair = 't 2_core' , -1 , 't2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''177 sn_humi = 'q 2_core' , -1 , 'q2' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''178 sn_prec = ' precip_core', -1 , 'precip' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''179 sn_snow = ' snow_core' , -1 , 'snow' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''180 sn_tdif = 'taudif_core' , 24 , 'taudif' , .true. , .true. , 'yearly' ,'bilinear_weights_orca2.nc', ''181 ! 205 &namsbc_core ! namsbc_core CORE bulk formulae 206 !----------------------------------------------------------------------- 207 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 208 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 209 sn_wndi = 'u_10.15JUNE2009_orca2' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , '' , 'Uwnd' 210 sn_wndj = 'v_10.15JUNE2009_orca2' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , '' , 'Vwnd' 211 sn_qsr = 'ncar_rad.15JUNE2009_orca2' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , '' , '' 212 sn_qlw = 'ncar_rad.15JUNE2009_orca2' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , '' , '' 213 sn_tair = 't_10.15JUNE2009_orca2' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , '' , '' 214 sn_humi = 'q_10.15JUNE2009_orca2' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , '' , '' 215 sn_prec = 'ncar_precip.15JUNE2009_orca2', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , '' , '' 216 sn_snow = 'ncar_precip.15JUNE2009_orca2', -1 , 'SNOW' , .false. , .true. , 'yearly' , '' , '' 217 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , '' , '' 218 182 219 cn_dir = './' ! root directory for the location of the bulk files 183 ln_2m = . true.! air temperature and humidity referenced at 2m (T) instead 10m (F)184 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data ?220 ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F) 221 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 185 222 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 186 223 / 187 224 !----------------------------------------------------------------------- 188 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 189 !----------------------------------------------------------------------- 190 ! send 191 cn_snd_temperature= 'weighted oce and ice' ! 'oce only' 'weighted oce and ice' 'mixed oce-ice' 192 cn_snd_albedo = 'weighted ice' ! 'none' 'weighted ice' 'mixed oce-ice' 193 cn_snd_thickness = 'none' ! 'none' 'weighted ice and snow' 194 cn_snd_crt_nature = 'mixed oce-ice' ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice' 195 cn_snd_crt_refere = 'cartesian' ! 'spherical' 'cartesian' 196 cn_snd_crt_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 197 cn_snd_crt_grid = 'T' ! 'T' 198 ! receive 199 cn_rcv_w10m = 'coupled' ! 'none' 'coupled' 200 cn_rcv_taumod = 'none' ! 'none' 'coupled' 201 cn_rcv_tau_nature = 'mixed oce-ice' ! 'oce only' 'oce and ice' 'mixed oce-ice' 202 cn_rcv_tau_refere = 'cartesian' ! 'spherical' 'cartesian' 203 cn_rcv_tau_orient = 'eastward-northward' ! 'eastward-northward' or 'local grid' 204 cn_rcv_tau_grid = 'U,V' ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V' 205 cn_rcv_dqnsdt = 'coupled' ! 'none' 'coupled' 206 cn_rcv_qsr = 'conservative' ! 'conservative' 'oce and ice' 'mixed oce-ice' 207 cn_rcv_qns = 'conservative' ! 'conservative' 'oce and ice' 'mixed oce-ice' 208 cn_rcv_emp = 'conservative' ! 'conservative' 'oce and ice' 'mixed oce-ice' 209 cn_rcv_rnf = 'coupled' ! 'coupled' 'climato' 'mixed' 210 cn_rcv_cal = 'coupled' ! 'none' 'coupled' 211 / 212 !----------------------------------------------------------------------- 213 &namsbc_cpl_co2 ! coupled ocean/biogeo/atmosphere model ("key_cpl_carbon_cycle") 214 !----------------------------------------------------------------------- 215 cn_snd_co2 = 'coupled' ! send : 'none' 'coupled' 216 cn_rcv_co2 = 'coupled' ! receive : 'none' 'coupled' 225 &namsbc_mfs ! namsbc_mfs MFS bulk formulae 226 !----------------------------------------------------------------------- 227 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 228 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 229 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 230 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 231 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false. , 'daily' ,'bilinear.nc', '' 232 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 233 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' 234 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' 235 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' 236 237 cn_dir = './ECMWF/' ! root directory for the location of the bulk files 238 / 239 !----------------------------------------------------------------------- 240 &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") 241 !----------------------------------------------------------------------- 242 ! ! description ! multiple ! vector ! vector ! vector ! 243 ! ! ! categories ! reference ! orientation ! grids ! 244 ! send 245 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 246 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 247 sn_snd_thick = 'none' , 'no' , '' , '' , '' 248 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 249 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 250 ! receive 251 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 252 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 253 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 254 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 255 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 256 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 257 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 258 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 259 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 260 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 217 261 / 218 262 !----------------------------------------------------------------------- 219 263 &namtra_qsr ! penetrative solar radiation 220 264 !----------------------------------------------------------------------- 221 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !222 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !223 sn_chl = 'chlorophyll', -1 , 'CHLA' , .true. , .true., 'yearly' , '' , ''224 265 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 266 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 267 sn_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' 268 225 269 cn_dir = './' ! root directory for the location of the runoff files 226 270 ln_traqsr = .true. ! Light penetration (T) or not (F) … … 228 272 ln_qsr_2bd = .false. ! 2 bands light penetration 229 273 ln_qsr_bio = .true. ! bio-model light penetration 230 nn_chldta = 0 ! RGB : Chl 2D data (=1), 3D data (=2) or cst value (=0)274 nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0) 231 275 rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1) 232 276 rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction 233 277 rn_si1 = 23.0 ! 2 bands: longest depth of extinction 234 rn_si2 = 62.0 ! 3 bands: longest depth of extinction (for blue waveband & 0.01 mg/m2 Chl)235 278 / 236 279 !----------------------------------------------------------------------- 237 280 &namsbc_rnf ! runoffs namelist surface boundary condition 238 281 !----------------------------------------------------------------------- 239 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 240 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 241 sn_rnf = 'runoff_1m_nomask' , -1 , 'sorunoff' , .true. , .true. , 'yearly' , '' , '' 242 sn_cnf = 'runoff_1m_nomask' , 0 , 'socoefr' , .false. , .true. , 'yearly' , '' , '' 243 282 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 283 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 284 sn_rnf = 'runoff_core_monthly', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' 285 sn_cnf = 'runoff_core_monthly', 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' 286 sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' 287 sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' 288 sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' 289 244 290 cn_dir = './' ! root directory for the location of the runoff files 245 ln_rnf_emp = .true.! runoffs included into precipitation field (T) or into a file (F)246 ln_rnf_mouth = .true.! specific treatment at rivers mouths247 rn_hrnf = 15.e0! depth over which enhanced vertical mixing is used291 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F) 292 ln_rnf_mouth = .true. ! specific treatment at rivers mouths 293 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used 248 294 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] 249 295 rn_rfact = 1.e0 ! multiplicative factor for runoff 296 ln_rnf_depth = .false. ! read in depth information for runoff 297 ln_rnf_tem = .false. ! read in temperature information for runoff 298 ln_rnf_sal = .false. ! read in salinity information for runoff 299 / 300 !----------------------------------------------------------------------- 301 &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk 302 !----------------------------------------------------------------------- 303 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 304 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 305 sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' 306 307 cn_dir = './' ! root directory for the location of the bulk files 308 rn_pref = 101000._wp ! reference atmospheric pressure [N/m2]/ 309 ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) 310 ln_apr_obc = .false. ! inverse barometer added to OBC ssh data 250 311 / 251 312 !----------------------------------------------------------------------- 252 313 &namsbc_ssr ! surface boundary condition : sea surface restoring 253 314 !----------------------------------------------------------------------- 254 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim! 'yearly'/ ! weights ! rotation !255 ! ! ! (if <0 months) ! name ! (logical) ! (T/F)! 'monthly' ! filename ! pairing !256 sn_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , ''257 sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true., 'yearly' , '' , ''258 315 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 316 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 317 sn_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , '' 318 sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' 319 259 320 cn_dir = './' ! root directory for the location of the runoff files 260 321 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) 261 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 322 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) 262 323 ! or to SSS only (=1) or no damping term (=0) 263 324 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 264 rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu]325 rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] 265 326 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 266 327 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 267 / 328 / 268 329 !----------------------------------------------------------------------- 269 330 &namsbc_alb ! albedo parameters 270 331 !----------------------------------------------------------------------- 271 rn_cloud = 0.0 ! cloud correction to snow and ice albedo272 rn_albice = 0.5 332 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 333 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 273 334 rn_alphd = 0.80 ! coefficients for linear interpolation used to 274 rn_alphc = 0.65 ! compute albedo between two extremes values 335 rn_alphc = 0.65 ! compute albedo between two extremes values 275 336 rn_alphdi = 0.72 ! (Pyane, 1972) 337 / 338 !----------------------------------------------------------------------- 339 &namberg ! iceberg parameters 340 !----------------------------------------------------------------------- 341 ln_icebergs = .false. 342 ln_bergdia = .true. ! Calculate budgets 343 nn_verbose_level = 1 ! Turn on more verbose output if level > 0 344 nn_verbose_write = 15 ! Timesteps between verbose messages 345 nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage 346 ! Initial mass required for an iceberg of each class 347 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 348 ! Proportion of calving mass to apportion to each class 349 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 350 ! Ratio between effective and real iceberg mass (non-dim) 351 ! i.e. number of icebergs represented at a point 352 rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 353 ! thickness of newly calved bergs (m) 354 rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 355 rn_rho_bergs = 850. ! Density of icebergs 356 rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs 357 ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics 358 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits 359 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 360 ln_passive_mode = .false. ! iceberg - ocean decoupling 361 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no) 362 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 363 rn_test_box = 108.0, 116.0, -66.0, -58.0 364 rn_speed_limit = 0. ! CFL speed limit for a berg 365 366 ! filename ! freq (hours) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! 367 ! ! (<0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 368 sn_icb = 'calving' , -1 , 'calvingmask', .true. , .true., 'yearly' , ' ' , ' ' 369 370 cn_dir = './' 276 371 / 277 372 … … 282 377 !! namcla cross land advection 283 378 !! namobc open boundaries parameters ("key_obc") 284 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 379 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 285 380 !! nambdy Unstructured open boundaries ("key_bdy") 286 381 !! namtide Tidal forcing at open boundaries ("key_bdy_tides") 287 382 !!====================================================================== 288 383 ! 289 384 !----------------------------------------------------------------------- 290 385 &namlbc ! lateral momentum boundary condition … … 292 387 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 293 388 ! free slip ! partial slip ! no slip ! strong slip 389 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 294 390 / 295 391 !----------------------------------------------------------------------- … … 301 397 &namobc ! open boundaries parameters ("key_obc") 302 398 !----------------------------------------------------------------------- 303 ln_obc_clim= .false. ! climatological obc data files (T) or not (F)304 ln_vol_cst= .true. ! impose the total volume conservation (T) or not (F)305 ln_obc_fla = .false. ! Flather open boundary condition306 nn_obcdta= 1 ! = 0 the obc data are equal to the initial state399 ln_obc_clim = .false. ! climatological obc data files (T) or not (F) 400 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 401 ln_obc_fla = .false. ! Flather open boundary condition 402 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 307 403 ! = 1 the obc data are read in 'obc.dta' files 308 cn_obcdta= 'annual' ! set to annual if obc datafile hold 1 year of data404 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 309 405 ! set to monthly if obc datafile hold 1 month of data 310 rn_dpein= 1. ! damping time scale for inflow at east open boundary311 rn_dpwin= 1. ! - - - west - -312 rn_dpnin= 1. ! - - - north - -313 rn_dpsin= 1. ! - - - south - -314 rn_dpeob= 3000. ! time relaxation (days) for the east open boundary315 rn_dpwob= 15. ! - - - west - -316 rn_dpnob= 3000. ! - - - north - -317 rn_dpsob= 15. ! - - - south - -318 rn_volemp = 1.! = 0 the total volume change with the surface flux (E-P-R)406 rn_dpein = 1. ! damping time scale for inflow at east open boundary 407 rn_dpwin = 1. ! - - - west - - 408 rn_dpnin = 1. ! - - - north - - 409 rn_dpsin = 1. ! - - - south - - 410 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 411 rn_dpwob = 15. ! - - - west - - 412 rn_dpnob = 3000. ! - - - north - - 413 rn_dpsob = 15. ! - - - south - - 414 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 319 415 ! = 1 the total volume remains constant 320 416 / … … 322 418 &namagrif ! AGRIF zoom ("key_agrif") 323 419 !----------------------------------------------------------------------- 324 nn_cln_update = 3 ! baroclinic update frequency 325 ln_spc_dyn = .true. ! use 0 as special value for dynamics 326 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [s] 327 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [s] 420 nn_cln_update = 3 ! baroclinic update frequency 421 ln_spc_dyn = .true. ! use 0 as special value for dynamics 422 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 423 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 424 / 425 !----------------------------------------------------------------------- 426 &nam_tide ! tide parameters (#ifdef key_tide) 427 !----------------------------------------------------------------------- 428 ln_tide_pot = .true. ! use tidal potential forcing 429 nb_harmo = 11 ! number of constituents used 430 clname(1) = 'M2' ! name of constituent 431 clname(2) = 'S2' 432 clname(3) = 'N2' 433 clname(4) = 'K1' 434 clname(5) = 'O1' 435 clname(6) = 'Q1' 436 clname(7) = 'M4' 437 clname(8) = 'K2' 438 clname(9) = 'P1' 439 clname(10) = 'Mf' 440 clname(11) = 'Mm' 328 441 / 329 442 !----------------------------------------------------------------------- 330 443 &nambdy ! unstructured open boundaries ("key_bdy") 331 444 !----------------------------------------------------------------------- 332 filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) 333 filbdy_data_T = 'bdydata_grid_T.nc' ! name of data file (T-points) 334 filbdy_data_U = 'bdydata_grid_U.nc' ! name of data file (U-points) 335 filbdy_data_V = 'bdydata_grid_V.nc' ! name of data file (V-points) 336 ln_bdy_clim = .false. ! contain 1 (T) or 12 (F) time dumps and be cyclic 337 ln_bdy_vol = .true. ! total volume correction (see volbdy parameter) 338 ln_bdy_mask = .false. ! boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F) 339 ln_bdy_tides = .true. ! Apply tidal harmonic forcing with Flather condition 340 ln_bdy_dyn_fla = .true. ! Apply Flather condition to velocities 341 ln_bdy_tra_frs = .false. ! Apply FRS condition to temperature and salinity 342 ln_bdy_dyn_frs = .false. ! Apply FRS condition to velocities 343 nbdy_dta = 1 ! = 0, bdy data are equal to the initial state 445 nb_bdy = 1 ! number of open boundary sets 446 ln_coords_file = .true. ! =T : read bdy coordinates from file 447 cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files 448 ln_mask_file = .false. ! =T : read mask from file 449 cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.) 450 nn_dyn2d = 2 ! boundary conditions for barotropic fields 451 nn_dyn2d_dta = 3 ! = 0, bdy data are equal to the initial state 344 452 ! = 1, bdy data are read in 'bdydata .nc' files 345 nb_rimwidth = 9 ! width of the relaxation zone 346 volbdy = 0 ! = 0, the total water flux across open boundaries is zero 347 ! = 1, the total volume of the system is conserved 348 / 349 !----------------------------------------------------------------------- 350 &nambdy_tide ! tidal forcing at unstructured boundaries 351 !----------------------------------------------------------------------- 352 filtide = 'bdytide_' ! file name root of tidal forcing files 353 tide_cpt = 'M2','S1' ! names of tidal components used 354 tide_speed = 28.984106, 15.000001 ! phase speeds of tidal components (deg/hour) 355 ln_tide_date = .false. ! adjust tidal harmonics for start date of run 356 / 357 453 ! = 2, use tidal harmonic forcing data from files 454 ! = 3, use external data AND tidal harmonic forcing 455 nn_dyn3d = 0 ! boundary conditions for baroclinic velocities 456 nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state 457 ! = 1, bdy data are read in 'bdydata .nc' files 458 nn_tra = 1 ! boundary conditions for T and S 459 nn_tra_dta = 1 ! = 0, bdy data are equal to the initial state 460 ! = 1, bdy data are read in 'bdydata .nc' files 461 nn_rimwidth = 10 ! width of the relaxation zone 462 ln_vol = .false. ! total volume correction (see nn_volctl parameter) 463 nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero 464 / 465 !----------------------------------------------------------------------- 466 &nambdy_dta ! open boundaries - external data ("key_bdy") 467 !----------------------------------------------------------------------- 468 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 469 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 470 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true. , .false. , 'daily' , '' , '' 471 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true. , .false. , 'daily' , '' , '' 472 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true. , .false. , 'daily' , '' , '' 473 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true. , .false. , 'daily' , '' , '' 474 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true. , .false. , 'daily' , '' , '' 475 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'daily' , '' , '' 476 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'daily' , '' , '' 477 cn_dir = 'bdydta/' 478 ln_full_vel = .false. 479 / 480 !----------------------------------------------------------------------- 481 &nambdy_tide ! tidal forcing at open boundaries 482 !----------------------------------------------------------------------- 483 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 484 tide_cpt(1) ='Q1' ! names of tidal components used 485 tide_cpt(2) ='O1' ! names of tidal components used 486 tide_cpt(3) ='P1' ! names of tidal components used 487 tide_cpt(4) ='S1' ! names of tidal components used 488 tide_cpt(5) ='K1' ! names of tidal components used 489 tide_cpt(6) ='2N2' ! names of tidal components used 490 tide_cpt(7) ='MU2' ! names of tidal components used 491 tide_cpt(8) ='N2' ! names of tidal components used 492 tide_cpt(9) ='NU2' ! names of tidal components used 493 tide_cpt(10) ='M2' ! names of tidal components used 494 tide_cpt(11) ='L2' ! names of tidal components used 495 tide_cpt(12) ='T2' ! names of tidal components used 496 tide_cpt(13) ='S2' ! names of tidal components used 497 tide_cpt(14) ='K2' ! names of tidal components used 498 tide_cpt(15) ='M4' ! names of tidal components used 499 tide_speed(1) = 13.398661 ! phase speeds of tidal components (deg/hour) 500 tide_speed(2) = 13.943036 ! phase speeds of tidal components (deg/hour) 501 tide_speed(3) = 14.958932 ! phase speeds of tidal components (deg/hour) 502 tide_speed(4) = 15.000001 ! phase speeds of tidal components (deg/hour) 503 tide_speed(5) = 15.041069 ! phase speeds of tidal components (deg/hour) 504 tide_speed(6) = 27.895355 ! phase speeds of tidal components (deg/hour) 505 tide_speed(7) = 27.968210 ! phase speeds of tidal components (deg/hour) 506 tide_speed(8) = 28.439730 ! phase speeds of tidal components (deg/hour) 507 tide_speed(9) = 28.512585 ! phase speeds of tidal components (deg/hour) 508 tide_speed(10) = 28.984106 ! phase speeds of tidal components (deg/hour) 509 tide_speed(11) = 29.528479 ! phase speeds of tidal components (deg/hour) 510 tide_speed(12) = 29.958935 ! phase speeds of tidal components (deg/hour) 511 tide_speed(13) = 30.000002 ! phase speeds of tidal components (deg/hour) 512 tide_speed(14) = 30.082138 ! phase speeds of tidal components (deg/hour) 513 tide_speed(15) = 57.968212 ! phase speeds of tidal components (deg/hour) 514 ln_tide_date = .true. ! adjust tidal harmonics for start date of run 515 / 358 516 !!====================================================================== 359 517 !! *** Bottom boundary condition *** 360 518 !!====================================================================== 361 519 !! nambfr bottom friction 362 !! nambbc bottom temperature boundary condition ("key_trabbc")363 !! nambbl bottom boundary layer scheme ("key_trabbl _dif","key_trabbl_adv")364 !!====================================================================== 365 520 !! nambbc bottom temperature boundary condition 521 !! nambbl bottom boundary layer scheme ("key_trabbl") 522 !!====================================================================== 523 ! 366 524 !----------------------------------------------------------------------- 367 525 &nambfr ! bottom friction … … 371 529 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 372 530 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 373 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) 374 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 375 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d = .true.) 531 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 532 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 533 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) 534 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 376 535 / 377 536 !----------------------------------------------------------------------- 378 537 &nambbc ! bottom temperature boundary condition 379 538 !----------------------------------------------------------------------- 380 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 539 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 540 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 381 541 ! = 1 constant flux 382 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 542 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 383 543 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 384 544 / … … 386 546 &nambbl ! bottom boundary layer scheme 387 547 !----------------------------------------------------------------------- 388 ! ! diffusive bbl ("key_trabbl") 389 ! ! advective bbl ("key_trabbl_adv") 390 rn_ahtbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] 391 / 548 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) 549 nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0) 550 rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] 551 rn_gambbl = 10. ! advective bbl coefficient [s] 552 / 553 392 554 !!====================================================================== 393 555 !! Tracer (T & S ) namelists … … 396 558 !! namtra_adv advection scheme 397 559 !! namtra_ldf lateral diffusion scheme 398 !! namtra_dmp T & S newtonian damping ("key_tradmp")399 !!====================================================================== 400 560 !! namtra_dmp T & S newtonian damping 561 !!====================================================================== 562 ! 401 563 !----------------------------------------------------------------------- 402 564 &nameos ! ocean physical parameters 403 565 !----------------------------------------------------------------------- 404 nn_eos = 0! type of equation of state and Brunt-Vaisala frequency566 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 405 567 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 406 568 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 407 569 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 408 rn_alpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) 409 rn_beta = 0.001 ! saline expension coefficient (neos= 2) 410 / 411 !----------------------------------------------------------------------- 412 &namtra_adv ! advection scheme for tracer 413 !----------------------------------------------------------------------- 414 ln_traadv_cen2 = .false. ! 2nd order centered scheme 415 ln_traadv_tvd = .true. ! TVD scheme 416 ln_traadv_muscl = .false. ! MUSCL scheme 417 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 418 ln_traadv_ubs = .false. ! UBS scheme 419 / 420 !----------------------------------------------------------------------- 421 &namtra_ldf ! lateral diffusion scheme for tracer 422 !----------------------------------------------------------------------- 423 ! Type of the operator : 424 ln_traldf_lap = .true. ! laplacian operator 425 ln_traldf_bilap = .false. ! bilaplacian operator 426 ! Direction of action : 427 ln_traldf_level = .false. ! iso-level 428 ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) 429 ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") 430 ! Coefficient 431 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 432 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 433 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") 434 / 435 !----------------------------------------------------------------------- 436 &namtra_dmp ! tracer: T & S newtonian damping ('key_tradmp') 437 !----------------------------------------------------------------------- 570 rn_alpha = 2.0e-4 ! thermal expension coefficient (nn_eos= 1 or 2) 571 rn_beta = 7.7e-4 ! saline expension coefficient (nn_eos= 2) 572 / 573 !----------------------------------------------------------------------- 574 &namtra_adv ! advection scheme for tracer 575 !----------------------------------------------------------------------- 576 ln_traadv_cen2 = .false. ! 2nd order centered scheme 577 ln_traadv_tvd = .true. ! TVD scheme 578 ln_traadv_muscl = .false. ! MUSCL scheme 579 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 580 ln_traadv_ubs = .false. ! UBS scheme 581 ln_traadv_qck = .false. ! QUICKEST scheme 582 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 583 / 584 !---------------------------------------------------------------------------------- 585 &namtra_ldf ! lateral diffusion scheme for tracers 586 !---------------------------------------------------------------------------------- 587 ! ! Operator type: 588 ln_traldf_lap = .true. ! laplacian operator 589 ln_traldf_bilap = .false. ! bilaplacian operator 590 ! ! Direction of action: 591 ln_traldf_level = .false. ! iso-level 592 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 593 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 594 ! ! Griffies parameters (all need "key_ldfslp") 595 ln_traldf_grif = .false. ! use griffies triads 596 ln_traldf_gdia = .false. ! output griffies eddy velocities 597 ln_triad_iso = .false. ! pure lateral mixing in ML 598 ln_botmix_grif = .false. ! lateral mixing on bottom 599 ! ! Coefficients 600 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 601 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 602 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 603 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 604 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 605 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 606 ! (normally=0; not used with Griffies) 607 / 608 !----------------------------------------------------------------------- 609 &namtra_dmp ! tracer: T & S newtonian damping 610 !----------------------------------------------------------------------- 611 ln_tradmp = .true. ! add a damping termn (T) or not (F) 438 612 nn_hdmp = -1 ! horizontal shape =-1, damping in Med and Red Seas only 439 613 ! =XX, damping poleward of XX degrees (XX>0) 440 614 ! + F(distance-to-coast) + Red and Med Seas 441 nn_zdmp = 1! vertical shape =0 damping throughout the water column615 nn_zdmp = 0 ! vertical shape =0 damping throughout the water column 442 616 ! =1 no damping in the mixing layer (kz criteria) 443 617 ! =2 no damping in the mixed layer (rho crieria) … … 445 619 rn_bot = 360. ! bottom time scale of damping [days] 446 620 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 447 nn_file = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) 448 / 621 nn_file = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 622 / 623 449 624 !!====================================================================== 450 625 !! *** Dynamics namelists *** … … 456 631 !! namdyn_ldf lateral diffusion scheme 457 632 !!====================================================================== 458 633 ! 459 634 !----------------------------------------------------------------------- 460 635 &namdyn_adv ! formulation of the momentum advection 461 636 !----------------------------------------------------------------------- 462 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 637 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 463 638 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 464 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 465 / 639 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 640 / 466 641 !----------------------------------------------------------------------- 467 642 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 468 643 !----------------------------------------------------------------------- 469 ln_dynvor_ene = .false. ! enstrophy conserving scheme 470 ln_dynvor_ens = .false. ! energy conserving scheme 471 ln_dynvor_mix = .false. ! mixed scheme 472 ln_dynvor_een = .true. ! energy & enstrophy scheme 644 ln_dynvor_ene = .false. ! enstrophy conserving scheme 645 ln_dynvor_ens = .false. ! energy conserving scheme 646 ln_dynvor_mix = .false. ! mixed scheme 647 ln_dynvor_een = .true. ! energy & enstrophy scheme 473 648 / 474 649 !----------------------------------------------------------------------- 475 650 &namdyn_hpg ! Hydrostatic pressure gradient option 476 651 !----------------------------------------------------------------------- 477 ln_hpg_zco = .false. ! z-coordinate - full steps 652 ln_hpg_zco = .false. ! z-coordinate - full steps 478 653 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) 479 654 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 480 ln_hpg_hel = .false. ! s-coordinate (helsinki modification)481 ln_hpg_wdj = .false. ! s-coordinate (weighted density jacobian)482 655 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 483 ln_hpg_rot = .false. ! s-coordinate (ROTated axes scheme) 484 rn_gamma = 0.e0 ! weighting coefficient (wdj scheme) 656 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 485 657 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 486 658 ! centered time scheme (F) 487 nn_dynhpg_rst = 0 ! =1 dynhpg restartable restart or not (=0)488 659 / 489 660 !----------------------------------------------------------------------- … … 497 668 &namdyn_ldf ! lateral diffusion on momentum 498 669 !----------------------------------------------------------------------- 499 ! Type of the operator :500 ln_dynldf_lap = .true. ! laplacian operator 501 ln_dynldf_bilap = .false. ! bilaplacian operator 502 ! Direction of action :503 ln_dynldf_level = .false. ! iso-level 670 ! ! Type of the operator : 671 ln_dynldf_lap = .true. ! laplacian operator 672 ln_dynldf_bilap = .false. ! bilaplacian operator 673 ! ! Direction of action : 674 ln_dynldf_level = .false. ! iso-level 504 675 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.) 505 676 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp") 506 ! Coefficient 507 rn_ahm_0 = 40000. ! horizontal eddy viscosity [m2/s] 508 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 509 / 677 ! ! Coefficient 678 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 679 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 680 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 681 / 682 510 683 !!====================================================================== 511 684 !! Tracers & Dynamics vertical physics namelists 512 685 !!====================================================================== 513 !! 514 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric")515 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke")516 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")517 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm")518 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx")519 !!====================================================================== 520 686 !! namzdf vertical physics 687 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 688 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 689 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp") 690 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 691 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") 692 !!====================================================================== 693 ! 521 694 !----------------------------------------------------------------------- 522 695 &namzdf ! vertical physics … … 529 702 nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1) 530 703 rn_avevd = 100. ! evd mixing coefficient [m2/s] 531 ln_zdfnpc = .false. ! Non-Penetrative algorithm (T) or not (F)704 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F) 532 705 nn_npc = 1 ! frequency of application of npc 533 706 nn_npcp = 365 ! npc control print frequency … … 541 714 rn_alp = 5. ! coefficient of the parameterization 542 715 nn_ric = 2 ! coefficient of the parameterization 716 rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation 717 rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m) 718 rn_mldmax =1000.0 ! maximum allowable mixed-layer depth estimate (m) 719 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer 720 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer 721 ln_mldw = .true. ! Flag to use or not the mized layer depth param. 543 722 / 544 723 !----------------------------------------------------------------------- … … 547 726 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 548 727 rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation 549 rn_ebb = 6 0. ! coef. of the surface input of tke728 rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 550 729 rn_emin = 1.e-6 ! minimum value of tke [m2/s2] 551 730 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] 552 rn_bshear = 1.e-20 ! background shear (>0)553 731 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom 554 732 ! = 1 bounded by the local vertical scale factor 555 733 ! = 2 first vertical derivative of mixing length bounded by 1 556 ! = 3 same criteria as case 2 but applied in a different way734 ! = 3 as =2 with distinct disspipative an mixing length scale 557 735 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 558 ln_mxl0 = .true. ! mixing length scale surface value as function of wind stress (T) or not (F) 559 rn_lmin = 0.001 ! interior buoyancy lenght scale minimum value 560 rn_lmin0 = 0.01 ! surface buoyancy lenght scale minimum value 561 nn_etau = 0 ! exponentially deceasing penetration of tke due to internal & intertial waves 562 ! = 0 no penetration ( O(2 km) resolution) 563 ! = 1 additional tke source (rn_efr * en) 564 ! = 2 additional tke source (rn_efr * en) applied only at the base of the mixed layer 565 ! = 3 additional tke source (HF contribution: mean of stress module - module of mean stress) 566 nn_htau = 1 ! type of exponential decrease of tke penetration 736 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 737 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value 738 ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002) 739 rn_lc = 0.15 ! coef. associated to Langmuir cells 740 nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves 741 ! = 0 no penetration 742 ! = 1 add a tke source below the ML 743 ! = 2 add a tke source just at the base of the ML 744 ! = 3 as = 1 applied on HF part of the stress ("key_coupled") 745 rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 746 nn_htau = 1 ! type of exponential decrease of tke penetration below the ML 567 747 ! = 0 constant 10 m length scale 568 ! = 1 0.5m at the equator to 30m at high latitudes 569 ! = 2 30 meters constant depth penetration 570 ! otion used only if nn_etau = 1 or 2: 571 rn_efr = 0.05 ! fraction of surface tke value which penetrates inside the ocean 572 ! otion used only if nn_etau = 3: 573 rn_addhft = -1.e-3 ! add offset applied to the "mean of stress module - module of mean stress" (always kept > 0) 574 rn_sclhft = 1. ! scale factor applied to the "mean of stress module - module of mean stress" 575 ln_lc = .true. ! Langmuir cell parameterisation 576 rn_lc = 0.15 ! coef. associated to Langmuir cells 748 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 577 749 / 578 750 !------------------------------------------------------------------------ 579 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and option nally:751 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally: 580 752 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 581 ln_kpprimix = .true. ! shear instability mixing 753 ln_kpprimix = .true. ! shear instability mixing 582 754 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s] 583 755 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s] 584 756 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability 585 757 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s] 586 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 587 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 758 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2] 759 rn_difcon = 1. ! maximum mixing in interior convection [m2/s] 588 760 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv 589 761 nn_ave = 1 ! constant (=0) or profile (=1) background on avt 762 / 763 !----------------------------------------------------------------------- 764 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") 765 !----------------------------------------------------------------------- 766 rn_emin = 1.e-6 ! minimum value of e [m2/s2] 767 rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] 768 ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) 769 rn_clim_galp = 0.53 ! galperin limit 770 ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation 771 ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case 772 rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux 773 rn_charn = 70000. ! Charnock constant for wb induced roughness length 774 nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 775 nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum) 776 nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 777 nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum) 778 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 779 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 590 780 / 591 781 !----------------------------------------------------------------------- … … 601 791 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1) 602 792 rn_tfe = 0.333 ! tidal dissipation efficiency 603 rn_me = 0.2 ! mixing efficiency 604 ln_tmx_itf = . TRUE. ! ITF specific parameterisation793 rn_me = 0.2 ! mixing efficiency 794 ln_tmx_itf = .true. ! ITF specific parameterisation 605 795 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 606 796 / 607 !!====================================================================== 608 !! *** Miscelaneous namelists *** 797 798 !!====================================================================== 799 !! *** Miscellaneous namelists *** 609 800 !!====================================================================== 610 801 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 611 !! nammpp_dyndist Massively Parallel domain decomposition ("key_agrif" && "key_mpp_dyndist")612 802 !! namctl Control prints & Benchmark 613 !! namsol elliptic solver / island / free surface 614 !!====================================================================== 615 616 !----------------------------------------------------------------------- 617 &namsol ! elliptic solver / island / free surface 803 !! namsol elliptic solver / island / free surface 804 !!====================================================================== 805 ! 806 !----------------------------------------------------------------------- 807 &namsol ! elliptic solver / island / free surface 618 808 !----------------------------------------------------------------------- 619 809 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) … … 633 823 ! buffer blocking send or immediate non-blocking sends, resp. 634 824 nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation 635 / 636 !----------------------------------------------------------------------- 637 &nammpp_dyndist ! Massively Parallel Distribution for AGRIF zoom ("key_agrif" && "key_mpp_dyndist") 638 !----------------------------------------------------------------------- 639 jpni = 1 ! jpni number of processors following i 640 jpnj = 1 ! jpnj number of processors following j 641 jpnij = 1 ! jpnij number of local domains 825 ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold 826 jpni = 0 ! jpni number of processors following i (set automatically if < 1) 827 jpnj = 0 ! jpnj number of processors following j (set automatically if < 1) 828 jpnij = 0 ! jpnij number of local domains (set automatically if < 1) 642 829 / 643 830 !----------------------------------------------------------------------- … … 654 841 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 655 842 ! (no physical validity of the results) 656 nn_ bit_cmp = 0 ! bit comparison mode (1/0): CAUTION use zero except for test657 ! of comparison between single and multiple processor runs 658 / 659 660 !! ======================================================================661 !! *** Diagnostics namelists ***662 !! ======================================================================843 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 844 / 845 846 !!====================================================================== 847 !! *** Diagnostics namelists *** 848 !!====================================================================== 849 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 663 850 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 664 !! namgap level mean model-data gap ("key_diagap")665 851 !! namflo float parameters ("key_float") 666 852 !! namptr Poleward Transport Diagnostics 667 !!====================================================================== 668 853 !! namhsb Heat and salt budgets 854 !!====================================================================== 855 ! 856 !----------------------------------------------------------------------- 857 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 858 !----------------------------------------------------------------------- 859 nn_nchunks_i= 4 ! number of chunks in i-dimension 860 nn_nchunks_j= 4 ! number of chunks in j-dimension 861 nn_nchunks_k= 31 ! number of chunks in k-dimension 862 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 863 ! is optimal for postprocessing which works exclusively with horizontal slabs 864 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 865 ! (F) ignore chunking information and produce netcdf3-compatible files 866 / 669 867 !----------------------------------------------------------------------- 670 868 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 671 ! ! or mixed-layer trends or barotropic vorticity ( 'key_trdmld' or"key_trdvor")869 ! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor") 672 870 !----------------------------------------------------------------------- 673 871 nn_trd = 365 ! time step frequency dynamics and tracers trends … … 680 878 / 681 879 !----------------------------------------------------------------------- 682 &namgap ! level mean model-data gap ('key_diagap')683 !-----------------------------------------------------------------------684 nn_gap = 15 ! time-step frequency of model-data gap computation685 nn_prg = 10 ! AUTO - time-step frequency of gap print in model output686 /687 !-----------------------------------------------------------------------688 880 &namflo ! float parameters ("key_float") 689 881 !----------------------------------------------------------------------- 690 ln_rstflo = .false. ! float restart (T) or not (F) 691 nn_writefl= 75 ! frequency of writing in float output file 692 nn_stockfl= 5475 ! frequency of creation of the float restart file 693 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 694 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 695 ! or computed with Blanke' scheme (F) 882 jpnfl = 1 ! total number of floats during the run 883 jpnnewflo = 0 ! number of floats for the restart 884 ln_rstflo = .false. ! float restart (T) or not (F) 885 nn_writefl = 75 ! frequency of writing in float output file 886 nn_stockfl = 5475 ! frequency of creation of the float restart file 887 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 888 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 889 ! or computed with Blanke' scheme (F) 890 ln_ariane = .true. ! Input with Ariane tool convention(T) 891 ln_flo_ascii = .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 696 892 / 697 893 !----------------------------------------------------------------------- 698 894 &namptr ! Poleward Transport Diagnostic 699 895 !----------------------------------------------------------------------- 700 ln_diaptr = .true. 896 ln_diaptr = .true. ! Poleward heat and salt transport (T) or not (F) 701 897 ln_diaznl = .true. ! Add zonal means and meridional stream functions 702 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 898 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 703 899 ! (orca configuration only, need input basins mask file named "subbasins.nc" 704 900 ln_ptrcomp = .true. ! Add decomposition : overturning 705 nf_ptr = 1 ! Frequency of ptr computation [time step] 706 nf_ptr_wri = 15 ! AUTO - Frequency of ptr outputs 707 / 901 nn_fptr = 1 ! Frequency of ptr computation [time step] 902 nn_fwri = 15 ! Frequency of ptr outputs [time step] 903 / 904 !----------------------------------------------------------------------- 905 &namhsb ! Heat and salt budgets 906 !----------------------------------------------------------------------- 907 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 908 / 909 !----------------------------------------------------------------------- 910 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 911 !----------------------------------------------------------------------- 912 nit000_han = 1 ! First time step used for harmonic analysis 913 nitend_han = 75 ! Last time step used for harmonic analysis 914 nstep_han = 15 ! Time step frequency for harmonic analysis 915 tname(1) = 'M2' ! Name of tidal constituents 916 tname(2) = 'K1' 917 / 918 !----------------------------------------------------------------------- 919 &namdct ! transports through sections 920 !----------------------------------------------------------------------- 921 nn_dct = 15 ! time step frequency for transports computing 922 nn_dctwri = 15 ! time step frequency for transports writing 923 nn_secdebug = 112 ! 0 : no section to debug 924 ! -1 : debug all section 925 ! 0 < n : debug section number n 926 / 927 928 !!====================================================================== 929 !! *** Observation & Assimilation namelists *** 930 !!====================================================================== 931 !! namobs observation and model comparison ('key_diaobs') 932 !! nam_asminc assimilation increments ('key_asminc') 933 !!====================================================================== 934 ! 935 !----------------------------------------------------------------------- 936 &namobs ! observation usage switch ('key_diaobs') 937 !----------------------------------------------------------------------- 938 ln_t3d = .false. ! Logical switch for T profile observations 939 ln_s3d = .false. ! Logical switch for S profile observations 940 ln_ena = .false. ! Logical switch for ENACT insitu data set 941 ! ! ln_cor Logical switch for Coriolis insitu data set 942 ln_profb = .false. ! Logical switch for feedback insitu data set 943 ln_sla = .false. ! Logical switch for SLA observations 944 945 ln_sladt = .false. ! Logical switch for AVISO SLA data 946 947 ln_slafb = .false. ! Logical switch for feedback SLA data 948 ! ln_ssh Logical switch for SSH observations 949 950 ln_sst = .true. ! Logical switch for SST observations 951 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 952 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 953 954 ln_sstfb = .false. ! Logical switch for feedback SST data 955 ! ln_sss Logical switch for SSS observations 956 ! ln_seaice Logical switch for Sea Ice observations 957 ! ln_vel3d Logical switch for velocity observations 958 ! ln_velavcur Logical switch for velocity daily av. cur. 959 ! ln_velhrcur Logical switch for velocity high freq. cur. 960 ! ln_velavadcp Logical switch for velocity daily av. ADCP 961 ! ln_velhradcp Logical switch for velocity high freq. ADCP 962 ! ln_velfb Logical switch for feedback velocity data 963 ! ln_grid_global Global distribtion of observations 964 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 965 ! grid_search_file Grid search lookup file header 966 ! enactfiles ENACT input observation file names 967 ! coriofiles Coriolis input observation file name 968 ! ! profbfiles: Profile feedback input observation file name 969 profbfiles = 'profiles_01.nc' 970 ! ln_profb_enatim Enact feedback input time setting switch 971 ! slafilesact Active SLA input observation file name 972 ! slafilespas Passive SLA input observation file name 973 ! ! slafbfiles: Feedback SLA input observation file name 974 slafbfiles = 'sla_01.nc' 975 ! sstfiles GHRSST input observation file name 976 ! ! sstfbfiles: Feedback SST input observation file name 977 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 978 ! seaicefiles Sea Ice input observation file name 979 ! velavcurfiles Vel. cur. daily av. input file name 980 ! velhvcurfiles Vel. cur. high freq. input file name 981 ! velavadcpfiles Vel. ADCP daily av. input file name 982 ! velhvadcpfiles Vel. ADCP high freq. input file name 983 ! velfbfiles Vel. feedback input observation file name 984 ! dobsini Initial date in window YYYYMMDD.HHMMSS 985 ! dobsend Final date in window YYYYMMDD.HHMMSS 986 ! n1dint Type of vertical interpolation method 987 ! n2dint Type of horizontal interpolation method 988 ! ln_nea Rejection of observations near land switch 989 nmsshc = 0 ! MSSH correction scheme 990 ! mdtcorr MDT correction 991 ! mdtcutoff MDT cutoff for computed correction 992 ln_altbias = .false. ! Logical switch for alt bias 993 ln_ignmis = .true. ! Logical switch for ignoring missing files 994 ! endailyavtypes ENACT daily average types 995 ln_grid_global = .true. 996 ln_grid_search_lookup = .false. 997 / 998 !----------------------------------------------------------------------- 999 &nam_asminc ! assimilation increments ('key_asminc') 1000 !----------------------------------------------------------------------- 1001 ln_bkgwri = .false. ! Logical switch for writing out background state 1002 ln_trainc = .false. ! Logical switch for applying tracer increments 1003 ln_dyninc = .false. ! Logical switch for applying velocity increments 1004 ln_sshinc = .false. ! Logical switch for applying SSH increments 1005 ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) 1006 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) 1007 nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] 1008 nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1] 1009 nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1] 1010 nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] 1011 niaufn = 0 ! Type of IAU weighting function 1012 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 1013 salfixmin = -9999 ! Minimum salinity after applying the increments 1014 nn_divdmp = 0 ! Number of iterations of divergence damping operator 1015 / 1016 !----------------------------------------------------------------------- 1017 &namsbc_wave ! External fields from wave model 1018 !----------------------------------------------------------------------- 1019 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 1020 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 1021 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' ,'' , '' 1022 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' ,'' , '' 1023 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' ,'' , '' 1024 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' ,'' , '' 1025 ! 1026 cn_dir_cdg = './' ! root directory for the location of drag coefficient files 1027 / 1028 !----------------------------------------------------------------------- 1029 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 1030 !----------------------------------------------------------------------- 1031 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 1032 ln_neptsimp = .false. ! yes/no use simplified neptune 1033 1034 ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep 1035 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 1036 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 1037 ! Specify whether to ramp down the Neptune velocity in shallow 1038 ! water, and if so the depth range controlling such ramping down 1039 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 1040 rn_htrmin = 100.0 ! min. depth of transition range 1041 rn_htrmax = 200.0 ! max. depth of transition range 1042 / -
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ice_ORCA2
r705 r2061 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/LIM 2 :1 - dynamics/advection/thermo (namicerun)2 !! NEMO/LIM-2 : 1 - dynamics/advection/thermo (namicerun) 3 3 !! namelists 2 - ice intialisation (namiceini) 4 4 !! 3 - ice dynamic (namicedyn) 5 !! 5 - ice advection (namicetrp) 6 !! 6 - thermodynamic (namicethd) 7 !! 7 - ice salinity (namicesal) 8 !! 8 - mechanical redistribution of ice (namiceitdme) 9 !! 3 - ice diagnostics (namicedia) 10 !! 9 - ice outputs (namiceout) 5 !! 4 - ice advection (namicetrp) 6 !! 5 - thermodynamic (namicethd) 7 !! 6 - ice damping (namice_dmp) 8 !! 7 - ice diagnostics (namicedia) 9 !! 8 - ice outputs (namiceout) 11 10 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 12 11 … … 16 15 cn_icerst_in = "restart_ice_in" ! suffix of ice restart name (input) 17 16 cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) 18 ln_limdyn = .true.! ice dynamics (T) or thermodynamics only (F)19 ln_limdmp = .false. ! restoring ice thickness and fraction leads (T) or not (F)20 acrit = 1.0e-06 , 1.0e-06 ! minimum fraction for leads in the Northern (Southern) Hemisphere21 hsndif = 0.0! computation of temperature in snow (=0.0) or not22 hicdif = 0.0! computation of temperature in ice (=0.0) or not (=9999.0)17 ln_limdyn = .true. ! ice dynamics (T) or thermodynamics only (F) 18 ln_limdmp = .false. ! restoring ice thickness and fraction leads (T => fill namice_dmp) 19 acrit = 1.0e-06 , 1.0e-06 ! minimum lead fraction in the Northern & Southern Hemispheres 20 hsndif = 0.e0 ! computation of temperature in snow (=0.0) or not 21 hicdif = 0.e0 ! computation of temperature in ice (=0.0) or not (=9999.0) 23 22 / 24 23 !----------------------------------------------------------------------- 25 24 &namiceini ! ice initialisation 26 25 !----------------------------------------------------------------------- 27 ln_limini = .false. ! read the i ce initial state in the file'Ice_initialization.nc' (T) or not (F)26 ln_limini = .false. ! read the initial state in 'Ice_initialization.nc' (T) or not (F) 28 27 ttest = 2.0 ! threshold water temperature for initial sea ice 29 28 hninn = 0.5 ! initial snow thickness in the north … … 49 48 c_rhg = 20.0 ! 2nd bulk-rhelogy parameter 50 49 etamn = 0.0e+07 ! minimun value for viscosity 51 creepl = 2.0e-08 ! creep limit50 creepl = 1.0e-08 ! creep limit 52 51 ecc = 2.0 ! eccentricity of the elliptical yield curve 53 52 ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice [m2/s] 53 nevp = 120 ! number of EVP subcycling iterations 54 telast = 9600 ! timescale for EVP elastic waves 55 alphaevp = 1.0 ! coefficient for the solution of EVP int. stresses 54 56 / 55 57 !----------------------------------------------------------------------- 56 58 &namicetrp ! ice transport 57 59 !----------------------------------------------------------------------- 58 bound = 0. ! boundary conditions (=0.0 no-slip, =1.0free-slip)60 bound = 0.0 ! boundary conditions (=0. no-slip, =1. free-slip) 59 61 / 60 62 !----------------------------------------------------------------------- 61 63 &namicethd ! ice thermodynamic 62 64 !----------------------------------------------------------------------- 63 ! hmelt : maximum melting at the bottom64 ! hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere65 ! caution 1.0, 1.0 best value to be used!!! (gilles G.)66 ! hicmin : ice thickness corr. to max. energy stored in brine pocket67 ! hiclim : minimum ice thickness68 ! amax : maximum lead fraction69 ! swiqst : energy stored in brine pocket (=1) or not (=0)70 ! sbeta : numerical caracteritic of the scheme for diffusion in ice71 ! Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)72 ! parlat : percentage of energy used for lateral ablation73 ! hakspl : slope of distr. for Hakkinen-Mellor's lateral melting74 ! hibspl : slope of distribution for Hibler's lateral melting75 ! exld : exponent for leads-closure rate76 ! hakdif : coefficient for diffusions of ice and snow77 ! thth : threshold thickness for comp. of eq. thermal conductivity78 ! hnzst : thickness of the surf. layer in temp. computation79 ! parsub : switch for snow sublimation or not80 ! alphs : coefficient for snow density when snow ice formation81 !82 65 hmelt = -0.15 ! maximum melting at the bottom 83 66 hiccrit = 0.3 , 0.3 ! ice thickness for lateral accretion in the Northern (Southern) Hemisphere 67 ! ! (caution 1.0, 1.0 best value to be used!!! (gilles G.)) 84 68 hicmin = 0.2 ! ice thickness corr. to max. energy stored in brine pocket 85 69 hiclim = 0.05 ! minimum ice thickness … … 95 79 thth = 0.2 ! threshold thickness for comp. of eq. thermal conductivity 96 80 hnzst = 0.1 ! thickness of the surf. layer in temp. computation 97 parsub = 0.0 ! switch for snow sublimation or not81 parsub = 1.0 ! switch for snow sublimation or not 98 82 alphs = 1.0 ! coefficient for snow density when snow ice formation 83 / 84 !----------------------------------------------------------------------- 85 &namice_dmp ! damping of sea ice alone open boundaries 86 ! ! (hard coded damping area: check if it fit your config) 87 !----------------------------------------------------------------------- 88 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 89 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 90 sn_hicif = 'ice_damping', -1. , 'hicif' , .true. , .true. , 'yearly' , '' , '' 91 sn_cnf = 'ice_damping', -1. , 'frld' , .true. , .true. , 'yearly' , '' , '' 92 ! 93 cn_dir = './' ! root directory for the location of the runoff files 99 94 / 100 95 !----------------------------------------------------------------------- 101 96 &namicedia ! ice diagnostics 102 97 !----------------------------------------------------------------------- 103 fmtinf = '1PE13.5' ! format of the output values104 nfrinf = 4! number of variables written in one line105 ntmoy = 1! instantaneous values of ice evolution or averaging106 ninfo = 1! frequency of ouputs on file ice_evolu in case of averaging98 fmtinf = '1PE13.5' ! format of the output values 99 nfrinf = 4 ! number of variables written in one line 100 ntmoy = 1 ! instantaneous values of ice evolution or averaging 101 ninfo = 1 ! frequency of ouputs on file ice_evolu in case of averaging 107 102 / 108 103 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 109 104 &namiceout ! parameters for outputs 110 105 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 106 ! ! title of the field ! name ! units ! save ! multipl. ! additive ! 107 ! ! ! ! ! or not ! factor ! factor ! 108 field_1 = 'Snow thickness ', 'isnowthi', 'm ', 1 , 1.0 , 0.0 109 field_2 = 'Ice thickness ', 'iicethic', 'm ', 1 , 1.0 , 0.0 110 field_3 = 'Ice produced ', 'iiceprod', 'm/kt ', 1 , 1.0 , 0.0 111 field_4 = 'Ice concentration ', 'ileadfra', '% ', 1 , -1.0 , 1.0 112 field_5 = 'Ice temperature ', 'iicetemp', 'C ', 1 , 1.0 , -273.15 113 field_6 = 'Oceanic flux at the ice base ', 'ioceflxb', 'w/m2 ', 1 , 1.0 , 0.0 114 field_7 = 'Ice velocity u ', 'iicevelu', 'm/s ', 1 , 1.0 , 0.0 115 field_8 = 'Ice velocity v ', 'iicevelv', 'm/s ', 1 , 1.0 , 0.0 116 field_9 = 'Sea surface temperature ', 'isstempe', 'C ', 1 , 1.0 , -273.15 117 field_10 = 'Sea surface salinity ', 'isssalin', 'PSU ', 1 , 1.0 , 0.0 118 field_11 = 'Total flux at ocean surface ', 'iocetflx', 'w/m2 ', 1 , 1.0 , 0.0 119 field_12 = 'Solar flux at ocean surface ', 'iocesflx', 'w/m2 ', 1 , 1.0 , 0.0 120 field_13 = 'Non-solar flux at ocean surface ', 'iocwnsfl', 'w/m2 ', 1 , 1.0 , 0.0 121 field_14 = 'Salt flux at ocean surface ', 'iocesafl', 'kg/m2/kt', 1 , 1.0 , 0.0 122 field_15 = 'Wind stress u ', 'iocestru', 'Pa ', 1 , 1.0 , 0.0 123 field_16 = 'Wind stress v ', 'iocestrv', 'Pa ', 1 , 1.0 , 0.0 124 field_17 = 'Solar flux at ice/ocean surface ', 'iicesflx', 'w/m2 ', 1 , 1.0 , 0.0 125 field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2 ', 1 , 1.0 , 0.0 126 field_19 = 'Snow precipitation ', 'isnowpre', 'kg/day ', 1 , 1.0 , 0.0 127 ! 111 128 noumef = 19 ! number of fields 112 !113 ! ! title of the field ! name ! units ! save ! multipl. ! additive !114 ! ! ! ! ! or not ! factor ! factor !115 field_1 = 'Snow thickness ', 'isnowthi', 'm ', 1 , 1.0 , 0.0116 field_2 = 'Ice thickness ', 'iicethic', 'm ', 1 , 1.0 , 0.0117 field_3 = 'Ice produced ', 'iiceprod', 'm/kt ', 1 , 1.0 , 0.0118 field_4 = 'Ice concentration ', 'ileadfra', '% ', 1 , -1.0 , 1.0119 field_5 = 'Ice temperature ', 'iicetemp', 'C ', 1 , 1.0 , -273.15120 field_6 = 'Oceanic flux at the ice base ', 'ioceflxb', 'w/m2 ', 1 , 1.0 , 0.0121 field_7 = 'Ice velocity u ', 'iicevelu', 'm/s ', 1 , 1.0 , 0.0122 field_8 = 'Ice velocity v ', 'iicevelv', 'm/s ', 1 , 1.0 , 0.0123 field_9 = 'Sea surface temperature ', 'isstempe', 'C ', 1 , 1.0 , -273.15124 field_10 = 'Sea surface salinity ', 'isssalin', 'PSU ', 1 , 1.0 , 0.0125 field_11 = 'Total flux at ocean surface ', 'iocetflx', 'w/m2 ', 1 , 1.0 , 0.0126 field_12 = 'Solar flux at ocean surface ', 'iocesflx', 'w/m2 ', 1 , 1.0 , 0.0127 field_13 = 'Non-solar flux at ocean surface ', 'iocwnsfl', 'w/m2 ', 1 , 1.0 , 0.0128 field_14 = 'Salt flux at ocean surface ', 'iocesafl', 'kg/m2/kt', 1 , 1.0 , 0.0129 field_15 = 'Wind stress u ', 'iocestru', 'Pa ', 1 , 1.0 , 0.0130 field_16 = 'Wind stress v ', 'iocestrv', 'Pa ', 1 , 1.0 , 0.0131 field_17 = 'Solar flux at ice/ocean surface ', 'iicesflx', 'w/m2 ', 1 , 1.0 , 0.0132 field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2 ', 1 , 1.0 , 0.0133 field_19 = 'Snow precipitation ', 'isnowpre', 'kg/day ', 1 , 1.0 , 0.0134 129 / -
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_pisces_ORCA2
r1775 r2061 15 15 &nampisext ! air-sea exchange 16 16 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 17 ln_co2int = .false. ! read atm pco2 from a file (T) or constant (F) 17 18 atcco2 = 0.36886500E+03 ! AUTO - atmospheric pCO2 19 clname = 'atcco2.txt' ! Name of atm pCO2 file - ln_co2int = T 20 nn_offset = 0 ! Offset model-data start year - ln_co2int = T 21 ! ! If your model year is iyy, nn_offset=(years(1)-iyy) 22 ! ! then the first atmospheric CO2 record read is at years(1) 23 / 24 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 25 &nampisatm ! Atmospheric prrssure 26 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 27 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 28 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 29 sn_patm = 'presatm' , -1 , 'patm' , .true. , .true. , 'yearly' , '' , '' 30 cn_dir = './' ! root directory for the location of the dynamical files 31 ! 32 ln_presatm = .true. ! constant atmopsheric pressure (F) or from a file (T) 18 33 / 19 34 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 20 35 &nampisbio ! biological parameters 21 36 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 22 part = 0.85 ! part of calcite not dissolved in guts 23 nrdttrc = 1 ! time step frequency for biology 24 wsbio = 2. ! POC sinking speed 25 xkmort = 1.E-7 ! half saturation constant for mortality 26 ferat3 = 3.E-6 ! Fe/C in zooplankton 27 wsbio2 = 50. ! Big particles sinking speed 37 nrdttrc = 1 ! time step frequency for biology 38 wsbio = 2. ! POC sinking speed 39 xkmort = 2.E-7 ! half saturation constant for mortality 40 ferat3 = 10.E-6 ! Fe/C in zooplankton 41 wsbio2 = 30. ! Big particles sinking speed 42 niter1max = 1 ! Maximum number of iterations for POC 43 niter2max = 1 ! Maximum number of iterations for GOC 28 44 / 29 45 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 30 46 &nampislim ! parameters for nutrient limitations 31 47 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 32 conc0 = 2.e-6 ! Phosphate half saturation 33 conc1 = 10E-6 ! Phosphate half saturation for diatoms 34 conc2 = 0.02E-9 ! Iron half saturation for phyto 35 conc2m = 0.08E-9 ! Max iron half saturation for phyto 36 conc3 = 0.1E-9 ! Iron half saturation for diatoms 37 conc3m = 0.4E-9 ! Maxi iron half saturation for diatoms 48 concnno3 = 1.e-6 ! Nitrate half saturation of nanophytoplankton 49 concdno3 = 3.E-6 ! Phosphate half saturation for diatoms 38 50 concnnh4 = 1.E-7 ! NH4 half saturation for phyto 39 concdnh4 = 5.E-7 ! NH4 half saturation for diatoms 51 concdnh4 = 3.E-7 ! NH4 half saturation for diatoms 52 concnfer = 1.E-9 ! Iron half saturation for phyto 53 concdfer = 3.E-9 ! Iron half saturation for diatoms 54 concbfe = 1.E-11 ! Half-saturation for Fe limitation of Bacteria 55 concbnh4 = 2.5E-8 ! NH4 half saturation for phyto 56 concbno3 = 2.5E-7 ! Phosphate half saturation for diatoms 57 xsizedia = 1.E-6 ! Minimum size criteria for diatoms 58 xsizephy = 1.E-6 ! Minimum size criteria for phyto 59 xsizern = 3.0 ! Size ratio for nanophytoplankton 60 xsizerd = 3.0 ! Size ratio for diatoms 40 61 xksi1 = 2.E-6 ! half saturation constant for Si uptake 41 xksi2 = 3.33E-6 ! half saturation constant for Si/C62 xksi2 = 20E-6 ! half saturation constant for Si/C 42 63 xkdoc = 417.E-6 ! half-saturation constant of DOC remineralization 43 caco3r = 0.15 ! mean rain ratio 44 / 64 qnfelim = 7.E-6 ! Optimal quota of phyto 65 qdfelim = 7.E-6 ! Optimal quota of diatoms 66 caco3r = 0.3 ! mean rain ratio 67 / 68 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 69 &nampisopt ! parameters for optics 70 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 71 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 72 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 73 sn_par = 'par.orca' , 24 , 'fr_par' , .true. , .true. , 'yearly' , '' , '' 74 cn_dir = './' ! root directory for the location of the dynamical files 75 ln_varpar = .true. ! boolean for PAR variable 76 parlux = 0.43 ! Fraction of shortwave as PAR 77 / 45 78 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 46 79 &nampisprod ! parameters for phytoplankton growth 47 80 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 48 pislope = 4. ! P-I slope49 pislope2 = 4. ! P-I slope for diatoms81 pislope = 2. ! P-I slope 82 pislope2 = 2. ! P-I slope for diatoms 50 83 excret = 0.05 ! excretion ratio of phytoplankton 51 84 excret2 = 0.05 ! excretion ratio of diatoms 85 ln_newprod = .true. ! Enable new parame. of production (T/F) 86 bresp = 0.00333 ! Basal respiration rate 52 87 chlcnm = 0.033 ! Minimum Chl/C in nanophytoplankton 53 88 chlcdm = 0.05 ! Minimum Chl/C in diatoms 54 fecnm = 10E-6 ! Maximum Fe/C in nanophytoplankton 55 fecdm = 15E-6 ! Minimum Fe/C in diatoms 56 grosip = 0.151 ! mean Si/C ratio 89 chlcmin = 0.004 ! Maximum Chl/c in phytoplankton 90 fecnm = 40E-6 ! Maximum Fe/C in nanophytoplankton 91 fecdm = 40E-6 ! Minimum Fe/C in diatoms 92 grosip = 0.159 ! mean Si/C ratio 57 93 / 58 94 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 59 95 &nampismort ! parameters for phytoplankton sinks 60 96 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 61 wchl = 0.001 ! quadratic mortality of phytoplankton 62 wchld = 0.02 ! maximum quadratic mortality of diatoms 97 wchl = 0.01 ! quadratic mortality of phytoplankton 98 wchld = 0.01 ! maximum quadratic mortality of diatoms 99 wchldm = 0.03 ! maximum quadratic mortality of diatoms 63 100 mprat = 0.01 ! phytoplankton mortality rate 64 101 mprat2 = 0.01 ! Diatoms mortality rate 65 mpratm = 0.01 ! Phytoplankton minimum mortality rate66 102 / 67 103 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 68 104 &nampismes ! parameters for mesozooplankton 69 105 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 70 grazrat2 = 0.7 ! maximal mesozoo grazing rate 106 part2 = 0.75 ! part of calcite not dissolved in mesozoo guts 107 grazrat2 = 0.75 ! maximal mesozoo grazing rate 71 108 resrat2 = 0.005 ! exsudation rate of mesozooplankton 72 109 mzrat2 = 0.03 ! mesozooplankton mortality rate 73 110 xprefc = 1. ! zoo preference for phyto 74 xprefp = 0. 2! zoo preference for POC111 xprefp = 0.3 ! zoo preference for POC 75 112 xprefz = 1. ! zoo preference for zoo 76 xprefpoc = 0.2 ! zoo preference for poc 113 xprefpoc = 0.3 ! zoo preference for poc 114 xthresh2zoo = 1E-8 ! zoo feeding threshold for mesozooplankton 115 xthresh2dia = 1E-8 ! diatoms feeding threshold for mesozooplankton 116 xthresh2phy = 1E-8 ! nanophyto feeding threshold for mesozooplankton 117 xthresh2poc = 1E-8 ! poc feeding threshold for mesozooplankton 118 xthresh2 = 3E-7 ! Food threshold for grazing 77 119 xkgraz2 = 20.E-6 ! half sturation constant for meso grazing 78 epsher2 = 0.3 3 ! Efficicency of Mesozoo growth120 epsher2 = 0.3 ! Efficicency of Mesozoo growth 79 121 sigma2 = 0.6 ! Fraction of mesozoo excretion as DOM 80 122 unass2 = 0.3 ! non assimilated fraction of P by mesozoo 81 grazflux = 5.e3 ! flux-feeding rate123 grazflux = 2.e3 ! flux-feeding rate 82 124 / 83 125 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 84 126 &nampiszoo ! parameters for microzooplankton 85 127 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 86 grazrat = 4.0 ! maximal zoo grazing rate 128 part = 0.5 ! part of calcite not dissolved in microzoo gutsa 129 grazrat = 3.0 ! maximal zoo grazing rate 87 130 resrat = 0.03 ! exsudation rate of zooplankton 88 mzrat = 0.0 89 xpref2c = 0. 0 ! Microzoo preference for POM90 xpref2p = 0.5! Microzoo preference for Nanophyto131 mzrat = 0.004 ! zooplankton mortality rate 132 xpref2c = 0.1 ! Microzoo preference for POM 133 xpref2p = 1. ! Microzoo preference for Nanophyto 91 134 xpref2d = 0.5 ! Microzoo preference for Diatoms 92 xkgraz = 20.E-6 ! half sturation constant for grazing 93 epsher = 0.33 ! Efficiency of microzoo growth 135 xthreshdia = 1.E-8 ! Diatoms feeding threshold for microzooplankton 136 xthreshphy = 1.E-8 ! Nanophyto feeding threshold for microzooplankton 137 xthreshpoc = 1.E-8 ! POC feeding threshold for microzooplankton 138 xthresh = 3.E-7 ! Food threshold for feeding 139 xkgraz = 20.E-6 ! half sturation constant for grazing 140 epsher = 0.3 ! Efficiency of microzoo growth 94 141 sigma1 = 0.6 ! Fraction of microzoo excretion as DOM 95 142 unass = 0.3 ! non assimilated fraction of phyto by zoo 96 143 / 97 144 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 145 &nampisfer ! parameters for iron chemistry 146 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 147 ln_fechem = .false. ! complex iron chemistry ( T/F ) 148 ln_ligvar = .true. ! variable ligand concentration 149 xlam1 = 0.005 ! scavenging rate of Iron 150 xlamdust = 150.0 ! Scavenging rate of dust 151 ligand = 0.6E-9 ! Ligands concentration 152 / 153 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 98 154 &nampisrem ! parameters for remineralization 99 155 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 100 xremik = 0.3 156 xremik = 0.35 ! remineralization rate of DOC 101 157 xremip = 0.025 ! remineralisation rate of POC 102 158 nitrif = 0.05 ! NH4 nitrification rate 103 xsirem = 0.015 ! remineralization rate of Si 104 xlam1 = 0.005 ! scavenging rate of Iron 105 oxymin = 1.E-6 ! Half-saturation constant for anoxia 159 xsirem = 0.003 ! remineralization rate of Si 160 xsiremlab = 0.03 ! fast remineralization rate of Si 161 xsilab = 0.5 ! Fraction of labile biogenic silica 162 oxymin = 1.E-6 ! Half-saturation constant for anoxia 106 163 / 107 164 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 108 165 &nampiscal ! parameters for Calcite chemistry 109 166 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 110 kdca = 0.327e3! calcite dissolution rate constant (1/time)167 kdca = 6. ! calcite dissolution rate constant (1/time) 111 168 nca = 1. ! order of dissolution reaction (dimensionless) 112 169 / 113 170 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 114 &nampissed ! parameters for inputs deposition 115 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 116 ln_dustfer = .true. ! boolean for dust input from the atmosphere 117 ln_river = .true. ! AUTO boolean for river input of nutrients 171 &nampissbc ! parameters for inputs deposition 172 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 173 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 174 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 175 sn_dust = 'dust.orca' , -1 , 'dust' , .true. , .true. , 'yearly' , '' , '' 176 sn_solub = 'solubility.orca' , -12 , 'solubility1' , .false. , .true. , 'yearly' , '' , '' 177 sn_riverdic = 'river.orca' , 120 , 'riverdic' , .true. , .true. , 'yearly' , '' , '' 178 sn_riverdoc = 'river.orca' , 120 , 'riverdoc' , .true. , .true. , 'yearly' , '' , '' 179 sn_riverdin = 'river.orca' , 120 , 'riverdin' , .true. , .true. , 'yearly' , '' , '' 180 sn_riverdon = 'river.orca' , 120 , 'riverdon' , .true. , .true. , 'yearly' , '' , '' 181 sn_riverdip = 'river.orca' , 120 , 'riverdip' , .true. , .true. , 'yearly' , '' , '' 182 sn_riverdop = 'river.orca' , 120 , 'riverdop' , .true. , .true. , 'yearly' , '' , '' 183 sn_riverdsi = 'river.orca' , 120 , 'riverdsi' , .true. , .true. , 'yearly' , '' , '' 184 sn_ndepo = 'ndeposition.orca', -12 , 'ndep' , .false. , .true. , 'yearly' , '' , '' 185 sn_ironsed = 'bathy.orca' , -12 , 'bathy' , .false. , .true. , 'yearly' , '' , '' 186 sn_hydrofe = 'hydrofe' , -12 , 'epsdb' , .false. , .true. , 'yearly' , '' , '' 187 ! 188 cn_dir = './' ! root directory for the location of the dynamical files 189 ln_dust = .true. ! boolean for dust input from the atmosphere 190 ln_solub = .true. ! boolean for variable solubility of atm. Iron 191 ln_river = .true. ! boolean for river input of nutrients 118 192 ln_ndepo = .true. ! boolean for atmospheric deposition of N 119 ln_sedinput = .true. ! boolean for Fe input from sediments 120 sedfeinput = 1E-9 ! Coastal release of Iron 121 dustsolub = 0.014 ! Solubility of the dust 193 ln_ironsed = .true. ! boolean for Fe input from sediments 194 ln_ironice = .true. ! boolean for Fe input from sea ice 195 ln_hydrofe = .false. ! boolean for from hydrothermal vents 196 sedfeinput = 1.e-9 ! Coastal release of Iron 197 dustsolub = 0.02 ! Solubility of the dust 198 wdust = 2.0 ! Dust sinking speed 199 icefeinput = 15.e-9 ! Iron concentration in sea ice 200 nitrfix = 1.e-7 ! Nitrogen fixation rate 201 diazolight = 50. ! Diazotrophs sensitivity to light (W/m2) 202 concfediaz = 1.e-10 ! Diazotrophs half-saturation Cste for Iron 203 hratio = 9.e+5 ! Fe to 3He ratio assumed for vent iron supply 122 204 / 123 205 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' … … 126 208 xkr_eta = 1.17 ! Sinking exponent 127 209 xkr_zeta = 2.28 ! N content exponent 210 xkr_ncontent = 5.7E-6 ! N content factor 128 211 xkr_mass_min = 0.0002 ! Minimum mass for Aggregates 129 212 xkr_mass_max = 1. ! Maximum mass for Aggregates … … 140 223 / 141 224 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 142 &nampisdia ! additional 2D/3D tracers diagnostics ("key_trc_diaadd") 143 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 144 nwritedia = 5475 ! time step frequency for tracers diagnostics 145 ! 225 &nampisdia ! additional 2D/3D tracers diagnostics 226 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 146 227 ! ! name ! title of the field ! units ! 147 228 ! ! ! ! ! … … 175 256 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 176 257 ln_pisdmp = .true. ! Relaxation fo some tracers to a mean value 177 ln_pisclo = .true. ! Restoring of tracer to initial value on closed sea ("key_dtatrc") 178 / 258 nn_pisdmp = 5475 ! Frequency of Relaxation 259 / 260 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 261 &nampismass ! Mass conservation 262 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 263 ln_check_mass = .true. ! Check mass conservation 264 / -
CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_top_ORCA2
r1775 r2061 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/TOP1 : 1 - tracer definition (namtoptrc) 3 !! namelists 2 - dynamical tracer trends (namtoptrd) 4 !! 6 - tracer advection (namtopadv) 5 !! 7 - tracer bottom boundary (namtopbbl) 6 !! 8 - tracer lateral diffusion (namtopldf) 7 !! 3 - tracer vertical physics (namtopzdf) 8 !! 9 - tracer newtonian damping (namtopdmp) 2 !! NEMO/TOP1 : 1 - tracer definition (namtrc ) 3 !! 2 - tracer data initialisation (namtrc_dta) 4 !! 3 - tracer advection (namtrc_adv) 5 !! 4 - tracer lateral diffusion (namtrc_ldf) 6 !! 5 - tracer vertical physics (namtrc_zdf) 7 !! 6 - tracer newtonian damping (namtrc_dmp) 8 !! 7 - dynamical tracer trends (namtrc_trd) 9 !! 8 - tracer output diagonstics (namtrc_dia) 9 10 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 10 11 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 11 &namt optrc ! tracers definition12 &namtrc ! tracers definition 12 13 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 13 ndttrc = 1 ! time step frequency for passive tracers 14 nwritetrc = 5475 ! time step frequency for tracer outputs 15 ln_rsttr = .false. ! AUTO start from a restart file (T) or not (F) 16 nrsttr = 0 ! AUTO restart control = 0 initial time step is not compared to the restart file value 17 ! = 1 do not use the value in the restart file 18 ! = 2 calendar parameters read in the restart file 19 cn_trcrst_in = "restart_trc" ! suffix of pass. tracer restart name (input) 20 cn_trcrst_out = "restart_trc" ! suffix of pass. tracer restart name (output) 14 nn_dttrc = 1 ! time step frequency for passive sn_tracers 15 nn_writetrc = 5475 ! time step frequency for sn_tracer outputs 16 ln_top_euler = .true. ! use Euler time-stepping for TOP 17 ln_rsttr = .false. ! AUTO start from a restart file (T) or not (F) 18 nn_rsttr = 0 ! AUTO restart control = 0 initial time step is not compared to the restart file value 19 ! = 1 do not use the value in the restart file 20 ! = 2 calendar parameters read in the restart file 21 cn_trcrst_in = "restart_trc" ! suffix of pass. sn_tracer restart name (input) 22 cn_trcrst_out = "restart_trc" ! suffix of pass. sn_tracer restart name (output) 23 ln_trcdta = .true. ! Initialisation from data input file (T) or not (F) 24 ln_trcdmp = .false. ! add a damping termn (T) or not (F) 25 ln_trcdmp_clo = .true. ! restoring on close sea + baltic sea (T) or not (F) 21 26 ! 22 ! ! name ! title of the field ! units! initial data ! save !23 ! ! ! !! from file ! or not !24 ! ! ! !! or not ! !25 tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'molC/L', .true. , .true.26 tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', 'eq/L ' , .true. , .true.27 tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', 'molO2/L' , .true. , .true.28 tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', 'molC/L', .false. , .true.29 tracer(5) = 'PO4 ' , 'Phosphate Concentration ', 'molC/L', .true. , .true.30 tracer(6) = 'POC ' , 'Small organic carbon Concentration ', 'molC/L', .false. , .true.31 tracer(7) = 'Si ' , 'Silicate Concentration ', 'molSi/L' , .true. , .true.32 tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', 'molC/L', .false. , .true.33 tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', 'molC/L', .false. , .true.34 tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', 'molC/L' , .false., .true.35 tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', 'molC/L', .false. , .true.36 tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', 'molC/L', .false. , .true.37 tracer(13) = 'BSi ' , 'Diatoms Silicate Concentration ', 'molC/L', .false. , .true.38 tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', 'molFe/L' , .true. , .true.39 tracer(15) = 'BFe ' , 'Big iron particles Concentration ', 'molFe/L' , .false. , .true.40 tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', 'molC/L', .false. , .true.41 tracer(17) = 'SFe ' , 'Small iron particles Concentration ', 'molFe/L' , .false. , .true.42 tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', 'molFe/L' , .false. , .true.43 tracer(19) = 'DSi ' , 'Sinking biogenic Silicate Concentration', 'molC/L', .false. , .true.44 tracer(20) = 'NFe ' , 'Nano iron Concentration ', 'molFe/L' , .false. , .true.45 tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', 'gChl/L', .false. , .true.46 tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', 'gChl/L', .false. , .true.47 tracer(23) = 'NO3 ' , 'Nitrates Concentration ', 'molC/L', .true. , .true.48 tracer(24) = 'NH4 ' , 'Ammonium Concentration ', 'molC/L', .false. , .true.27 ! ! name ! title of the field ! initial data ! initial data ! save ! 28 ! ! ! ! units ! from file ! or not ! 29 ! ! ! ! ! or not ! ! 30 sn_tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'mol-C/L' , .true. , .true. 31 sn_tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', 'eq/L ' , .true. , .true. 32 sn_tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', 'mol-C/L' , .true. , .true. 33 sn_tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', 'mol-C/L' , .false. , .true. 34 sn_tracer(5) = 'PO4 ' , 'Phosphate Concentration ', 'mol-C/L' , .true. , .true. 35 sn_tracer(6) = 'POC ' , 'Small organic carbon Concentration ', 'mol-C/L' , .false. , .true. 36 sn_tracer(7) = 'Si ' , 'Silicate Concentration ', 'mol-C/L' , .true. , .true. 37 sn_tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', 'mol-C/L' , .false. , .true. 38 sn_tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', 'mol-C/L' , .false. , .true. 39 sn_tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', 'mol-C/L' , .true. , .true. 40 sn_tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', 'mol-C/L' , .false. , .true. 41 sn_tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', 'mol-C/L' , .false. , .true. 42 sn_tracer(13) = 'DSi ' , 'Diatoms Silicate Concentration ', 'mol-C/L' , .false. , .true. 43 sn_tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', 'mol-C/L' , .true. , .true. 44 sn_tracer(15) = 'BFe ' , 'Big iron particles Concentration ', 'mol-C/L' , .false. , .true. 45 sn_tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', 'mol-C/L' , .false. , .true. 46 sn_tracer(17) = 'SFe ' , 'Small iron particles Concentration ', 'mol-C/L' , .false. , .true. 47 sn_tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', 'mol-C/L' , .false. , .true. 48 sn_tracer(19) = 'GSi ' , 'Sinking biogenic Silicate Concentration', 'mol-C/L' , .false. , .true. 49 sn_tracer(20) = 'NFe ' , 'Nano iron Concentration ', 'mol-C/L' , .false. , .true. 50 sn_tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', 'mol-C/L' , .false. , .true. 51 sn_tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', 'mol-C/L' , .false. , .true. 52 sn_tracer(23) = 'NO3 ' , 'Nitrates Concentration ', 'mol-C/L' , .true. , .true. 53 sn_tracer(24) = 'NH4 ' , 'Ammonium Concentration ', 'mol-C/L' , .false. , .true. 49 54 / 50 55 !----------------------------------------------------------------------- 51 &namt opadv ! advection scheme for passive tracer56 &namtrc_dta ! Initialisation from data input file 52 57 !----------------------------------------------------------------------- 53 ln_trcadv_cen2 = .false. ! 2nd order centered scheme 54 ln_trcadv_tvd = .false. ! TVD scheme 55 ln_trcadv_muscl = .true. ! MUSCL scheme 56 ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 57 ln_trcadv_smolar = .false. ! SMOLAR scheme 58 rsc = 1. ! tuning coefficient for Smol-Car. scheme 59 ncortrc = 1 ! number of corrective phases for Smol-Car. scheme 60 crosster = .false. ! computes Smol-Car crossterms (T) or not (F) 58 ! 59 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! 60 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 61 sn_trcdta(1) = 'data_DIC_nomask' , -12 , 'DIC' , .false. , .true. , 'yearly' , '' , '' 62 sn_trcdta(2) = 'data_Alkalini_nomask' , -12 , 'Alkalini', .false. , .true. , 'yearly' , '' , '' 63 sn_trcdta(3) = 'data_O2_nomask' , -1 , 'O2' , .true. , .true. , 'yearly' , '' , '' 64 sn_trcdta(5) = 'data_PO4_nomask' , -1 , 'PO4' , .true. , .true. , 'yearly' , '' , '' 65 sn_trcdta(7) = 'data_Si_nomask' , -1 , 'Si' , .true. , .true. , 'yearly' , '' , '' 66 sn_trcdta(10) = 'data_DOC_nomask' , -12 , 'DOC' , .false. , .true. , 'yearly' , '' , '' 67 sn_trcdta(14) = 'data_Fer_nomask' , -12 , 'Fer' , .false. , .true. , 'yearly' , '' , '' 68 sn_trcdta(23) = 'data_NO3_nomask' , -1 , 'NO3' , .true. , .true. , 'yearly' , '' , '' 69 ! 70 cn_dir = './' ! root directory for the location of the data files 71 rn_trfac(1) = 1.0e-06 ! multiplicative factor 72 rn_trfac(2) = 1.0e-06 ! - - - - 73 rn_trfac(3) = 44.6e-06 ! - - - - 74 rn_trfac(5) = 122.0e-06 ! - - - - 75 rn_trfac(7) = 1.0e-06 ! - - - - 76 rn_trfac(10) = 1.0 ! - - - - 77 rn_trfac(14) = 1.0 ! - - - - 78 rn_trfac(23) = 7.6e-06 ! - - - - 61 79 / 62 80 !----------------------------------------------------------------------- 63 &namt opbbl ! bottom boundary layerscheme for passive tracer81 &namtrc_adv ! advection scheme for passive tracer 64 82 !----------------------------------------------------------------------- 65 atrcbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] 83 ln_trcadv_cen2 = .false. ! 2nd order centered scheme 84 ln_trcadv_tvd = .false. ! TVD scheme 85 ln_trcadv_muscl = .true. ! MUSCL scheme 86 ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 87 ln_trcadv_ubs = .false. ! UBS scheme 88 ln_trcadv_qck = .false. ! QUICKEST scheme 89 ln_trcadv_msc_ups = .false. ! use upstream scheme within muscl 66 90 / 67 91 !----------------------------------------------------------------------- 68 &namt opldf ! lateral diffusion scheme for passive tracer92 &namtrc_ldf ! lateral diffusion scheme for passive tracer 69 93 !----------------------------------------------------------------------- 70 94 ln_trcldf_diff = .true. ! performs lateral diffusion (T) or not (F) … … 77 101 ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") 78 102 ! ! Coefficient 79 ahtrc0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 80 ahtrb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 81 aeivtr0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_trcldf_eiv") 82 trcrat = 1. ! ratio betweeen passive and active tracer diffusion coeff 103 rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 104 rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 83 105 / 84 106 !----------------------------------------------------------------------- 85 &namt opzdf ! vertical physics107 &namtrc_zdf ! vertical physics 86 108 !----------------------------------------------------------------------- 87 109 ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping 88 n _trczdf_exp= 3 ! number of sub-timestep for ln_trczdfexp=T110 nn_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T 89 111 / 90 112 !----------------------------------------------------------------------- 91 &namt oprad ! treatment of negative concentrations113 &namtrc_rad ! treatment of negative concentrations 92 114 !----------------------------------------------------------------------- 93 115 ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F) 94 116 / 95 117 !----------------------------------------------------------------------- 96 &namt opdmp ! passive tracer newtonian damping ('key_trcdmp')118 &namtrc_dmp ! passive tracer newtonian damping 97 119 !----------------------------------------------------------------------- 98 ndmptr = 20 ! type of damping in passive tracers 99 ! ='latitude', damping poleward of 'ndmp' degrees and function 100 ! of the distance-to-coast. Red and Med Seas as ndmptr=-1 101 ! =-1 damping only in Med and Red Seas 102 ndmpftr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 103 nmldmptr = 1 ! type of damping: =0 damping throughout the water column 104 ! =1 no damping in the mixed layer defined by avt >5cm2/s ) 105 ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 106 sdmptr = 50. ! surface time scale for internal damping (days) 107 bdmptr = 360. ! bottom time scale for internal damping (days) 108 hdmptr = 800. ! depth of transition between sdmptr and bdmptr (meters) 120 nn_hdmp_tr = -1 ! horizontal shape =-1, damping in Med and Red Seas only 121 ! =XX, damping poleward of XX degrees (XX>0) 122 ! + F(distance-to-coast) + Red and Med Seas 123 nn_zdmp_tr = 1 ! vertical shape =0 damping throughout the water column 124 ! =1 no damping in the mixing layer (kz criteria) 125 ! =2 no damping in the mixed layer (rho crieria) 126 rn_surf_tr = 50. ! surface time scale of damping [days] 127 rn_bot_tr = 360. ! bottom time scale of damping [days] 128 rn_dep_tr = 800. ! depth of transition between rn_surf and rn_bot [meters] 129 nn_file_tr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 109 130 / 110 131 !----------------------------------------------------------------------- 111 &namt optrd ! diagnostics on tracer trends132 &namtrc_trd ! diagnostics on tracer trends ('key_trdtrc') 112 133 ! or mixed-layer trends ('key_trdmld_trc') 113 134 !---------------------------------------------------------------------- 114 n trd_trc= 5475 ! time step frequency and tracers trends115 n ctls_trc= 0 ! control surface type in mixed-layer trends (0,1 or n<jpk)116 ucf_trc= 1 ! unit conversion factor (=1 -> /seconds ; =86400. -> /day)135 nn_trd_trc = 5475 ! time step frequency and tracers trends 136 nn_ctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 137 rn_ucf_trc = 1 ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 117 138 ln_trdmld_trc_restart = .false. ! restart for ML diagnostics 118 139 ln_trdmld_trc_instant = .true. ! flag to diagnose trends of instantantaneous or mean ML T/S 119 luttrd(1) = .true. 120 luttrd(2) = .true. 121 luttrd(3) = .false. 122 luttrd(4) = .false. 123 luttrd(5) = .false. 124 luttrd(6) = .false. 125 luttrd(7) = .false. 126 luttrd(8) = .false. 127 luttrd(9) = .false. 128 luttrd(10) = .false. 129 luttrd(11) = .false. 130 luttrd(12) = .false. 131 luttrd(13) = .false. 132 luttrd(14) = .false. 133 luttrd(15) = .false. 134 luttrd(16) = .false. 135 luttrd(17) = .false. 136 luttrd(18) = .false. 137 luttrd(19) = .false. 138 luttrd(20) = .false. 139 luttrd(21) = .false. 140 luttrd(22) = .false. 141 luttrd(23) = .true. 142 luttrd(24) = .false. 140 ln_trdtrc(1) = .true. 141 ln_trdtrc(2) = .true. 142 ln_trdtrc(23) = .true. 143 143 / 144 !----------------------------------------------------------------------- 145 &namtrc_dia ! parameters for passive tracer additional diagnostics 146 !---------------------------------------------------------------------- 147 ln_diatrc = .true. ! save additional diag. (T) or not (F) 148 nn_writedia = 5475 ! time step frequency for diagnostics 149 /
Note: See TracChangeset
for help on using the changeset viewer.