Changes between Version 4 and Version 5 of ticket/0663_MLF


Ignore:
Timestamp:
2010-04-22T19:42:48+02:00 (11 years ago)
Author:
gm
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0663_MLF

    v4 v5  
    2424This update is part of the LOCEAN.8 task (Update of the ocean physics) due by June 25th 2010. 
    2525 
    26 The main remaining modifications concern :  
     26The 3 main modifications concern :  
    2727 
    28 (1) forcing terms as a mean between time step n-1/2 and n+1/2 (sbcmod.F90, sbc.F90, trasbc.F90, dynzdf....F90) 
     28'''(1) The forcing terms'''. They are defined as a mean between time step n-1/2 and n+1/2 (sbcmod.F90, sbc.F90, trasbc.F90, dynzdf....F90) 
    2929  
    3030[[Image(NEW_LF_RA.jpg, 50%)]][[BR]] 
     
    3838      Beginning of sbs routine, introduce the swap, but only for kt >= nit000+1. If no restart, the '_b' fields is known only at the end of the sbc calculation. Therefore we have choosen to regroup the setting of '_b' fields at nit000 at the end of sbc routine. CAUTION : here the swap must be done at each kt, what ever the nn_fsbc value is. Otherwise the ocean forcing will be wrong. 
    3939 
    40 (2) Asselin filter: remove the forcing from the filter (dynnxt (?), tranxt, traqsr) 
     40'''(2) The Asselin filter''' remove the forcing from the filter (dynnxt (?), tranxt, traqsr) 
    4141 
    42 (3) tranxt.F90 : semi-implicit hydrostatic pressure gradient : change the coefficients to take into account the Asselin filter parameter 
     42'''(3) The semi-implicit hpg'''. In tranxt.F90, for the semi-implicit hydrostatic pressure gradient,  the coefficients are changed in order to take into account the Asselin filter parameter 
    4343 
    4444   time mean T and S for the semi implicit hpg computation was (ta+2*tn + tb) / 4 . It is now replaced by 
     
    5050   The change extend further the stability limit of the time stepping scheme when ln_dynhpg_imp = T 
    5151 
     52'''N.B.''' : Three key elements of the LF-RA scheme have already been implemented in former version of NEMO (v3.0 and v3.1). As a remainder, there are: 
    5253 
    53 === Open questions === 
     54   Forcing evaluation : the momentum and tracers forcing are evaluated at kt+1/2, no-more at kt. This is done in the fldread routine since v3.0 
     55   Flow chart : general structure has changed for key_vvl (variable volume or equivalently non-linear free surface) (v3.1). In particular the ssh equation is solved first together followed immediately by the computation of the now vertical velocity, so that wn is available to solve the tracer and the dynamics equations.  
    5456 
    55 (2) Do we need to remove the stress component from the filter in dynnxt?  
     57   The enhanced vertical diffusion (evd) is no more a source of small divergence of two consecutive time step as the static instability is now tested on both before and now fields. See ticket #401. 
    5658 
    57 (3) in vvl case : thickness weighted tracer in the time mean, do we need to compute and use a time mean thickness in the hpg computation? 
     59 
     60=== Pending issues === 
     61 
     62• Do we need to remove the stress component from the filter in dynnxt?  
     63 
     64• in vvl case : thickness weighted tracer in the time mean, do we need to compute and use a time mean thickness in the hpg computation? 
     65 
     66• The associated changes in TOP are still unclear. A priori, in PISCES p4zopt.F90 should be modified. Nevertheless, the TOP time stepping strategy seems to be a forward in time scheme. If yes, then the evaluation of the forcing at kt+1/2 is an improvement (if not to say a bug fix!). If not, action is required...  
     67 
     68 
     69• Check the consistency with the change planned in the emps field (switch to fsalt, a true salt flux).  
     70 
     71 
     72• modification required in limsbc, limsbc_2 and CICE interface : conversion of emps into a salt flux (new name 'fsalt' ?)  
     73 
     74 
     75• the modifications in trasbc will interfere with the new runoff scheme (see 2009WP/2009Stream3#Substream3.2:fromNEMOteamshort-term) 
    5876 
    5977----