/[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.7 by gforget, Tue Sep 18 21:30:36 2012 UTC revision 1.8 by gforget, Wed Sep 19 18:48:18 2012 UTC
# Line 45  C     == Local variables == Line 45  C     == Local variables ==
45        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
46        EXTERNAL DIFFERENT_MULTIPLE        EXTERNAL DIFFERENT_MULTIPLE
47        CHARACTER*(MAX_LEN_MBUF) suff        CHARACTER*(MAX_LEN_MBUF) suff
48          INTEGER iLa
49  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
50        INTEGER bi, bj        INTEGER bi, bj
51  #endif  #endif
# Line 59  C---+----1----+----2----+----3----+----4 Line 60  C---+----1----+----2----+----3----+----4
60       &  CALL LAYERS_DIAGNOSTICS_FILL( myThid )       &  CALL LAYERS_DIAGNOSTICS_FILL( myThid )
61  #endif  #endif
62    
63    #ifdef ALLOW_LAYERS_OUTPUT
64    
65    cgf layers_maxNum loop and dimension would be needed for
66    cgf the following and tave output to work beyond iLa.EQ.1
67          DO iLa=1,layers_maxNum
68          IF ( layers_num(iLa).EQ.1 ) THEN      
69    
70  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
71        IF ( writeBinaryPrec .EQ. precFloat64 ) THEN        IF ( writeBinaryPrec .EQ. precFloat64 ) THEN
72          pf(1:1) = 'D'          pf(1:1) = 'D'
# Line 70  C---+----1----+----2----+----3----+----4 Line 78  C---+----1----+----2----+----3----+----4
78        IF ( DIFFERENT_MULTIPLE(layers_diagFreq,myTime,deltaTClock)        IF ( DIFFERENT_MULTIPLE(layers_diagFreq,myTime,deltaTClock)
79       &  .AND. myIter .GT. 0 ) THEN       &  .AND. myIter .GT. 0 ) THEN
80    
81    cgf this needs layers_maxNum loop
82          IF ( layers_MDSIO ) THEN          IF ( layers_MDSIO ) THEN
83            WRITE(suff,'(I10.10)') myIter            WRITE(suff,'(I10.10)') myIter
84  #ifdef LAYERS_UFLUX  #ifdef LAYERS_UFLUX
# Line 89  C---+----1----+----2----+----3----+----4 Line 98  C---+----1----+----2----+----3----+----4
98  #endif /* LAYERS_THICKNESS */  #endif /* LAYERS_THICKNESS */
99  #endif /* LAYERS_VFLUX */  #endif /* LAYERS_VFLUX */
100  #ifdef LAYERS_PRHO_REF  #ifdef LAYERS_PRHO_REF
101         IF (LAYER_nb .EQ. 3) THEN         IF (layers_num(1) .EQ. 3) THEN
102              CALL WRITE_FLD_3D_RL( 'layers_prho.', suff, Nr,              CALL WRITE_FLD_3D_RL( 'layers_prho.', suff, Nr,
103       &                             prho, myIter, myThid )       &                             prho, myIter, myThid )
104         ENDIF         ENDIF
# Line 158  C      Normalize by integrated time Line 167  C      Normalize by integrated time
167  #endif /* LAYERS_VFLUX */  #endif /* LAYERS_VFLUX */
168    
169  #ifdef LAYERS_PRHO_REF  #ifdef LAYERS_PRHO_REF
170         IF (LAYER_nb .EQ. 3) THEN         IF (layers_num(1) .EQ. 3) THEN
171            CALL TIMEAVE_NORMALIZE( prho_tave, layers_timeave,            CALL TIMEAVE_NORMALIZE( prho_tave, layers_timeave,
172       &                            Nr, bi, bj, myThid )       &                            Nr, bi, bj, myThid )
173         ENDIF         ENDIF
# Line 187  C      Normalize by integrated time Line 196  C      Normalize by integrated time
196  #endif /* LAYERS_VFLUX */  #endif /* LAYERS_VFLUX */
197    
198  #ifdef LAYERS_PRHO_REF  #ifdef LAYERS_PRHO_REF
199        IF (LAYER_nb .EQ. 3) THEN        IF (layers_num(1) .EQ. 3) THEN
200           CALL WRITE_FLD_3D_RL( 'layers_prho-tave.', suff, Nr,           CALL WRITE_FLD_3D_RL( 'layers_prho-tave.', suff, Nr,
201       &                          prho_tave, myIter, myThid )       &                          prho_tave, myIter, myThid )
202        ENDIF        ENDIF
# Line 217  C      Reset averages to zero Line 226  C      Reset averages to zero
226  #endif /* LAYERS_VFLUX */  #endif /* LAYERS_VFLUX */
227    
228  #ifdef LAYERS_PRHO_REF  #ifdef LAYERS_PRHO_REF
229            IF (LAYER_nb .EQ. 3) THEN            IF (layers_num(1) .EQ. 3) THEN
230            CALL TIMEAVE_RESET(prho_tave,Nr,bi,bj,myThid)            CALL TIMEAVE_RESET(prho_tave,Nr,bi,bj,myThid)
231            ENDIF            ENDIF
232  #endif /* LAYERS_PRHO_REF */  #endif /* LAYERS_PRHO_REF */
# Line 232  C--   end of bloc: if time is a multiple Line 241  C--   end of bloc: if time is a multiple
241        ENDIF        ENDIF
242  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
243    
244          ENDIF !IF ( layers_num(iLa).EQ.1 ) THEN      
245          ENDDO !DO iLa=1,layers_maxNum
246    
247    #endif /* ALLOW_LAYERS_OUTPUT */
248    
249  #endif /* ALLOW_LAYERS */  #endif /* ALLOW_LAYERS */
250    
251        RETURN        RETURN
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22