/[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.6 by jmc, Mon Feb 7 03:07:49 2005 UTC revision 1.7 by edhill, Tue Feb 8 20:39:58 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 60  C     !LOCAL VARIABLES: Line 61  C     !LOCAL VARIABLES:
61  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
62        INTEGER ii        INTEGER ii
63        CHARACTER*(MAX_LEN_FNAM) diag_mnc_bn        CHARACTER*(MAX_LEN_FNAM) diag_mnc_bn
64          CHARACTER*(5) ctmp
65        INTEGER CW_DIMS, NLEN        INTEGER CW_DIMS, NLEN
66        PARAMETER ( CW_DIMS = 10 )        PARAMETER ( CW_DIMS = 10 )
67        PARAMETER ( NLEN    = 80 )        PARAMETER ( NLEN    = 80 )
68        INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)        INTEGER dim(CW_DIMS), ib(CW_DIMS), ie(CW_DIMS)
69        CHARACTER*(NLEN) dn(CW_DIMS)        CHARACTER*(NLEN) dn(CW_DIMS)
70          CHARACTER*(NLEN) d_cw_name
71        CHARACTER*(NLEN) dn_blnk        CHARACTER*(NLEN) dn_blnk
72          _RS ztmp(Nr+Nrphys)
73  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
74    
75  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
# Line 96  C       Update the record dimension by w Line 99  C       Update the record dimension by w
99          CALL MNC_CW_SET_UDIM(diag_mnc_bn, 0, myThid)          CALL MNC_CW_SET_UDIM(diag_mnc_bn, 0, myThid)
100    
101          dn(1)(1:NLEN) = dn_blnk(1:NLEN)          dn(1)(1:NLEN) = dn_blnk(1:NLEN)
102          WRITE(dn(1),'(a,i6.6)') 'Zd', nlevels(listnum)          WRITE(dn(1),'(a,i6.6)') 'Zmd', nlevels(listnum)
103          dim(1) = nlevels(listnum)          dim(1) = nlevels(listnum)
104          ib(1)  = 1          ib(1)  = 1
105          ie(1)  = nlevels(listnum)          ie(1)  = nlevels(listnum)
106    
107          CALL MNC_CW_ADD_GNAME('diag_levels', 1,          CALL MNC_CW_ADD_GNAME('diag_levels', 1,
108       &       dim, dn, ib, ie, myThid)       &       dim, dn, ib, ie, myThid)
109          CALL MNC_CW_ADD_VNAME(dn(1), 'diag_levels',          CALL MNC_CW_ADD_VNAME('diag_levels', 'diag_levels',
110       &       0,0, myThid)       &       0,0, myThid)
111          CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',          CALL MNC_CW_ADD_VATTR_TEXT('diag_levels','description',
112       &    'Idicies of vertical levels within the data source arrays',       &       'Idicies of vertical levels within the source arrays',
113       &       myThid)       &       myThid)
114                    
115          CALL MNC_CW_RL_W('I',diag_mnc_bn,0,0,          CALL MNC_CW_RL_W('I',diag_mnc_bn,0,0,
116       &       dn(1), levs(1,listnum), myThid)       &       'diag_levels', levs(1,listnum), myThid)
117    
118          CALL MNC_CW_DEL_VNAME(dn(1), myThid)          CALL MNC_CW_DEL_VNAME('diag_levels', myThid)
119          CALL MNC_CW_DEL_GNAME('diag_levels', myThid)          CALL MNC_CW_DEL_GNAME('diag_levels', myThid)
120    
121    C       Now define:  Zmdxxxxxx, Zudxxxxxx, Zldxxxxxx
122            ctmp(1:5) = 'mul  '
123            DO i = 1,3
124              dn(1)(1:NLEN) = dn_blnk(1:NLEN)
125              WRITE(dn(1),'(3a,i6.6)') 'Z',ctmp(i:i),'d',nlevels(listnum)
126              CALL MNC_CW_ADD_GNAME(dn(1), 1, dim, dn, ib, ie, myThid)
127              CALL MNC_CW_ADD_VNAME(dn(1), dn(1), 0,0, myThid)
128              DO j = 1,nlevels(listnum)
129                IF (i .EQ. 1) THEN
130                  ztmp(j) = rC(levs(j,listnum))
131                  CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
132         &             'Dimensional coordinate value at the mid point',
133         &             myThid)
134                ELSEIF (i .EQ. 2) THEN
135                  ztmp(j) = rF(levs(j,listnum))
136                  CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
137         &             'Dimensional coordinate value at the upper point',
138         &             myThid)
139                ELSEIF (i .EQ. 3) THEN
140                  ztmp(j) = rF(levs(j,listnum) + 1)
141                  CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
142         &             'Dimensional coordinate value at the lower point',
143         &             myThid)
144                ENDIF
145              ENDDO
146              CALL MNC_CW_RS_W('D',diag_mnc_bn,0,0, dn(1), ztmp, myThid)
147              CALL MNC_CW_DEL_VNAME(dn(1), myThid)
148              CALL MNC_CW_DEL_GNAME(dn(1), myThid)
149            ENDDO
150    
151        ENDIF        ENDIF
152  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
153    
# Line 216  C         Prepare for mdsio optionality Line 250  C         Prepare for mdsio optionality
250              _BEGIN_MASTER( myThid )              _BEGIN_MASTER( myThid )
251    
252              DO ii = 1,CW_DIMS              DO ii = 1,CW_DIMS
253                  d_cw_name(1:NLEN) = dn_blnk(1:NLEN)
254                dn(ii)(1:NLEN) = dn_blnk(1:NLEN)                dn(ii)(1:NLEN) = dn_blnk(1:NLEN)
255              ENDDO              ENDDO
256    
257    C           Note that the "d_cw_name" variable is a hack that hides a
258    C           subtlety within MNC.  Basically, each MNC-wrapped file is
259    C           caching its own concept of what each "grid name" (that is, a
260    C           dimension group name) means.  So one cannot re-use the same
261    C           "grid" name for different collections of dimensions within a
262    C           given file.  By appending the "m" values to each name, we
263    C           guarantee uniqueness within each MNC-produced file.
264                WRITE(d_cw_name,'(a,i6.6)') 'd_cw_',m
265    
266  C           XY dimensions  C           XY dimensions
267              dim(1)       = sNx + 2*OLx              dim(1)       = sNx + 2*OLx
268              dim(2)       = sNy + 2*OLy              dim(2)       = sNy + 2*OLy
269              ib(1)        = OLx + 1              ib(1)        = OLx + 1
270              ib(2)        = OLy + 1              ib(2)        = OLy + 1
271              IF (gdiag(n)(2:2) .EQ. 'M') THEN              IF (gdiag(m)(2:2) .EQ. 'M') THEN
272                dn(1)(1:2) = 'X'                dn(1)(1:2) = 'X'
273                ie(1)      = OLx + sNx                ie(1)      = OLx + sNx
274                dn(2)(1:2) = 'Y'                dn(2)(1:2) = 'Y'
275                ie(2)      = OLy + sNy                ie(2)      = OLy + sNy
276              ELSEIF (gdiag(n)(2:2) .EQ. 'U') THEN              ELSEIF (gdiag(m)(2:2) .EQ. 'U') THEN
277                dn(1)(1:3) = 'Xp1'                dn(1)(1:3) = 'Xp1'
278                ie(1)      = OLx + sNx + 1                ie(1)      = OLx + sNx + 1
279                dn(2)(1:2) = 'Y'                dn(2)(1:2) = 'Y'
280                ie(2)      = OLy + sNy                ie(2)      = OLy + sNy
281              ELSEIF (gdiag(n)(2:2) .EQ. 'V') THEN              ELSEIF (gdiag(m)(2:2) .EQ. 'V') THEN
282                dn(1)(1:2) = 'X'                dn(1)(1:2) = 'X'
283                ie(1)      = OLx + sNx                ie(1)      = OLx + sNx
284                dn(2)(1:3) = 'Yp1'                dn(2)(1:3) = 'Yp1'
285                ie(2)      = OLy + sNy + 1                ie(2)      = OLy + sNy + 1
286              ELSEIF (gdiag(n)(2:2) .EQ. 'Z') THEN              ELSEIF (gdiag(m)(2:2) .EQ. 'Z') THEN
287                dn(1)(1:3) = 'Xp1'                dn(1)(1:3) = 'Xp1'
288                ie(1)      = OLx + sNx + 1                ie(1)      = OLx + sNx + 1
289                dn(2)(1:3) = 'Yp1'                dn(2)(1:3) = 'Yp1'
# Line 248  C           XY dimensions Line 292  C           XY dimensions
292                            
293  C           Z is special since it varies  C           Z is special since it varies
294              WRITE(dn(3),'(a,i6.6)') 'Zd', nlevels(listnum)              WRITE(dn(3),'(a,i6.6)') 'Zd', nlevels(listnum)
295                IF ( (gdiag(m)(10:10) .EQ. 'R')
296         &           .AND. (gdiag(m)(9:9) .EQ. 'M') ) THEN
297                  WRITE(dn(3),'(a,i6.6)') 'Zmd', nlevels(listnum)
298                ENDIF
299                IF ( (gdiag(m)(10:10) .EQ. 'R')
300         &           .AND. (gdiag(m)(9:9) .EQ. 'L') ) THEN
301                  WRITE(dn(3),'(a,i6.6)') 'Zld', nlevels(listnum)
302                ENDIF
303                IF ( (gdiag(m)(10:10) .EQ. 'R')
304         &           .AND. (gdiag(m)(9:9) .EQ. 'U') ) THEN
305                  WRITE(dn(3),'(a,i6.6)') 'Zud', nlevels(listnum)
306                ENDIF
307              dim(3) = Nr+Nrphys              dim(3) = Nr+Nrphys
308              ib(3)  = 1              ib(3)  = 1
309              ie(3)  = nlevels(listnum)              ie(3)  = nlevels(listnum)
# Line 258  C           Time dimension Line 314  C           Time dimension
314              ib(4)  = 1              ib(4)  = 1
315              ie(4)  = 1              ie(4)  = 1
316    
317              CALL MNC_CW_ADD_GNAME('diag_cw_temp', 4,              CALL MNC_CW_ADD_GNAME(d_cw_name, 4,
318       &             dim, dn, ib, ie, myThid)       &             dim, dn, ib, ie, myThid)
319              CALL MNC_CW_ADD_VNAME(cdiag(m), 'diag_cw_temp',              CALL MNC_CW_ADD_VNAME(cdiag(m), d_cw_name,
320       &             4,5, myThid)       &             4,5, myThid)
321              CALL MNC_CW_ADD_VATTR_TEXT(cdiag(m),'description',              CALL MNC_CW_ADD_VATTR_TEXT(cdiag(m),'description',
322       &             tdiag(m),myThid)       &             tdiag(m),myThid)
# Line 271  C           Time dimension Line 327  C           Time dimension
327       &             cdiag(m), qtmp1, myThid)       &             cdiag(m), qtmp1, myThid)
328    
329              CALL MNC_CW_DEL_VNAME(cdiag(m), myThid)              CALL MNC_CW_DEL_VNAME(cdiag(m), myThid)
330              CALL MNC_CW_DEL_GNAME('diag_cw_temp', myThid)              CALL MNC_CW_DEL_GNAME(d_cw_name, myThid)
331    
332              _END_MASTER( myThid )              _END_MASTER( myThid )
333    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22