Changeset 3142
- Timestamp:
- 2011-11-17T15:54:45+01:00 (12 years ago)
- 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 607 607 Pressure gradient formulations in an $s$-coordinate have been the subject of a vast 608 608 number of papers ($e.g.$, \citet{Song1998, Shchepetkin_McWilliams_OM05}). 609 A number of different pressure gradient options are coded , but they are not yet fully610 documented or tested. 609 A number of different pressure gradient options are coded but the ROMS-like, density Jacobian with 610 cubic polynomial method is currently disabled whilst known bugs are under investigation. 611 611 612 612 $\bullet$ Traditional coding (see for example \citet{Madec_al_JPO96}: (\np{ln\_dynhpg\_sco}=true) … … 625 625 ($e_{3w}$). 626 626 627 $\bullet$ Pressure Jacobian scheme (prj) (a research paper in preparation) (\np{ln\_dynhpg\_prj}=true) 628 627 629 $\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 632 Note that expression \eqref{Eq_dynhpg_sco} is commonly used when the variable volume formulation is 633 activated (\key{vvl}) because in that case, even with a flat bottom, the coordinate surfaces are not 634 horizontal 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 637 the density profile across the water column. This method maintains the monotonicity between the 638 density nodes The pressure can be calculated by analytical integration of the density profile and a 639 pressure Jacobian method is used to solve the horizontal pressure gradient. This method can provide 640 a more accurate calculation of the horizontal pressure gradient than the standard scheme. 642 641 643 642 %-------------------------------------------------------------------------------------------------------------- -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r3118 r3142 139 139 WRITE(numout,*) ' time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = ', ln_dynhpg_imp 140 140 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') 141 146 ! 142 147 IF( lk_vvl .AND. .NOT. (ln_hpg_sco.OR.ln_hpg_prj) ) &
Note: See TracChangeset
for help on using the changeset viewer.