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

Timestamp:

2009-06-14T09:02:37+02:00 (15 years ago)

Author:

gm

Comment:

--

2009WP/2009Stream3/LeapFrogC

@@ -5 +5 @@
 [[PageOutline]]
 
-'''Object''': 
-   Introduce in the trunk the modified Leap-Frog + Robert-Asselin filter time stepping scheme (thereafter LF-RA) described in (Leclair and Madec Ocean Modelling 2009, in press). The scheme allows a exact conservation of the heat and salt contents when using variable volume thickness (key_vvl). The non-conservation of the standard LF-RA due to the time diffusion of the forcing through the filter. In the modified scheme the exact conservation is achieved by removing the forcing from the filter. This requires 3 modifications, on the forcing evaluation, on its time discretisation, and on the filter, resp. : [[BR]] 
+'''Object''': [[BR]]
+Introduce in the trunk the modified Leap-Frog + Robert-Asselin filter time stepping scheme (thereafter LF-RA) described in (Leclair and Madec Ocean Modelling 2009, in press). The scheme allows a exact conservation of the heat and salt contents when using variable volume thickness (key_vvl). The non-conservation of the standard LF-RA due to the time diffusion of the forcing through the filter. In the modified scheme the exact conservation is achieved by removing the forcing from the filter. This requires 3 modifications regarding the forcing evaluation, its time discretisation, and the filter : [[BR]] 
  - Forcing evaluation : the momentum and tracers forcing are evaluated at kt+1/2, no-more at kt. [[BR]]
  - Forcing time discretisation : 0.5 * ( Q(k+1/2) + Q(kt-1/2) ), instead of simply Q(kt). [[BR]]
@@ -15 +15 @@
 
 '''Module involved''': 
-    fldread.F90   
-    sbcmod.F90   ;   trasbc.F90 , traqsr.F90 , dynzdf_exp.F90 , dynzdf_imp.F90   
-    sshnxt.F90 , tranxt.F90 , dynnxt.F90
+    fldread.F90   [[BR]]
+    sbcmod.F90   ;   trasbc.F90 , traqsr.F90 , dynzdf_exp.F90 , dynzdf_imp.F90   [[BR]]
+    sshnxt.F90 , tranxt.F90 , dynnxt.F90[[BR]]
 
 ----
@@ -28 +28 @@
 === 2.  Forcing time discretisation ===
 
-[[Image(NEW_LF_RA.jpg, 50%)]]
+[[Image(NEW_LF_RA.jpg, 50%)]][[BR]]
+Illustration of forcing integration methods. Old formulation (top) and new formulation (bottom).
 
 Blah blah blah...