Changeset 3142


Ignore:
Timestamp:
2011-11-17T15:54:45+01:00 (10 years ago)
Author:
acc
Message:

Branch dev_NEMO_MERGE_2011. Minor documentation changes and check in dynhpg.F90 to prevent use of buggy horizontal pressure gradient option

Location:
branches/2011/dev_NEMO_MERGE_2011
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_NEMO_MERGE_2011/DOC/TexFiles/Chapters/Chap_DYN.tex

    r3116 r3142  
    607607Pressure gradient formulations in an $s$-coordinate have been the subject of a vast  
    608608number of papers ($e.g.$, \citet{Song1998, Shchepetkin_McWilliams_OM05}).  
    609 A number of different pressure gradient options are coded, but they are not yet fully  
    610 documented or tested.  
     609A number of different pressure gradient options are coded but the ROMS-like, density Jacobian with  
     610cubic polynomial method is currently disabled whilst known bugs are under investigation. 
    611611 
    612612$\bullet$ Traditional coding (see for example \citet{Madec_al_JPO96}: (\np{ln\_dynhpg\_sco}=true) 
     
    625625($e_{3w}$).  
    626626 
     627$\bullet$ Pressure Jacobian scheme (prj) (a research paper in preparation) (\np{ln\_dynhpg\_prj}=true) 
     628 
    627629$\bullet$ Density Jacobian with cubic polynomial scheme (DJC) \citep{Shchepetkin_McWilliams_OM05}  
    628 (\np{ln\_dynhpg\_djc}=true) 
    629  
    630 $\bullet$ Pressure Jacobian scheme (prj) \citep{Thiem_Berntsen_OM06} (\np{ln\_dynhpg\_prj}=true) 
    631  
    632 Note that expression \eqref{Eq_dynhpg_sco} is commonly used when the variable volume  
    633 formulation is activated (\key{vvl}) because in that case, even with a flat bottom,  
    634 the coordinate surfaces are not horizontal but follow the free surface  
    635 \citep{Levier2007}. Only the pressure jacobian scheme (\np{ln\_dynhpg\_prj}=true) is available as an  
    636 alternative to the default \np{ln\_dynhpg\_sco}=true when \key{vvl} is active.  The pressure Jacobian scheme uses  
    637 a constrained cubic spline to reconstruct the density profile across the water column. This method 
    638 maintains the monotonicity between the density nodes and is of a higher order than the linear 
    639 interpolation method. The pressure can be calculated by analytical integration of the density profile and 
    640 a pressure Jacobian method is used to solve the horizontal pressure gradient. This method should 
    641 provide a more accurate calculation of the horizontal pressure gradient than the standard scheme. 
     630(\np{ln\_dynhpg\_djc}=true) (currently disabled; under development) 
     631 
     632Note that expression \eqref{Eq_dynhpg_sco} is commonly used when the variable volume formulation is 
     633activated (\key{vvl}) because in that case, even with a flat bottom, the coordinate surfaces are not 
     634horizontal but follow the free surface \citep{Levier2007}. The pressure jacobian scheme 
     635(\np{ln\_dynhpg\_prj}=true) is available as an improved option to \np{ln\_dynhpg\_sco}=true when 
     636\key{vvl} is active.  The pressure Jacobian scheme uses a constrained cubic spline to reconstruct 
     637the density profile across the water column. This method maintains the monotonicity between the 
     638density nodes  The pressure can be calculated by analytical integration of the density profile and a 
     639pressure Jacobian method is used to solve the horizontal pressure gradient. This method can provide 
     640a more accurate calculation of the horizontal pressure gradient than the standard scheme. 
    642641 
    643642%-------------------------------------------------------------------------------------------------------------- 
  • branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

    r3118 r3142  
    139139         WRITE(numout,*) '      time stepping: centered (F) or semi-implicit (T)  ln_dynhpg_imp = ', ln_dynhpg_imp 
    140140      ENDIF 
     141      ! 
     142      IF( ln_hpg_djc )   & 
     143         &   CALL ctl_stop('dyn_hpg_init : Density Jacobian: Cubic polynominal method & 
     144                           & currently disabled (bugs under investigation). Please select & 
     145                           & either  ln_hpg_sco or  ln_hpg_prj instead') 
    141146      ! 
    142147      IF( lk_vvl .AND. .NOT. (ln_hpg_sco.OR.ln_hpg_prj) )   & 
Note: See TracChangeset for help on using the changeset viewer.