/[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.15 by edhill, Wed Aug 10 19:18:52 2005 UTC revision 1.16 by edhill, Thu Aug 11 02:50:12 2005 UTC
# Line 39  C     == Local variables == Line 39  C     == Local variables ==
39                    
40        IF (SEAICEwriteState) THEN        IF (SEAICEwriteState) THEN
41    
42           IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,deltaTClock)          IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,deltaTClock)
43       &      ) THEN       &       ) THEN
44                            
45  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
46             IF (useMNC .AND. SEAICE_dump_mnc) THEN            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
47               CALL MNC_CW_SET_UDIM('sice', -1, myThid)              CALL MNC_CW_SET_UDIM('sice', -1, myThid)
48               CALL MNC_CW_RL_W_S('D','sice',0,0,'T', myTime, myThid)              CALL MNC_CW_RL_W_S('D','sice',0,0,'T', myTime, myThid)
49               CALL MNC_CW_SET_UDIM('sice', 0, myThid)              CALL MNC_CW_SET_UDIM('sice', 0, myThid)
50               CALL MNC_CW_I_W_S('I','sice',0,0,'iter', myIter, myThid)              CALL MNC_CW_I_W_S('I','sice',0,0,'iter', myIter, myThid)
51               CALL MNC_CW_RL_W_S('D','sice',0,0,'model_time',              CALL MNC_CW_RL_W_S('D','sice',0,0,'model_time',
52       &            myTime,myThid)       &           myTime,myThid)
53               CALL MNC_CW_RL_W('D','sice',0,0,'si_UWIND',uwind,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'si_UWIND',uwind,myThid)
54               CALL MNC_CW_RL_W('D','sice',0,0,'si_VWIND',vwind,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'si_VWIND',vwind,myThid)
55               CALL MNC_CW_RS_W('D','sice',0,0,'fu',fu,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'fu',fu,myThid)
56               CALL MNC_CW_RS_W('D','sice',0,0,'fv',fv,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'fv',fv,myThid)
57               CALL MNC_CW_RS_W('D','sice',0,0,'EmPmR',EmPmR,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'EmPmR',EmPmR,myThid)
58               CALL MNC_CW_RS_W('D','sice',0,0,'Qnet',Qnet,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'Qnet',Qnet,myThid)
59               CALL MNC_CW_RS_W('D','sice',0,0,'Qsw',Qsw,myThid)              CALL MNC_CW_RS_W('D','sice',0,0,'Qsw',Qsw,myThid)
60             ENDIF            ENDIF
61  #endif  #endif
62             IF (SEAICE_dump_mdsio) THEN            IF (SEAICE_dump_mdsio) THEN
63               WRITE(suff,'(I10.10)') myIter              WRITE(suff,'(I10.10)') myIter
64               _BARRIER              _BARRIER
65               _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
66               CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)
67               CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)
68               CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)
69               CALL WRITE_FLD_XY_RS( 'FV.',suff,fv,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'FV.',suff,fv,myIter,myThid)
70               CALL WRITE_FLD_XY_RS( 'EmPmR.',suff,EmPmR,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'EmPmR.',suff,EmPmR,myIter,myThid)
71               CALL WRITE_FLD_XY_RS( 'Qnet.',suff,Qnet,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'Qnet.',suff,Qnet,myIter,myThid)
72               CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)
73               _END_MASTER( myThid )              _END_MASTER( myThid )
74               _BARRIER              _BARRIER
75             ENDIF            ENDIF
76              
77  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
78         CALL PLOT_FIELD_XYRS( uwind , 'Current uwind ', myIter, myThid )         CALL PLOT_FIELD_XYRS( uwind , 'Current uwind ', myIter, myThid )
79         CALL PLOT_FIELD_XYRS( vwind , 'Current vwind ', myIter, myThid )         CALL PLOT_FIELD_XYRS( vwind , 'Current vwind ', myIter, myThid )
# Line 93  C     == Local variables == Line 93  C     == Local variables ==
93         CALL PLOT_FIELD_XYRL( Qsw   , 'Current Qsw   ', myIter, myThid )         CALL PLOT_FIELD_XYRL( Qsw   , 'Current Qsw   ', myIter, myThid )
94  #endif  #endif
95    
96              DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
97                 DO bi=myBxLo(myThid),myBxHi(myThid)              DO bi=myBxLo(myThid),myBxHi(myThid)
98                    DO j=1,sNy                DO j=1,sNy
99                       DO i=1,sNx                  DO i=1,sNx
100                          arr(i,j,bi,bj)=UICE(i,j,1,bi,bj)                    arr(i,j,bi,bj)=UICE(i,j,1,bi,bj)
101                       ENDDO                  ENDDO
102                    ENDDO                ENDDO
                ENDDO  
