Changes between Version 1 and Version 2 of 2014WP/2014_SystemSimplification
- Timestamp:
- 2013-11-11T16:40:38+01:00 (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
2014WP/2014_SystemSimplification
v1 v2 1 = = '''Note on possible system simplification''' ==1 = '''Note on possible system simplification''' = 2 2 We need to simplify the code for both maintenance and future evolution. This means fewer options, and a consensus on the default (and thus well validated) set of numerics. 3 3 4 The idea is to change the used choices from : 4 The idea is to change the used choices from : 5 5 6 linear filtered free surface ocean with TVD, "Madec" isoneutral diffusion, and standard equation of state ; 6 linear filtered free surface ocean with TVD, "Madec" isoneutral diffusion, and standard equation of state ; levitating sea-ice with virtual salt fluxes 7 7 8 8 to : 9 9 10 Non linear (vvl) split-explicit (ts) free surface with 3rd order advection (UBS), 10 Non linear (vvl) split-explicit (ts) free surface with 3rd order advection (UBS), "Griffies" isoneutral diffusion and TEOS10 ; embedded multi-category sea-ice (exact fluxes). 11 11 12 '''Free surface & vertical coordinate :''' 12 == '''Free surface & vertical coordinate :''' == 13 • remove filtered free surface and all associated solvers 13 14 14 • remove filtered free surface and all associated solvers 15 • default : non-linear free surface 15 16 16 • default : non-linear free surface 17 • linear option set without optimization (i.e. starting from the non-linear case with minimal changes, not optimal changes), or linear case suppressed ?? 17 18 18 • linear option set without optimization (i.e. starting from the non-linear case with minimal changes, not optimal changes), or linear case suppressed ?? 19 • default : variable volume level (vvl), either z*/s*, or z/s-tilde 19 20 20 • default : variable volume level (vvl), either z*/s*, or z/s-tilde 21 Needs : 21 22 22 Needs : 23 * split-explicit free surface adapted to AGRIF refinements and TAM 23 24 24 - split-explicit free surface adapted to AGRIF refinements and TAM 25 * tracer solved systematically as thickness weighted tracer: D,,t,,(e,,3t,, T) 25 26 26 - tracer solved systematically as thickness weighted tracer: D,,t,,(e,,3t,, T) 27 * no more substitute (domzgr_substitute.h90): all e3 always defined at before, now and after time steps. 27 28 28 - no more substitute (domzgr_substitute.h90): all e3 always defined at before, now and after time steps. 29 * sea-ice : embedded sea-ice only, so management of land-fast ice and especially sea-ice approaching the ocean-bottom.[[BR]] 29 30 30 - sea-ice : embedded sea-ice only, so management of land-fast ice and especially sea-ice approaching the ocean-bottom. 31 == '''Equation of State (eos) :''' == 32 • EOS10 as default + 1 simplified eos. Suppression of the now used EOS. 31 33 32 '''Equation of State (eos) :''' 34 • EOS provides systematically alpha and beta (used in ldf_iso, bbl, mld, npc…) 33 35 34 • EOS10 as default + 1 simplified eos. Suppression of the now used EOS. 36 • no more potential density (mixed layer depth computed from alpha & beta) 35 37 36 • EOS provides systematically alpha and beta (used in ldf_iso, bbl, mld, npc…) 38 • efficient transfer function from conservative temperature to potential temperature and vice versa. [[BR]] 37 39 38 • no more potential density (mixed layer depth computed from alpha & beta) 40 == '''Coupled interface to the atmosphere :''' == 41 • OASIS-MCT only: remove OASIS 3 interface. 39 42 40 • efficient transfer function from conservative temperature to potential temperature and vice versa. 43 • Multi-category sea-ice interface [[BR]] 41 44 42 '''Coupled interface to the atmosphere :''' 45 == '''Physics :''' == 46 • avm defined at t-point ==> suppress avmu, avmv and ~10 lbc_lnk calls 43 47 44 • OASIS-MCT only: remove OASIS 3 interface. 48 • LDF: remove 1D and 2D cases, keep only 3D ; laplacian and bilaplacian coefficients defined at the same points (T and F) ; time varying coefficients (Smagorinsky, function of growth rate of baroclinic instabilities): better integration rather than configuration dependent cases 45 49 46 • Multi-category sea-ice interface 50 • TKE remove (use GLS instead?) 47 51 48 '''Physics :''' 52 • KPP remove unless if-less KPP is implemented (Lemarié et al 2012). Keeping a validated KPP scheme in the system make sense for comparison purposes. 49 53 50 • avm defined at t-point ==> suppress avmu, avmv and ~10 lbc_lnk calls 54 • Non-penetrative convection : either correct it (using alpha & beta) or remove it[[BR]] 51 55 52 • LDF: remove 1D and 2D cases, keep only 3D ; laplacian and bilaplacian coefficients defined at the same points (T and F) ; time varying coefficients (Smagorinsky, function of growth rate of baroclinic instabilities): better integration rather than configuration dependent cases 56 == '''Configurations :''' == 57 • remove all hard coded specific alterations associated with configurations (ORCA,…) while introducing associated input files 53 58 54 • TKE remove (use GLS instead?) 59 • same for the definition of coordinate system (s, mixed z-s): only a simple case hard coded. For more sophisticated cases, the required information should be read in an input file.[[BR]] 55 60 56 • KPP remove unless if-less KPP is implemented (Lemarié et al 2012). Keeping a validated KPP scheme in the system make sense for comparison purposes. 61 == '''Obsolescent features : '''to be suppressed == 62 • remove LIM2 : only LIM3 remains in the system 57 63 58 • Non-penetrative convection : either correct it (using alpha & beta) or remove it 64 • remove OBC : only BDY remains in the system 59 65 60 '''Configurations :''' 66 • acceleration of convergence (in particular, rdttra function of k ==> unique rdt) 61 67 62 • remove all hard coded specific alterations associated with configurations (ORCA,…) while introducing associated input files 68 • remove Neptune effect : ORCA1 run shows its inefficiency. 63 69 64 • same for the definition of coordinate system (s, mixed z-s): only a simple case hard coded. For more sophisticated cases, the required information should be read in an input file. 70 • duplicate isoneutral operator: keep only the Griffies' one 65 71 66 '''Obsolescent features : '''to be suppressed 72 • Smagorinsky (?) : use 3rd order operator (UBS) instead and remove it ? 67 73 68 • remove LIM2 : only LIM3 remains in the system 74 • tracer advection schemes: diminish their number. Typically one module for 2^ nd^ /4^ th^ order and one for a 3^ rd^ order scheme, both associated with an optional flux limiter to ensure positivity.[[BR]] 69 75 70 • remove OBC : only BDY remains in the system 76 == '''Outputs :''' == 77 • remove Dimg format ? (agreement needed from DRAKKAR group) 71 78 72 • acceleration of convergence (in particular, rdttra function of k ==> unique rdt) 73 74 • remove Neptune effect : ORCA1 run shows its inefficiency. 75 76 • duplicate isoneutral operator: keep only the Griffies' one 77 78 • Smagorinsky (?) : use 3rd order operator (UBS) instead and remove it ? 79 80 • tracer advection schemes: diminish their number. Typically one module for 2^nd^/4^th^ order and one for a 3^rd^ order scheme, both associated with an optional flux limiter to ensure positivity. 81 82 '''Outputs :''' 83 84 • remove Dimg format ? (agreement needed from DRAKKAR group) 85 86 • remove IOIPSL: the key issue here is do we believe XIOs is able to work on all existing and foreseen platform? 79 • remove IOIPSL: the key issue here is do we believe XIOs is able to work on all existing and foreseen platform?
