86 |
#include "PARAMS.h" |
#include "PARAMS.h" |
87 |
#include "CG2D.h" |
#include "CG2D.h" |
88 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
89 |
|
#include "AVER.h" |
90 |
|
|
91 |
C == Routine arguments == |
C == Routine arguments == |
92 |
C myThid - Thread number for this instance of the routine. |
C myThid - Thread number for this instance of the routine. |
109 |
myCurrentTime = startTime |
myCurrentTime = startTime |
110 |
myCurrentIter = nIter0 |
myCurrentIter = nIter0 |
111 |
|
|
112 |
|
C-- Dump for start state |
113 |
|
CALL WRITE_STATE( .TRUE., myCurrentTime, myCurrentIter, myThid ) |
114 |
|
|
115 |
C-- Begin time stepping loop |
C-- Begin time stepping loop |
116 |
DO I=1, nTimeSteps |
DO I=1, nTimeSteps |
117 |
|
|
121 |
C-- Step forward fields and calculate time tendency terms |
C-- Step forward fields and calculate time tendency terms |
122 |
CALL DYNAMICS( myCurrentTime, myCurrentIter, myThid ) |
CALL DYNAMICS( myCurrentTime, myCurrentIter, myThid ) |
123 |
|
|
124 |
|
C-- Do time averages |
125 |
|
IF (taveFreq.GT.0.) THEN |
126 |
|
CALL DO_TIME_AVERAGES( myCurrentTime, myCurrentIter, myThid ) |
127 |
|
ENDIF |
128 |
|
|
129 |
C-- Do IO if needed. |
C-- Do IO if needed. |
130 |
C Note: |
C Note: |
131 |
C ===== |
C ===== |
168 |
C-- Step-forward U/V/Theta/Salt for purposes of final I/O dump |
C-- Step-forward U/V/Theta/Salt for purposes of final I/O dump |
169 |
CALL DYNAMICS( myCurrentTime, myCurrentIter, myThid ) |
CALL DYNAMICS( myCurrentTime, myCurrentIter, myThid ) |
170 |
|
|
171 |
|
C-- Do time averages |
172 |
|
IF (taveFreq.GT.0.) THEN |
173 |
|
CALL DO_TIME_AVERAGES( myCurrentTime, myCurrentIter, myThid ) |
174 |
|
ENDIF |
175 |
|
|
176 |
C-- Dump for end state |
C-- Dump for end state |
177 |
CALL WRITE_STATE( myCurrentTime, myCurrentIter, myThid ) |
CALL WRITE_STATE( .FALSE., myCurrentTime, myCurrentIter, myThid ) |
178 |
|
|
179 |
|
|
180 |
RETURN |
RETURN |