1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
|
cswdice modified for ice -- 25 Oct 01 |
|
|
cswddic modified for dic -- 16 July 02 |
|
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
5 |
#include "GCHEM_OPTIONS.h" |
#include "GCHEM_OPTIONS.h" |
17 |
#include "SIZE.h" |
#include "SIZE.h" |
18 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
19 |
#include "PARAMS.h" |
#include "PARAMS.h" |
|
cswdice --- add -- |
|
|
#ifdef DIC_BIOTIC |
|
20 |
#include "DIC_BIOTIC.h" |
#include "DIC_BIOTIC.h" |
|
#endif |
|
|
cswdice -- end add --- |
|
21 |
|
|
22 |
C == Routine arguments == |
C == Routine arguments == |
23 |
C myTime - Current time of simulation ( s ) |
C myTime - Current time of simulation ( s ) |
40 |
|
|
41 |
#ifdef ALLOW_TIMEAVE |
#ifdef ALLOW_TIMEAVE |
42 |
|
|
|
c set arrays to zero if first timestep |
|
|
IF ( myIter .EQ. nIter0+1 ) THEN |
|
|
DO bj = myByLo(myThid), myByHi(myThid) |
|
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
|
|
CALL TIMEAVE_RESET(BIOave, Nr, bi, bj, myThid) |
|
|
CALL TIMEAVE_RESET(CARave, Nr, bi, bj, myThid) |
|
|
CALL TIMEAVE_RESET(SURave, 1, bi, bj, myThid) |
|
|
CALL TIMEAVE_RESET(SUROave, 1, bi, bj, myThid) |
|
|
CALL TIMEAVE_RESET(pCO2ave, 1, bi, bj, myThid) |
|
|
CALL TIMEAVE_RESET(pHave, 1, bi, bj, myThid) |
|
|
CALL TIMEAVE_RESET(fluxCO2ave, 1, bi, bj, myThid) |
|
|
print*,'QQ start timeave' |
|
|
do k=1,Nr |
|
|
DIC_TimeAve(bi,bj,k)=0.d0 |
|
|
enddo |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
|
|
|
43 |
C Dump files and restart average computation if needed |
C Dump files and restart average computation if needed |
44 |
IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock) |
IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock) |
45 |
& ) THEN |
& ) THEN |