/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_out.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_out.F

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

revision 1.4 by edhill, Sat Dec 18 02:18:55 2004 UTC revision 1.13 by edhill, Mon May 2 21:24:12 2005 UTC
# Line 21  C     !USES: Line 21  C     !USES:
21  #include "SIZE.h"  #include "SIZE.h"
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24    #include "GRID.h"
25  #include "DIAGNOSTICS_SIZE.h"  #include "DIAGNOSTICS_SIZE.h"
26  #include "DIAGNOSTICS.h"  #include "DIAGNOSTICS.h"
27    
# Line 43  C     !LOCAL VARIABLES: Line 44  C     !LOCAL VARIABLES:
44        INTEGER i, j, k, m, n, bi, bj        INTEGER i, j, k, m, n, bi, bj
45        CHARACTER*8 parms1        CHARACTER*8 parms1
46        CHARACTER*3 mate_index        CHARACTER*3 mate_index
47        INTEGER mate        INTEGER mate, mVec
48        _RL qtmp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+Nrphys,nSx,nSy)        _RL qtmp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+Nrphys,nSx,nSy)
       _RL qtmp2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+Nrphys,nSx,nSy)  
49        _RL undef, getcon        _RL undef, getcon
50        EXTERNAL getcon        EXTERNAL getcon
51        INTEGER ILNBLNK        INTEGER ILNBLNK
52        EXTERNAL ILNBLNK        EXTERNAL ILNBLNK
53        INTEGER ilen        INTEGER ilen
54    
55        CHARACTER*(MAX_LEN_FNAM) pref        INTEGER ioUnit
56          CHARACTER*(MAX_LEN_FNAM) fn
57        CHARACTER*(MAX_LEN_MBUF) suff        CHARACTER*(MAX_LEN_MBUF) suff
58        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
       CHARACTER*(80) fn  
59        LOGICAL glf        LOGICAL glf
60  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
61        INTEGER ii        INTEGER ii
62        CHARACTER*(MAX_LEN_FNAM) diag_mnc_bn        CHARACTER*(MAX_LEN_FNAM) diag_mnc_bn
63          CHARACTER*(5) ctmp
64        INTEGER CW_DIMS, NLEN        INTEGER CW_DIMS, NLEN
65        PARAMETER ( CW_DIMS = 10 )        PARAMETER ( CW_DIMS = 10 )
66        PARAMETER ( NLEN    = 80 )        PARAMETER ( NLEN    = 80 )
67        INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)        INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)
68        CHARACTER*(NLEN) dn(CW_DIMS)        CHARACTER*(NLEN) dn(CW_DIMS)
69          CHARACTER*(NLEN) d_cw_name
70        CHARACTER*(NLEN) dn_blnk        CHARACTER*(NLEN) dn_blnk
71          _RS ztmp(Nr+Nrphys)
72  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
73    
74  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
75    
76          ioUnit= standardMessageUnit
77        undef = getcon('UNDEF')        undef = getcon('UNDEF')
78        glf = globalFiles        glf = globalFiles
79        WRITE(suff,'(I10.10)') myIter        WRITE(suff,'(I10.10)') myIter
80        pref = fnames(listnum)        ilen = ILNBLNK(fnames(listnum))
81        ilen=ILNBLNK( pref )        WRITE( fn, '(A,A,A)' ) fnames(listnum)(1:ilen),'.',suff(1:10)
       WRITE( fn, '(A,A,A)' ) pref(1:ilen),'.',suff(1:10)  
82    
83  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
84        IF (useMNC .AND. diag_mnc) THEN        IF (useMNC .AND. diag_mnc) THEN
# Line 86  C---+----1----+----2----+----3----+----4 Line 88  C---+----1----+----2----+----3----+----4
88          DO i = 1,NLEN          DO i = 1,NLEN
89            dn_blnk(i:i) = ' '            dn_blnk(i:i) = ' '
90          ENDDO          ENDDO
91  c       WRITE( diag_mnc_bn, '(A,A)' ) 'diag.', pref(1:ilen)          WRITE( diag_mnc_bn, '(A)' ) fnames(listnum)(1:ilen)
         WRITE( diag_mnc_bn, '(A)' ) pref(1:ilen)  
