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 */ |