Changeset 12063 for NEMO/branches/2019/dev_ASINTER-01-05_merged/tests
- Timestamp:
- 2019-12-05T11:46:38+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_ASINTER-01-05_merged/tests
- Files:
-
- 9 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_cen-ahm1000_cfg
r10075 r12063 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_ubs_cfg
r10075 r12063 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_ubs_cfg
r10075 r12063 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/OVERFLOW/EXPREF/namelist_zps_FCT2_flux_ubs_cfg
r10075 r12063 211 211 ! ! coefficients 212 212 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 213 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)213 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 214 214 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 215 215 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg
r11586 r12063 34 34 cn_exp = "OVF_zps_FCT4_flux_ubs" ! experience name 35 35 nn_it000 = 1 ! first time step 36 !nn_itend = 6120 ! here 17h of simulation (=6120 time-step) 37 nn_itend = 5760 ! here 16h of simulation (=5760 time-step) abort after 5802 for zps: pb of physics conditions 36 nn_itend = 6120 ! here 17h of simulation (=6120 time-step) 38 37 nn_istate = 0 ! output the initial state (1) or not (0) 39 nn_stock = 1080 ! frequency of creation of a restart file (modulo referenced to 1)40 nn_write = 1080 ! frequency of write in the output file (modulo referenced to nn_it000)38 nn_stock = 6120 ! frequency of creation of a restart file (modulo referenced to 1) 39 nn_write = 6120 ! frequency of write in the output file (modulo referenced to nn_it000) 41 40 / 42 41 !----------------------------------------------------------------------- … … 255 254 &namzdf ! vertical physics manager (default: NO selection) 256 255 !----------------------------------------------------------------------- 256 ! ! adaptive-implicit vertical advection 257 ln_zad_Aimp = .true. ! Courant number dependent scheme (Shchepetkin 2015) 258 ! 257 259 ! ! type of vertical closure (required) 258 260 ln_zdfcst = .true. ! constant mixing … … 279 281 ! ! coefficients 280 282 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 281 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)283 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 282 284 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 283 285 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_vect_een_cfg
r10075 r12063 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/README.rst
r10605 r12063 3 3 ********************** 4 4 5 The complete and up-to-date set of test cases is available on 6 `NEMO test cases Github repository <http://github.com/NEMO-ocean/NEMO-examples>`_. 7 8 Download it directly in the ``./tests`` root directory with 5 .. todo:: 6 7 CANAL animated gif is missing 8 9 .. contents:: 10 :local: 11 :depth: 1 12 13 Installation 14 ============ 15 16 Download 17 -------- 18 19 | The complete and up-to-date set of test cases is available on 20 :github:`NEMO test cases repository <NEMO-examples>`. 21 | Download it directly into the :file:`./tests` root directory with 9 22 10 23 .. code-block:: console … … 12 25 $ git clone http://github.com/NEMO-ocean/NEMO-examples 13 26 14 .. contents:: 15 :local: 16 17 Procedure 18 ========= 19 20 Compile test cases 21 ------------------ 22 23 The compilation of the test cases is very similar to the manner the reference configurations are compiled. 24 If you are not familiar on how to compile NEMO, it is first recomended to read :doc:`the instructions <install>` 25 26 | In the same manner as the ref. cfg are compiled with '-r' option, test cases can be compile by the use of makenemo with '-a' option. 27 | Here an example to compile a copy named WAD2 of the wetting and drying test case (WAD) on the macport_osx architecture on 4 cores: 27 Compilation 28 ----------- 29 30 The compilation of the test cases is very similar to 31 the manner the reference configurations are compiled. 32 If you are not familiar on how to compile NEMO, 33 it is first recomended to read :doc:`the instructions <install>`. 34 35 | As the reference configurations are compiled with ``-r`` option, 36 test cases can be compiled by the use of :file:`makenemo` with ``-a`` option. 37 | Here an example to compile a copy named WAD2 of the wetting and drying test case (WAD): 28 38 29 39 .. code-block:: console 30 40 31 $ ./makenemo -n 'WAD2' -a 'WAD' -m 'macport_osx' -j '4' 32 33 Run and analyse the test cases 34 ------------------------------ 35 36 There no requirement of specific input file for the test_cases presented here. The XIOS xml input files and namelist are already setup correctly. 41 $ ./makenemo -n 'WAD2' -a 'WAD' -m 'my_arch' -j '4' 42 43 Run and analysis 44 ---------------- 45 46 There no requirement of specific input file for the test_cases presented here. 47 The XIOS xml input files and namelist are already setup correctly. 37 48 For detailed description and Jupyter notebook, the reader is directed on 38 the `NEMO test cases repository <http://github.com/NEMO-ocean/NEMO-examples>`_49 the :github:`NEMO test cases repository <NEMO-examples>` 39 50 40 51 The description below is a brief advertisement of some test cases. 41 52 53 List of test cases 54 ================== 55 42 56 ICE_AGRIF 43 ========= 44 45 This test case illustrates the advection of an ice patch across an East/West and North/South periodic channel 46 over a slab ocean (i.e. one ocean layer), and with an AGRIF zoom (1:3) in the center 47 The purpose of this configuration is to test the advection of the ice patch in 48 and across the AGRIF boundary 49 One can either impose ice velocities or ice-atm. stresses and let rheology define velocities 50 (see README for details) 51 52 .. image:: _static/ICE_AGRIF_UDIAG_43days_UM5.gif 57 --------- 58 59 .. figure:: _static/ICE_AGRIF_UDIAG_43days_UM5.gif 60 :width: 200px 61 :align: left 62 63 .. 64 65 | This test case illustrates the advection of an ice patch across 66 an East/West and North/South periodic channel over a slab ocean (i.e. one ocean layer), 67 and with an AGRIF zoom (1:3) in the center. 68 | The purpose of this configuration is to 69 test the advection of the ice patch in and across the AGRIF boundary. 70 One can either impose ice velocities or ice-atm. 71 Stresses and let rheology define velocities (see :file:`README` for details) 53 72 54 73 VORTEX 55 ====== 56 57 This test case illustrates the propagation of an anticyclonic eddy over a Beta plan and a flat bottom. 58 It is implemented here with an online refined subdomain (1:3) out of which the vortex propagates. 59 It serves as a benchmark for quantitative estimates of nesting errors as in Debreu et al. (2012) :cite:`DEBREU2012`, 60 Penven et al. (2006) :cite:`PENVEN2006` or Spall and Holland (1991) :cite:`SPALL1991`. 61 62 The animation below (sea level anomaly in meters) illustrates with two 1:2 successively nested grids how 63 the vortex smoothly propagates out of the refined grids. 64 65 .. image:: _static/VORTEX_anim.gif 74 ------ 75 76 .. figure:: _static/VORTEX_anim.gif 77 :width: 200px 78 :align: right 79 80 .. 81 82 This test case illustrates the propagation of an anticyclonic eddy over a Beta plan and a flat bottom. 83 It is implemented here with an online refined subdomain (1:3) out of which the vortex propagates. 84 It serves as a benchmark for quantitative estimates of nesting errors as in :cite:`DEBREU2012`, 85 :cite:`PENVEN2006` or :cite:`SPALL1991`. 86 87 The animation (sea level anomaly in meters) illustrates with 88 two 1:2 successively nested grids how the vortex smoothly propagates out of the refined grids. 66 89 67 90 ISOMIP 68 ====== 69 70 The purpose of this test case is to evaluate the impact of various schemes and new development with the iceshelf cavities circulation and melt. 71 This configuration served as initial assesment of the ice shelf module in Losh et al. (2008) :cite:`LOSCH2008` and Mathiot et al. (2017) :cite:`MATHIOT2017`. 72 The default setup is the one described `here <http://staff.acecrc.org.au/~bkgalton/ISOMIP/test_cavities.pdf>`_. 73 74 The figure below (meridional overturning circulation) illustrates the circulation generated after 10000 days by the ice shelf melting (ice pump). 75 76 .. image:: _static/ISOMIP_moc.png 91 ------ 92 93 .. figure:: _static/ISOMIP_moc.png 94 :width: 200px 95 :align: left 96 97 .. 98 99 | The purpose of this test case is to evaluate the impact of various schemes and new development with 100 the iceshelf cavities circulation and melt. 101 This configuration served as initial assesment of the ice shelf module in :cite:`LOSCH2008` and 102 :cite:`MATHIOT2017`. 103 The default setup is the one described |ISOMIP|_. 104 | The figure (meridional overturning circulation) illustrates 105 the circulation generated after 10000 days by the ice shelf melting (ice pump). 106 107 .. |ISOMIP| replace:: here 77 108 78 109 LOCK_EXCHANGE 79 ============= 80 81 The LOCK EXCHANGE experiment is a classical fluid dynamics experiment that has been adapted 82 by Haidvogel and Beckmann (1999) :cite:`HAIDVOGEL1999` for testing advection schemes in ocean circulation models. 83 It has been used by several authors including Burchard and Bolding (2002) :cite:`BURCHARD2002` and Ilicak et al. (2012) :cite:`ILICAK2012`. 84 The LOCK EXCHANGE experiment can in particular illustrate the impact of different choices of numerical schemes 85 and/or subgrid closures on spurious interior mixing. 86 87 Below the animation of the LOCK_EXCHANGE test case using the advection scheme FCT4 (forth order) for tracer and ubs for dynamics. 88 89 .. image:: _static/LOCK-FCT4_flux_ubs.gif 110 ------------- 111 112 .. figure:: _static/LOCK-FCT4_flux_ubs.gif 113 :width: 200px 114 :align: right 115 116 .. 117 118 | The LOCK EXCHANGE experiment is a classical fluid dynamics experiment that has been adapted 119 by :cite:`HAIDVOGEL1999` for testing advection schemes in ocean circulation models. 120 It has been used by several authors including :cite:`BURCHARD2002` and :cite:`ILICAK2012`. 121 The LOCK EXCHANGE experiment can in particular illustrate 122 the impact of different choices of numerical schemes and/or subgrid closures on 123 spurious interior mixing. 124 | Here the animation of the LOCK_EXCHANGE test case using 125 the advection scheme FCT4 (forth order) for tracer and ubs for dynamics. 90 126 91 127 OVERFLOW 92 ======== 93 94 The OVERFLOW experiment illustrates the impact of different choices of numerical schemes 95 and/or subgrid closures on spurious interior mixing close to bottom topography. 96 The OVERFLOW experiment is adapted from the non-rotating overflow configuration described 97 in Haidvogel and Beckmann (1999) :cite:`HAIDVOGEL1999` and further used by Ilicak et al. (2012) :cite:`ILICAK2012`. 98 Here we can assess the behaviour of the second-order tracer advection scheme FCT2 and fortht-order FCT4, z-coordinate and sigma coordinate (...). 99 100 Below the animation of the OVERFLOW test case in sigma coordinate with the forth-order advection scheme FCT4. 101 102 .. image:: _static/OVF-sco_FCT4_flux_cen-ahm1000.gif 128 -------- 129 130 .. figure:: _static/OVF-sco_FCT4_flux_cen-ahm1000.gif 131 :width: 200px 132 :align: left 133 134 .. 135 136 | The OVERFLOW experiment illustrates the impact of different choices of numerical schemes and/or 137 subgrid closures on spurious interior mixing close to bottom topography. 138 The OVERFLOW experiment is adapted from the non-rotating overflow configuration described in 139 :cite:`HAIDVOGEL1999` and further used by :cite:`ILICAK2012`. 140 Here we can assess the behaviour of the second-order tracer advection scheme FCT2 and 141 forth-order FCT4, z-coordinate and sigma coordinate (...). 142 | Here the animation of the OVERFLOW test case in sigma coordinate with 143 the forth-order advection scheme FCT4. 103 144 104 145 WAD 105 === 106 107 A set of simple closed basin geometries for testing the Wetting and drying capabilities. 108 Examples range from a closed channel with EW linear bottom slope to a parabolic EW channel with a Gaussian ridge. 109 110 Below the animation of the test case 7. This test case is a simple linear slope with a mid-depth shelf with an open boundary forced with a sinusoidally varying ssh. 111 This test case has been introduced to emulate a typical coastal application with a tidally forced open boundary with an adverse SSH gradient that, when released, creates a surge up the slope. 112 The parameters are chosen such that the surge rises above sea-level before falling back and oscillating towards an equilibrium position 113 114 .. image:: _static/wad_testcase_7.gif 146 --- 147 148 .. figure:: _static/wad_testcase_7.gif 149 :width: 200px 150 :align: right 151 152 .. 153 154 | A set of simple closed basin geometries for testing the Wetting and drying capabilities. 155 Examples range from a closed channel with EW linear bottom slope to 156 a parabolic EW channel with a Gaussian ridge. 157 | Here the animation of the test case 7. 158 This test case is a simple linear slope with a mid-depth shelf with 159 an open boundary forced with a sinusoidally varying ssh. 160 This test case has been introduced to emulate a typical coastal application with 161 a tidally forced open boundary with an adverse SSH gradient that, 162 when released, creates a surge up the slope. 163 The parameters are chosen such that 164 the surge rises above sea-level before falling back and oscillating towards an equilibrium position. 115 165 116 166 CANAL 117 ===== 118 119 East-west periodic canal of variable size with several initial states and associated geostrophic currents (zonal jets or vortex). 120 121 .. image::_static/CANAL_image.gif 167 ----- 168 169 .. figure:: _static/CANAL_image.gif 170 :width: 200px 171 :align: left 172 173 .. 174 175 East-west periodic canal of variable size with several initial states and 176 associated geostrophic currents (zonal jets or vortex). 122 177 123 178 ICE_ADV2D 124 ========= 125 126 This test case illustrates the advection of an ice patch across an East/West and North/South periodic channel 127 over a slab ocean (i.e. one ocean layer).179 --------- 180 181 | This test case illustrates the advection of an ice patch across 182 an East/West and North/South periodic channel over a slab ocean (i.e. one ocean layer). 128 183 The configuration is similar to ICE_AGRIF, except for the AGRIF zoom. 129 184 | The purpose of this configuration is to test the advection schemes available in the sea-ice code 130 185 (for now, Prather and Ultimate-Macho from 1st to 5th order), 131 186 especially the occurence of overshoots in ice thickness 132 133 187 134 188 ICE_ADV1D 135 ========= 136 137 This experiment is the classical Schar & Smolarkiewicz (1996) test case :cite:`SCHAR1996`,138 which has been used in :cite:`LIPSCOMB2004`, 139 and in which very specific shapes of ice concentration,thickness and volume converge toward the center of a basin.189 --------- 190 191 | This experiment is the classical :cite:`SCHAR1996` test case , 192 which has been used in :cite:`LIPSCOMB2004`, and in which very specific shapes of ice concentration, 193 thickness and volume converge toward the center of a basin. 140 194 Convergence is unidirectional (in x) while fields are homogeneous in y. 141 The purpose of this configuration is to test the caracteristics of advection schemes available in the sea-ice code 195 | The purpose of this configuration is to 196 test the caracteristics of advection schemes available in the sea-ice code 142 197 (for now, Prather and Ultimate-Macho from 1st to 5th order), 143 especially the constitency between concentration, thickness and volume, and the preservation of initial shapes.144 145 References 146 ========== 147 148 .. bibliography:: test _cases.bib149 150 151 198 especially the constitency between concentration, thickness and volume, 199 and the preservation of initial shapes. 200 201 .. rubric:: References 202 203 .. bibliography:: tests.bib 204 :all: 205 :style: unsrt 206 :labelprefix: T -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/demo_cfgs.txt
r10516 r12063 9 9 WAD OCE 10 10 BENCH OCE ICE TOP 11 STATION_ASF OCE -
NEMO/branches/2019/dev_ASINTER-01-05_merged/tests/test_cases.bib
r10627 r12063 115 115 abstract = {Many problems of fluid dynamics involve the coupled transport of several, density-like, dependent variables (for instance, densities of mass and momenta in elastic flows). In this paper, a conservative and synchronous flux-corrected transport (FCT) formalism is developed which aims at a consistent transport of such variables. The technique differs from traditional FCT algorithms in two respects. First, the limiting of transportive fluxes of the primary variables (e.g., mass and momentum) does not derive from smooth estimates of the variables, but it derives from analytic constraints implied by the Lagrangian form of the governing continuity equations, which are imposed on the specific mixing ratios of the variables (e.g., velocity components). Second, the traditional FCT limiting based on sufficiency conditions is augmented by an iterative procedure which approaches the necessity requirements. This procedure can also be used in the framework of traditional FCT schemes, and a demonstration is provided that it can significantly reduce some of the pathological behaviors of FCT algorithms. Although the approach derived is applicable to the transport of arbitrary conserved quantities, it is particularly useful for the synchronous transport of mass and momenta in elastic flows, where it assures intrinsic stability of the algorithm regardless of the magnitude of the mass-density variable. This latter property becomes especially important in fluids with large density variations, or in models with a material “vertical” coordinate (e.g., geophysical hydrostatic stratified flows in isopycnic/isentropic coordinates), where material surfaces can collapse to zero-mass layers admitting, therefore, arbitrarily large local Courant numbers.} 116 116 } 117 118 @article{Brodeau_al_2017, 119 author={Laurent Brodeau and Bernard Barnier and Sergey Gulev and Cian Woods}, 120 title={Climatologically significant effects of some approximations in the bulk parameterizations of turbulent air-sea fluxes}, 121 journal={J. Phys. Oceanogr.}, 122 doi={10.1175/JPO-D-16-0169.1}, 123 year={2017}, 124 pages = {5-28}, 125 volume={47}, 126 number={1} 127 }
Note: See TracChangeset
for help on using the changeset viewer.