103              ENDDO              ENDDO
104              ENDDO
105    
106              IF (SEAICE_dump_mdsio) THEN
107              _BARRIER              _BARRIER
108              _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
109              CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid)
110              _END_MASTER( myThid )              _END_MASTER( myThid )
111              _BARRIER              _BARRIER
112              ENDIF
113    #ifdef ALLOW_MNC
114              IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
115                CALL MNC_CW_RS_W('D','sice',0,0,'si_UICE',arr,myThid)
116              ENDIF
117    #endif
118  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
119         _EXCH_XY_R4( arr, myThid )            _EXCH_XY_R4( arr, myThid )
120         CALL PLOT_FIELD_XYRS( arr   , 'Current uice  ', myIter, myThid )            CALL PLOT_FIELD_XYRS( arr   , 'Current uice  ',
121         &         myIter, myThid )
122  #endif  #endif
123    
124              DO bj=myByLo(myThid),myByHi(myThid)            DO bj=myByLo(myThid),myByHi(myThid)
125                 DO bi=myBxLo(myThid),myBxHi(myThid)              DO bi=myBxLo(myThid),myBxHi(myThid)
126                    DO j=1,sNy                DO j=1,sNy
127                       DO i=1,sNx                  DO i=1,sNx
128                          arr(i,j,bi,bj)=VICE(i,j,1,bi,bj)                    arr(i,j,bi,bj)=VICE(i,j,1,bi,bj)
129                       ENDDO                  ENDDO
130                    ENDDO                ENDDO
                ENDDO  
131              ENDDO              ENDDO
132              ENDDO
133              IF (SEAICE_dump_mdsio) THEN
134              _BARRIER              _BARRIER
135              _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
136              CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid)
137              _END_MASTER( myThid )              _END_MASTER( myThid )
138              _BARRIER              _BARRIER
139  #ifdef SEAICE_DEBUG            ENDIF
140         _EXCH_XY_R4( arr, myThid )  #ifdef ALLOW_MNC
141         CALL PLOT_FIELD_XYRS( arr   , 'Current vice  ', myIter, myThid )            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
142                CALL MNC_CW_RS_W('D','sice',0,0,'si_VICE',arr,myThid)
143              ENDIF
144  #endif  #endif
145                #ifdef SEAICE_DEBUG
146              DO bj=myByLo(myThid),myByHi(myThid)            _EXCH_XY_R4( arr, myThid )
147                 DO bi=myBxLo(myThid),myBxHi(myThid)            CALL PLOT_FIELD_XYRS( arr   , 'Current vice  ',
148                    DO j=1,sNy       &         myIter, myThid )
149                       DO i=1,sNx  #endif
150                          arr(i,j,bi,bj)=HEFF(i,j,1,bi,bj)            
151                       ENDDO            DO bj=myByLo(myThid),myByHi(myThid)
152                    ENDDO              DO bi=myBxLo(myThid),myBxHi(myThid)
153                 ENDDO                DO j=1,sNy
154                    DO i=1,sNx
155                      arr(i,j,bi,bj)=HEFF(i,j,1,bi,bj)
156                    ENDDO
157                  ENDDO
158              ENDDO              ENDDO
159              ENDDO
160              IF (SEAICE_dump_mdsio) THEN
161              _BARRIER              _BARRIER
162              _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
163              CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid)
164              _END_MASTER( myThid )              _END_MASTER( myThid )
165              _BARRIER              _BARRIER
166  #ifdef SEAICE_DEBUG            ENDIF
167         _EXCH_XY_R4( arr, myThid )  #ifdef ALLOW_MNC
168         CALL PLOT_FIELD_XYRS( arr   , 'Current heff  ', myIter, myThid )            IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
169                CALL MNC_CW_RS_W('D','sice',0,0,'si_HEFF',arr,myThid)
170              ENDIF
171  #endif  #endif
172                #ifdef SEAICE_DEBUG
173              DO bj=myByLo(myThid),myByHi(myThid)            _EXCH_XY_R4( arr, myThid )
174                 DO bi=myBxLo(myThid),myBxHi(myThid)            CALL PLOT_FIELD_XYRS( arr   , 'Current heff  ',
175                    DO j=1,sNy       &         myIter, myThid )
176                       DO i=1,sNx  #endif
177                          arr(i,j,bi,bj)=AREA(i,j,1,bi,bj)            
178                       ENDDO            DO bj=myByLo(myThid),myByHi(myThid)
179                    ENDDO              DO bi=myBxLo(myThid),myBxHi(myThid)
180                 ENDDO                DO j=1,sNy
181                    DO i=1,sNx
182                      arr(i,j,bi,bj)=AREA(i,j,1,bi,bj)
183                    ENDDO
184                  ENDDO
185              ENDDO              ENDDO
186              ENDDO
187              IF (SEAICE_dump_mdsio) THEN
188              _BARRIER              _BARRIER
189              _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
190              CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid)              CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid)
191              _END_MASTER( myThid )              _END_MASTER( myThid )
192              _BARRIER              _BARRIER
193              ENDIF
194    #ifdef ALLOW_MNC
195              IF ( useMNC .AND. SEAICE_dump_mnc ) THEN
196                CALL MNC_CW_RS_W('D','sice',0,0,'si_AREA',arr,myThid)
197              ENDIF
198    #endif
199  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
200         _EXCH_XY_R4( arr, myThid )            _EXCH_XY_R4( arr, myThid )
201         CALL PLOT_FIELD_XYRS( arr   , 'Current area  ', myIter, myThid )            CALL PLOT_FIELD_XYRS( arr   , 'Current area  ',
202         &         myIter, myThid )
203  #endif  #endif
204    
205           ENDIF           ENDIF
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22