13 |
c check conservation in the model |
c check conservation in the model |
14 |
C============================================================== |
C============================================================== |
15 |
SUBROUTINE MONOD_CONS( |
SUBROUTINE MONOD_CONS( |
16 |
& myIter,myTime,myThid) |
& myTime,myIter,myThid) |
17 |
#include "SIZE.h" |
#include "SIZE.h" |
18 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
19 |
#include "PARAMS.h" |
#include "PARAMS.h" |
163 |
c print*,'tmptotO',tmptotO |
c print*,'tmptotO',tmptotO |
164 |
#endif |
#endif |
165 |
|
|
166 |
|
_END_MASTER(myThid) |
167 |
|
|
168 |
|
IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN |
169 |
write(DAR_cons_unit1,*), myIter, tmptotP, tmptotP/voltot |
write(DAR_cons_unit1,*), myIter, tmptotP, tmptotP/voltot |
170 |
write(DAR_cons_unit2,*), myIter, tmptotN, tmptotN/voltot |
write(DAR_cons_unit2,*), myIter, tmptotN, tmptotN/voltot |
171 |
write(DAR_cons_unit3,*), myIter, tmptotFe, tmptotFe/voltot |
write(DAR_cons_unit3,*), myIter, tmptotFe, tmptotFe/voltot |
175 |
write(DAR_cons_unit6,*), myIter, tmptotA, tmptotA/voltot |
write(DAR_cons_unit6,*), myIter, tmptotA, tmptotA/voltot |
176 |
write(DAR_cons_unit7,*), myIter, tmptotO, tmptotO/voltot |
write(DAR_cons_unit7,*), myIter, tmptotO, tmptotO/voltot |
177 |
#endif |
#endif |
178 |
|
ENDIF |
|
_END_MASTER(myThid) |
|
179 |
|
|
180 |
C-- Everyone else must wait |
C-- Everyone else must wait |
181 |
_BARRIER |
_BARRIER |