Changes between Version 2 and Version 3 of 2009WP/2009Stream3/VVL

Timestamp:

2009-04-20T16:48:32+02:00 (15 years ago)

Author:

rblod

Comment:

--

2009WP/2009Stream3/VVL

@@ -6 +6 @@
 [[PageOutline]]
 
+This page gives a non exhaustive  description of the work done to restore VVL option from NEMO version 3.1[[BR]]
+At this stage (april 2009)  :
+ - variable volume option can be used with either z, zps and s coordinate
+ - time splitting issue has not been really investigated
+
 ----
-== Getting the code ==
+== 1. Getting the code ==
 [[BR]]
 
 At this time (04-20-09), the work has been done on a branch starting from svn rev1359.[[BR]]
  - Outside modipsl :[[BR]]
-svn co --username nemo_user http://forge.ipsl.jussieu.fr/nemo/svn/branches/dev_004_VVL/NEMO[[BR]]
- - Inside modipsl :[[BR]]
+''svn co --username nemo_user !http://forge.ipsl.jussieu.fr/nemo/svn/branches/dev_004_VVL/NEMO
+  - Using modipsl :[[BR]]
 Edit modipsl/util/mod.def and replace [[BR]]
-tags/nemo_v3_1[[BR]]
+t''ags/nemo_v3_1''[[BR]]
 by [[BR]]
-branches/dev_004_VVL in NEMO section[[BR]]
-[[BR]]
+''branches/dev_004_VVL in NEMO section''[[BR]]
 
 
@@ -35 +39 @@
 
 
-== Code organisation ==
+== 2.  Code organisation ==
 [[BR]]
-
-Flow chart (original for Mathieu Leclair) :
+=== 2.1 Flow chart ===
+(original from Mathieu Leclair)
+[[Image(NEW_flowchart.png, 50%)]]
 
 
+=== 2.2 Sea surface height repartition  ===
 
-=== TITLE2.1  ===
+Repartition on the whole water column, so key_sigma_vll has been suppressed. Reference coordinate is referred as e3t_0 for instance, and pre-processing (in domzgr_substitute.h90) is used to define before, now, after scale factors, for instance :[[BR]]
+#   define  fsdept_n(i,j,k)  (fsdept_0(i,j,k)*(1+sshn(i,j)*mut(i,j,k)))[[BR]]
+
+Note: we shouldn't need 3d array mut (muu, muv) , but the code blows up without.[[BR]]
+
+
+== 3. Modified routines ==
+
+ - domvvl is now used only for initialisation (called from istate)
+ - dynspg_flt, dynspg_ts, dynspg_exp : no more update of ssh, just barotropic contribution for velocities
+ - trazdf_exp, trazdf_imp : ponderation by the correct scale factors
+ - wzvmod : suprresion of wzv_mod routine, addition of ssh_wzv to compute ssh after from the ssh equation, update the vertical coordinate and compute vertical velocity, addition of ssh_nxt for ssh time stepping
+ - dynnxt, tranxt : rewritting, take vvl case into account
+ - zdfevd : for stability issue, test to apply enhanced diffusion has to be performed on now Brunt-Vaisala frequency
+
+== 4. Cautions ==
+ - use s-coordinate pressure gradient when using variable volume[[BR]]
+ - the time splitting is minimalistic [[BR]]
+ - sea ice and vvl
+ - vertical velocity in output is shifted by 1 time step 
+
+== 5. To be done ==
+ - embedded sea-ice (Gurvan and Yevgeni)
+ - in vector invariant formulation, temporal derivatif of h is not treated (for both baroclinic and barotropic mode)
+ - a simplification for code and diagnostic could be to work on the tracer content everywhere
+ - for global ocean, resolve the damping force using time splitting routine
+ - use a predictor-corrector shceme for time splitting
+ - suppress mu* arrays
+ - we should used correct before, now, after everywhere in the code
+ - develop a pressure gradient algorithm for zps-vvl 
+
+== 6. Tests ==
 [[BR]]
-
-=== TITLE2.2  ===
-[[BR]]
-
-== Modified routines ==
-== Cautions ==
-pressure gradient[[BR]]
-time splitting[[BR]]
-== To be done ==
-
-== Tests ==
-
+=== 6.1 Paris ===
+ - ORCA2_LIM zps fixed volume and filtred free surface
+ - ORCA2_LIM zps variable volume and filtred free surface
+ - ORCA2_LIM zps fixed volume and time splitting free surface
+ - ORCA2_LIM zps variable volume and time splitting free surface
+ - restartability
+ - Agrif test case
+ - parallelisation
+ - conservation (with Asselin filter to 0 and ice tracer fluxes to 0)
 ----