Changes between Version 3 and Version 4 of schwarz_nemo

Timestamp:

07/23/19 15:07:30 (4 years ago)

Author:

snguyen

Comment:

--

schwarz_nemo

@@ -11 +11 @@
 With the Schwarz algorithm, the total amount of time steps is divided in coupling windows (called Schwarz loops in the code), each coupling window is repeated the number of schwarz iterations given as a parameter (mswr), and for each schwarz iteration a number of time steps (the length of a coupling window) is done.
 
-For a computation of 5 days with 15 time steps per day, a coupling window of 1 day and 6 schwarz iterations (mswr=6) there will be 3 loops.
-
-* One global loop over the coupling windows : 5 to iterate. Counter `iswloop`, limits `1 : nsloops`.
-* One inner loop over the number of schwarz iteration for each coupling window : 6. Counter `kswr`, limits `1 : mswr`.
-* One most inner loop over the time steps done in one Schwarz coupling window : 15 time steps. Counter `istp`, limits `nit000 + (iswloop-1)*ntsinswr : nit000+iswloop*ntsinswr-1`
+For a computation of 5 days with 15 time steps per day, a coupling window of 1 day and 6 Schwarz iterations there will be 3 loops.
+
+* One global loop over the coupling windows : 5 to iterate. Counter `iswloop`, limits `1 to nsloops`.
+* One inner loop over the number of schwarz iteration for each coupling window : 6. Counter `kswr`, limits `1 to mswr`.
+* One most inner loop over the time steps done in one Schwarz coupling window : 15 time steps. Counter `istp`, limits `nit000 + (iswloop-1)*ntsinswr to nit000+iswloop*ntsinswr-1`
 
 The counter given to subroutine `stp(istp)` is rewind to its value before entering the Schwarz loop when the Schwarz coupling window is repeated.
@@ -112 +112 @@
 === The Schwarz algorithm: storing and restoring ===
 
-The Schwarz algorithm is implemented by storing and restoring the state of the ocean and ice at the proper times.
+The Schwarz algorithm is implemented by storing and restoring the state of the ocean and ice at the proper position inside the three loops structure.
 
 The routines used to store and restore a state are all in `schwarz.F90` which is a new module.
@@ -125 +125 @@
 IF ( kswr > 1 ) CALL swz_reinit(istp) ! reset Ocean state to first time step for new schwarz iteration
 }}}
+
+The routine `swz_reinit` is reproducing the initialisation structure from `nemo_init` in `nemogcm.F90` by selecting only the relevant subroutines and coding the corresponding restoring of variables. 
+
+(In case you wish to use PISCES you will have to add the relevant routines in `swz_reinit` and module `schwarz.F90`)
 
 `Unfortunately this is not the only place variables are stored and restored!`