Changes between Version 1 and Version 2 of schwarz_nemo

Timestamp:

07/22/19 17:32:01 (4 years ago)

Author:

snguyen

Comment:

--

schwarz_nemo

@@ -114 +114 @@
 The Schwarz algorithm is implemented by storing and restoring the state of the ocean and ice at the proper times.
 
-The routines used to store and restore a state are all in `schwarz.F90` which is a new subroutine.
-
-The storing is done just before the schwarz looping is done: inside `iswloop` loop, before `kswr` loop.
-Like this we are sure to store the initial condition of a schwarz loop.
-
+The routines used to store and restore a state are all in `schwarz.F90` which is a new module.
+
+The storing is done just before the schwarz looping is done: inside the `iswloop` loop, before the `kswr` loop.
+
+Like this we are certain to store the initial condition for a Schwarz loop.
+
+The restoring is done at the beginning of the `kswr` loop before the `istp` loop.
+If we are not in the first iteration the variables state is restored:
+{{{
+IF ( kswr > 1 ) CALL swz_reinit(istp) ! reset Ocean state to first time step for new schwarz iteration
+}}}
+
+`Unfortunately this is not the only place variables are stored and restored!`
+
+Because NEMO is modular some implementation details are repeated in the modules some modules have their own storing and restoring details.
+
+The files concerned are:
+{{{
+sbcmod.F90
+sbcrnf.F90
+traqsr.F90
+trasbc.F90
+}}}
+
+== Coding details: modification of other files ==
+
+The routines:
+{{{
+iom.F90
+step.F90
+}}}
+have been modified to output only one time serie for a set of Schwarz loops. The parameter `ksout` selects which Schwarz loop is given for outputs. 
+
+Some routines have a particular behavior when `istp=nit000` or/and `istp=nitend`. They where modified to behave correctly when doing Schwarz iterations:
+{{{
+closea.F90
+diaar5.F90
+diafwb.F90
+diahth.F90
+dynspg_ts.F90
+dynvor.F90
+fldread.F90
+limdyn.F90
+limitd_me.F90
+limrst.F90
+limthd.F90
+limtrp.F90
+limupdate1.F90
+limupdate2.F90
+restart.F90
+sbcblk_core.F90
+sbcice_lim.F90
+sbcssm.F90
+stpctl.F90
+traadv_eiv.F90
+traadv_tvd.F90
+trazdf_imp.F90
+zdftke.F90
+zdftmx.F90
+}}}
+
+Some other routines which are not compiled have also been modified for the sake of consistency. These are the routines relative to PISCES present in `ORCA1_LIM3_PISCES/MY_SRC` which get copied in `ORCA2_LIM3_PISCES`.
+
+They are:
+{{{
+p4zflx.F90
+p4zsed.F90
+p4zsms.F90
+trcrst.F90
+trcsms_age.F90
+trcsms_cfc.F90
+trcstp.F90
+trcwri.F90
+}}}
+
+You will need these updates if you want to add PISCES support to the Schwarz algorithm.
+
+Lastly, but much more important, some files have been modified to declare the variables needed to store the parameters of the Schwarz algorithm and the storing/restoring variables.
+
+The parameters are in:
+{{{
+in_out_manager.F90
+}}}
+and read in
+{{{
+domain.F90
+}}}
+
+
+The "Schwarz fields" named with `_s` appended at the end of the original variable name are defined in these subroutines:
+{{{
+dom_oce.F90
+dynspg_oce.F90
+ice.F90
+limdiahsb.F90
+oce.F90
+sbc_oce.F90
+sbcrnf.F90
+zdf_oce.F90
+zdftke.F90
+}}}