New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
2021WP/KNL-01_Sibylle_RK3_stage1 (diff) – NEMO

Changes between Version 4 and Version 5 of 2021WP/KNL-01_Sibylle_RK3_stage1


Ignore:
Timestamp:
2021-03-01T18:18:30+01:00 (3 years ago)
Author:
techene
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • 2021WP/KNL-01_Sibylle_RK3_stage1

    v4 v5  
    4040give for each the chosen name and description wrt coding rules. 
    4141}}} 
    42 RK3 implementation is split up into :  
    43 - dynamic and active tracers  
    44 - vertical physics  
    45 - barotropic mode  
    46 - mass forcing 
     42RK3 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  
    4748- passive tracers 
    4849 
     50''' Code preparation ''' 
     51In order to preserve constancy property velocity for momentum and active tracers must be the same.  
     52Advection routines in flux form are modified to take (u,v,w) as an input argument. 
     53In order to use advection routines for the barotropic mode we need the possibility to de-activate vertical advection computation.  
     54Advection routines in flux and vector form are modified to take an optional argument (no_zad) to do so.  
    4955 
    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 ''' 
     58For 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). 
     59Forcing 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 
    5562  * Change eos divhor and sshwzv interface. 
    5663* Add RK3 time stepping routines 
     
    5865  * stprk3 orchestrates  
    5966 
    60 2) Barotrope part 
    61 TBD 
     67''' Barotrope part ''' 
     68In order to validate 2D mode implementation we remove above zero forcing for barotropic variables mass forcing remains to zero.  
    6269 
    63  
     70*  
    6471 
    6572''...''