/[MITgcm]/MITgcm/model/src/write_grid.F
ViewVC logotype

Diff of /MITgcm/model/src/write_grid.F

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

revision 1.25 by mlosch, Sun Oct 26 11:09:46 2008 UTC revision 1.26 by jmc, Tue Dec 30 02:14:05 2008 UTC
# Line 103  C     Write 3D geometry arrays Line 103  C     Write 3D geometry arrays
103    
104  #ifdef ALLOW_MDSIO  #ifdef ALLOW_MDSIO
105  C     Write vertical grid arrays  C     Write vertical grid arrays
106          CALL MDSWRITEVEC_LOC_RS('RC', writeBinaryPrec, Nr, rC,          CALL MDS_WRITEVEC_LOC('RC', writeBinaryPrec, 'RS', Nr, rC,
107       &                                0,0, 1, -1, myThid)       &                                0,0, 1, -1, myThid)
108          CALL MDSWRITEVEC_LOC_RS('RF', writeBinaryPrec,1+Nr, rF,          CALL MDS_WRITEVEC_LOC('RF', writeBinaryPrec, 'RS',1+Nr,rF,
109       &                                0,0, 1, -1, myThid)       &                                0,0, 1, -1, myThid)
110          CALL MDSWRITEVEC_LOC_RS('DRC',writeBinaryPrec, Nr, drC,          CALL MDS_WRITEVEC_LOC('DRC',writeBinaryPrec, 'RS', Nr, drC,
111       &                                0,0, 1, -1, myThid)       &                                0,0, 1, -1, myThid)
112          CALL MDSWRITEVEC_LOC_RS('DRF',writeBinaryPrec, Nr, drF,          CALL MDS_WRITEVEC_LOC('DRF',writeBinaryPrec, 'RS', Nr, drF,
113       &                                0,0, 1, -1, myThid)       &                                0,0, 1, -1, myThid)
114  C--  comment out the "if": also usefull for the Ocean (to get the full Presure)  C--  comment out the "if": also usefull for the Ocean (to get the full Presure)
115  c       IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN  c       IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN
116          DO k=1,Nr+1          DO k=1,Nr+1
117            tmpVar(k) = phiRef(2*k-1)            tmpVar(k) = phiRef(2*k-1)
118          ENDDO          ENDDO
119          CALL MDSWRITEVEC_LOC_RS('PHrefF', writeBinaryPrec, 1+Nr,          CALL MDS_WRITEVEC_LOC('PHrefF', writeBinaryPrec, 'RS',1+Nr,
120       &                            tmpVar, 0, 0, 1, -1, myThid)       &                            tmpVar, 0, 0, 1, -1, myThid)
121          DO k=1,Nr          DO k=1,Nr
122            tmpVar(k) = phiRef(2*k)            tmpVar(k) = phiRef(2*k)
123          ENDDO          ENDDO
124          CALL MDSWRITEVEC_LOC_RS('PHrefC', writeBinaryPrec, Nr,          CALL MDS_WRITEVEC_LOC('PHrefC', writeBinaryPrec, 'RS', Nr,
125       &                            tmpVar, 0, 0, 1, -1, myThid)       &                            tmpVar, 0, 0, 1, -1, myThid)
126  c       ENDIF  c       ENDIF
127  #endif /* ALLOW_MDSIO */  #endif /* ALLOW_MDSIO */
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22