Changeset 7467 for branches/2016/dev_merge_2016
- Timestamp:
- 2016-12-06T19:07:53+01:00 (8 years ago)
- Location:
- branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES
- Files:
-
- 12 added
- 4 deleted
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/EXP00/field_def.xml
r7412 r7467 42 42 <field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" /> 43 43 <field id="sbt" long_name="sea bottom temperature" unit="degC" /> 44 <field id="tosmint" long_name="vertical integral of temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" /> 44 45 <field id="sst_wl" long_name="Delta SST of warm layer" unit="degC" /> 45 46 <field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" /> … … 50 51 <field id="sssmax" long_name="max of sea surface salinity" field_ref="sss" operation="maximum" /> 51 52 <field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum" /> 52 <field id="sbs" long_name="sea bottom salinity" unit="1e-3" /> 53 <field id="sbs" long_name="sea bottom salinity" unit="0.001" /> 54 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 53 55 54 56 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> … … 95 97 96 98 <!-- variables available with key_diaar5 --> 97 <field id="botpres" long_name="Pressure at sea floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" /> 99 <field id="botpres" long_name="Sea Water Pressure at Sea Floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" /> 100 <field id="sshdyn" long_name="dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid" unit="m" /> 101 <field id="sshdyn2" long_name="square of dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid_squared" unit="m2" > sshdyn * sshdyn </field> 102 <field id="tnpeo" long_name="Tendency of ocean potential energy content" unit="W/m2" /> 98 103 99 104 <!-- variables available with key_vvl --> … … 185 190 <field id="empmr" long_name="Net Upward Water Flux" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" /> 186 191 <field id="empbmr" long_name="Net Upward Water Flux at pre. tstep" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" /> 192 <field id="emp_oce" long_name="Evap minus Precip over ocean" standard_name="evap_minus_precip_over_sea_water" unit="kg/m2/s" /> 193 <field id="emp_ice" long_name="Evap minus Precip over ice" standard_name="evap_minus_precip_over_sea_ice" unit="kg/m2/s" /> 187 194 <field id="saltflx" long_name="Downward salt flux" unit="1e-3/m2/s" /> 188 195 <field id="fmmflx" long_name="Water flux due to freezing/melting" unit="kg/m2/s" /> … … 205 212 <!-- * variable related to ice shelf forcing * --> 206 213 <field id="fwfisf" long_name="Ice shelf melting" unit="Kg/m2/s" /> 207 <field id="qisf" long_name="Ice Shelf Heat Flux" unit="W/m2" /> 214 <field id="fwfisf3d" long_name="Ice shelf melting" unit="kg/m2/s" grid_ref="grid_T_3D" /> 215 <field id="qlatisf" long_name="Ice shelf latent heat flux" unit="W/m2" /> 216 <field id="qlatisf3d" long_name="Ice shelf latent heat flux" unit="W/m2" grid_ref="grid_T_3D" /> 217 <field id="qhcisf" long_name="Ice shelf heat content flux" unit="W/m2" /> 218 <field id="qhcisf3d" long_name="Ice shelf heat content flux" unit="W/m2" grid_ref="grid_T_3D" /> 208 219 <field id="isfgammat" long_name="transfert coefficient for isf (temperature) " unit="m/s" /> 209 220 <field id="isfgammas" long_name="transfert coefficient for isf (salinity) " unit="m/s" /> … … 239 250 <field id="hflx_rain_cea" long_name="heat flux due to rainfall" standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> 240 251 <field id="hflx_evap_cea" long_name="heat flux due to evaporation" standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water" unit="W/m2" /> 241 <field id="hflx_snow_cea" long_name="heat flux due to snow falling over ice-free ocean" standard_name="heat_flux_into_sea_water_due_to_snow_thermodynamics" unit="W/m2" /> 252 <field id="hflx_snow_cea" long_name="heat flux due to snow falling" standard_name="heat_flux_onto_ocean_and_ice_due_to_snow_thermodynamics" unit="W/m2" /> 253 <field id="hflx_snow_ai_cea" long_name="heat flux due to snow falling over ice" standard_name="heat_flux_onto_ice_due_to_snow_thermodynamics" unit="W/m2" /> 254 <field id="hflx_snow_ao_cea" long_name="heat flux due to snow falling over ice-free ocean" standard_name="heat_flux_onto_sea_water_due_to_snow_thermodynamics" unit="W/m2" /> 242 255 <field id="hflx_ice_cea" long_name="heat flux due to ice thermodynamics" standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="W/m2" /> 243 256 <field id="hflx_rnf_cea" long_name="heat flux due to runoffs" standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> … … 277 290 <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" /> 278 291 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 292 <field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" /> 293 <field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" /> 279 294 280 295 <field id="iceconc" long_name="ice concentration" standard_name="sea_ice_area_fraction" unit="%" /> … … 291 306 <field id="micesalt" long_name="Mean ice salinity" unit="1e-3" /> 292 307 <field id="miceage" long_name="Mean ice age" unit="years" /> 308 <field id="alb_ice" long_name="Mean albedo over sea ice" unit="" /> 309 <field id="albedo" long_name="Mean albedo over sea ice and ocean" unit="" /> 293 310 294 311 <field id="iceage_cat" long_name="Ice age for categories" unit="days" axis_ref="ncatice" /> … … 313 330 <field id="icevolu" long_name="ice volume" unit="m" /> 314 331 <field id="snowvol" long_name="snow volume" unit="m" /> 332 <field id="tau_icebfr" long_name="ice friction on ocean bottom for landfast ice" unit="N/m2" /> 315 333 316 334 <field id="icetrp" long_name="ice volume transport" unit="m/day" /> … … 326 344 <field id="sfxbom" long_name="salt flux from bot melt" unit="1e-3*kg/m2/day" /> 327 345 <field id="sfxsum" long_name="salt flux from surf melt" unit="1e-3*kg/m2/day" /> 346 <field id="sfxlam" long_name="salt flux from lateral melt" unit="1e-3*kg/m2/day" /> 328 347 <field id="sfxsni" long_name="salt flux from snow-ice formation" unit="1e-3*kg/m2/day" /> 329 348 <field id="sfxopw" long_name="salt flux from open water ice formation" unit="1e-3*kg/m2/day" /> 349 <field id="sfxsub" long_name="salt flux from sublimation" unit="1e-3*kg/m2/day" /> 330 350 <field id="sfx" long_name="salt flux total" unit="1e-3*kg/m2/day" /> 331 351 … … 335 355 <field id="vfxsni" long_name="daily snowice ice prod." unit="m/day" /> 336 356 <field id="vfxsum" long_name="surface melt" unit="m/day" /> 357 <field id="vfxlam" long_name="lateral melt" unit="m/day" /> 337 358 <field id="vfxbom" long_name="bottom melt" unit="m/day" /> 338 359 <field id="vfxres" long_name="daily resultant ice prod./melting from limupdate" unit="m/day" /> … … 340 361 <field id="vfxsnw" long_name="snw melt/growth" unit="m/day" /> 341 362 <field id="vfxsub" long_name="snw sublimation" unit="m/day" /> 363 <field id="vfxsub_err" long_name="excess of snw sublimation sent to ocean" unit="m/day" /> 342 364 <field id="vfxspr" long_name="snw precipitation on ice" unit="m/day" /> 365 <field id="vfxthin" long_name="daily thermo ice prod. for thin ice(20cm) + open water" unit="m/day" /> 343 366 344 367 <field id="afxtot" long_name="area tendency (total)" unit="day-1" /> … … 384 407 <field id="e3u_0" long_name="Initial U-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_U_3D"/> 385 408 <field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" /> 386 <field id="ucli" long_name="uncorrected ocean current along i-axis" standard_name="uc_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" />387 <field id="ucli2" long_name="uncorrected2 ocean current along i-axis" standard_name="uc2_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" />388 <field id="un_b" long_name="un_b ocean barotropic current along i-axis" unit="m/s" />389 409 <field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" /> 390 410 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 391 411 <field id="ssu" long_name="ocean surface current along i-axis" unit="m/s" /> 392 <field id="wdlmtu" long_name="ocean current limiter along i-axis" unit="#" />393 412 <field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" /> 394 413 <field id="ubar" long_name="ocean barotropic current along i-axis" unit="m/s" /> … … 416 435 <field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" /> 417 436 418 < !-- variables available with key_diaar5 -->419 <field id="u_masstr " long_name="ocean eulerian mass transport along i-axis" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" />437 <field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" /> 438 <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport" unit="kg/s" /> 420 439 <field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" /> 421 440 <field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="1e-3*kg/s" /> 441 <field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W" /> 442 <field id="uadv_salttr" long_name="ocean advective salt transport along i-axis" standard_name="advectice_ocean_salt_x_transport" unit="1e-3*kg/s" /> 422 443 <field id="ueiv_heattr" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" /> 444 <field id="ueiv_salttr" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" /> 445 <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" grid_ref="grid_U_3D" /> 446 <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_U_3D" /> 423 447 <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W" /> 448 <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" /> 424 449 </field_group> 425 450 … … 462 487 <field id="v_heattr" long_name="ocean eulerian heat transport along j-axis" standard_name="ocean_heat_y_transport" unit="W" /> 463 488 <field id="v_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_y_transport" unit="1e-3*kg/s" /> 489 <field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W" /> 490 <field id="vadv_salttr" long_name="ocean advective salt transport along j-axis" standard_name="advectice_ocean_salt_y_transport" unit="1e-3*kg/s" /> 464 491 <field id="veiv_heattr" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" /> 492 <field id="veiv_salttr" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" /> 493 <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" grid_ref="grid_V_3D" /> 494 <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_V_3D" /> 465 495 <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W" /> 496 <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" /> 466 497 </field_group> 467 498 … … 519 550 520 551 <!-- available with ln_diahsb --> 521 <field id="bgtemper" long_name="drift in global mean temperature wrt timestep 1" standard_name="change_over_time_in_sea_water_potential_temperature" unit="degC" />552 <field id="bgtemper" long_name="drift in global mean temperature wrt timestep 1" standard_name="change_over_time_in_sea_water_potential_temperature" unit="degC" /> 522 553 <field id="bgsaline" long_name="drift in global mean salinity wrt timestep 1" standard_name="change_over_time_in_sea_water_practical_salinity" unit="1e-3" /> 523 <field id="bgheatco" long_name="drift in global mean heat content wrt timestep 1" unit="10^9J" /> 524 <field id="bgsaltco" long_name="drift in global mean salt content wrt timestep 1" unit="1e-3*m3" /> 554 <field id="bgheatco" long_name="drift in global mean heat content wrt timestep 1" unit="1.e20J" /> 555 <field id="bgheatfx" long_name="drift in global mean heat flux wrt timestep 1" unit="W/m2" /> 556 <field id="bgsaltco" long_name="drift in global mean salt content wrt timestep 1" unit="1e-3*km3" /> 525 557 <field id="bgvolssh" long_name="drift in global mean ssh volume wrt timestep 1" unit="km3" /> 526 558 <field id="bgvole3t" long_name="drift in global mean volume variation (e3t) wrt timestep 1" unit="km3" /> 527 <field id="bgvoltot" long_name="drift in global mean volume wrt timestep 1" unit="km3" /> 528 <!-- NOTE: No matching iom_put call --> 529 <field id="bgsshtot" long_name="drift in global mean ssh wrt timestep 1" standard_name="global_average_sea_level_change" unit="m" /> 530 <field id="bgfrcvol" long_name="drift in global mean volume from forcing wrt timestep 1" unit="km3" /> 531 <field id="bgfrctem" long_name="drift in global mean heat content from forcing wrt timestep 1" unit="10^9J" /> 532 <field id="bgfrcsal" long_name="drift in global mean salt content from forcing wrt timestep 1" unit="1e-3*km3" /> 533 <field id="bgmistem" long_name="global mean temperature error due to free surface" unit="degC" /> 534 <field id="bgmissal" long_name="global mean salinity error due to free surface" unit="1e-3" /> 535 </field_group> 559 <field id="bgfrcvol" long_name="global mean volume from forcing" unit="km3" /> 560 <field id="bgfrctem" long_name="global mean heat content from forcing" unit="1.e20J" /> 561 <field id="bgfrchfx" long_name="global mean heat flux from forcing" unit="W/m2" /> 562 <field id="bgfrcsal" long_name="global mean salt content from forcing" unit="1e-3*km3" /> 563 <field id="bgmistem" long_name="global mean temperature error due to free surface (no vvl)" unit="degC" /> 564 <field id="bgmissal" long_name="global mean salinity error due to free surface (no vvl)" unit="1e-3" /> 565 </field_group> 536 566 537 567 <!-- LIM3 scalar variables --> … … 539 569 <field_group id="SBC_scalar" domain_ref="1point" > 540 570 <!-- available with ln_limdiaout --> 541 <field id="ibgvoltot" long_name="global mean ice volume" unit="km3" /> 542 <field id="sbgvoltot" long_name="global mean snow volume" unit="km3" /> 543 <field id="ibgarea" long_name="global mean ice area" unit="km2" /> 544 <field id="ibgsaline" long_name="global mean ice salinity" unit="1e-3" /> 545 <field id="ibgtemper" long_name="global mean ice temperature" unit="degC" /> 546 <field id="ibgheatco" long_name="global mean ice heat content" unit="10^20J" /> 547 <field id="sbgheatco" long_name="global mean snow heat content" unit="10^20J" /> 548 <field id="ibgsaltco" long_name="global mean ice salt content" unit="1e-3*km3" /> 549 550 <field id="ibgvfx" long_name="global mean volume flux (emp)" unit="m/day" /> 551 <field id="ibgvfxbog" long_name="global mean volume flux (bottom growth)" unit="m/day" /> 552 <field id="ibgvfxopw" long_name="global mean volume flux (open water growth)" unit="m/day" /> 553 <field id="ibgvfxsni" long_name="global mean volume flux (snow-ice growth)" unit="m/day" /> 554 <field id="ibgvfxdyn" long_name="global mean volume flux (dynamic growth)" unit="m/day" /> 555 <field id="ibgvfxbom" long_name="global mean volume flux (bottom melt)" unit="m/day" /> 556 <field id="ibgvfxsum" long_name="global mean volume flux (surface melt)" unit="m/day" /> 557 <field id="ibgvfxres" long_name="global mean volume flux (resultant)" unit="m/day" /> 558 <field id="ibgvfxspr" long_name="global mean volume flux (snow precip)" unit="m/day" /> 559 <field id="ibgvfxsnw" long_name="global mean volume flux (snow melt)" unit="m/day" /> 560 <field id="ibgvfxsub" long_name="global mean volume flux (snow sublimation)" unit="m/day" /> 561 562 <field id="ibgsfx" long_name="global mean salt flux (total)" unit="1e-3*m/day" /> 563 <field id="ibgsfxbri" long_name="global mean salt flux (brines)" unit="1e-3*m/day" /> 564 <field id="ibgsfxdyn" long_name="global mean salt flux (dynamic)" unit="1e-3*m/day" /> 565 <field id="ibgsfxres" long_name="global mean salt flux (resultant)" unit="1e-3*m/day" /> 566 <field id="ibgsfxbog" long_name="global mean salt flux (thermo)" unit="1e-3*m/day" /> 567 <field id="ibgsfxopw" long_name="global mean salt flux (thermo)" unit="1e-3*m/day" /> 568 <field id="ibgsfxsni" long_name="global mean salt flux (thermo)" unit="1e-3*m/day" /> 569 <field id="ibgsfxbom" long_name="global mean salt flux (thermo)" unit="1e-3*m/day" /> 570 <field id="ibgsfxsum" long_name="global mean salt flux (thermo)" unit="1e-3*m/day" /> 571 572 <field id="ibghfxdhc" long_name="Heat content variation in snow and ice" unit="W" /> 573 <field id="ibghfxspr" long_name="Heat content of snow precip" unit="W" /> 574 575 <field id="ibghfxthd" long_name="heat fluxes from ice-ocean exchange during thermo" unit="W" /> 576 <field id="ibghfxsum" long_name="heat fluxes causing surface ice melt" unit="W" /> 577 <field id="ibghfxbom" long_name="heat fluxes causing bottom ice melt" unit="W" /> 578 <field id="ibghfxbog" long_name="heat fluxes causing bottom ice growth" unit="W" /> 579 <field id="ibghfxdif" long_name="heat fluxes causing ice temperature change" unit="W" /> 580 <field id="ibghfxopw" long_name="heat fluxes causing open water ice formation" unit="W" /> 581 <field id="ibghfxdyn" long_name="heat fluxes from ice-ocean exchange during dynamic" unit="W" /> 582 <field id="ibghfxres" long_name="heat fluxes from ice-ocean exchange during resultant" unit="W" /> 583 <field id="ibghfxsub" long_name="heat fluxes from sublimation" unit="W" /> 584 <field id="ibghfxsnw" long_name="heat fluxes from snow-ocean exchange" unit="W" /> 585 <field id="ibghfxout" long_name="non solar heat fluxes received by the ocean" unit="W" /> 586 <field id="ibghfxin" long_name="total heat fluxes at the ice surface" unit="W" /> 587 588 <field id="ibgfrcvol" long_name="global mean forcing volume (emp)" unit="km3" /> 589 <field id="ibgfrcsfx" long_name="global mean forcing salt (sfx)" unit="1e-3*km3" /> 590 <field id="ibgvolgrm" long_name="global mean ice growth+melt volume" unit="km3" /> 571 <field id="ibgfrcvoltop" long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)" unit="km3" /> 572 <field id="ibgfrcvolbot" long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3" /> 573 <field id="ibgfrctemtop" long_name="global mean heat on top of ice/snw/ocean-atm " unit="1e20J" /> 574 <field id="ibgfrctembot" long_name="global mean heat below ice (on top of ocean) " unit="1e20J" /> 575 <field id="ibgfrcsal" long_name="global mean ice/snow forcing (salt equivalent ocean volume)" unit="pss*km3" /> 576 <field id="ibgfrchfxtop" long_name="global mean heat flux on top of ice/snw/ocean-atm " unit="W/m2" /> 577 <field id="ibgfrchfxbot" long_name="global mean heat flux below ice (on top of ocean) " unit="W/m2" /> 578 579 <field id="ibgvolume" long_name="drift in ice/snow volume (equivalent ocean volume)" unit="km3" /> 580 <field id="ibgsaltco" long_name="drift in ice salt content (equivalent ocean volume)" unit="pss*km3" /> 581 <field id="ibgheatco" long_name="drift in ice/snow heat content" unit="1e20J" /> 582 <field id="ibgheatfx" long_name="drift in ice/snow heat flux" unit="W/m2" /> 583 <field id="sbgvol_tot" long_name="global mean snow volume" unit="km3" /> 584 <field id="ibgarea_tot" long_name="global mean ice area" unit="km2" /> 585 <field id="ibgsalt_tot" long_name="global mean ice salt content" unit="1e-3*km3" /> 586 <field id="ibgheat_tot" long_name="global mean ice heat content" unit="1e20J" /> 587 <field id="sbgheat_tot" long_name="global mean snow heat content" unit="1e20J" /> 591 588 </field_group> 592 589 … … 622 619 623 620 <!-- Poleward transport : ptr --> 624 <field_group id="diaptr" domain_ref="ptr" > 621 <field_group id="diaptr" domain_ref="ptr" > 625 622 <field id="zomsfglo" long_name="Meridional Stream-Function: Global" unit="Sv" grid_ref="gznl_W_3D" /> 626 623 <field id="zomsfatl" long_name="Meridional Stream-Function: Atlantic" unit="Sv" grid_ref="gznl_W_3D" /> … … 628 625 <field id="zomsfind" long_name="Meridional Stream-Function: Indian" unit="Sv" grid_ref="gznl_W_3D" /> 629 626 <field id="zomsfipc" long_name="Meridional Stream-Function: Pacific+Indian" unit="Sv" grid_ref="gznl_W_3D" /> 630 <field id="zotemglo" long_name="Zonal Mean Temperature : Global" unit="deg C" grid_ref="gznl_T_3D" />631 <field id="zotematl" long_name="Zonal Mean Temperature : Atlantic" unit="deg C" grid_ref="gznl_T_3D" />632 <field id="zotempac" long_name="Zonal Mean Temperature : Pacific" unit="deg C" grid_ref="gznl_T_3D" />633 <field id="zotemind" long_name="Zonal Mean Temperature : Indian" unit="deg C" grid_ref="gznl_T_3D" />634 <field id="zotemipc" long_name="Zonal Mean Temperature : Pacific+Indian" unit="deg C" grid_ref="gznl_T_3D" />635 <field id="zosalglo" long_name="Zonal Mean Salinity : Global" unit=" 1e-3" grid_ref="gznl_T_3D" />636 <field id="zosalatl" long_name="Zonal Mean Salinity : Atlantic" unit=" 1e-3" grid_ref="gznl_T_3D" />637 <field id="zosalpac" long_name="Zonal Mean Salinity : Pacific" unit=" 1e-3" grid_ref="gznl_T_3D" />638 <field id="zosalind" long_name="Zonal Mean Salinity : Indian" unit=" 1e-3" grid_ref="gznl_T_3D" />639 <field id="zosalipc" long_name="Zonal Mean Salinity : Pacific+Indian" unit=" 1e-3" grid_ref="gznl_T_3D" />627 <field id="zotemglo" long_name="Zonal Mean Temperature : Global" unit="degree_C" grid_ref="gznl_T_3D" /> 628 <field id="zotematl" long_name="Zonal Mean Temperature : Atlantic" unit="degree_C" grid_ref="gznl_T_3D" /> 629 <field id="zotempac" long_name="Zonal Mean Temperature : Pacific" unit="degree_C" grid_ref="gznl_T_3D" /> 630 <field id="zotemind" long_name="Zonal Mean Temperature : Indian" unit="degree_C" grid_ref="gznl_T_3D" /> 631 <field id="zotemipc" long_name="Zonal Mean Temperature : Pacific+Indian" unit="degree_C" grid_ref="gznl_T_3D" /> 632 <field id="zosalglo" long_name="Zonal Mean Salinity : Global" unit="0.001" grid_ref="gznl_T_3D" /> 633 <field id="zosalatl" long_name="Zonal Mean Salinity : Atlantic" unit="0.001" grid_ref="gznl_T_3D" /> 634 <field id="zosalpac" long_name="Zonal Mean Salinity : Pacific" unit="0.001" grid_ref="gznl_T_3D" /> 635 <field id="zosalind" long_name="Zonal Mean Salinity : Indian" unit="0.001" grid_ref="gznl_T_3D" /> 636 <field id="zosalipc" long_name="Zonal Mean Salinity : Pacific+Indian" unit="0.001" grid_ref="gznl_T_3D" /> 640 637 <field id="zosrfglo" long_name="Zonal Mean Surface" unit="m2" grid_ref="gznl_T_3D" /> 641 638 <field id="zosrfatl" long_name="Zonal Mean Surface : Atlantic" unit="m2" grid_ref="gznl_T_3D" /> … … 644 641 <field id="zosrfipc" long_name="Zonal Mean Surface : Pacific+Indian" unit="m2" grid_ref="gznl_T_3D" /> 645 642 <field id="sophtadv" long_name="Advective Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 643 <field id="sophtadv_atl" long_name="Advective Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 644 <field id="sophtadv_pac" long_name="Advective Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 645 <field id="sophtadv_ind" long_name="Advective Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 646 <field id="sophtadv_ipc" long_name="Advective Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 646 647 <field id="sophtldf" long_name="Diffusive Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 648 <field id="sophtldf_atl" long_name="Diffusive Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 649 <field id="sophtldf_pac" long_name="Diffusive Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 650 <field id="sophtldf_ind" long_name="Diffusive Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 651 <field id="sophtldf_ipc" long_name="Diffusive Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 652 <field id="sophtove" long_name="Overturning Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 653 <field id="sophtove_atl" long_name="Overturning Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 654 <field id="sophtove_pac" long_name="Overturning Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 655 <field id="sophtove_ind" long_name="Overturning Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 656 <field id="sophtove_ipc" long_name="Overturning Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 657 <field id="sophtbtr" long_name="Barotropic Heat Transport" unit="PW" grid_ref="gznl_T_2D" /> 658 <field id="sophtbtr_atl" long_name="Barotropic Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 659 <field id="sophtbtr_pac" long_name="Barotropic Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 660 <field id="sophtbtr_ind" long_name="Barotropic Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" /> 661 <field id="sophtbtr_ipc" long_name="Barotropic Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 662 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection" unit="PW" grid_ref="gznl_T_2D" /> 663 <field id="sophteiv_atl" long_name="Heat Transport from mesoscale eddy advection: Atlantic" unit="PW" grid_ref="gznl_T_2D" /> 664 <field id="sophteiv_pac" long_name="Heat Transport from mesoscale eddy advection: Pacific" unit="PW" grid_ref="gznl_T_2D" /> 665 <field id="sophteiv_ind" long_name="Heat Transport from mesoscale eddy advection: Indian" unit="PW" grid_ref="gznl_T_2D" /> 666 <field id="sophteiv_ipc" long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" /> 647 667 <field id="sopstadv" long_name="Advective Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 668 <field id="sopstadv_atl" long_name="Advective Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 669 <field id="sopstadv_pac" long_name="Advective Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 670 <field id="sopstadv_ind" long_name="Advective Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 671 <field id="sopstadv_ipc" long_name="Advective Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 672 <field id="sopstove" long_name="Overturning Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 673 <field id="sopstove_atl" long_name="Overturning Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 674 <field id="sopstove_pac" long_name="Overturning Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 675 <field id="sopstove_ind" long_name="Overturning Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 676 <field id="sopstove_ipc" long_name="Overturning Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 677 <field id="sopstbtr" long_name="Barotropic Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 678 <field id="sopstbtr_atl" long_name="Barotropic Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 679 <field id="sopstbtr_pac" long_name="Barotropic Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 680 <field id="sopstbtr_ind" long_name="Barotropic Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 681 <field id="sopstbtr_ipc" long_name="Barotropic Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 648 682 <field id="sopstldf" long_name="Diffusive Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" /> 683 <field id="sopstldf_atl" long_name="Diffusive Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 684 <field id="sopstldf_pac" long_name="Diffusive Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 685 <field id="sopstldf_ind" long_name="Diffusive Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 686 <field id="sopstldf_ipc" long_name="Diffusive Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 687 <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection" unit="Giga g/s" grid_ref="gznl_T_2D" /> 688 <field id="sopsteiv_atl" long_name="Salt Transport from mesoscale eddy advection: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" /> 689 <field id="sopsteiv_pac" long_name="Salt Transport from mesoscale eddy advection: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" /> 690 <field id="sopsteiv_ind" long_name="Salt Transport from mesoscale eddy advection: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 691 <field id="sopsteiv_ipc" long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" /> 649 692 </field_group> 650 693 … … 665 708 <field id="ttrd_ad" long_name="temperature-trend: advection" standard_name="tendency_of_sea_water_temperature_due_to_advection" unit="degC/s" > sqrt( ttrd_xad^2 + ttrd_yad^2 + ttrd_zad^2 ) </field> 666 709 <field id="strd_ad" long_name="salinity -trend: advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field> 710 <field id="ttrd_totad" long_name="temperature-trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="degC/s" /> 711 <field id="strd_totad" long_name="salinity -trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" /> 667 712 <field id="ttrd_sad" long_name="temperature-trend: surface adv. (no-vvl)" unit="degC/s" grid_ref="grid_T_2D" /> 668 713 <field id="strd_sad" long_name="salinity -trend: surface adv. (no-vvl)" unit="1e-3/s" grid_ref="grid_T_2D" /> … … 671 716 <field id="ttrd_zdf" long_name="temperature-trend: vertical diffusion" standard_name="tendency_of_sea_water_temperature_due_to_vertical_mixing" unit="degC/s" /> 672 717 <field id="strd_zdf" long_name="salinity -trend: vertical diffusion" standard_name="tendency_of_sea_water_salinity_due_to_vertical_mixing" unit="1e-3/s" /> 718 <field id="ttrd_evd" long_name="temperature-trend: EVD convection" unit="degC/s" /> 719 <field id="strd_evd" long_name="salinity -trend: EVD convection" unit="1e-3/s" /> 673 720 674 721 <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 722 <field id="ttrd_iso" long_name="temperature-trend: isopycnal diffusion" unit="degC/s" > ttrd_ldf + ttrd_zdf - ttrd_zdfp </field> 723 <field id="strd_iso" long_name="salinity -trend: isopycnal diffusion" unit="1e-3/s" > strd_ldf + strd_zdf - strd_zdfp </field> 675 724 <field id="ttrd_zdfp" long_name="temperature-trend: pure vert. diffusion" unit="degC/s" /> 676 725 <field id="strd_zdfp" long_name="salinity -trend: pure vert. diffusion" unit="1e-3/s" /> … … 689 738 <field id="ttrd_atf" long_name="temperature-trend: asselin time filter" unit="degC/s" /> 690 739 <field id="strd_atf" long_name="salinity -trend: asselin time filter" unit="1e-3/s" /> 740 <field id="ttrd_tot" long_name="temperature-trend: total model trend" unit="degC/s" /> 741 <field id="strd_tot" long_name="salinity -trend: total model trend" unit="1e-3/s" /> 742 743 <!-- Thickness weighted versions: --> 744 <field id="ttrd_xad_e3t" unit="degC/s * m" > ttrd_xad * e3t </field> 745 <field id="strd_xad_e3t" unit="1e-3/s * m" > strd_xad * e3t </field> 746 <field id="ttrd_yad_e3t" unit="degC/s * m" > ttrd_yad * e3t </field> 747 <field id="strd_yad_e3t" unit="1e-3/s * m" > strd_yad * e3t </field> 748 <field id="ttrd_zad_e3t" unit="degC/s * m" > ttrd_zad * e3t </field> 749 <field id="strd_zad_e3t" unit="1e-3/s * m" > strd_zad * e3t </field> 750 <field id="ttrd_ad_e3t" unit="degC/s * m" > ttrd_ad * e3t </field> 751 <field id="strd_ad_e3t" unit="1e-3/s * m" > strd_ad * e3t </field> 752 <field id="ttrd_totad_e3t" unit="degC/s * m" > ttrd_totad * e3t </field> 753 <field id="strd_totad_e3t" unit="1e-3/s * m" > strd_totad * e3t </field> 754 <field id="ttrd_ldf_e3t" unit="degC/s * m" > ttrd_ldf * e3t </field> 755 <field id="strd_ldf_e3t" unit="1e-3/s * m" > strd_ldf * e3t </field> 756 <field id="ttrd_zdf_e3t" unit="degC/s * m" > ttrd_zdf * e3t </field> 757 <field id="strd_zdf_e3t" unit="1e-3/s * m" > strd_zdf * e3t </field> 758 <field id="ttrd_evd_e3t" unit="degC/s * m" > ttrd_evd * e3t </field> 759 <field id="strd_evd_e3t" unit="1e-3/s * m" > strd_evd * e3t </field> 760 761 <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 762 <field id="ttrd_iso_e3t" unit="degC/s * m" > ttrd_iso * e3t </field> 763 <field id="strd_iso_e3t" unit="1e-3/s * m" > strd_iso * e3t </field> 764 <field id="ttrd_zdfp_e3t" unit="degC/s * m" > ttrd_zdfp * e3t </field> 765 <field id="strd_zdfp_e3t" unit="1e-3/s * m" > strd_zdfp * e3t </field> 766 767 <!-- --> 768 <field id="ttrd_dmp_e3t" unit="degC/s * m" > ttrd_dmp * e3t </field> 769 <field id="strd_dmp_e3t" unit="1e-3/s * m" > strd_dmp * e3t </field> 770 <field id="ttrd_bbl_e3t" unit="degC/s * m" > ttrd_bbl * e3t </field> 771 <field id="strd_bbl_e3t" unit="1e-3/s * m" > strd_bbl * e3t </field> 772 <field id="ttrd_npc_e3t" unit="degC/s * m" > ttrd_npc * e3t </field> 773 <field id="strd_npc_e3t" unit="1e-3/s * m" > strd_npc * e3t </field> 774 <field id="ttrd_qns_e3t" unit="degC/s * m" > ttrd_qns * e3t_surf </field> 775 <field id="strd_cdt_e3t" unit="degC/s * m" > strd_cdt * e3t_surf </field> 776 <field id="ttrd_qsr_e3t" unit="degC/s * m" > ttrd_qsr * e3t </field> 777 <field id="ttrd_bbc_e3t" unit="degC/s * m" > ttrd_bbc * e3t </field> 778 <field id="ttrd_atf_e3t" unit="degC/s * m" > ttrd_atf * e3t </field> 779 <field id="strd_atf_e3t" unit="1e-3/s * m" > strd_atf * e3t </field> 780 <field id="ttrd_tot_e3t" unit="degC/s * m" > ttrd_tot * e3t </field> 781 <field id="strd_tot_e3t" unit="1e-3/s * m" > strd_tot * e3t </field> 782 691 783 692 784 <!-- variables available with ln_KE_trd --> … … 782 874 783 875 <field_group id="ptrc_T" grid_ref="grid_T_3D"> 876 <!-- PISCES standard : variables available with ln_p4z --> 784 877 <field id="DIC" long_name="Dissolved inorganic Concentration" unit="mmol/m3" /> 785 <field id="DIC_ E3T" long_name="DIC * E3T" unit="mmol/m2" > DIC * e3t </field >878 <field id="DIC_e3t" long_name="DIC * e3t" unit="mmol/m2" > DIC * e3t </field > 786 879 <field id="Alkalini" long_name="Total Alkalinity Concentration" unit="mmol/m3" /> 787 <field id="Alkalini_ E3T" long_name="Alkalini * E3T" unit="mmol/m2" > Alkalini * e3t </field >880 <field id="Alkalini_e3t" long_name="Alkalini * e3t" unit="mmol/m2" > Alkalini * e3t </field > 788 881 <field id="O2" long_name="Oxygen Concentration" unit="mmol/m3" /> 789 <field id="O2_ E3T" long_name="O2 * E3T" unit="mmol/m2" > O2 * e3t </field >882 <field id="O2_e3t" long_name="O2 * e3t" unit="mmol/m2" > O2 * e3t </field > 790 883 <field id="CaCO3" long_name="Calcite Concentration" unit="mmol/m3" /> 791 <field id="CaCO3_ E3T" long_name="CaCO3 * E3T" unit="mmol/m2" > CaCO3 * e3t </field >884 <field id="CaCO3_e3t" long_name="CaCO3 * e3t" unit="mmol/m2" > CaCO3 * e3t </field > 792 885 <field id="PO4" long_name="Phosphate Concentration" unit="mmol/m3" /> 793 <field id="PO4_ E3T" long_name="PO4 * E3T" unit="mmol/m2" > PO4 * e3t </field >886 <field id="PO4_e3t" long_name="PO4 * e3t" unit="mmol/m2" > PO4 * e3t </field > 794 887 <field id="POC" long_name="Small organic carbon Concentration" unit="mmol/m3" /> 795 <field id="POC_ E3T" long_name="POC * E3T" unit="mmol/m2" > POC * e3t </field >888 <field id="POC_e3t" long_name="POC * e3t" unit="mmol/m2" > POC * e3t </field > 796 889 <field id="Si" long_name="Silicate Concentration" unit="mmol/m3" /> 797 <field id="Si_ E3T" long_name="Si * E3T" unit="mmol/m2" > Si * e3t </field >890 <field id="Si_e3t" long_name="Si * e3t" unit="mmol/m2" > Si * e3t </field > 798 891 <field id="PHY" long_name="(Nano)Phytoplankton Concentration" unit="mmol/m3" /> 799 <field id="PHY_ E3T" long_name="PHY * E3T" unit="mmol/m2" > PHY * e3t </field >892 <field id="PHY_e3t" long_name="PHY * e3t" unit="mmol/m2" > PHY * e3t </field > 800 893 <field id="ZOO" long_name="(Micro)Zooplankton Concentration" unit="mmol/m3" /> 801 <field id="ZOO_ E3T" long_name="ZOO2 * E3T" unit="mmol/m2" > ZOO * e3t </field >894 <field id="ZOO_e3t" long_name="ZOO2 * e3t" unit="mmol/m2" > ZOO * e3t </field > 802 895 <field id="DOC" long_name="Dissolved organic Concentration" unit="mmol/m3" /> 803 <field id="DOC_ E3T" long_name="DOC * E3T" unit="mmol/m2" > DOC * e3t </field >896 <field id="DOC_e3t" long_name="DOC * e3t" unit="mmol/m2" > DOC * e3t </field > 804 897 <field id="PHY2" long_name="Diatoms Concentration" unit="mmol/m3" /> 805 <field id="PHY2_ E3T" long_name="PHY2 * E3T" unit="mmol/m2" > PHY2 * e3t </field >898 <field id="PHY2_e3t" long_name="PHY2 * e3t" unit="mmol/m2" > PHY2 * e3t </field > 806 899 <field id="ZOO2" long_name="Mesozooplankton Concentration" unit="mmol/m3" /> 807 <field id="ZOO2_ E3T" long_name="ZOO2 * E3T" unit="mmol/m2" > ZOO2 * e3t </field >900 <field id="ZOO2_e3t" long_name="ZOO2 * e3t" unit="mmol/m2" > ZOO2 * e3t </field > 808 901 <field id="DSi" long_name="Diatoms Silicate Concentration" unit="mmol/m3" /> 809 <field id="DSi_ E3T" long_name="Dsi * E3T" unit="mmol/m2" > DSi * e3t </field >902 <field id="DSi_e3t" long_name="Dsi * e3t" unit="mmol/m2" > DSi * e3t </field > 810 903 <field id="Fer" long_name="Dissolved Iron Concentration" unit="mmol/m3" /> 811 <field id="Fer_ E3T" long_name="Fer * E3T" unit="mmol/m2" > Fer * e3t </field >904 <field id="Fer_e3t" long_name="Fer * e3t" unit="mmol/m2" > Fer * e3t </field > 812 905 <field id="BFe" long_name="Big iron particles Concentration" unit="mmol/m3" /> 813 <field id="BFe_ E3T" long_name="BFe * E3T" unit="mmol/m2" > BFe * e3t </field >906 <field id="BFe_e3t" long_name="BFe * e3t" unit="mmol/m2" > BFe * e3t </field > 814 907 <field id="GOC" long_name="Big organic carbon Concentration" unit="mmol/m3" /> 815 <field id="GOC_ E3T" long_name="GOC * E3T" unit="mmol/m2" > GOC * e3t </field >908 <field id="GOC_e3t" long_name="GOC * e3t" unit="mmol/m2" > GOC * e3t </field > 816 909 <field id="SFe" long_name="Small iron particles Concentration" unit="mmol/m3" /> 817 <field id="SFe_ E3T" long_name="SFe * E3T" unit="mmol/m2" > SFe * e3t </field >910 <field id="SFe_e3t" long_name="SFe * e3t" unit="mmol/m2" > SFe * e3t </field > 818 911 <field id="DFe" long_name="Diatoms iron Concentration" unit="mmol/m3" /> 819 <field id="DFe_ E3T" long_name="DFe * E3T" unit="mmol/m2" > DFe * e3t </field >912 <field id="DFe_e3t" long_name="DFe * e3t" unit="mmol/m2" > DFe * e3t </field > 820 913 <field id="GSi" long_name="Sinking biogenic Silicate Concentration" unit="mmol/m3" /> 821 <field id="GSi_ E3T" long_name="GSi * E3T" unit="mmol/m2" > GSi * e3t </field >914 <field id="GSi_e3t" long_name="GSi * e3t" unit="mmol/m2" > GSi * e3t </field > 822 915 <field id="NFe" long_name="Nano iron Concentration" unit="mmol/m3" /> 823 <field id="NFe_ E3T" long_name="NFe * E3T" unit="mmol/m2" > NFe * e3t </field >916 <field id="NFe_e3t" long_name="NFe * e3t" unit="mmol/m2" > NFe * e3t </field > 824 917 <field id="NCHL" long_name="Nano chlorophyl Concentration" unit="mg/m3" /> 825 <field id="NCHL_ E3T" long_name="NCHL * E3T" unit="mmol/m2" > NCHL * e3t </field >918 <field id="NCHL_e3t" long_name="NCHL * e3t" unit="mmol/m2" > NCHL * e3t </field > 826 919 <field id="DCHL" long_name="Diatoms chlorophyl Concentration" unit="mg/m3" /> 827 <field id="DCHL_ E3T" long_name="DCHL * E3T" unit="mmol/m2" > DCHL * e3t </field >920 <field id="DCHL_e3t" long_name="DCHL * e3t" unit="mmol/m2" > DCHL * e3t </field > 828 921 <field id="NO3" long_name="Nitrate Concentration" unit="mmol/m3" /> 829 <field id="NO3_ E3T" long_name="NO3 * E3T" unit="mmol/m2" > NO3 * e3t </field >922 <field id="NO3_e3t" long_name="NO3 * e3t" unit="mmol/m2" > NO3 * e3t </field > 830 923 <field id="NH4" long_name="Ammonium Concentration" unit="mmol/m3" /> 831 <field id="NH4_E3T" long_name="NH4 * E3T" unit="mmol/m2" > NH4 * e3t </field > 832 833 <!-- PISCES with Kriest parametisation : variables available with key_kriest --> 834 <field id="Num" long_name="Number of organic particles" unit="1" /> 835 <field id="Num_E3T" long_name="Num * E3T" unit="m" > Num * e3t </field > 836 837 <!-- PISCES light : variables available with key_pisces_reduced --> 924 <field id="NH4_e3t" long_name="NH4 * e3t" unit="mmol/m2" > NH4 * e3t </field > 925 926 <!-- PISCES quota : variables available with ln_p5z --> 927 928 <field id="DON" long_name="Dissolved organic N Concentration" unit="mmol/m3" /> 929 <field id="DON_e3t" long_name="DON * e3t" unit="mmol/m2" > DON * e3t </field > 930 <field id="DOP" long_name="Dissolved organic P Concentration" unit="mmol/m3" /> 931 <field id="DOP_e3t" long_name="DOP * e3t" unit="mmol/m2" > DOP * e3t </field > 932 <field id="PON" long_name="Small PON Concentration" unit="mmol/m3" /> 933 <field id="PON_e3t" long_name="PON * e3t" unit="mmol/m2" > PON * e3t </field > 934 <field id="POP" long_name="Small POP Concentration" unit="mmol/m3" /> 935 <field id="POP_e3t" long_name="POP * e3t" unit="mmol/m2" > POP * e3t </field > 936 <field id="GON" long_name="Big PON Concentration" unit="mmol/m3" /> 937 <field id="GON_e3t" long_name="GON * e3t" unit="mmol/m2" > GON * e3t </field > 938 <field id="GOP" long_name="Big POP Concentration" unit="mmol/m3" /> 939 <field id="GOP_e3t" long_name="GOP * e3t" unit="mmol/m2" > GOP * e3t </field > 940 <field id="PHYN" long_name="Nanophytoplankton N biomass" unit="mmol/m3" /> 941 <field id="PHYN_e3t" long_name="PHYN * e3t" unit="mmol/m2" > PHYN * e3t </field > 942 <field id="PHYP" long_name="Nanophytoplankton P biomass" unit="mmol/m3" /> 943 <field id="PHYP_e3t" long_name="PHYP * e3t" unit="mmol/m2" > PHYP * e3t </field > 944 <field id="DIAN" long_name="Diatoms N biomass" unit="mmol/m3" /> 945 <field id="DIAN_e3t" long_name="DIAN * e3t" unit="mmol/m2" > DIAN * e3t </field > 946 <field id="DIAP" long_name="Diatoms P biomass" unit="mmol/m3" /> 947 <field id="DIAP_e3t" long_name="DIAP * e3t" unit="mmol/m2" > DIAP * e3t </field > 948 <field id="PIC" long_name="Picophytoplankton C biomass" unit="mmol/m3" /> 949 <field id="PIC_e3t" long_name="PIC * e3t" unit="mmol/m2" > PIC * e3t </field > 950 <field id="PICN" long_name="Picophytoplankton N biomass" unit="mmol/m3" /> 951 <field id="PICN_e3t" long_name="PICN * e3t" unit="mmol/m2" > PICN * e3t </field > 952 <field id="PICP" long_name="Picophytoplankton P biomass" unit="mmol/m3" /> 953 <field id="PICP_e3t" long_name="PICP * e3t" unit="mmol/m2" > PICP * e3t </field > 954 <field id="PFe" long_name="Picophytoplankton Fe biomass" unit="mmol/m3" /> 955 <field id="PFe_e3t" long_name="PFe * e3t" unit="mmol/m2" > PFe * e3t </field > 956 <field id="PCHL" long_name="Picophytoplankton Chl biomass" unit="mg/m3" /> 957 <field id="PCHL_e3t" long_name="PCHL * e3t" unit="mmol/m2" > PCHL * e3t </field > 958 959 <!-- PISCES with ligand parametisation : variables available namelist paramter ln_ligand --> 960 <field id="LGW" long_name="Weak ligands concentration" unit="mmol/m3" /> 961 <field id="LGW_e3t" long_name="LGW * e3t" unit="mmol/m2" > LGW * e3t </field > 962 <field id="LFe" long_name="Lithogenic iron concentration" unit="mmol/m3" /> 963 <field id="LFe_e3t" long_name="LFe * e3t" unit="mmol/m2" > LFe * e3t </field > 964 965 <!-- PISCES light : variables available with ln_p2z --> 838 966 <field id="DET" long_name="Detritus" unit="mmol-N/m3" /> 839 <field id="DET_ E3T" long_name="DET * E3T" unit="mmol-N/m2" > DET * e3t </field >967 <field id="DET_e3t" long_name="DET * e3t" unit="mmol-N/m2" > DET * e3t </field > 840 968 <field id="DOM" long_name="Dissolved Organic Matter" unit="mmol-N/m3" /> 841 <field id="DOM_E3T" long_name="DOM * E3T" unit="mmol-N/m2" > DOM * e3t </field > 842 843 <!-- CFC11 : variables available with key_cfc --> 844 <field id="CFC11" long_name="CFC-11 Concentration" unit="umol/m3" /> 845 <field id="CFC11_E3T" long_name="CFC11 * E3T" unit="umol/m2" > CFC11 * e3t </field > 846 <!-- Bomb C14 : variables available with key_c14b --> 847 <field id="C14B" long_name="Bomb C14 Concentration" unit="1" /> 848 <field id="C14B_E3T" long_name="C14B * E3T" unit="m" > C14B * e3t </field > 969 <field id="DOM_e3t" long_name="DOM * e3t" unit="mmol-N/m2" > DOM * e3t </field > 970 971 <!-- CFC11 : variables available with ln_cfc11 --> 972 <field id="CFC11" long_name="Chlorofluoro carbon11 Concentration" unit="umol/m3" /> 973 <field id="CFC11_e3t" long_name="CFC11 * e3t" unit="umol/m2" > CFC11 * e3t </field > 974 975 <!-- CFC12 : variables available with ln_cfc12 --> 976 <field id="CFC12" long_name="Chlorofluoro carbon12 Concentration" unit="umol/m3" /> 977 <field id="CFC12_e3t" long_name="CFC12 * e3t" unit="umol/m2" > CFC12 * e3t </field > 978 979 <!-- SF6 : variables available with ln_sf6 --> 980 <field id="SF6" long_name="Sulfur hexafluoride Concentration" unit="umol/m3" /> 981 <field id="SF6_e3t" long_name="SF6 * e3t" unit="umol/m2" > SF6 * e3t </field > 982 983 <!-- C14 : variables available with ln_c14 --> 984 <field id="RC14" long_name="Radiocarbon ratio" unit="-" /> 985 <field id="RC14_e3t" long_name="RC14 * e3t" unit="m" > RC14 * e3t </field > 986 987 <!-- AGE : variables available with ln_age --> 988 <field id="Age" long_name="Sea water age since surface contact" unit="yr" /> 989 <field id="Age_e3t" long_name="Age * e3t" unit="yr * m" > Age * e3t </field > 990 849 991 </field_group> 850 992 … … 857 999 <field id="PAR" long_name="Photosynthetically Available Radiation" unit="W/m2" grid_ref="grid_T_3D" /> 858 1000 <field id="PARDM" long_name="Daily mean PAR" unit="W/m2" grid_ref="grid_T_3D" /> 859 <field id="PPPHY" long_name="Primary production of nanophyto" unit="mol/m3/s" grid_ref="grid_T_3D" /> 860 <field id="PPPHY2" long_name="Primary production of diatoms" unit="mol/m3/s" grid_ref="grid_T_3D" /> 861 <field id="PPNEWN" long_name="New Primary production of nanophyto" unit="mol/m3/s" grid_ref="grid_T_3D" /> 862 <field id="PPNEWD" long_name="New Primary production of diatoms" unit="mol/m3/s" grid_ref="grid_T_3D" /> 863 <field id="PBSi" long_name="Primary production of Si diatoms" unit="mol/m3/s" grid_ref="grid_T_3D" /> 864 <field id="PFeN" long_name="Primary production of nano iron" unit="mol/m3/s" grid_ref="grid_T_3D" /> 1001 <field id="PPPHYN" long_name="Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1002 <field id="PPPHYP" long_name="Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1003 <field id="PPPHYD" long_name="Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1004 <field id="PPNEWN" long_name="New Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1005 <field id="PPNEWP" long_name="New Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1006 <field id="PPNEWD" long_name="New Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1007 <field id="PBSi" long_name="Primary production of Si diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1008 <field id="PFeN" long_name="Primary production of nano iron" unit="molC/m3/s" grid_ref="grid_T_3D" /> 1009 <field id="PFeP" long_name="Primary production of pico iron" unit="molC/m3/s" grid_ref="grid_T_3D" /> 865 1010 <field id="PFeD" long_name="Primary production of diatoms iron" unit="mol/m3/s" grid_ref="grid_T_3D" /> 866 1011 <field id="xfracal" long_name="Calcifying fraction" unit="1" grid_ref="grid_T_3D" /> … … 871 1016 <field id="REMIN" long_name="Oxic remineralization of OM" unit="mol/m3/s" grid_ref="grid_T_3D" /> 872 1017 <field id="DENIT" long_name="Anoxic remineralization of OM" unit="mol/m3/s" grid_ref="grid_T_3D" /> 1018 <field id="REMINP" long_name="Oxic remineralization rate of POC" unit="d-1" grid_ref="grid_T_3D" /> 1019 <field id="REMING" long_name="Oxic remineralization rate of GOC" unit="d-1" grid_ref="grid_T_3D" /> 873 1020 <field id="Nfix" long_name="Nitrogen fixation" unit="mol/m3/s" grid_ref="grid_T_3D" /> 874 1021 <field id="Mumax" long_name="Maximum growth rate" unit="s-1" grid_ref="grid_T_3D" /> 875 1022 <field id="MuN" long_name="Realized growth rate for nanophyto" unit="s-1" grid_ref="grid_T_3D" /> 1023 <field id="MuP" long_name="Realized growth rate for picophyto" unit="s-1" grid_ref="grid_T_3D" /> 876 1024 <field id="MuD" long_name="Realized growth rate for diatomes" unit="s-1" grid_ref="grid_T_3D" /> 1025 <field id="MunetN" long_name="Net growth rate for nanophyto" unit="s-1" grid_ref="grid_T_3D" /> 1026 <field id="MunetP" long_name="Net growth rate for picophyto" unit="s-1" grid_ref="grid_T_3D" /> 1027 <field id="MunetD" long_name="Net growth rate for diatomes" unit="s-1" grid_ref="grid_T_3D" /> 877 1028 <field id="LNnut" long_name="Nutrient limitation term in Nanophyto" unit="" grid_ref="grid_T_3D" /> 1029 <field id="LPnut" long_name="Nutrient limitation term in Picophyto" unit="-" grid_ref="grid_T_3D" /> 878 1030 <field id="LDnut" long_name="Nutrient limitation term in Diatoms" unit="" grid_ref="grid_T_3D" /> 879 1031 <field id="LNFe" long_name="Iron limitation term in Nanophyto" unit="" grid_ref="grid_T_3D" /> 1032 <field id="LPFe" long_name="Iron limitation term in Picophyto" unit="-" grid_ref="grid_T_3D" /> 880 1033 <field id="LDFe" long_name="Iron limitation term in Diatoms" unit="" grid_ref="grid_T_3D" /> 881 1034 <field id="LNlight" long_name="Light limitation term in Nanophyto" unit="" grid_ref="grid_T_3D" /> 1035 <field id="LPlight" long_name="Light limitation term in Picophyto" unit="-" grid_ref="grid_T_3D" /> 882 1036 <field id="LDlight" long_name="Light limitation term in Diatoms" unit="" grid_ref="grid_T_3D" /> 1037 <field id="SIZEN" long_name="Mean relative size of nanophyto." unit="-" grid_ref="grid_T_3D" /> 1038 <field id="SIZEP" long_name="Mean relative size of picophyto." unit="-" grid_ref="grid_T_3D" /> 1039 <field id="SIZED" long_name="Mean relative size of diatoms" unit="-" grid_ref="grid_T_3D" /> 883 1040 <field id="Fe2" long_name="Iron II concentration" unit="nmol/m3" grid_ref="grid_T_3D" /> 884 1041 <field id="Fe3" long_name="Iron III concentration" unit="nmol/m3" grid_ref="grid_T_3D" /> … … 893 1050 <field id="Sdenit" long_name="Nitrate reduction in the sediments" unit="mol/m2/s" /> 894 1051 <field id="Ironice" long_name="Iron input/uptake due to sea ice" unit="mol/m2/s" /> 1052 <field id="SedCal" long_name="Calcite burial in the sediments" unit="molC/m2/s" /> 1053 <field id="SedSi" long_name="Silicon burial in the sediments" unit="molSi/m2/s" /> 1054 <field id="SedC" long_name="Organic C burial in the sediments" unit="molC/m2/s" /> 895 1055 <field id="HYDR" long_name="Iron input from hydrothemal vents" unit="mol/m2/s" grid_ref="grid_T_3D" /> 896 1056 <field id="EPC100" long_name="Export of carbon particles at 100 m" unit="mol/m2/s" /> … … 911 1071 <field id="Ironsed" long_name="Iron deposition from sediment" unit="mol/m2/s" grid_ref="grid_T_3D" /> 912 1072 913 914 <!-- PISCES with Kriest parametisation : variables available with key_kriest -->915 <field id="EPN100" long_name="Particulate number flux at 100 m" unit="mol/m2/s" />916 <field id="EXPN" long_name="Particulate number flux" unit="mol/m2/s" grid_ref="grid_T_3D" />917 <field id="XNUM" long_name="Number of particles in aggregats" unit="1" grid_ref="grid_T_3D" />918 <field id="WSC" long_name="sinking speed of mass flux" unit="m2/s" grid_ref="grid_T_3D" />919 <field id="WSN" long_name="sinking speed of number flux" unit="m2/s" grid_ref="grid_T_3D" />920 921 1073 <!-- dbio_T on T grid : variables available with key_diaar5 --> 922 1074 <field id="TPP" long_name="Total Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D" /> … … 927 1079 <field id="ZO2MIN" long_name="Depth of oxygen minimum concentration" unit="m" /> 928 1080 <field id="INTNFIX" long_name="Nitrogen fixation rate : vert. integrated" unit="mol/m2/s" /> 929 <field id="INTPPPHY " long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s" />930 <field id="INTPPPHY 2" long_name="Vertically integrated primary production by diatom" unit="mol/m2/s" />1081 <field id="INTPPPHYN" long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s" /> 1082 <field id="INTPPPHYD" long_name="Vertically integrated primary production by diatom" unit="mol/m2/s" /> 931 1083 <field id="INTPP" long_name="Vertically integrated primary production by phyto" unit="mol/m2/s" /> 932 1084 <field id="INTPNEW" long_name="Vertically integrated new primary production" unit="mol/m2/s" /> … … 958 1110 <field id="TDETSED" long_name="TDETSED" unit="" /> 959 1111 960 <!-- CFC11 : variables available with key_cfc --> 961 <field id="qtrCFC11" long_name="Air-sea flux of CFC-11" unit="mol/m2/s" /> 962 <field id="qintCFC11" long_name="Cumulative air-sea flux of CFC-11" unit="mol/m2" /> 963 964 <!-- Bomb C14 : variables available with key_c14b --> 965 <field id="qtrC14b" long_name="Air-sea flux of Bomb C14" unit="mol/m2/s" /> 966 <field id="qintC14b" long_name="Cumulative air-sea flux of Bomb C14" unit="mol/m2" /> 967 <field id="fdecay" long_name="Radiactive decay of Bomb C14" unit="mol/m3" grid_ref="grid_T_3D" /> 1112 <!-- CFC11 : variables available with ln_cfc11 --> 1113 <field id="qtr_CFC11" long_name="Air-sea flux of CFC-11" unit="mol/m2/s" /> 1114 <field id="qint_CFC11" long_name="Cumulative air-sea flux of CFC-11" unit="mol/m2" /> 1115 1116 <!-- CFC12 : variables available with ln_cfc12 --> 1117 <field id="qtr_CFC12" long_name="Air-sea flux of CFC12" unit="mol/m2/s" /> 1118 <field id="qint_CFC12" long_name="Cumulative air-sea flux of CFC12" unit="mol/m2" /> 1119 1120 <!-- SF6 : variables available with ln_sf6 --> 1121 <field id="qtr_SF6" long_name="Air-sea flux of SF6" unit="mol/m2/s" /> 1122 <field id="qint_SF6" long_name="Cumulative air-sea flux of SF6" unit="mol/m2" /> 1123 1124 <!-- C14 : variables available with ln_c14 --> 1125 <field id="DeltaC14" long_name="Delta C14" unit="permil" grid_ref="grid_T_3D" /> 1126 <field id="C14Age" long_name="Radiocarbon age" unit="yr" grid_ref="grid_T_3D" /> 1127 <field id="RAge" long_name="Reservoir Age" unit="yr" /> 1128 <field id="qtr_C14" long_name="Air-sea flux of C14" unit="1/m2/s" /> 1129 <field id="qint_C14" long_name="Cumulative air-sea flux of C14" unit="1/m2" /> 968 1130 </field_group> 969 1131 970 <field_group id="PISCES_scalar" domain_ref="1point" > 971 <field id="pno3tot" long_name="global mean nitrate concentration" unit="mol/m3" /> 1132 <field_group id="tracer_scalar" domain_ref="1point" > 1133 <!-- PISCES scalar --> 1134 <field id="pno3tot" long_name="Global mean nitrate concentration" unit="mol/m3" /> 972 1135 <field id="ppo4tot" long_name="global mean phosphorus concentration" unit="mol/m3" /> 973 <field id="psiltot" long_name=" global mean silicate concentration" unit="mol/m3" />974 <field id="palktot" long_name=" global mean alkalinity concentration" unit="mol/m3" />975 <field id="pfertot" long_name=" global mean iron concentration" unit="mol/m3" />976 <field id="tcflx" long_name=" total Flux of Carbon out of the ocean" unit="mol/s" />977 <field id="tcflxcum" long_name=" cumulative total Flux of Carbon out of the ocean" unit="mol/s" />978 <field id="tcexp" long_name=" total Carbon export at 100m" unit="mol/s" />979 <field id="tintpp" long_name=" global total integrated primary production" unit="mol/s" />980 <field id="tnfix" long_name=" global total nitrogen fixation" unit="mol/s" />1136 <field id="psiltot" long_name="Global mean silicate concentration" unit="mol/m3" /> 1137 <field id="palktot" long_name="Global mean alkalinity concentration" unit="mol/m3" /> 1138 <field id="pfertot" long_name="Global mean iron concentration" unit="mol/m3" /> 1139 <field id="tcflx" long_name="Total Flux of Carbon out of the ocean" unit="mol/s" /> 1140 <field id="tcflxcum" long_name="Cumulative total Flux of Carbon out of the ocean" unit="mol/s" /> 1141 <field id="tcexp" long_name="Total Carbon export at 100m" unit="mol/s" /> 1142 <field id="tintpp" long_name="Global total integrated primary production" unit="mol/s" /> 1143 <field id="tnfix" long_name="Global total nitrogen fixation" unit="mol/s" /> 981 1144 <field id="tdenit" long_name="Total denitrification" unit="mol/s" /> 1145 <!-- C14 scalar --> 1146 <field id="AtmCO2" long_name="Global atmospheric CO2" unit="ppm" /> 1147 <field id="AtmC14" long_name="Global atmospheric DeltaC14" unit="permil" /> 1148 <field id="K_C14" long_name="Global 14C/C exchange velocity" unit="m/yr" /> 1149 <field id="K_CO2" long_name="Global CO2 piston velocity" unit="cm/h" /> 1150 <field id="C14Inv" long_name="global Radiocarbon ocean inventory" unit="10^26 atoms" /> 982 1151 </field_group> 983 1152 -
branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/EXP00/namelist_cfg
r7412 r7467 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : GYRE_PISCES Configuration namelist to overwrite reference dynamical namelist2 !! NEMO/OPA Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref 3 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 4 ! 5 !----------------------------------------------------------------------- 6 &namusr_def 7 !----------------------------------------------------------------------- 8 rn_dx = 1000.0 9 rn_dz = 1.0 10 nn_wad_test = 1 11 / 4 12 !----------------------------------------------------------------------- 5 13 &namrun ! parameters of the run … … 7 15 cn_exp = "WAD" ! experience name 8 16 nn_it000 = 1 ! first time step 9 nn_itend = 5760 ! last time step17 nn_itend = 3840 ! last time step 10 18 nn_leapy = 30 ! Leap year calendar (1) or not (0) 11 19 nn_stock = 48000 ! frequency of creation of a restart file (modulo referenced to 1) … … 16 24 / 17 25 !----------------------------------------------------------------------- 18 &namcfg ! parameters of the configuration 19 !----------------------------------------------------------------------- 20 cp_cfg = "wad" ! name of the configuration 21 jp_cfg = 1 ! resolution of the configuration 22 jpidta = 51 ! 1st lateral dimension ( >= jpi ) = 30*jp_cfg+2 23 jpjdta = 34 ! 2nd " " ( >= jpj ) = 20*jp_cfg+2 24 jpkdta = 10 ! number of levels ( >= jpk ) 25 jpiglo = 51 ! 1st dimension of global domain --> i = jpidta 26 jpjglo = 34 ! 2nd - - --> j = jpjdta 27 jpizoom = 1 ! left bottom (i,j) indices of the zoom 28 jpjzoom = 1 ! in data domain indices 29 jperio = 0 ! lateral cond. type (between 0 and 6) 26 &namcfg ! parameters of the configuration 27 !----------------------------------------------------------------------- 28 ln_read_cfg = .false. ! (=T) read the domain configuration file 29 ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules 30 30 / 31 31 !----------------------------------------------------------------------- 32 32 &namzgr ! vertical coordinate 33 33 !----------------------------------------------------------------------- 34 ln_sco = .true. ! s- or hybrid z-s-coordinate 35 ln_linssh = .false. ! linear free surface 36 / 37 !----------------------------------------------------------------------- 38 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 39 !----------------------------------------------------------------------- 40 ln_s_sh94 = .false. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 41 ln_s_sf12 = .true. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 42 ln_sigcrit = .true. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 43 ! stretching coefficients for all functions 44 rn_sbot_min = 0.01 ! minimum depth of s-bottom surface (>0) (m) 45 rn_sbot_max = 15.0 ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 46 rn_hc = 3.0 ! critical depth for transition to stretched coordinates 34 ln_zps = .false. ! z-coordinate - partial steps 35 ln_sco = .true. ! s-coordinate 47 36 / 48 37 !----------------------------------------------------------------------- 49 38 &namdom ! space and time domain (bathymetry, mesh, timestep) 50 39 !----------------------------------------------------------------------- 51 nn_msh = 1 ! create (=1) a mesh file or not (=0) 52 nn_bathy = 0 ! compute (=0) or read (=1) the bathymetry file 53 rn_bathy = 10. ! value of the bathymetry. if (=0) bottom flat at jpkm1 54 rn_rdt = 12. ! time step for the dynamics 55 jphgr_msh = 1 ! type of horizontal mesh 56 ppglam0 = 0.0 ! longitude of first raw and column T-point (jphgr_msh = 1) 57 ppgphi0 = 10.0 ! latitude of first raw and column T-point (jphgr_msh = 1) 58 ppe1_deg = 0.01 ! zonal grid-spacing (degrees) 59 ppe2_deg = 0.01 ! meridional grid-spacing (degrees) 60 ppe1_m = 999999.0 ! zonal grid-spacing (degrees) 61 ppe2_m = 999999.0 ! meridional grid-spacing (degrees) 62 ppsur = 999999.0 ! ORCA r4, r2 and r05 coefficients 63 ppa0 = 999999.0 ! (default coefficients) 64 ppa1 = 999999.0 ! 65 ppkth = 999999.0 ! 66 ppacr = 999999.0 ! 67 ppdzmin = 0.2 ! Minimum vertical spacing 68 pphmax = 10.0 ! Maximum depth 69 ldbletanh = .FALSE. ! Use/do not use double tanf function for vertical coordinates 70 ppa2 = 999999.0 ! Double tanh function parameters 71 ppkth2 = 999999.0 ! 72 ppacr2 = 999999.0 ! 40 ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time 41 ! 42 nn_msh = 1 ! create (>0) a mesh file or not (=0) 43 rn_rdt = 18. ! time step for the dynamics 73 44 / 74 45 !----------------------------------------------------------------------- … … 89 60 nn_fsbc = 1 ! frequency of surface boundary condition computation 90 61 ! ! (also = the frequency of sea-ice model call) 91 ln_ ana= .true. ! analytical formulation (T => fill namsbc_ana )92 ln_blk _core= .false. ! CORE bulk formulation (T => fill namsbc_core)62 ln_usr = .true. ! analytical formulation (T => fill namsbc_ana ) 63 ln_blk = .false. ! CORE bulk formulation (T => fill namsbc_core) 93 64 nn_ice = 0 ! =0 no ice boundary condition , 94 65 ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) … … 257 228 &nameos ! ocean physical parameters 258 229 !----------------------------------------------------------------------- 259 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 260 ! =-1, TEOS-10 261 ! = 0, EOS-80 262 ! = 1, S-EOS (simplified eos) 263 ln_useCT = .false. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 264 ! ! 265 ! ! S-EOS coefficients : 230 ln_teos10 = .false. ! = Use TEOS-10 equation of state 231 ln_eos80 = .false. ! = Use EOS80 equation of state 232 ln_seos = .true. ! = Use simplified equation of state (S-EOS) 233 ! 234 ! ! S-EOS coefficients (ln_seos=T): 266 235 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 267 236 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) … … 339 308 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 340 309 !----------------------------------------------------------------------- 341 ln_dynvor_ene = . true.! enstrophy conserving scheme310 ln_dynvor_ene = .false. ! enstrophy conserving scheme 342 311 ln_dynvor_ens = .false. ! energy conserving scheme 343 312 ln_dynvor_mix = .false. ! mixed scheme 344 ln_dynvor_een = . false.! energy & enstrophy scheme345 nn_een_e3f = 1! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)313 ln_dynvor_een = .true. ! energy & enstrophy scheme 314 nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 346 315 / 347 316 !----------------------------------------------------------------------- … … 356 325 !----------------------------------------------------------------------- 357 326 ln_dynspg_ts = .true. ! split-explicit free surface 327 ln_bt_auto = .false. ! Number of sub-step defined from: 328 nn_baro = 12 ! =F : the number of sub-step in rn_rdt seconds 358 329 / 359 330 !----------------------------------------------------------------------- … … 476 447 !----------------------------------------------------------------------- 477 448 ln_wd = .true. ! T/F activation of wetting and drying 478 !rn_wdmin1 = 0.25 ! Minimum wet depth on dried cells 479 rn_wdmin1 = 0.4 ! Minimum wet depth on dried cells 480 rn_wdmin2 = 0.00001 ! Tolerance of min wet depth on dried cells 481 rn_wdld = 10.0 ! Land elevation below which wetting/drying is allowed 449 rn_wdmin1 = 0.25 ! Minimum wet depth on dried cells 450 rn_wdmin2 = 0.000 ! Tolerance of min wet depth on dried cells 451 rn_wdld = 2.5 ! Land elevation below which wetting/drying is allowed 482 452 nn_wdit = 50 ! Max iterations for W/D limiter 483 453 / -
branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/MY_SRC/domwri.F90
r7412 r7467 8 8 !! NEMO 1.0 ! 2002-08 (G. Madec) F90 and several file 9 9 !! 3.0 ! 2008-01 (S. Masson) add dom_uniq 10 !! 4.0 ! 2016-01 (G. Madec) simplified mesh_mask.nc file 10 11 !!---------------------------------------------------------------------- 11 12 … … 13 14 !! dom_wri : create and write mesh and mask file(s) 14 15 !! dom_uniq : identify unique point of a grid (TUVF) 16 !! dom_stiff : diagnose maximum grid stiffness/hydrostatic consistency (s-coordinate) 15 17 !!---------------------------------------------------------------------- 16 18 USE dom_oce ! ocean space and time domain 19 USE phycst , ONLY : rsmall 20 USE wet_dry, ONLY : ln_wd, ht_wd 21 ! 17 22 USE in_out_manager ! I/O manager 18 23 USE iom ! I/O library … … 26 31 27 32 PUBLIC dom_wri ! routine called by inidom.F90 28 PUBLIC dom_wri_coordinate ! routine called by domhgr.F90 33 PUBLIC dom_stiff ! routine called by inidom.F90 34 29 35 !! * Substitutions 30 36 # include "vectopt_loop_substitute.h90" 31 37 !!---------------------------------------------------------------------- 32 !! NEMO/OPA 3.3 , NEMO Consortium (2010)33 !! $Id: domwri.F90 5836 2015-10-26 14:49:40Z cetlod$38 !! NEMO/OPA 4.0 , NEMO Consortium (2016) 39 !! $Id: domwri.F90 7421 2016-12-01 17:10:41Z flavoni $ 34 40 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 35 41 !!---------------------------------------------------------------------- 36 42 CONTAINS 37 38 SUBROUTINE dom_wri_coordinate39 !!----------------------------------------------------------------------40 !! *** ROUTINE dom_wri_coordinate ***41 !!42 !! ** Purpose : Create the NetCDF file which contains all the43 !! standard coordinate information plus the surface,44 !! e1e2u and e1e2v. By doing so, those surface will45 !! not be changed by the reduction of e1u or e2v scale46 !! factors in some straits.47 !! NB: call just after the read of standard coordinate48 !! and the reduction of scale factors in some straits49 !!50 !! ** output file : coordinate_e1e2u_v.nc51 !!----------------------------------------------------------------------52 INTEGER :: inum0 ! temprary units for 'coordinate_e1e2u_v.nc' file53 CHARACTER(len=21) :: clnam0 ! filename (mesh and mask informations)54 ! ! workspaces55 REAL(wp), POINTER, DIMENSION(:,: ) :: zprt, zprw56 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepu, zdepv57 !!----------------------------------------------------------------------58 !59 IF( nn_timing == 1 ) CALL timing_start('dom_wri_coordinate')60 !61 IF(lwp) WRITE(numout,*)62 IF(lwp) WRITE(numout,*) 'dom_wri_coordinate : create NetCDF coordinate file'63 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~~'64 65 clnam0 = 'coordinate_e1e2u_v' ! filename (mesh and mask informations)66 67 ! create 'coordinate_e1e2u_v.nc' file68 ! ============================69 !70 CALL iom_open( TRIM(clnam0), inum0, ldwrt = .TRUE., kiolib = jprstlib )71 !72 ! ! horizontal mesh (inum3)73 CALL iom_rstput( 0, 0, inum0, 'glamt', glamt, ktype = jp_r4 ) ! ! latitude74 CALL iom_rstput( 0, 0, inum0, 'glamu', glamu, ktype = jp_r4 )75 CALL iom_rstput( 0, 0, inum0, 'glamv', glamv, ktype = jp_r4 )76 CALL iom_rstput( 0, 0, inum0, 'glamf', glamf, ktype = jp_r4 )77 78 CALL iom_rstput( 0, 0, inum0, 'gphit', gphit, ktype = jp_r4 ) ! ! longitude79 CALL iom_rstput( 0, 0, inum0, 'gphiu', gphiu, ktype = jp_r4 )80 CALL iom_rstput( 0, 0, inum0, 'gphiv', gphiv, ktype = jp_r4 )81 CALL iom_rstput( 0, 0, inum0, 'gphif', gphif, ktype = jp_r4 )82 83 CALL iom_rstput( 0, 0, inum0, 'e1t', e1t, ktype = jp_r8 ) ! ! e1 scale factors84 CALL iom_rstput( 0, 0, inum0, 'e1u', e1u, ktype = jp_r8 )85 CALL iom_rstput( 0, 0, inum0, 'e1v', e1v, ktype = jp_r8 )86 CALL iom_rstput( 0, 0, inum0, 'e1f', e1f, ktype = jp_r8 )87 88 CALL iom_rstput( 0, 0, inum0, 'e2t', e2t, ktype = jp_r8 ) ! ! e2 scale factors89 CALL iom_rstput( 0, 0, inum0, 'e2u', e2u, ktype = jp_r8 )90 CALL iom_rstput( 0, 0, inum0, 'e2v', e2v, ktype = jp_r8 )91 CALL iom_rstput( 0, 0, inum0, 'e2f', e2f, ktype = jp_r8 )92 93 CALL iom_rstput( 0, 0, inum0, 'e1e2u', e1e2u, ktype = jp_r8 )94 CALL iom_rstput( 0, 0, inum0, 'e1e2v', e1e2v, ktype = jp_r8 )95 96 CALL iom_close( inum0 )97 !98 IF( nn_timing == 1 ) CALL timing_stop('dom_wri_coordinate')99 !100 END SUBROUTINE dom_wri_coordinate101 102 43 103 44 SUBROUTINE dom_wri … … 113 54 !! domhgr, domzgr, and dommsk. Note: the file contain depends on 114 55 !! the vertical coord. used (z-coord, partial steps, s-coord) 115 !! MOD(n msh, 3) = 1 : 'mesh_mask.nc' file56 !! MOD(nn_msh, 3) = 1 : 'mesh_mask.nc' file 116 57 !! = 2 : 'mesh.nc' and mask.nc' files 117 58 !! = 0 : 'mesh_hgr.nc', 'mesh_zgr.nc' and … … 120 61 !! vertical coordinate. 121 62 !! 122 !! if n msh <= 3: write full 3D arrays for e3[tuvw] and gdep[tuvw]123 !! if 3 < n msh <= 6: write full 3D arrays for e3[tuvw] and 2D arrays63 !! if nn_msh <= 3: write full 3D arrays for e3[tuvw] and gdep[tuvw] 64 !! if 3 < nn_msh <= 6: write full 3D arrays for e3[tuvw] and 2D arrays 124 65 !! corresponding to the depth of the bottom t- and w-points 125 !! if 6 < n msh <= 9: write 2D arrays corresponding to the depth and the66 !! if 6 < nn_msh <= 9: write 2D arrays corresponding to the depth and the 126 67 !! thickness (e3[tw]_ps) of the bottom points 127 68 !! … … 129 70 !! masks, depth and vertical scale factors 130 71 !!---------------------------------------------------------------------- 131 !! 132 INTEGER :: inum0 ! temprary units for 'mesh_mask.nc' file 133 INTEGER :: inum1 ! temprary units for 'mesh.nc' file 134 INTEGER :: inum2 ! temprary units for 'mask.nc' file 135 INTEGER :: inum3 ! temprary units for 'mesh_hgr.nc' file 136 INTEGER :: inum4 ! temprary units for 'mesh_zgr.nc' file 137 CHARACTER(len=21) :: clnam0 ! filename (mesh and mask informations) 138 CHARACTER(len=21) :: clnam1 ! filename (mesh informations) 139 CHARACTER(len=21) :: clnam2 ! filename (mask informations) 140 CHARACTER(len=21) :: clnam3 ! filename (horizontal mesh informations) 141 CHARACTER(len=21) :: clnam4 ! filename (vertical mesh informations) 72 INTEGER :: inum ! temprary units for 'mesh_mask.nc' file 73 CHARACTER(len=21) :: clnam ! filename (mesh and mask informations) 142 74 INTEGER :: ji, jj, jk ! dummy loop indices 143 ! ! workspaces 144 REAL(wp), POINTER, DIMENSION(:,: ) :: zprt, zprw 145 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepu, zdepv 75 INTEGER :: izco, izps, isco, icav 76 ! 77 REAL(wp), POINTER, DIMENSION(:,:) :: zprt, zprw ! 2D workspace 78 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepu, zdepv ! 3D workspace 146 79 !!---------------------------------------------------------------------- 147 80 ! 148 81 IF( nn_timing == 1 ) CALL timing_start('dom_wri') 149 82 ! 150 CALL wrk_alloc( jpi, jpj, zprt, zprw)151 CALL wrk_alloc( jpi, jpj, jpk,zdepu, zdepv )83 CALL wrk_alloc( jpi,jpj, zprt , zprw ) 84 CALL wrk_alloc( jpi,jpj,jpk, zdepu, zdepv ) 152 85 ! 153 86 IF(lwp) WRITE(numout,*) … … 155 88 IF(lwp) WRITE(numout,*) '~~~~~~~' 156 89 157 clnam0 = 'mesh_mask' ! filename (mesh and mask informations) 158 clnam1 = 'mesh' ! filename (mesh informations) 159 clnam2 = 'mask' ! filename (mask informations) 160 clnam3 = 'mesh_hgr' ! filename (horizontal mesh informations) 161 clnam4 = 'mesh_zgr' ! filename (vertical mesh informations) 162 163 SELECT CASE ( MOD(nmsh, 3) ) 164 ! ! ============================ 165 CASE ( 1 ) ! create 'mesh_mask.nc' file 166 ! ! ============================ 167 CALL iom_open( TRIM(clnam0), inum0, ldwrt = .TRUE., kiolib = jprstlib ) 168 inum2 = inum0 ! put all the informations 169 inum3 = inum0 ! in unit inum0 170 inum4 = inum0 171 172 ! ! ============================ 173 CASE ( 2 ) ! create 'mesh.nc' and 174 ! ! 'mask.nc' files 175 ! ! ============================ 176 CALL iom_open( TRIM(clnam1), inum1, ldwrt = .TRUE., kiolib = jprstlib ) 177 CALL iom_open( TRIM(clnam2), inum2, ldwrt = .TRUE., kiolib = jprstlib ) 178 inum3 = inum1 ! put mesh informations 179 inum4 = inum1 ! in unit inum1 180 ! ! ============================ 181 CASE ( 0 ) ! create 'mesh_hgr.nc' 182 ! ! 'mesh_zgr.nc' and 183 ! ! 'mask.nc' files 184 ! ! ============================ 185 CALL iom_open( TRIM(clnam2), inum2, ldwrt = .TRUE., kiolib = jprstlib ) 186 CALL iom_open( TRIM(clnam3), inum3, ldwrt = .TRUE., kiolib = jprstlib ) 187 CALL iom_open( TRIM(clnam4), inum4, ldwrt = .TRUE., kiolib = jprstlib ) 188 ! 189 END SELECT 190 191 ! ! masks (inum2) 192 CALL iom_rstput( 0, 0, inum2, 'tmask', tmask, ktype = jp_i1 ) ! ! land-sea mask 193 CALL iom_rstput( 0, 0, inum2, 'umask', umask, ktype = jp_i1 ) 194 CALL iom_rstput( 0, 0, inum2, 'vmask', vmask, ktype = jp_i1 ) 195 CALL iom_rstput( 0, 0, inum2, 'fmask', fmask, ktype = jp_i1 ) 90 clnam = 'mesh_mask' ! filename (mesh and mask informations) 91 92 ! ! ============================ 93 ! ! create 'mesh_mask.nc' file 94 ! ! ============================ 95 CALL iom_open( TRIM(clnam), inum, ldwrt = .TRUE., kiolib = jprstlib ) 96 ! 97 ! ! global domain size 98 CALL iom_rstput( 0, 0, inum, 'jpiglo', REAL( jpiglo, wp), ktype = jp_i4 ) 99 CALL iom_rstput( 0, 0, inum, 'jpjglo', REAL( jpjglo, wp), ktype = jp_i4 ) 100 CALL iom_rstput( 0, 0, inum, 'jpkglo', REAL( jpkglo, wp), ktype = jp_i4 ) 101 102 ! ! domain characteristics 103 CALL iom_rstput( 0, 0, inum, 'jperio', REAL( jperio, wp), ktype = jp_i4 ) 104 ! ! type of vertical coordinate 105 IF( ln_zco ) THEN ; izco = 1 ; ELSE ; izco = 0 ; ENDIF 106 IF( ln_zps ) THEN ; izps = 1 ; ELSE ; izps = 0 ; ENDIF 107 IF( ln_sco ) THEN ; isco = 1 ; ELSE ; isco = 0 ; ENDIF 108 CALL iom_rstput( 0, 0, inum, 'ln_zco' , REAL( izco, wp), ktype = jp_i4 ) 109 CALL iom_rstput( 0, 0, inum, 'ln_zps' , REAL( izps, wp), ktype = jp_i4 ) 110 CALL iom_rstput( 0, 0, inum, 'ln_sco' , REAL( isco, wp), ktype = jp_i4 ) 111 ! ! ocean cavities under iceshelves 112 IF( ln_isfcav ) THEN ; icav = 1 ; ELSE ; icav = 0 ; ENDIF 113 CALL iom_rstput( 0, 0, inum, 'ln_isfcav', REAL( icav, wp), ktype = jp_i4 ) 114 115 ! ! masks 116 CALL iom_rstput( 0, 0, inum, 'tmask', tmask, ktype = jp_i1 ) ! ! land-sea mask 117 CALL iom_rstput( 0, 0, inum, 'umask', umask, ktype = jp_i1 ) 118 CALL iom_rstput( 0, 0, inum, 'vmask', vmask, ktype = jp_i1 ) 119 CALL iom_rstput( 0, 0, inum, 'fmask', fmask, ktype = jp_i1 ) 196 120 197 121 CALL dom_uniq( zprw, 'T' ) 198 122 DO jj = 1, jpj 199 123 DO ji = 1, jpi 200 jk=mikt(ji,jj) 201 zprt(ji,jj) = tmask(ji,jj,jk) * zprw(ji,jj) ! ! unique point mask 124 zprt(ji,jj) = ssmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 202 125 END DO 203 126 END DO ! ! unique point mask 204 CALL iom_rstput( 0, 0, inum 2, 'tmaskutil', zprt, ktype = jp_i1 )127 CALL iom_rstput( 0, 0, inum, 'tmaskutil', zprt, ktype = jp_i1 ) 205 128 CALL dom_uniq( zprw, 'U' ) 206 129 DO jj = 1, jpj 207 130 DO ji = 1, jpi 208 jk=miku(ji,jj) 209 zprt(ji,jj) = umask(ji,jj,jk) * zprw(ji,jj) ! ! unique point mask 131 zprt(ji,jj) = ssumask(ji,jj) * zprw(ji,jj) ! ! unique point mask 210 132 END DO 211 133 END DO 212 CALL iom_rstput( 0, 0, inum 2, 'umaskutil', zprt, ktype = jp_i1 )134 CALL iom_rstput( 0, 0, inum, 'umaskutil', zprt, ktype = jp_i1 ) 213 135 CALL dom_uniq( zprw, 'V' ) 214 136 DO jj = 1, jpj 215 137 DO ji = 1, jpi 216 jk=mikv(ji,jj) 217 zprt(ji,jj) = vmask(ji,jj,jk) * zprw(ji,jj) ! ! unique point mask 138 zprt(ji,jj) = ssvmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 218 139 END DO 219 140 END DO 220 CALL iom_rstput( 0, 0, inum2, 'vmaskutil', zprt, ktype = jp_i1 ) 221 CALL dom_uniq( zprw, 'F' ) 222 DO jj = 1, jpj 223 DO ji = 1, jpi 224 jk=mikf(ji,jj) 225 zprt(ji,jj) = fmask(ji,jj,jk) * zprw(ji,jj) ! ! unique point mask 226 END DO 227 END DO 228 CALL iom_rstput( 0, 0, inum2, 'fmaskutil', zprt, ktype = jp_i1 ) 141 CALL iom_rstput( 0, 0, inum, 'vmaskutil', zprt, ktype = jp_i1 ) 142 !!gm ssfmask has been removed ==>> find another solution to defined fmaskutil 143 !! Here we just remove the output of fmaskutil. 144 ! CALL dom_uniq( zprw, 'F' ) 145 ! DO jj = 1, jpj 146 ! DO ji = 1, jpi 147 ! zprt(ji,jj) = ssfmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 148 ! END DO 149 ! END DO 150 ! CALL iom_rstput( 0, 0, inum, 'fmaskutil', zprt, ktype = jp_i1 ) 151 !!gm 229 152 230 153 ! ! horizontal mesh (inum3) 231 CALL iom_rstput( 0, 0, inum3, 'glamt', glamt, ktype = jp_r4 ) ! ! latitude 232 CALL iom_rstput( 0, 0, inum3, 'glamu', glamu, ktype = jp_r4 ) 233 CALL iom_rstput( 0, 0, inum3, 'glamv', glamv, ktype = jp_r4 ) 234 CALL iom_rstput( 0, 0, inum3, 'glamf', glamf, ktype = jp_r4 ) 235 236 CALL iom_rstput( 0, 0, inum3, 'gphit', gphit, ktype = jp_r4 ) ! ! longitude 237 CALL iom_rstput( 0, 0, inum3, 'gphiu', gphiu, ktype = jp_r4 ) 238 CALL iom_rstput( 0, 0, inum3, 'gphiv', gphiv, ktype = jp_r4 ) 239 CALL iom_rstput( 0, 0, inum3, 'gphif', gphif, ktype = jp_r4 ) 240 241 CALL iom_rstput( 0, 0, inum3, 'e1t', e1t, ktype = jp_r8 ) ! ! e1 scale factors 242 CALL iom_rstput( 0, 0, inum3, 'e1u', e1u, ktype = jp_r8 ) 243 CALL iom_rstput( 0, 0, inum3, 'e1v', e1v, ktype = jp_r8 ) 244 CALL iom_rstput( 0, 0, inum3, 'e1f', e1f, ktype = jp_r8 ) 245 246 CALL iom_rstput( 0, 0, inum3, 'e2t', e2t, ktype = jp_r8 ) ! ! e2 scale factors 247 CALL iom_rstput( 0, 0, inum3, 'e2u', e2u, ktype = jp_r8 ) 248 CALL iom_rstput( 0, 0, inum3, 'e2v', e2v, ktype = jp_r8 ) 249 CALL iom_rstput( 0, 0, inum3, 'e2f', e2f, ktype = jp_r8 ) 250 251 CALL iom_rstput( 0, 0, inum3, 'ff', ff, ktype = jp_r8 ) ! ! coriolis factor 154 CALL iom_rstput( 0, 0, inum, 'glamt', glamt, ktype = jp_r8 ) ! ! latitude 155 CALL iom_rstput( 0, 0, inum, 'glamu', glamu, ktype = jp_r8 ) 156 CALL iom_rstput( 0, 0, inum, 'glamv', glamv, ktype = jp_r8 ) 157 CALL iom_rstput( 0, 0, inum, 'glamf', glamf, ktype = jp_r8 ) 158 159 CALL iom_rstput( 0, 0, inum, 'gphit', gphit, ktype = jp_r8 ) ! ! longitude 160 CALL iom_rstput( 0, 0, inum, 'gphiu', gphiu, ktype = jp_r8 ) 161 CALL iom_rstput( 0, 0, inum, 'gphiv', gphiv, ktype = jp_r8 ) 162 CALL iom_rstput( 0, 0, inum, 'gphif', gphif, ktype = jp_r8 ) 163 164 CALL iom_rstput( 0, 0, inum, 'e1t', e1t, ktype = jp_r8 ) ! ! e1 scale factors 165 CALL iom_rstput( 0, 0, inum, 'e1u', e1u, ktype = jp_r8 ) 166 CALL iom_rstput( 0, 0, inum, 'e1v', e1v, ktype = jp_r8 ) 167 CALL iom_rstput( 0, 0, inum, 'e1f', e1f, ktype = jp_r8 ) 168 169 CALL iom_rstput( 0, 0, inum, 'e2t', e2t, ktype = jp_r8 ) ! ! e2 scale factors 170 CALL iom_rstput( 0, 0, inum, 'e2u', e2u, ktype = jp_r8 ) 171 CALL iom_rstput( 0, 0, inum, 'e2v', e2v, ktype = jp_r8 ) 172 CALL iom_rstput( 0, 0, inum, 'e2f', e2f, ktype = jp_r8 ) 173 174 CALL iom_rstput( 0, 0, inum, 'ff_f', ff_f, ktype = jp_r8 ) ! ! coriolis factor 175 CALL iom_rstput( 0, 0, inum, 'ff_t', ff_t, ktype = jp_r8 ) 252 176 253 177 ! note that mbkt is set to 1 over land ==> use surface tmask 254 178 zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) , wp ) 255 CALL iom_rstput( 0, 0, inum4, 'mbathy', zprt, ktype = jp_i2 ) ! ! nb of ocean T-points 256 CALL iom_rstput( 0, 0, inum4, 'bathy', bathy, ktype = jp_r8 ) ! ! nb of ocean T-points 179 CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 ) ! ! nb of ocean T-points 257 180 zprt(:,:) = ssmask(:,:) * REAL( mikt(:,:) , wp ) 258 CALL iom_rstput( 0, 0, inum 4, 'misf', zprt, ktype = jp_i2) ! ! nb of ocean T-points181 CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 ) ! ! nb of ocean T-points 259 182 zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp ) 260 CALL iom_rstput( 0, 0, inum4, 'isfdraft', zprt, ktype = jp_r4 ) ! ! nb of ocean T-points 261 262 IF( ln_sco ) THEN ! s-coordinate 263 CALL iom_rstput( 0, 0, inum4, 'hbatt', hbatt ) 264 CALL iom_rstput( 0, 0, inum4, 'hbatu', hbatu ) 265 CALL iom_rstput( 0, 0, inum4, 'hbatv', hbatv ) 266 CALL iom_rstput( 0, 0, inum4, 'hbatf', hbatf ) 267 ! 268 CALL iom_rstput( 0, 0, inum4, 'gsigt', gsigt ) ! ! scaling coef. 269 CALL iom_rstput( 0, 0, inum4, 'gsigw', gsigw ) 270 CALL iom_rstput( 0, 0, inum4, 'gsi3w', gsi3w ) 271 CALL iom_rstput( 0, 0, inum4, 'esigt', esigt ) 272 CALL iom_rstput( 0, 0, inum4, 'esigw', esigw ) 273 ! 274 CALL iom_rstput( 0, 0, inum4, 'e3t_0', e3t_0 ) ! ! scale factors 275 CALL iom_rstput( 0, 0, inum4, 'e3u_0', e3u_0 ) 276 CALL iom_rstput( 0, 0, inum4, 'e3v_0', e3v_0 ) 277 CALL iom_rstput( 0, 0, inum4, 'e3w_0', e3w_0 ) 278 CALL iom_rstput( 0, 0, inum4, 'rx1', rx1 ) ! ! Max. grid stiffness ratio 279 ! 280 CALL iom_rstput( 0, 0, inum4, 'gdept_1d' , gdept_1d ) ! ! stretched system 281 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d' , gdepw_1d ) 282 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r4 ) 283 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r4 ) 284 ENDIF 285 286 IF( ln_zps ) THEN ! z-coordinate - partial steps 287 ! 288 IF( nmsh <= 6 ) THEN ! ! 3D vertical scale factors 289 CALL iom_rstput( 0, 0, inum4, 'e3t_0', e3t_0 ) 290 CALL iom_rstput( 0, 0, inum4, 'e3u_0', e3u_0 ) 291 CALL iom_rstput( 0, 0, inum4, 'e3v_0', e3v_0 ) 292 CALL iom_rstput( 0, 0, inum4, 'e3w_0', e3w_0 ) 293 ELSE ! ! 2D masked bottom ocean scale factors 294 DO jj = 1,jpj 295 DO ji = 1,jpi 296 e3tp(ji,jj) = e3t_0(ji,jj,mbkt(ji,jj)) * ssmask(ji,jj) 297 e3wp(ji,jj) = e3w_0(ji,jj,mbkt(ji,jj)) * ssmask(ji,jj) 298 END DO 299 END DO 300 CALL iom_rstput( 0, 0, inum4, 'e3t_ps', e3tp ) 301 CALL iom_rstput( 0, 0, inum4, 'e3w_ps', e3wp ) 302 END IF 303 ! 304 IF( nmsh <= 3 ) THEN ! ! 3D depth 305 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r4 ) 306 DO jk = 1,jpk 307 DO jj = 1, jpjm1 308 DO ji = 1, fs_jpim1 ! vector opt. 309 zdepu(ji,jj,jk) = MIN( gdept_0(ji,jj,jk) , gdept_0(ji+1,jj ,jk) ) 310 zdepv(ji,jj,jk) = MIN( gdept_0(ji,jj,jk) , gdept_0(ji ,jj+1,jk) ) 311 END DO 312 END DO 313 END DO 314 CALL lbc_lnk( zdepu, 'U', 1. ) ; CALL lbc_lnk( zdepv, 'V', 1. ) 315 CALL iom_rstput( 0, 0, inum4, 'gdepu', zdepu, ktype = jp_r4 ) 316 CALL iom_rstput( 0, 0, inum4, 'gdepv', zdepv, ktype = jp_r4 ) 317 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r4 ) 318 ELSE ! ! 2D bottom depth 319 DO jj = 1,jpj 320 DO ji = 1,jpi 321 zprt(ji,jj) = gdept_0(ji,jj,mbkt(ji,jj) ) * ssmask(ji,jj) 322 zprw(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1) * ssmask(ji,jj) 323 END DO 324 END DO 325 CALL iom_rstput( 0, 0, inum4, 'hdept', zprt, ktype = jp_r4 ) 326 CALL iom_rstput( 0, 0, inum4, 'hdepw', zprw, ktype = jp_r4 ) 327 ENDIF 328 ! 329 CALL iom_rstput( 0, 0, inum4, 'gdept_1d', gdept_1d ) ! ! reference z-coord. 330 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d', gdepw_1d ) 331 CALL iom_rstput( 0, 0, inum4, 'e3t_1d' , e3t_1d ) 332 CALL iom_rstput( 0, 0, inum4, 'e3w_1d' , e3w_1d ) 333 ENDIF 334 335 IF( ln_zco ) THEN 336 ! ! z-coordinate - full steps 337 CALL iom_rstput( 0, 0, inum4, 'gdept_1d', gdept_1d ) ! ! depth 338 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d', gdepw_1d ) 339 CALL iom_rstput( 0, 0, inum4, 'e3t_1d' , e3t_1d ) ! ! scale factors 340 CALL iom_rstput( 0, 0, inum4, 'e3w_1d' , e3w_1d ) 183 CALL iom_rstput( 0, 0, inum, 'isfdraft', zprt, ktype = jp_r8 ) ! ! nb of ocean T-points 184 ! ! vertical mesh 185 CALL iom_rstput( 0, 0, inum, 'e3t_0', e3t_0, ktype = jp_r8 ) ! ! scale factors 186 CALL iom_rstput( 0, 0, inum, 'e3u_0', e3u_0, ktype = jp_r8 ) 187 CALL iom_rstput( 0, 0, inum, 'e3v_0', e3v_0, ktype = jp_r8 ) 188 CALL iom_rstput( 0, 0, inum, 'e3w_0', e3w_0, ktype = jp_r8 ) 189 ! 190 CALL iom_rstput( 0, 0, inum, 'gdept_1d' , gdept_1d , ktype = jp_r8 ) ! stretched system 191 CALL iom_rstput( 0, 0, inum, 'gdepw_1d' , gdepw_1d , ktype = jp_r8 ) 192 CALL iom_rstput( 0, 0, inum, 'gdept_0' , gdept_0 , ktype = jp_r8 ) 193 CALL iom_rstput( 0, 0, inum, 'gdepw_0' , gdepw_0 , ktype = jp_r8 ) 194 ! 195 IF( ln_sco ) THEN ! s-coordinate stiffness 196 CALL dom_stiff( zprt ) 197 CALL iom_rstput( 0, 0, inum, 'stiffness', zprt ) ! ! Max. grid stiffness ratio 341 198 ENDIF 342 199 ! ! ============================ 343 !! close the files200 CALL iom_close( inum ) ! close the files 344 201 ! ! ============================ 345 SELECT CASE ( MOD(nmsh, 3) )346 CASE ( 1 )347 CALL iom_close( inum0 )348 CASE ( 2 )349 CALL iom_close( inum1 )350 CALL iom_close( inum2 )351 CASE ( 0 )352 CALL iom_close( inum2 )353 CALL iom_close( inum3 )354 CALL iom_close( inum4 )355 END SELECT356 202 ! 357 203 CALL wrk_dealloc( jpi, jpj, zprt, zprw ) … … 372 218 !! 2) check which elements have been changed 373 219 !!---------------------------------------------------------------------- 374 !375 220 CHARACTER(len=1) , INTENT(in ) :: cdgrd ! 376 221 REAL(wp), DIMENSION(:,:), INTENT(inout) :: puniq ! … … 406 251 END SUBROUTINE dom_uniq 407 252 253 254 SUBROUTINE dom_stiff( px1 ) 255 !!---------------------------------------------------------------------- 256 !! *** ROUTINE dom_stiff *** 257 !! 258 !! ** Purpose : Diagnose maximum grid stiffness/hydrostatic consistency 259 !! 260 !! ** Method : Compute Haney (1991) hydrostatic condition ratio 261 !! Save the maximum in the vertical direction 262 !! (this number is only relevant in s-coordinates) 263 !! 264 !! Haney, 1991, J. Phys. Oceanogr., 21, 610-619. 265 !!---------------------------------------------------------------------- 266 REAL(wp), DIMENSION(:,:), INTENT(out), OPTIONAL :: px1 ! stiffness 267 ! 268 INTEGER :: ji, jj, jk 269 REAL(wp) :: zrxmax 270 REAL(wp), DIMENSION(4) :: zr1 271 REAL(wp), DIMENSION(jpi,jpj) :: zx1 272 !!---------------------------------------------------------------------- 273 zx1(:,:) = 0._wp 274 zrxmax = 0._wp 275 zr1(:) = 0._wp 276 ! 277 DO ji = 2, jpim1 278 DO jj = 2, jpjm1 279 DO jk = 1, jpkm1 280 !!gm remark: dk(gdepw) = e3t ===>>> possible simplification of the following calculation.... 281 !! especially since it is gde3w which is used to compute the pressure gradient 282 !! furthermore, I think gdept_0 should be used below instead of w point in the numerator 283 !! so that the ratio is computed at the same point (i.e. uw and vw) .... 284 zr1(1) = ABS( ( gdepw_0(ji ,jj,jk )-gdepw_0(ji-1,jj,jk ) & 285 & +gdepw_0(ji ,jj,jk+1)-gdepw_0(ji-1,jj,jk+1) ) & 286 & / ( gdepw_0(ji ,jj,jk )+gdepw_0(ji-1,jj,jk ) & 287 & -gdepw_0(ji ,jj,jk+1)-gdepw_0(ji-1,jj,jk+1) + rsmall ) ) * umask(ji-1,jj,jk) 288 zr1(2) = ABS( ( gdepw_0(ji+1,jj,jk )-gdepw_0(ji ,jj,jk ) & 289 & +gdepw_0(ji+1,jj,jk+1)-gdepw_0(ji ,jj,jk+1) ) & 290 & / ( gdepw_0(ji+1,jj,jk )+gdepw_0(ji ,jj,jk ) & 291 & -gdepw_0(ji+1,jj,jk+1)-gdepw_0(ji ,jj,jk+1) + rsmall ) ) * umask(ji ,jj,jk) 292 zr1(3) = ABS( ( gdepw_0(ji,jj+1,jk )-gdepw_0(ji,jj ,jk ) & 293 & +gdepw_0(ji,jj+1,jk+1)-gdepw_0(ji,jj ,jk+1) ) & 294 & / ( gdepw_0(ji,jj+1,jk )+gdepw_0(ji,jj ,jk ) & 295 & -gdepw_0(ji,jj+1,jk+1)-gdepw_0(ji,jj ,jk+1) + rsmall ) ) * vmask(ji,jj ,jk) 296 zr1(4) = ABS( ( gdepw_0(ji,jj ,jk )-gdepw_0(ji,jj-1,jk ) & 297 & +gdepw_0(ji,jj ,jk+1)-gdepw_0(ji,jj-1,jk+1) ) & 298 & / ( gdepw_0(ji,jj ,jk )+gdepw_0(ji,jj-1,jk ) & 299 & -gdepw_0(ji,jj ,jk+1)-gdepw_0(ji,jj-1,jk+1) + rsmall ) ) * vmask(ji,jj-1,jk) 300 zrxmax = MAXVAL( zr1(1:4) ) 301 zx1(ji,jj) = MAX( zx1(ji,jj) , zrxmax ) 302 END DO 303 END DO 304 END DO 305 CALL lbc_lnk( zx1, 'T', 1. ) 306 ! 307 IF( PRESENT( px1 ) ) px1 = zx1 308 ! 309 zrxmax = MAXVAL( zx1 ) 310 ! 311 IF( lk_mpp ) CALL mpp_max( zrxmax ) ! max over the global domain 312 ! 313 IF(lwp) THEN 314 WRITE(numout,*) 315 WRITE(numout,*) 'dom_stiff : maximum grid stiffness ratio: ', zrxmax 316 WRITE(numout,*) '~~~~~~~~~' 317 ENDIF 318 ! 319 END SUBROUTINE dom_stiff 320 408 321 !!====================================================================== 409 322 END MODULE domwri
Note: See TracChangeset
for help on using the changeset viewer.