/[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.34 by gforget, Tue Feb 28 00:35:24 2012 UTC revision 1.35 by jmc, Fri Jun 1 03:19:53 2012 UTC
# Line 58  C     model.  This may differ from what Line 58  C     model.  This may differ from what
58  C     partial cells, etc.  C     partial cells, etc.
59        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
60         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
61           DO j=1-Oly,sNy+Oly           DO j=1-OLy,sNy+OLy
62            DO i=1-Olx,sNx+Olx            DO i=1-OLx,sNx+OLx
63  C     Total fluid column thickness (r_unit) :  C     Total fluid column thickness (r_unit) :
64              tmpfld(i,j,bi,bj) = Ro_surf(i,j,bi,bj) - R_low(i,j,bi,bj)              tmpfld(i,j,bi,bj) = Ro_surf(i,j,bi,bj) - R_low(i,j,bi,bj)
65            ENDDO            ENDDO
# Line 96  C     Write horizontal grid arrays Line 96  C     Write horizontal grid arrays
96          ENDIF          ENDIF
97  C     Write 3D geometry arrays  C     Write 3D geometry arrays
98          CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,-1,myThid)          CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,-1,myThid)
99            IF ( debugLevel.GE.debLevC ) THEN
100             CALL WRITE_FLD_XY_RS( 'rLowC', ' ', R_low, -1,myThid)
101             CALL WRITE_FLD_XY_RS( 'rLowW', ' ', rLowW, -1,myThid)
102             CALL WRITE_FLD_XY_RS( 'rLowS', ' ', rLowS, -1,myThid)
103             CALL WRITE_FLD_XY_RS( 'rSurfC',' ',Ro_surf,-1,myThid)
104             CALL WRITE_FLD_XY_RS( 'rSurfW',' ',rSurfW, -1,myThid)
105             CALL WRITE_FLD_XY_RS( 'rSurfS',' ',rSurfS, -1,myThid)
106            ENDIF
107          CALL WRITE_FLD_XYZ_RS( 'hFacC',' ',hFacC, 0,myThid)          CALL WRITE_FLD_XYZ_RS( 'hFacC',' ',hFacC, 0,myThid)
108          CALL WRITE_FLD_XYZ_RS( 'hFacW',' ',hFacW, 0,myThid)          CALL WRITE_FLD_XYZ_RS( 'hFacW',' ',hFacW, 0,myThid)
109          CALL WRITE_FLD_XYZ_RS( 'hFacS',' ',hFacS, 0,myThid)          CALL WRITE_FLD_XYZ_RS( 'hFacS',' ',hFacS, 0,myThid)
110          IF (buoyancyRelation .EQ. 'ATMOSPHERIC')          IF ( fluidIsAir )
111       &     CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,-1,myThid)       &    CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,-1,myThid)
112          IF ( useOBCS ) THEN          IF ( useOBCS ) THEN
113            CALL WRITE_FLD_XY_RS( 'maskInC',' ', maskInC, -1, myThid )            CALL WRITE_FLD_XY_RS( 'maskInC',' ', maskInC, -1, myThid )
114            CALL WRITE_FLD_XY_RS( 'maskInW',' ', maskInW, -1, myThid )            CALL WRITE_FLD_XY_RS( 'maskInW',' ', maskInW, -1, myThid )
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22