Changes between Version 2 and Version 3 of 2020WP/KERNEL-06_techene_better_e3_management

Timestamp:

2020-01-22T10:27:31+01:00 (5 years ago)

Author:

techene

Comment:

--

2020WP/KERNEL-06_techene_better_e3_management

@@ -11 +11 @@
 == Summary
 
-||=Action       || improve e3* management                        ||
+||=Action       || optimisation of the vertical scale factor e3 computation          ||
 ||=PI(S)        || Techene, Madec                                                 ||
-||=Digest       || get rid of the very large number of e3* arrays and improve some of the e3* computation  ||
+||=Digest       || compute e3 on the fly from e3_0(:,:,:,Klt) * ( 1 + ssh(:,:,Klt) / h_0(:,:) * mask(:,:,:) instead of storing e3t/u/v/w/f... ||
 ||=Dependencies || If any                                                ||
 ||=Branch       || source:/NEMO/branches/2020/dev_r{REV}_{ACTION_NAME} ||
@@ -23 +23 @@
 
 
-use e3_0 and eta to recompute all e3* when needed.
+The current e3[P] at P-point computation uses interpolation of the r3t 4D table at P = {u-, v-, w-, f-, uw-, vw-} points. This means 7 4D tables stored in memory. 
+The proposed optimisation consists in computing e3[P](ji,jj,jk,Ktl) on the fly using the r3[P] = ssh[P] / h_0 and the e3[P]_0. r3[P] is a 2D table, then this means only 4 2D tables stored in memory.
+z-tilde management is done through e3[P]_0 that may varies with time in the z-tilde case.
+Asselin filter management is done recomputing r3[P] directly with the filtered ssh.
 
 
@@ -34 +37 @@
 give for each the chosen name and description wrt coding rules.
 }}}
+
+1/ REFERENCE version change the way to compute the e3 using ssh instead of interpolating \\
+For this we need to modify : \\
+- cfgs/GYRE_PISCES/cpp_GYRE_PISCES.fcm (top de-activation) \\
+- src/OCE/DOM/dom_oce.F90 (declaration and allocation of new variable) \\
+- src/OCE/DOM/domain.F90 (initialisation of h new variables) \\
+- src/OCE/DOM/dommsk.F90 (initialisation of mask new variable) \\
+- src/OCE/DOM/domvvl.F90 (dom_vvl_interpol and dom_vvl interpol_3d) \\
+- src/OCE/DYN/dynatf.F90 (change but no effect since ln_dynadv_vec = .true.)
+- src/OCE/TRA/traatf.F90 (change the e3t update with Asselin filter)
+- src/OCE/step.F90 (change the ordering ssh >> tra >> dyn) \\
+
+
+2/ 
+
 
 ''...''