/[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.26 by jmc, Tue Dec 30 02:14:05 2008 UTC revision 1.27 by jmc, Tue Feb 3 22:59:39 2009 UTC
# Line 22  C     INITIALIZE_FIXED Line 22  C     INITIALIZE_FIXED
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"
# Line 41  C               loaded in multi-threaded Line 43  C               loaded in multi-threaded
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
# Line 101  C     Write 3D geometry arrays Line 100  C     Write 3D geometry arrays
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    
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22