Changes between Version 4 and Version 5 of 2021WP/KNL-01_Sibylle_RK3_stage1
- Timestamp:
- 2021-03-01T18:18:30+01:00 (3 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
2021WP/KNL-01_Sibylle_RK3_stage1
v4 v5 40 40 give for each the chosen name and description wrt coding rules. 41 41 }}} 42 RK3 implementation is split up into : 43 - dynamic and active tracers 44 - vertical physics 45 - barotropic mode 46 - mass forcing 42 RK3 implementation is splitted up into : 43 - code preparation 44 - dynamic and active tracers (barocline) 45 - vertical physics (TKE) 46 - barotropic mode (barotrope) 47 - mass forcing 47 48 - passive tracers 48 49 50 ''' Code preparation ''' 51 In order to preserve constancy property velocity for momentum and active tracers must be the same. 52 Advection routines in flux form are modified to take (u,v,w) as an input argument. 53 In order to use advection routines for the barotropic mode we need the possibility to de-activate vertical advection computation. 54 Advection routines in flux and vector form are modified to take an optional argument (no_zad) to do so. 49 55 50 1) Barocline part 51 * Set up a reference : GYRE only configuration with ssh, uu_b, un_adv are set to zero at each time step, no mass forcing, constant vertical viscosity. 52 * Prepare advection routines 53 * Add advective velocities as arguments for dynadv_cen2, dynadv_ubs (flux form) and traadv routines 54 * Allow vertical advection de-activation dynadv_cen2, dynadv_ubs (flux form) and dynzad routines 56 57 ''' Barocline part ''' 58 For sake of simplicity we started to implement RK3 regarding a GYRE configuration validation with no barotrope mode (ssh, uu_b, un_adv are set to zero at each time step). 59 Forcing have been removed except winds and heat flux. key_qco is active and vertical physics is modeled as constant with high viscosity coefficients. 60 61 * Prepare routines 55 62 * Change eos divhor and sshwzv interface. 56 63 * Add RK3 time stepping routines … … 58 65 * stprk3 orchestrates 59 66 60 2) Barotrope part 61 TBD 67 ''' Barotrope part ''' 68 In order to validate 2D mode implementation we remove above zero forcing for barotropic variables mass forcing remains to zero. 62 69 63 70 * 64 71 65 72 ''...''