/[MITgcm]/MITgcm/pkg/layers/layers_output.F
ViewVC logotype

Diff of /MITgcm/pkg/layers/layers_output.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by rpa, Sun Sep 20 20:54:48 2009 UTC revision 1.3 by dfer, Wed Sep 30 15:58:29 2009 UTC
# Line 50  C     == Local variables == Line 50  C     == Local variables ==
50          pf(1:1) = 'R'          pf(1:1) = 'R'
51        ENDIF        ENDIF
52    
53        IF ( DIFFERENT_MULTIPLE(dumpFreq,myTime,deltaTClock)        IF ( DIFFERENT_MULTIPLE(layers_diagFreq,myTime,deltaTClock)
54       &  .AND. myIter .GT. 0 ) THEN       &  .AND. myIter .GT. 0 ) THEN
55    
56          IF ( layers_MDSIO ) THEN          IF ( layers_MDSIO ) THEN
# Line 85  C           But how? Line 85  C           But how?
85    
86  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
87    
88    c set arrays to zero if first timestep
89          IF ( myIter.EQ.nIter0 ) THEN
90           DO bj = myByLo(myThid), myByHi(myThid)
91            DO bi = myBxLo(myThid), myBxHi(myThid)
92    #ifdef LAYERS_UFLUX
93              CALL TIMEAVE_RESET(layers_UFlux_T,Nlayers,bi,bj,myThid)
94    #ifdef LAYERS_THICKNESS
95              CALL TIMEAVE_RESET(layers_HU_T,Nlayers,bi,bj,myThid)
96    #endif /* LAYERS_THICKNESS */
97    #endif /* LAYERS_UFLUX */
98    
99    #ifdef LAYERS_VFLUX
100              CALL TIMEAVE_RESET(layers_VFlux_T,Nlayers,bi,bj,myThid)
101    #ifdef LAYERS_THICKNESS
102              CALL TIMEAVE_RESET(layers_HV_T,Nlayers,bi,bj,myThid)
103    #endif /* LAYERS_THICKNESS */
104    #endif /* LAYERS_VFLUX */
105             DO k=1,Nlayers
106              layers_TimeAve(k,bi,bj)=0.
107             ENDDO
108            ENDDO
109           ENDDO
110    
111  C     Dump files and restart average computation if needed  C     Dump files and restart average computation if needed
112        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,deltaTClock)        ELSEIF (
113       & .AND. myIter .GT. 0   ) THEN       &  DIFFERENT_MULTIPLE(layers_taveFreq,myTime,deltaTClock)
114         &       ) THEN
115    
116  C      Normalize by integrated time  C      Normalize by integrated time
117         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22