22 |
C | |
C | |
23 |
C |-- WRITE_GRID |
C |-- WRITE_GRID |
24 |
|
|
25 |
|
C !USES: |
26 |
IMPLICIT NONE |
IMPLICIT NONE |
27 |
|
|
28 |
#include "SIZE.h" |
#include "SIZE.h" |
29 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
30 |
#include "PARAMS.h" |
#include "PARAMS.h" |
43 |
C thread => needs to be in common to be shared by all threads |
C thread => needs to be in common to be shared by all threads |
44 |
COMMON / LOCAL_WRITE_GRID / tmpfld |
COMMON / LOCAL_WRITE_GRID / tmpfld |
45 |
_RS tmpfld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS tmpfld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
INTEGER i,j,bi,bj |
|
|
LOGICAL writegrid_mdsio |
|
|
#ifdef ALLOW_MDSIO |
|
|
INTEGER k |
|
46 |
_RS tmpVar(Nr+1) |
_RS tmpVar(Nr+1) |
47 |
#endif |
INTEGER i,j,k,bi,bj |
48 |
|
LOGICAL writegrid_mdsio |
49 |
#ifdef ALLOW_MNC |
#ifdef ALLOW_MNC |
50 |
CHARACTER*(4) bfname |
CHARACTER*(4) bfname |
51 |
#endif |
#endif |
100 |
IF (buoyancyRelation .EQ. 'ATMOSPHERIC') |
IF (buoyancyRelation .EQ. 'ATMOSPHERIC') |
101 |
& CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,-1,myThid) |
& CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,-1,myThid) |
102 |
|
|
|
#ifdef ALLOW_MDSIO |
|
103 |
C Write vertical grid arrays |
C Write vertical grid arrays |
104 |
CALL MDS_WRITEVEC_LOC('RC', writeBinaryPrec, 'RS', Nr, rC, |
CALL WRITE_GLVEC_RS( 'RC', ' ', rC, Nr, -1, myThid ) |
105 |
& 0,0, 1, -1, myThid) |
CALL WRITE_GLVEC_RS( 'RF', ' ', rF, 1+Nr,-1, myThid ) |
106 |
CALL MDS_WRITEVEC_LOC('RF', writeBinaryPrec, 'RS',1+Nr,rF, |
CALL WRITE_GLVEC_RS( 'DRC',' ', drC, Nr, -1, myThid ) |
107 |
& 0,0, 1, -1, myThid) |
CALL WRITE_GLVEC_RS( 'DRF',' ', drF, Nr, -1, myThid ) |
|
CALL MDS_WRITEVEC_LOC('DRC',writeBinaryPrec, 'RS', Nr, drC, |
|
|
& 0,0, 1, -1, myThid) |
|
|
CALL MDS_WRITEVEC_LOC('DRF',writeBinaryPrec, 'RS', Nr, drF, |
|
|
& 0,0, 1, -1, myThid) |
|
108 |
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) |
109 |
c IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN |
c IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN |
110 |
DO k=1,Nr+1 |
DO k=1,Nr+1 |
111 |
tmpVar(k) = phiRef(2*k-1) |
tmpVar(k) = phiRef(2*k-1) |
112 |
ENDDO |
ENDDO |
113 |
CALL MDS_WRITEVEC_LOC('PHrefF', writeBinaryPrec, 'RS',1+Nr, |
CALL WRITE_GLVEC_RS( 'PHrefF',' ',tmpVar,1+Nr,-1, myThid ) |
|
& tmpVar, 0, 0, 1, -1, myThid) |
|
114 |
DO k=1,Nr |
DO k=1,Nr |
115 |
tmpVar(k) = phiRef(2*k) |
tmpVar(k) = phiRef(2*k) |
116 |
ENDDO |
ENDDO |
117 |
CALL MDS_WRITEVEC_LOC('PHrefC', writeBinaryPrec, 'RS', Nr, |
CALL WRITE_GLVEC_RS( 'PHrefC',' ',tmpVar, Nr, -1, myThid ) |
|
& tmpVar, 0, 0, 1, -1, myThid) |
|
118 |
c ENDIF |
c ENDIF |
|
#endif /* ALLOW_MDSIO */ |
|
119 |
|
|
120 |
ENDIF |
ENDIF |
121 |
|
|