92    
93  C       Update the record dimension by writing the iteration number  C       Update the record dimension by writing the iteration number
94          CALL MNC_CW_SET_UDIM(diag_mnc_bn, -1, myThid)          CALL MNC_CW_SET_UDIM(diag_mnc_bn, -1, myThid)
# Line 95  C       Update the record dimension by w Line 96  C       Update the record dimension by w
96          CALL MNC_CW_SET_UDIM(diag_mnc_bn, 0, myThid)          CALL MNC_CW_SET_UDIM(diag_mnc_bn, 0, myThid)
97    
98          dn(1)(1:NLEN) = dn_blnk(1:NLEN)          dn(1)(1:NLEN) = dn_blnk(1:NLEN)
99          WRITE(dn(1),'(a,i6.6)') 'Zd', nlevels(listnum)          WRITE(dn(1),'(a,i6.6)') 'Zmd', nlevels(listnum)
100          dim(1) = nlevels(listnum)          dim(1) = nlevels(listnum)
101          ib(1)  = 1          ib(1)  = 1
102          ie(1)  = nlevels(listnum)          ie(1)  = nlevels(listnum)
103    
104          CALL MNC_CW_ADD_GNAME('diag_levels', 1,          CALL MNC_CW_ADD_GNAME('diag_levels', 1,
105       &       dim, dn, ib, ie, myThid)       &       dim, dn, ib, ie, myThid)
106          CALL MNC_CW_ADD_VNAME('level_indicies', 'diag_levels',          CALL MNC_CW_ADD_VNAME('diag_levels', 'diag_levels',
107       &       0,0, myThid)       &       0,0, myThid)
108          CALL MNC_CW_ADD_VATTR_TEXT('level_indicies','description',          CALL MNC_CW_ADD_VATTR_TEXT('diag_levels','description',
109       &    'Idicies of vertical levels within the data source arrays',       &       'Idicies of vertical levels within the source arrays',
110       &       myThid)       &       myThid)
111                    
112          CALL MNC_CW_RL_W('I',diag_mnc_bn,0,0,          CALL MNC_CW_RL_W('D',diag_mnc_bn,0,0,
113       &       'level_indicies', levs(1,listnum), myThid)       &       'diag_levels', levs(1,listnum), myThid)
114    
115          CALL MNC_CW_DEL_VNAME('level_indicies', myThid)          CALL MNC_CW_DEL_VNAME('diag_levels', myThid)
116          CALL MNC_CW_DEL_GNAME('diag_levels', myThid)          CALL MNC_CW_DEL_GNAME('diag_levels', myThid)
117    
118    C       Now define:  Zmdxxxxxx, Zudxxxxxx, Zldxxxxxx
119            ctmp(1:5) = 'mul  '
120            DO i = 1,3
121              dn(1)(1:NLEN) = dn_blnk(1:NLEN)
122              WRITE(dn(1),'(3a,i6.6)') 'Z',ctmp(i:i),'d',nlevels(listnum)
123              CALL MNC_CW_ADD_GNAME(dn(1), 1, dim, dn, ib, ie, myThid)
124              CALL MNC_CW_ADD_VNAME(dn(1), dn(1), 0,0, myThid)
125    
126    C         The following three ztmp() loops should eventually be modified
127    C         to reflect the fractional nature of levs(j,l) -- they should
128    C         do something like:
129    C            ztmp(j) = rC(INT(FLOOR(levs(j,l))))
130    C                      + ( rC(INT(FLOOR(levs(j,l))))
131    C                          + rC(INT(CEIL(levs(j,l)))) )
132    C                        / ( levs(j,l) - FLOOR(levs(j,l)) )
133    C         for averaged levels.
134              IF (i .EQ. 1) THEN
135                DO j = 1,nlevels(listnum)
136                  ztmp(j) = rC(NINT(levs(j,listnum)))
137                ENDDO
138                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
139         &           'Dimensional coordinate value at the mid point',
140         &           myThid)
141              ELSEIF (i .EQ. 2) THEN
142                DO j = 1,nlevels(listnum)
143                  ztmp(j) = rF(NINT(levs(j,listnum)) + 1)
144                ENDDO
145                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
146         &           'Dimensional coordinate value at the upper point',
147         &           myThid)
148              ELSEIF (i .EQ. 3) THEN
149                DO j = 1,nlevels(listnum)
150                  ztmp(j) = rF(NINT(levs(j,listnum)))
151                ENDDO
152                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
153         &           'Dimensional coordinate value at the lower point',
154         &           myThid)
155              ENDIF
156              CALL MNC_CW_RS_W('D',diag_mnc_bn,0,0, dn(1), ztmp, myThid)
157              CALL MNC_CW_DEL_VNAME(dn(1), myThid)
158              CALL MNC_CW_DEL_GNAME(dn(1), myThid)
159            ENDDO
160    
161        ENDIF        ENDIF
162  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
163    
# Line 160  C-        Empty diagnostics case : Line 205  C-        Empty diagnostics case :
205  C-        diagnostics is not empty :  C-        diagnostics is not empty :
206    
207              IF ( myThid.EQ.1 )              IF ( myThid.EQ.1 )
208       &             WRITE(6,2000) m,cdiag(m),ndiag(m),gdiag(m)       &             WRITE(ioUnit,2000) m,cdiag(m),ndiag(m),gdiag(m)
             IF ( parms1(5:5).NE.'C' ) THEN  
