/[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.9 by edhill, Mon Feb 28 19:38:30 2005 UTC revision 1.10 by edhill, Tue Mar 1 15:48:03 2005 UTC
# Line 125  C       Now define:  Zmdxxxxxx, Zudxxxxx Line 125  C       Now define:  Zmdxxxxxx, Zudxxxxx
125            WRITE(dn(1),'(3a,i6.6)') 'Z',ctmp(i:i),'d',nlevels(listnum)            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)            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)            CALL MNC_CW_ADD_VNAME(dn(1), dn(1), 0,0, myThid)
128            DO j = 1,nlevels(listnum)  
129              IF (i .EQ. 1) THEN  C         The following three ztmp() loops should eventually be modified
130                ztmp(j) = rC(levs(j,listnum))  C         to reflect the fractional nature of levs(j,l) -- they should
131                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',  C         do something like:
132       &             'Dimensional coordinate value at the mid point',  C            ztmp(j) = rC(INT(FLOOR(levs(j,l))))
133       &             myThid)  C                      + ( rC(INT(FLOOR(levs(j,l))))
134              ELSEIF (i .EQ. 2) THEN  C                          + rC(INT(CEIL(levs(j,l)))) )
135                ztmp(j) = rF(levs(j,listnum) + 1)  C                        / ( levs(j,l) - FLOOR(levs(j,l)) )
136                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',  C         for averaged levels.
137       &             'Dimensional coordinate value at the upper point',            IF (i .EQ. 1) THEN
138       &             myThid)              DO j = 1,nlevels(listnum)
139              ELSEIF (i .EQ. 3) THEN                ztmp(j) = rC(NINT(levs(j,listnum)))
140                ztmp(j) = rF(levs(j,listnum))              ENDDO
141                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',              CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
142       &             'Dimensional coordinate value at the lower point',       &           'Dimensional coordinate value at the mid point',
143       &             myThid)       &           myThid)
144              ENDIF            ELSEIF (i .EQ. 2) THEN
145            ENDDO              DO j = 1,nlevels(listnum)
146                  ztmp(j) = rF(NINT(levs(j,listnum)) + 1)
147                ENDDO
148                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
149         &           'Dimensional coordinate value at the upper point',
150         &           myThid)
151              ELSEIF (i .EQ. 3) THEN
152                DO j = 1,nlevels(listnum)
153                  ztmp(j) = rF(NINT(levs(j,listnum)))
154                ENDDO
155                CALL MNC_CW_ADD_VATTR_TEXT(dn(1),'description',
156         &           'Dimensional coordinate value at the lower point',
157         &           myThid)
158              ENDIF
159            CALL MNC_CW_RS_W('D',diag_mnc_bn,0,0, dn(1), ztmp, myThid)            CALL MNC_CW_RS_W('D',diag_mnc_bn,0,0, dn(1), ztmp, myThid)
160            CALL MNC_CW_DEL_VNAME(dn(1), myThid)            CALL MNC_CW_DEL_VNAME(dn(1), myThid)
161            CALL MNC_CW_DEL_GNAME(dn(1), myThid)            CALL MNC_CW_DEL_GNAME(dn(1), myThid)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22