Changes between Version 3 and Version 4 of schwarz_nemo


Ignore:
Timestamp:
07/23/19 15:07:30 (2 years ago)
Author:
snguyen
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • schwarz_nemo

    v3 v4  
    1111With 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. 
    1212 
    13 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. 
    14  
    15 * One global loop over the coupling windows : 5 to iterate. Counter `iswloop`, limits `1 : nsloops`. 
    16 * One inner loop over the number of schwarz iteration for each coupling window : 6. Counter `kswr`, limits `1 : mswr`. 
    17 * 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` 
     13For 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. 
     14 
     15* One global loop over the coupling windows : 5 to iterate. Counter `iswloop`, limits `1 to nsloops`. 
     16* One inner loop over the number of schwarz iteration for each coupling window : 6. Counter `kswr`, limits `1 to mswr`. 
     17* 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` 
    1818 
    1919The counter given to subroutine `stp(istp)` is rewind to its value before entering the Schwarz loop when the Schwarz coupling window is repeated. 
     
    112112=== The Schwarz algorithm: storing and restoring === 
    113113 
    114 The Schwarz algorithm is implemented by storing and restoring the state of the ocean and ice at the proper times. 
     114The Schwarz algorithm is implemented by storing and restoring the state of the ocean and ice at the proper position inside the three loops structure. 
    115115 
    116116The routines used to store and restore a state are all in `schwarz.F90` which is a new module. 
     
    125125IF ( kswr > 1 ) CALL swz_reinit(istp) ! reset Ocean state to first time step for new schwarz iteration 
    126126}}} 
     127 
     128The 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.  
     129 
     130(In case you wish to use PISCES you will have to add the relevant routines in `swz_reinit` and module `schwarz.F90`) 
    127131 
    128132`Unfortunately this is not the only place variables are stored and restored!`