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

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

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

revision 1.33 by jmc, Tue Mar 16 00:08:27 2010 UTC revision 1.34 by jmc, Sun Jun 17 02:20:46 2012 UTC
# Line 40  C     !USES: Line 40  C     !USES:
40  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
41  C     myThid  :: my Thread Id number  C     myThid  :: my Thread Id number
42        INTEGER myThid        INTEGER myThid
 CEOP  
43    
44  C     === Functions ====  C     !FUNCTIONS:
45        LOGICAL  MASTER_CPU_IO        LOGICAL  MASTER_CPU_IO
46        EXTERNAL MASTER_CPU_IO        EXTERNAL MASTER_CPU_IO
47    
48  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
49    C     bi, bj  :: tile indices
50    C     i, j    :: Loop counters
51  C     msgBuf  :: Informational/error message buffer  C     msgBuf  :: Informational/error message buffer
52          INTEGER bi, bj
53          INTEGER i, j
54        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
55  #ifdef ALLOW_MNC  CEOP
       INTEGER i  
 #endif  
56    
57  C     load grid spacing (vector) from files  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
58    
59    C     Load grid spacing (vector) from files
60        CALL LOAD_GRID_SPACING( myThid )        CALL LOAD_GRID_SPACING( myThid )
61    
62  C     Set up vertical grid and coordinate system  C     Set up vertical grid and coordinate system
# Line 89  C     Set up horizontal grid and coordin Line 92  C     Set up horizontal grid and coordin
92          WRITE(msgBuf,'(2A)') 'S/R INI_GRID: ',          WRITE(msgBuf,'(2A)') 'S/R INI_GRID: ',
93       &       'No grid coordinate system has been selected'       &       'No grid coordinate system has been selected'
94          CALL PRINT_ERROR( msgBuf , myThid)          CALL PRINT_ERROR( msgBuf , myThid)
95            CALL ALL_PROC_DIE( 0 )
96          STOP 'ABNORMAL END: S/R INI_GRID'          STOP 'ABNORMAL END: S/R INI_GRID'
97          _END_MASTER(myThid)          _END_MASTER(myThid)
98        ENDIF        ENDIF
99    
100    C---  Calculate reciprocals grid lengths (formerly part of INI_MASKS_ETC)
101          DO bj = myByLo(myThid), myByHi(myThid)
102           DO bi = myBxLo(myThid), myBxHi(myThid)
103            DO j=1-OLy,sNy+OLy
104             DO i=1-OLx,sNx+OLx
105              IF ( dxG(i,j,bi,bj) .NE. 0. )
106         &    recip_dxG(i,j,bi,bj) = 1. _d 0/dxG(i,j,bi,bj)
107              IF ( dyG(i,j,bi,bj) .NE. 0. )
108         &    recip_dyG(i,j,bi,bj) = 1. _d 0/dyG(i,j,bi,bj)
109              IF ( dxC(i,j,bi,bj) .NE. 0. )
110         &    recip_dxC(i,j,bi,bj) = 1. _d 0/dxC(i,j,bi,bj)
111              IF ( dyC(i,j,bi,bj) .NE. 0. )
112         &    recip_dyC(i,j,bi,bj) = 1. _d 0/dyC(i,j,bi,bj)
113              IF ( dxF(i,j,bi,bj) .NE. 0. )
114         &    recip_dxF(i,j,bi,bj) = 1. _d 0/dxF(i,j,bi,bj)
115              IF ( dyF(i,j,bi,bj) .NE. 0. )
116         &    recip_dyF(i,j,bi,bj) = 1. _d 0/dyF(i,j,bi,bj)
117              IF ( dxV(i,j,bi,bj) .NE. 0. )
118         &    recip_dxV(i,j,bi,bj) = 1. _d 0/dxV(i,j,bi,bj)
119              IF ( dyU(i,j,bi,bj) .NE. 0. )
120         &    recip_dyU(i,j,bi,bj) = 1. _d 0/dyU(i,j,bi,bj)
121              IF ( rA (i,j,bi,bj) .NE. 0. )
122         &    recip_rA (i,j,bi,bj) = 1. _d 0/rA (i,j,bi,bj)
123              IF ( rAs(i,j,bi,bj) .NE. 0. )
124         &    recip_rAs(i,j,bi,bj) = 1. _d 0/rAs(i,j,bi,bj)
125              IF ( rAw(i,j,bi,bj) .NE. 0. )
126         &    recip_rAw(i,j,bi,bj) = 1. _d 0/rAw(i,j,bi,bj)
127              IF ( rAz(i,j,bi,bj) .NE. 0. )
128         &    recip_rAz(i,j,bi,bj) = 1. _d 0/rAz(i,j,bi,bj)
129             ENDDO
130            ENDDO
131           ENDDO
132          ENDDO
133    
134    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
135    
136  #ifdef ALLOW_MONITOR  #ifdef ALLOW_MONITOR
137        IF ( MASTER_CPU_IO(myThid) ) THEN        IF ( MASTER_CPU_IO(myThid) ) THEN
138  C--   only the master thread is allowed to switch On/Off mon_write_stdout  C--   only the master thread is allowed to switch On/Off mon_write_stdout
# Line 117  C     & mon_write_mnc (since it is the o Line 157  C     & mon_write_mnc (since it is the o
157    
158        ENDIF        ENDIF
159    
160  C     Print out statistics of each horizontal grid array (helps when  C     Print out statistics of each horizontal grid array (helps when debugging)
 C     debugging)  
161        CALL MON_PRINTSTATS_RS(1,xC,'XC',myThid)        CALL MON_PRINTSTATS_RS(1,xC,'XC',myThid)
162        CALL MON_PRINTSTATS_RS(1,xG,'XG',myThid)        CALL MON_PRINTSTATS_RS(1,xG,'XG',myThid)
163        CALL MON_PRINTSTATS_RS(1,dxC,'DXC',myThid)        CALL MON_PRINTSTATS_RS(1,dxC,'DXC',myThid)
Legend:
Removed from v.1.33  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22