/[MITgcm]/MITgcm/pkg/seaice/seaice_do_diags.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_do_diags.F

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

revision 1.21 by jmc, Tue Jun 20 20:57:06 2006 UTC revision 1.22 by jmc, Wed Aug 9 04:33:25 2006 UTC
# Line 47  C     == Local variables == Line 47  C     == Local variables ==
47    
48          IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,deltaTClock)          IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,deltaTClock)
49       &       ) THEN       &       ) THEN
50                
51              _BARRIER
52  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
53            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
54              CALL MNC_CW_SET_UDIM('sice', -1, myThid)              CALL MNC_CW_SET_UDIM('sice', -1, myThid)
# Line 77  C     == Local variables == Line 78  C     == Local variables ==
78  #endif  #endif
79            IF (SEAICE_dump_mdsio) THEN            IF (SEAICE_dump_mdsio) THEN
80              WRITE(suff,'(I10.10)') myIter              WRITE(suff,'(I10.10)') myIter
             _BARRIER  
             _BEGIN_MASTER( myThid )  
81              CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)
82              CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)
83              CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)
# Line 88  C     == Local variables == Line 87  C     == Local variables ==
87              CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)
88  #if ( defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP) )  #if ( defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP) )
89              IF ( SEAICEuseEVP ) THEN              IF ( SEAICEuseEVP ) THEN
90               CALL WRITE_FLD_XY_Rs('SIGMA1.',suff,seaice_sigma1,               CALL WRITE_FLD_XY_RS('SIGMA1.',suff,seaice_sigma1,
91       &            myIter,myThid)       &            myIter,myThid)
92               CALL WRITE_FLD_XY_Rs('SIGMA2.',suff,seaice_sigma2,               CALL WRITE_FLD_XY_RS('SIGMA2.',suff,seaice_sigma2,
93       &            myIter,myThid)       &            myIter,myThid)
94               CALL WRITE_FLD_XY_Rs('SIGMA12.',suff,seaice_sigma12,               CALL WRITE_FLD_XY_RS('SIGMA12.',suff,seaice_sigma12,
95       &            myIter,myThid)       &            myIter,myThid)
96              ENDIF              ENDIF
97  #endif /* SEAICE_CGRID and SEAICE_ALLOW_EVP */  #endif /* SEAICE_CGRID and SEAICE_ALLOW_EVP */
             _END_MASTER( myThid )  
             _BARRIER  
98            ENDIF            ENDIF
99              
100  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
101         CALL PLOT_FIELD_XYRS( uwind , 'Current uwind ', myIter, myThid )         CALL PLOT_FIELD_XYRS( uwind , 'Current uwind ', myIter, myThid )
102         CALL PLOT_FIELD_XYRS( vwind , 'Current vwind ', myIter, myThid )         CALL PLOT_FIELD_XYRS( vwind , 'Current vwind ', myIter, myThid )
# Line 128  C     == Local variables == Line 125  C     == Local variables ==
125                ENDDO                ENDDO
126              ENDDO              ENDDO
127            ENDDO            ENDDO
128              _BARRIER
129            IF (SEAICE_dump_mdsio) THEN            IF (SEAICE_dump_mdsio) THEN
             _BARRIER  
             _BEGIN_MASTER( myThid )  
130              CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid)
             _END_MASTER( myThid )  
             _BARRIER  
131            ENDIF            ENDIF
132  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
133            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
# Line 143  C     == Local variables == Line 136  C     == Local variables ==
136  #endif  #endif
137  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
138            _EXCH_XY_R4( arr, myThid )            _EXCH_XY_R4( arr, myThid )
139            CALL PLOT_FIELD_XYRS( arr   , 'Current uice  ',            CALL PLOT_FIELD_XYRS( arr   , 'Current uice  ',
140       &         myIter, myThid )       &         myIter, myThid )
141  #endif  #endif
142              _BARRIER
143    
144            DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
145              DO bi=myBxLo(myThid),myBxHi(myThid)              DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 156  C     == Local variables == Line 150  C     == Local variables ==
150                ENDDO                ENDDO
151              ENDDO              ENDDO
152            ENDDO            ENDDO
153              _BARRIER
154            IF (SEAICE_dump_mdsio) THEN            IF (SEAICE_dump_mdsio) THEN
             _BARRIER  
             _BEGIN_MASTER( myThid )  
155              CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid)
             _END_MASTER( myThid )  
             _BARRIER  
156            ENDIF            ENDIF
157  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
158            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
# Line 170  C     == Local variables == Line 161  C     == Local variables ==
161  #endif  #endif
162  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
163            _EXCH_XY_R4( arr, myThid )            _EXCH_XY_R4( arr, myThid )
164            CALL PLOT_FIELD_XYRS( arr   , 'Current vice  ',            CALL PLOT_FIELD_XYRS( arr   , 'Current vice  ',
165       &         myIter, myThid )       &         myIter, myThid )
166  #endif  #endif
167                        _BARRIER
168    
169            DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
170              DO bi=myBxLo(myThid),myBxHi(myThid)              DO bi=myBxLo(myThid),myBxHi(myThid)
171                DO j=1-Oly,sNy+Oly                DO j=1-Oly,sNy+Oly
# Line 183  C     == Local variables == Line 175  C     == Local variables ==
175                ENDDO                ENDDO
176              ENDDO              ENDDO
177            ENDDO            ENDDO
178              _BARRIER
179            IF (SEAICE_dump_mdsio) THEN            IF (SEAICE_dump_mdsio) THEN
             _BARRIER  
             _BEGIN_MASTER( myThid )  
