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

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

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

revision 1.49 by jmc, Tue Dec 8 00:10:26 2009 UTC revision 1.50 by jmc, Wed May 18 01:15:18 2011 UTC
# Line 28  C     === Global variables === Line 28  C     === Global variables ===
28  #include "GRID.h"  #include "GRID.h"
29  #include "SURFACE.h"  #include "SURFACE.h"
30  #include "CG2D.h"  #include "CG2D.h"
 #ifdef ALLOW_OBCS  
 #include "OBCS.h"  
 #endif  
31    
32  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
33  C     === Routine arguments ===  C     === Routine arguments ===
# Line 41  C     !LOCAL VARIABLES: Line 38  C     !LOCAL VARIABLES:
38  C     === Local variables ===  C     === Local variables ===
39  C     bi,bj  :: tile indices  C     bi,bj  :: tile indices
40  C     i,j,k  :: Loop counters  C     i,j,k  :: Loop counters
41  C     faceArea - Temporary used to hold cell face areas.  C     faceArea :: Temporary used to hold cell face areas.
42  C     myNorm - Work variable used in calculating normalisation factor  C     myNorm   :: Work variable used in calculating normalisation factor
 C     sumArea - Work variable used to compute the total Domain Area  
43        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
44        INTEGER bi, bj        INTEGER bi, bj
45        INTEGER i, j, k, ks        INTEGER i, j, k, ks
# Line 112  C  deep-model: *deepFacC (faceArea), /de Line 108  C  deep-model: *deepFacC (faceArea), /de
108            ENDDO            ENDDO
109           ENDDO           ENDDO
110          ENDDO          ENDDO
 #ifdef ALLOW_OBCS  
         IF (useOBCS) THEN  
          DO i=1,sNx  
           IF (OB_Jn(i,bi,bj).NE.0) aS2d(i, OB_Jn(i,bi,bj), bi,bj)=0.  
           IF (OB_Jn(i,bi,bj).NE.0) aS2d(i,OB_Jn(i,bi,bj)+1,bi,bj)=0.  
           IF (OB_Js(i,bi,bj).NE.0) aS2d(i,OB_Js(i,bi,bj)+1,bi,bj)=0.  
           IF (OB_Js(i,bi,bj).NE.0) aS2d(i, OB_Js(i,bi,bj), bi,bj)=0.  
          ENDDO  
          DO j=1,sNy  
           IF (OB_Ie(j,bi,bj).NE.0) aW2d(OB_Ie(j,bi,bj), j, bi,bj)=0.  
           IF (OB_Ie(j,bi,bj).NE.0) aW2d(OB_Ie(j,bi,bj)+1,j,bi,bj)=0.  
           IF (OB_Iw(j,bi,bj).NE.0) aW2d(OB_Iw(j,bi,bj)+1,j,bi,bj)=0.  
           IF (OB_Iw(j,bi,bj).NE.0) aW2d(OB_Iw(j,bi,bj), j, bi,bj)=0.  
          ENDDO  
         ENDIF  
 #endif  
111          DO j=1,sNy          DO j=1,sNy
112           DO i=1,sNx           DO i=1,sNx
113            myNorm = MAX(ABS(aW2d(i,j,bi,bj)),myNorm)  #ifdef ALLOW_OBCS
114            myNorm = MAX(ABS(aS2d(i,j,bi,bj)),myNorm)             aW2d(i,j,bi,bj) = aW2d(i,j,bi,bj)
115         &                   *maskInC(i,j,bi,bj)*maskInC(i-1,j,bi,bj)
116               aS2d(i,j,bi,bj) = aS2d(i,j,bi,bj)
117         &                   *maskInC(i,j,bi,bj)*maskInC(i,j-1,bi,bj)
118    #endif /* ALLOW_OBCS */
119               myNorm = MAX(ABS(aW2d(i,j,bi,bj)),myNorm)
120               myNorm = MAX(ABS(aS2d(i,j,bi,bj)),myNorm)
121           ENDDO           ENDDO
122          ENDDO          ENDDO
123         ENDDO         ENDDO
Legend:
Removed from v.1.49  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22