209    
210                DO k = 1,nlevels(listnum)              IF ( parms1(5:5).EQ.'C' ) THEN
211                  CALL GETDIAG (levs(k,listnum),m,undef,qtmp1,myThid)  C             Check for Mate of a Counter Diagnostic
212                ENDDO  C             --------------------------------------
213                  mate_index = parms1(6:8)
214                  READ (mate_index,'(I3)') mate
215                  IF ( myThid.EQ.1 )
216         &                 WRITE(ioUnit,2003) cdiag(m),mate,cdiag(mate)
217                ELSE
218                  mate = 0
219    
220  C             Check for Mate of a Vector Diagnostic  C             Check for Mate of a Vector Diagnostic
221  C             -------------------------------------  C             -------------------------------------
222                IF ( parms1(1:1).EQ.'U' .OR. parms1(1:1).EQ.'V' ) THEN                IF ( parms1(1:1).EQ.'U' .OR. parms1(1:1).EQ.'V' ) THEN
223                  mate_index = parms1(6:8)                  mate_index = parms1(6:8)
224                  READ (mate_index,'(I3)') mate                  READ (mate_index,'(I3)') mVec
225                  IF ( idiag(mate).NE.0 ) THEN                  IF ( idiag(mVec).NE.0 ) THEN
226                    IF ( myThid.EQ.1 )                    IF ( myThid.EQ.1 )
227       &                 WRITE(6,2001) cdiag(m),mate,cdiag(mate)       &                 WRITE(ioUnit,2001) cdiag(m),mVec,cdiag(mVec)
228                  ELSE                  ELSE
229                    IF ( myThid.EQ.1 )                    IF ( myThid.EQ.1 )
230       &                 WRITE(6,2002) cdiag(m),mate,cdiag(mate)       &                 WRITE(ioUnit,2002) cdiag(m),mVec,cdiag(mVec)
231                  ENDIF                  ENDIF
232                ENDIF                ENDIF
233                                ENDIF
             ELSE  
234    
235  C             Check for Mate of a Counter Diagnostic              DO bj = myByLo(myThid), myByHi(myThid)
236  C             --------------------------------------               DO bi = myBxLo(myThid), myBxHi(myThid)
               mate_index = parms1(6:8)  
               READ (mate_index,'(I3)') mate  
               IF ( myThid.EQ.1 )  
      &                 WRITE(6,2003) cdiag(m),mate,cdiag(mate)  
237                DO k = 1,nlevels(listnum)                DO k = 1,nlevels(listnum)
238                  CALL getdiag2(levs(k,listnum),m,undef,qtmp1,myThid)                  CALL GETDIAG(
239                  CALL getdiag2(levs(k,listnum),mate,undef,qtmp2,myThid)       I                       levs(k,listnum),undef,
240                  DO bj = myByLo(myThid), myByHi(myThid)       O                       qtmp1(1-OLx,1-OLy,k,bi,bj),
241                    DO bi = myBxLo(myThid), myBxHi(myThid)       I                       m,mate,bi,bj,myThid)
                     DO j = 1,sNy  
                       DO i = 1,sNx  
                         IF (qtmp2(i,j,k,bi,bj).NE.0.) THEN  
                           qtmp1(i,j,k,bi,bj) =  
      &                         qtmp1(i,j,k,bi,bj) / qtmp2(i,j,k,bi,bj)  
                         ELSE  
                           qtmp1(i,j,k,bi,bj) = undef  
                         ENDIF  
                       ENDDO  
                     ENDDO  
                   ENDDO  
                 ENDDO  
