18 |
#include "SIZE.h" |
#include "SIZE.h" |
19 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
20 |
#include "PARAMS.h" |
#include "PARAMS.h" |
21 |
#ifdef ALLOW_MNC |
#ifdef ALLOW_MNC |
22 |
#include "MNC_PARAMS.h" |
#include "MNC_PARAMS.h" |
23 |
#endif |
#endif |
24 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
54 |
C Final Time Averages and Dump Files if needed |
C Final Time Averages and Dump Files if needed |
55 |
IF (DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock)) THEN |
IF (DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock)) THEN |
56 |
|
|
57 |
useVariableK = useKPP .OR. useGMredi .OR. ivdc_kappa.NE.0. |
useVariableK = useKPP .OR. usePP81 .OR. useMY82 .OR. useGGL90 |
58 |
|
& .OR. useGMredi .OR. ivdc_kappa.NE.0. |
59 |
WRITE(suff,'(I10.10)') myIter |
WRITE(suff,'(I10.10)') myIter |
60 |
|
|
61 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
164 |
CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave, |
CALL WRITE_FLD_XYZ_RL('PhHytave.',suff,phiHydtave, |
165 |
& myIter,myThid) |
& myIter,myThid) |
166 |
#endif /* ndef HRCUBE */ |
#endif /* ndef HRCUBE */ |
167 |
IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN |
IF ( fluidIsWater ) THEN |
168 |
#ifndef HRCUBE |
#ifndef HRCUBE |
169 |
CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave, |
CALL WRITE_FLD_XY_RL('PHLtave.',suff,phiHydLowtave, |
170 |
& myIter,myThid) |
& myIter,myThid) |
221 |
IF (useMNC .AND. timeave_mnc) THEN |
IF (useMNC .AND. timeave_mnc) THEN |
222 |
|
|
223 |
CALL MNC_CW_SET_UDIM('tave', -1, myThid) |
CALL MNC_CW_SET_UDIM('tave', -1, myThid) |
224 |
CALL MNC_CW_I_W('I','tave',0,0,'iter', myIter, myThid) |
CALL MNC_CW_I_W_S('I','tave',0,0,'iter', myIter, myThid) |
225 |
CALL MNC_CW_SET_UDIM('tave', 0, myThid) |
CALL MNC_CW_SET_UDIM('tave', 0, myThid) |
226 |
CALL MNC_CW_RL_W('D','tave',0,0,'model_time',myTime,myThid) |
CALL MNC_CW_RL_W_S('D','tave',0,0,'model_time',myTime,myThid) |
|
|
|
227 |
CALL MNC_CW_RL_W('D','tave',0,0,'ETAtave',etatave,myThid) |
CALL MNC_CW_RL_W('D','tave',0,0,'ETAtave',etatave,myThid) |
228 |
#ifndef HRCUBE |
#ifndef HRCUBE |
229 |
C EH3: something is wriong here... |
C EH3: something is wriong here... |
259 |
CALL MNC_CW_RL_W('D','tave',0,0,'PhHytave', |
CALL MNC_CW_RL_W('D','tave',0,0,'PhHytave', |
260 |
& phiHydtave,myThid) |
& phiHydtave,myThid) |
261 |
#endif /* ndef HRCUBE */ |
#endif /* ndef HRCUBE */ |
262 |
IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN |
IF ( fluidIsWater ) THEN |
263 |
#ifndef HRCUBE |
#ifndef HRCUBE |
264 |
CALL MNC_CW_RL_W('D','tave',0,0,'PHLtave', |
CALL MNC_CW_RL_W('D','tave',0,0,'PHLtave', |
265 |
& phiHydLowtave,myThid) |
& phiHydLowtave,myThid) |