Changes between Version 7 and Version 8 of 2014WP/2014_SystemSimplification_v2
- Timestamp:
- 2014-03-26T10:04:46+01:00 (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
2014WP/2014_SystemSimplification_v2
v7 v8 3 3 '''Motivation''': 4 4 5 We need to simplify the code for maintenance, robustness and future evolution (change in time stepping, …etc). This means fewer options, and a consensus on the default (and thus well validated) set of numerics. The idea is to change the current most used choices from : [[BR]]5 We need to simplify the code for maintenance, robustness and future evolution of the dynamical core. This means fewer options, and a consensus on the default (and thus well validated) set of numerics. The idea is to change the current most used choices from : [[BR]] 6 6 7 linear filtered free surface zps ocean with TVD, "Madec" iso-neutral diffusion, and standard equation of state; levitating sea-ice with virtual salt fluxes7 linear filtered free surface zps ocean with TVD, and standard iso-neutral diffusion ; levitating sea-ice with virtual salt fluxes 8 8 9 9 possibly to :[[BR]] 10 10 11 Non linear (vvl), split-explicit (ts) free surface, mixte zps/s-coordinate ocean with 3rd order advection (UBS with rotation of the diffusive part along iso-neutral surfaces), "Griffies" iso-neutral diffusion; embedded multi-category sea-ice (exact fluxes).11 Non linear (vvl), split-explicit (ts) free surface, mixte zps/s-coordinate ocean with high order advection schemes, triad iso-neutral diffusion (Griffies operator) ; embedded multi-category sea-ice (exact fluxes). 12 12 13 13 And this with a removing of obsolescent features and a cleaning of the code. 14 14 15 Timing: achievement within 2 years from now for NEMO v4.0 (to be ready for CMIP6, and on time for the foreseen evolution of the time stepping)15 Timing: achievement within 2 years from now for NEMO v4.0 (to be ready for CMIP6, and on time for the foreseen evolution of the dynamical core) 16 16 17 17 Notations: … … 33 33 '''*''' '''vertical coordinate''' : default option mixte zps/s-tilde coordinate (thus with variable volume level), with mixte zps*/s*, mixte z/s, zps, z and s coordinate system as sub-option (same code used) (2014-2015 ?) 34 34 35 !- tracer solved as thickness weighted tracer: D,, t,, (e,, 3t,, T) ('''WP 2014''' 35 !- tracer solved as thickness weighted tracer: D,, t,, (e,, 3t,, T) ('''WP 2014''') 36 36 37 !- thickness weighted quantities in the outputs (''' WP2014''')37 !- thickness weighted quantities in the outputs ('''WP2014''') 38 38 39 !- no more substitute (domzgr_substitute.h90): all e3 always defined at before, now and after time steps (''' partly done''' in 2013, '''to be achieved''' in 2014).39 !- no more substitute (domzgr_substitute.h90): all e3 always defined at before, now and after time steps ('''partly done''' in 2013, '''to be achieved''' in 2014). 40 40 41 41 !- adapte OFF line for tracer to variable volume (2015?) … … 57 57 !- 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. To be discussed 58 58 59 !- Addition of the rotation of the diffusive part of 3^ rd^ order scheme (2015 ?)59 !- Addition of the rotation of the diffusive part of 3^rd^ order scheme (2014-2015) 60 60 61 61 '''-''' '''vertical physics''': simplification versus diversity … … 65 65 !- 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. 66 66 67 !- suppress or keep the explicit vertical diffusion 67 !- suppress or keep the explicit vertical diffusion? 68 68 69 69 [[BR]] [[BR]]'''II.OBSOLESCENT features''' [[BR]] … … 71 71 '''mains''': 72 72 73 ''' *''' remove ''' 73 ''' *''' remove '''LIM2''' : only LIM3 remains in the system ('''WP 2014''') 74 74 75 ''' *''' remove ''' 75 ''' *''' remove '''OBC''' : only BDY remains in the system ('''==>> done in v3.6''') 76 76 77 ''' *''' remove ''' OASIS''' 3 and 4 : OASIS-MCT('''WP 2014''')77 ''' *''' remove '''OASIS''' 3 and 4 : OASIS-MCT instead ('''WP 2014''') 78 78 79 ''' -''' duplicate ''' isoneutral operator''' : keep only the Griffies' one (tests needed before decision, tests planed in 2014 NOCS+LOCEAN) 80 79 ''' -''' duplicate ''' isoneutral operator''' : keep only the triad operator (Griffies' one) ? Tests needed before decision. Those tests are planed in 2014 NOCS+LOCEAN) 81 80 82 81 '''smaller ones''' … … 86 85 ''' *''' remove '''Neptune effect''' : ORCA1 runs shows its inefficiency (''' WP 2014''' ). 87 86 88 ''' *''' remove '''no-slip accurate''' option (key_noslip_accurate). Module divcur.F90 (==>> done in branch [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2014 2014/dev_CNRS0_NOC1_LDF] 87 ''' *''' remove '''no-slip accurate''' option (key_noslip_accurate). Module divcur.F90 (==>> done in branch [https://forge.ipsl.jussieu.fr/nemo/browser/branches/2014 2014/dev_CNRS0_NOC1_LDF] ) 89 88 90 89 ''' *''' suppress '''cross-land advection''' : this option is hard coded for ORCA_R2_L31 and is nomore used. … … 95 94 96 95 ''' -''' '''Non-penetrative convection''' : with its correction and improvement using alpha & beta (work done by L. Brodeau, Stockholm), it seems now vbetter to keep it in the system. 96 97 ''' -''' '''runoffs''' : only introduced as lateral flow: remove the option of runoff as surface boundary condition as well as specific treatment in vicinity of river mouths. 98 97 99 98 100 … … 107 109 '''CODE cleaning''' 108 110 109 ''' *''' 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. This is also associated with the re-writing of bilaplacian operator on tracer and on dynamics as re-entrant laplacian operator. As a results, curl and div are no more required in the restart file. On going work, see #1260 and the associated wiki page [https://forge.ipsl.jussieu.fr/nemo/wiki/ticket/1260_CNRS0_NOC1_LDF 1260_CNRS0_NOC1_LDF]111 ''' *''' LDF: remove 1D and 2D cases, keep only 3D ; laplacian and bilaplacian coefficients defined at the same points (T and F) ; time varying coefficients (function of growth rate of baroclinic instabilities, of the local velocity...): better integration rather than configuration dependent cases. This is also associated with the re-writing of bilaplacian operator on tracer and on dynamics as re-entrant laplacian operator. As a results, curl and div are no more required in the restart file. On going work, see #1260 and the associated wiki page [https://forge.ipsl.jussieu.fr/nemo/wiki/ticket/1260_CNRS0_NOC1_LDF 1260_CNRS0_NOC1_LDF] 110 112 111 ''' *''' OBS operator : clean the stand-alone implementation ('''WP 2014''' 113 ''' *''' OBS operator : clean the stand-alone implementation ('''WP 2014''') 112 114 113 115 ''' *''' avm defined at t-point ==> suppress avmu, avmv and ~10 lbc_lnk calls (2015 ?)