180              CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid)
             _END_MASTER( myThid )  
             _BARRIER  
181            ENDIF            ENDIF
182  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
183            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
# Line 200  C     == Local variables == Line 189  C     == Local variables ==
189            CALL PLOT_FIELD_XYRS( arr   , 'Current heff  ',            CALL PLOT_FIELD_XYRS( arr   , 'Current heff  ',
190       &         myIter, myThid )       &         myIter, myThid )
191  #endif  #endif
192                        _BARRIER
193    
194            DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
195              DO bi=myBxLo(myThid),myBxHi(myThid)              DO bi=myBxLo(myThid),myBxHi(myThid)
196                DO j=1-Oly,sNy+Oly                DO j=1-Oly,sNy+Oly
# Line 210  C     == Local variables == Line 200  C     == Local variables ==
200                ENDDO                ENDDO
201              ENDDO              ENDDO
202            ENDDO            ENDDO
203              _BARRIER
204            IF (SEAICE_dump_mdsio) THEN            IF (SEAICE_dump_mdsio) THEN
             _BARRIER  
             _BEGIN_MASTER( myThid )  
205              CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid)
             _END_MASTER( myThid )  
             _BARRIER  
206            ENDIF            ENDIF
207  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
208            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
# Line 227  C     == Local variables == Line 214  C     == Local variables ==
214            CALL PLOT_FIELD_XYRS( arr   , 'Current area  ',            CALL PLOT_FIELD_XYRS( arr   , 'Current area  ',
215       &         myIter, myThid )       &         myIter, myThid )
216  #endif  #endif
217              _BARRIER
218    
219           ENDIF           ENDIF
220        ENDIF        ENDIF
# Line 306  C      Normalize by integrated time Line 294  C      Normalize by integrated time
294       &                 bi,bj,myThid)       &                 bi,bj,myThid)
295          ENDDO          ENDDO
296         ENDDO         ENDDO
297    c      IF (myIter.EQ.10) WRITE(0,*) myThid, dumpFiles
298    
299           _BARRIER
300  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
301         IF (useMNC .AND. SEAICE_tave_mnc) THEN         IF (useMNC .AND. SEAICE_tave_mnc) THEN
302           CALL MNC_CW_SET_UDIM('sice_tave', -1, myThid)           CALL MNC_CW_SET_UDIM('sice_tave', -1, myThid)
# Line 337  C    &        myTime,myThid) Line 327  C    &        myTime,myThid)
327  #endif  #endif
328         IF (SEAICE_tave_mdsio) THEN         IF (SEAICE_tave_mdsio) THEN
329           WRITE(suff,'(I10.10)') myIter           WRITE(suff,'(I10.10)') myIter
          _BARRIER  
          _BEGIN_MASTER( myThid )  
330           CALL WRITE_FLD_XY_RL('FUtave.'   ,suff,FUtave   ,myIter,myThid)           CALL WRITE_FLD_XY_RL('FUtave.'   ,suff,FUtave   ,myIter,myThid)
331           CALL WRITE_FLD_XY_RL('FVtave.'   ,suff,FVtave   ,myIter,myThid)           CALL WRITE_FLD_XY_RL('FVtave.'   ,suff,FVtave   ,myIter,myThid)
332           CALL WRITE_FLD_XY_RL('EmPmRtave.',suff,EmPmRtave,myIter,myThid)           CALL WRITE_FLD_XY_RL('EmPmRtave.',suff,EmPmRtave,myIter,myThid)
# Line 348  C    &        myTime,myThid) Line 336  C    &        myTime,myThid)
336           CALL WRITE_FLD_XY_RL('VICEtave.' ,suff,VICEtave ,myIter,myThid)           CALL WRITE_FLD_XY_RL('VICEtave.' ,suff,VICEtave ,myIter,myThid)
337           CALL WRITE_FLD_XY_RL('HEFFtave.' ,suff,HEFFtave ,myIter,myThid)           CALL WRITE_FLD_XY_RL('HEFFtave.' ,suff,HEFFtave ,myIter,myThid)
338           CALL WRITE_FLD_XY_RL('AREAtave.' ,suff,AREAtave ,myIter,myThid)           CALL WRITE_FLD_XY_RL('AREAtave.' ,suff,AREAtave ,myIter,myThid)
          _END_MASTER( myThid )  
          _BARRIER  
339         ENDIF         ENDIF
340                 _BARRIER
341    
342  C      Reset averages to zero  C      Reset averages to zero
343         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
344          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 369  C      Reset averages to zero Line 356  C      Reset averages to zero
356           ENDDO           ENDDO
357          ENDDO          ENDDO
358         ENDDO         ENDDO
359          
360        ENDIF        ENDIF
361          
362  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
363    
364  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22