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

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

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

revision 1.12 by jmc, Sun Aug 19 21:58:41 2007 UTC revision 1.13 by jmc, Mon Oct 29 18:15:29 2007 UTC
# Line 70  C---+----1----+----2----+----3----+----4 Line 70  C---+----1----+----2----+----3----+----4
70         DO bi=myBxLo(myThid), myBxHi(myThid)         DO bi=myBxLo(myThid), myBxHi(myThid)
71  C-    1rst bi,bj loop :  C-    1rst bi,bj loop :
72    
         IF (myIter.EQ.-1) THEN  
 C-- Initialise arrays :  
          DO j=1-Oly,sNy+Oly  
           DO i=1-Olx,sNx+Olx  
             rStarFacC(i,j,bi,bj) = 1.  
             rStarFacW(i,j,bi,bj) = 1.  
             rStarFacS(i,j,bi,bj) = 1.  
             rStarExpC(i,j,bi,bj) = 1.  
             rStarExpW(i,j,bi,bj) = 1.  
             rStarExpS(i,j,bi,bj) = 1.  
             rStarDhCDt(i,j,bi,bj) = 0.  
             rStarDhWDt(i,j,bi,bj) = 0.  
             rStarDhSDt(i,j,bi,bj) = 0.  
             PmEpR(i,j,bi,bj) = 0.  
           ENDDO  
          ENDDO  
          DO k=1,Nr  
           DO j=1-Oly,sNy+Oly  
            DO i=1-Olx,sNx+Olx  
             h0FacC(i,j,k,bi,bj) = _hFacC(i,j,k,bi,bj)  
             h0FacW(i,j,k,bi,bj) = _hFacW(i,j,k,bi,bj)  
             h0FacS(i,j,k,bi,bj) = _hFacS(i,j,k,bi,bj)  
            ENDDO  
           ENDDO  
          ENDDO  
         ELSE  
73  C-- copy rStarFacX -> rStarExpX  C-- copy rStarFacX -> rStarExpX
74           DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
75            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
76              rStarExpC(i,j,bi,bj) = rStarFacC(i,j,bi,bj)              rStarExpC(i,j,bi,bj) = rStarFacC(i,j,bi,bj)
77              rStarExpW(i,j,bi,bj) = rStarFacW(i,j,bi,bj)              rStarExpW(i,j,bi,bj) = rStarFacW(i,j,bi,bj)
78              rStarExpS(i,j,bi,bj) = rStarFacS(i,j,bi,bj)              rStarExpS(i,j,bi,bj) = rStarFacS(i,j,bi,bj)
79            ENDDO            ENDDO
80           ENDDO          ENDDO
         ENDIF  
81    
82  C-- Compute the new column thikness :  C-- Compute the new column thikness :
83          DO j=0,sNy+1          DO j=0,sNy+1
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22