Changes between Version 12 and Version 13 of ticket/0863

Timestamp:

2011-10-13T17:40:44+02:00 (13 years ago)

Author:

mlelod

Comment:

--

ticket/0863

@@ -32 +32 @@
 }}}
 == Core description ==
+=== Newly created and suppressed variables ===
 === Step 1 : Initialisation ===
 `dom_vvl_init` is called by `domain.F90:`
@@ -48 +49 @@
 
 === Step 2 : Next sea surface height ===
+Next sea surface height is computed by `ssh_nxt` routine in module `sshwzv`.
+
 === Step 3 : Next vertical scale factor ===
+Next vertical scale factors are computed by `dom_vvl_nxt` routine in `domvvl` module. 2 cases are distinguished:
+
+ * '''z_star coordinate''' [[BR]] Vertical scale factors are not obtained with the repartition of the sea surface heigt over the vertical levels anymore. The reason is that the conservation correction applied to the leap-fog time stepping scheme for tracers and volume only concerns the first ocean level. We thus have to treat the entire 3D vertical scale factors: [[BR]]
+{{{
+z_scale(:,:) = ( ssha(:,:) - sshb(:,:) ) * tmask(:,:,1) / ( ht_0(:,:) + sshn(:,:) + 1. - tmask(:,:,1) )
+DO jk = 1, jpkm1
+   fse3t_a(:,:,jk) = fse3t_b(:,:,jk) + fse3t_n(:,:,jk) * z_scale(:,:)
+END DO
+}}}
+ * '''z_tilde or layer coordinate''' [[BR]] 1. Baroclinic Vertical scale factor anomalies (`e3t_t_a`) are obtained thanks to 3 tendency terms:[[BR]] * High frequency baroclinic horizontal divergence (which is the entire divergence `hdivn` in the layer case).[[BR]] * Restoring towards 0 (only for z_tilde).[[BR]] * Thickness diffusion term. [[BR]] 2. Barotropic part of the scale factor is incremented by the repartition of the sea surface height variation (as in the z_star case)
+
 === Step 4 : Now vertical velocity (Eulerian cross-level velocity) ===
+Now Eulerian vertical velocity (`wn`) is computed by `wzv` routine in `sshwzv` module. There is no particular reason for letting this routine in this module but no real reason to put it elsewhere ... The computation of wn takes the thickness diffusion transports (`un_td` and `vn_td`) into account in order to avoid a compensating vertical velocity.
+
 === Step 5 : Thickness diffusion velocity in tracers advection ===
+In order to be consistent between volume and tracer equations the thickness diffusion transport is also taken into account in the lateral tracer advection terms.
+
 === Step 6 : Swap of sea surface height ===
+This part was previously done in `ssh_nxt` routine which could cause a bit of confusion here. the routine is now called ssh_swp, still in the `sshwzv` module.
+
 === Step 7 : Swap of vertical scale factors ===
+Swap of vertical scale factors is done in two different routines.
+
+ * The first part concerns `fse3t_b` and fse3t_nand is done in `dynnxt` routine because they are needed there.
+ * The second part concerns `e3t_t_b` and `e3t_t_n`. This is done by `dom_vvl_swp` routine in `domvvl` module.
+
+`dom_vvl_swp` routine also interpolates scale factors from t-, u- or v-points to other grid points and computes depths and water column heights.
+
 == Side modifications ==
+=== Time splitting external mode ===
 === Coupling between ocean and sea ice ===
 === Tracer lateral diffusion ===