242                ENDDO                ENDDO
243                 ENDDO
244              ENDIF              ENDDO
245    
246  C-        end of empty diag / not empty block  C-        end of empty diag / not empty block
247            ENDIF            ENDIF
# Line 216  C-        end of empty diag / not empty Line 249  C-        end of empty diag / not empty
249  #ifdef ALLOW_MDSIO  #ifdef ALLOW_MDSIO
250  C         Prepare for mdsio optionality  C         Prepare for mdsio optionality
251            IF (diag_mdsio) THEN            IF (diag_mdsio) THEN
252              CALL mdswritefield_new(fn,writeBinaryPrec,glf,'RL',              IF (fflags(listnum)(1:1) .EQ. ' ') THEN
253       &           Nr+Nrphys,nlevels(listnum),qtmp1,n,myIter,myThid)  C             This is the old default behavior
254                  CALL mdswritefield_new(fn,writeBinaryPrec,glf,'RL',
255         &             Nr+Nrphys,nlevels(listnum),qtmp1,n,myIter,myThid)
256                ELSEIF (fflags(listnum)(1:1) .EQ. 'R') THEN
257    C             Force it to be 32-bit precision
258                  CALL mdswritefield_new(fn,precFloat32,glf,'RL',
259         &             Nr+Nrphys,nlevels(listnum),qtmp1,n,myIter,myThid)
260                ELSEIF (fflags(listnum)(1:1) .EQ. 'D') THEN
261    C             Force it to be 64-bit precision
262                  CALL mdswritefield_new(fn,precFloat64,glf,'RL',
263         &             Nr+Nrphys,nlevels(listnum),qtmp1,n,myIter,myThid)
264                ENDIF
265            ENDIF            ENDIF
266  #endif /*  ALLOW_MDSIO  */  #endif /*  ALLOW_MDSIO  */
267    
# Line 227  C         Prepare for mdsio optionality Line 271  C         Prepare for mdsio optionality
271              _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
272    
273              DO ii = 1,CW_DIMS              DO ii = 1,CW_DIMS
274                  d_cw_name(1:NLEN) = dn_blnk(1:NLEN)
275                dn(ii)(1:NLEN) = dn_blnk(1:NLEN)                dn(ii)(1:NLEN) = dn_blnk(1:NLEN)
276              ENDDO              ENDDO
277    
278              dn(1)(1:2) = 'Xd'  C           Note that the "d_cw_name" variable is a hack that hides a
279              dim(1) = sNx + 2*OLx  C           subtlety within MNC.  Basically, each MNC-wrapped file is
280              ib(1)  = OLx + 1  C           caching its own concept of what each "grid name" (that is, a
281              ie(1)  = OLx + sNx  C           dimension group name) means.  So one cannot re-use the same
282              dn(2)(1:2) = 'Yd'  C           "grid" name for different collections of dimensions within a
283              dim(2) = sNy + 2*OLy  C           given file.  By appending the "m" values to each name, we
284              ib(2)  = OLy + 1  C           guarantee uniqueness within each MNC-produced file.
285              ie(2)  = OLy + sNy              WRITE(d_cw_name,'(a,i6.6)') 'd_cw_',m
286    
287    C           XY dimensions
288                dim(1)       = sNx + 2*OLx
289                dim(2)       = sNy + 2*OLy
290                ib(1)        = OLx + 1
291                ib(2)        = OLy + 1
292                IF (gdiag(m)(2:2) .EQ. 'M') THEN
293                  dn(1)(1:2) = 'X'
294                  ie(1)      = OLx + sNx
295                  dn(2)(1:2) = 'Y'
296                  ie(2)      = OLy + sNy
297                ELSEIF (gdiag(m)(2:2) .EQ. 'U') THEN
298                  dn(1)(1:3) = 'Xp1'
299                  ie(1)      = OLx + sNx + 1
300                  dn(2)(1:2) = 'Y'
301                  ie(2)      = OLy + sNy
302                ELSEIF (gdiag(m)(2:2) .EQ. 'V') THEN
303                  dn(1)(1:2) = 'X'
304                  ie(1)      = OLx + sNx
305                  dn(2)(1:3) = 'Yp1'
306                  ie(2)      = OLy + sNy + 1
307                ELSEIF (gdiag(m)(2:2) .EQ. 'Z') THEN
308                  dn(1)(1:3) = 'Xp1'
309                  ie(1)      = OLx + sNx + 1
310                  dn(2)(1:3) = 'Yp1'
311                  ie(2)      = OLy + sNy + 1
312                ENDIF
313                
314  C           Z is special since it varies  C           Z is special since it varies
315              WRITE(dn(3),'(a,i6.6)') 'Zd', nlevels(listnum)              WRITE(dn(3),'(a,i6.6)') 'Zd', nlevels(listnum)
316                IF ( (gdiag(m)(10:10) .EQ. 'R')
317         &           .AND. (gdiag(m)(9:9) .EQ. 'M') ) THEN
318                  WRITE(dn(3),'(a,i6.6)') 'Zmd', nlevels(listnum)
319                ENDIF
320                IF ( (gdiag(m)(10:10) .EQ. 'R')
321         &           .AND. (gdiag(m)(9:9) .EQ. 'L') ) THEN
322                  WRITE(dn(3),'(a,i6.6)') 'Zld', nlevels(listnum)
323                ENDIF
324                IF ( (gdiag(m)(10:10) .EQ. 'R')
325         &           .AND. (gdiag(m)(9:9) .EQ. 'U') ) THEN
326                  WRITE(dn(3),'(a,i6.6)') 'Zud', nlevels(listnum)
327                ENDIF
328              dim(3) = Nr+Nrphys              dim(3) = Nr+Nrphys
329              ib(3)  = 1              ib(3)  = 1
330              ie(3)  = nlevels(listnum)              ie(3)  = nlevels(listnum)
331    
332              CALL MNC_CW_ADD_GNAME('diag_cw_temp', 3,  C           Time dimension
333                dn(4)(1:1) = 'T'
334                dim(4) = -1
335                ib(4)  = 1
336                ie(4)  = 1
337    
338                CALL MNC_CW_ADD_GNAME(d_cw_name, 4,
339       &             dim, dn, ib, ie, myThid)       &             dim, dn, ib, ie, myThid)
340              CALL MNC_CW_ADD_VNAME(cdiag(m), 'diag_cw_temp',              CALL MNC_CW_ADD_VNAME(cdiag(m), d_cw_name,
341       &             4,5, myThid)       &             4,5, myThid)
342              CALL MNC_CW_ADD_VATTR_TEXT(cdiag(m),'description',              CALL MNC_CW_ADD_VATTR_TEXT(cdiag(m),'description',
343       &             tdiag(m),myThid)       &             tdiag(m),myThid)
344              CALL MNC_CW_ADD_VATTR_TEXT(cdiag(m),'units',              CALL MNC_CW_ADD_VATTR_TEXT(cdiag(m),'units',
345       &             udiag(m),myThid)       &             udiag(m),myThid)
346    
347              CALL MNC_CW_RL_W('D',diag_mnc_bn,0,0,              IF ((fflags(listnum)(1:1) .EQ. ' ')
348         &           .OR. (fflags(listnum)(1:1) .EQ. 'R')) THEN
349                  CALL MNC_CW_RL_W('R',diag_mnc_bn,0,0,
350       &             cdiag(m), qtmp1, myThid)       &             cdiag(m), qtmp1, myThid)
351                ELSEIF (fflags(listnum)(1:1) .EQ. 'D') THEN
352                  CALL MNC_CW_RL_W('D',diag_mnc_bn,0,0,
353         &             cdiag(m), qtmp1, myThid)
354                ENDIF
355                
356              CALL MNC_CW_DEL_VNAME(cdiag(m), myThid)              CALL MNC_CW_DEL_VNAME(cdiag(m), myThid)
357              CALL MNC_CW_DEL_GNAME('diag_cw_temp', myThid)              CALL MNC_CW_DEL_GNAME(d_cw_name, myThid)
358    
359              _END_MASTER( myThid )              _END_MASTER( myThid )
360    
# Line 269  C--     end of Processing Fld # n Line 365  C--     end of Processing Fld # n
365          ENDIF          ENDIF
366        ENDDO        ENDDO
367    
368   2000 format(1x,'Computing Diagnostic # ',i3,2x,a8,5x,'Counter: ',   2000 format(1x,'Computing Diagnostic # ',i3,2x,a8,5x,
369       &     i4,6x,'Parms: ',a16)       &     'Counter:',i8,3x,'Parms: ',a16)
370   2001 format(1x,'          Vector  Mate for  ',a8,5x,   2001 format(1x,'          Vector  Mate for  ',a8,5x,
371       &     'Diagnostic # ',i3,2x,a8,' exists ')       &     'Diagnostic # ',i3,2x,a8,' exists ')
372   2002 format(1x,'        Vector  Mate for  ',a8,5x,   2002 format(1x,'          Vector  Mate for  ',a8,5x,
373       &     'Diagnostic # ',i3,2x,a8,' not enabled')       &     'Diagnostic # ',i3,2x,a8,' not enabled')
374   2003 format(1x,'      use Counter Mate for  ',a8,5x,   2003 format(1x,'      use Counter Mate for  ',a8,5x,
375       &     'Diagnostic # ',i3,2x,a8)       &     'Diagnostic # ',i3,2x,a8)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22