/[MITgcm]/MITgcm/pkg/seaice/diffus.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/diffus.F

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

revision 1.4 by dimitri, Sat Dec 28 10:11:11 2002 UTC revision 1.5 by dimitri, Tue Feb 18 05:33:55 2003 UTC
# Line 32  C     i,j,bi,bj - Loop counters Line 32  C     i,j,bi,bj - Loop counters
32    
33        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
34        _RL DELTXX1, DELTYY1, DELTXX, DELTYY        _RL DELTXX1, DELTYY1, DELTXX, DELTYY
35        _RL HEFF1  (1-OLx:sNx+OLx, 1-OLy:sNy+OLy)        _RL HEFF1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
36    
37        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
38         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
39    
40          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
41           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
42            HEFF1(I,J)=0.0 _d 0            HEFF1(I,J,bi,bj)=0.0 _d 0
43           ENDDO           ENDDO
44          ENDDO          ENDDO
45    
# Line 52  C NOW DO DIFFUSION WITH NUIT CONVERSION Line 52  C NOW DO DIFFUSION WITH NUIT CONVERSION
52       &          *CSTICE(I,J,bi,bj)*CSTICE(I,J,bi,bj))       &          *CSTICE(I,J,bi,bj)*CSTICE(I,J,bi,bj))
53            DELTYY=DELTYY1/(DYTICE(I,J,bi,bj)*DYTICE(I,J,bi,bj)            DELTYY=DELTYY1/(DYTICE(I,J,bi,bj)*DYTICE(I,J,bi,bj)
54       &          *CSTICE(I,J,bi,bj))       &          *CSTICE(I,J,bi,bj))
55            HEFF1(I,J)=DELTXX*((HEFF(I+1,J,3,bi,bj)-HEFF(I,J,3,bi,bj))            HEFF1(I,J,bi,bj)=DELTXX*((HEFF(I+1,J,3,bi,bj)
56       &     *HEFFM(I+1,J,bi,bj)       &     -HEFF(I,J,3,bi,bj))*HEFFM(I+1,J,bi,bj)
57       &     -(HEFF(I,J,3,bi,bj)-HEFF(I-1,J,3,bi,bj))*HEFFM(I-1,J,bi,bj))       &     -(HEFF(I,J,3,bi,bj)-HEFF(I-1,J,3,bi,bj))*HEFFM(I-1,J,bi,bj))
58       &     +DELTYY*((HEFF(I,J+1,3,bi,bj)-HEFF(I,J,3,bi,bj))       &     +DELTYY*((HEFF(I,J+1,3,bi,bj)-HEFF(I,J,3,bi,bj))
59       &     *CSUICE(I,J,bi,bj)*HEFFM(I,J+1,bi,bj)       &     *CSUICE(I,J,bi,bj)*HEFFM(I,J+1,bi,bj)
# Line 71  C NOW DO DIFFUSION WITH NUIT CONVERSION Line 71  C NOW DO DIFFUSION WITH NUIT CONVERSION
71         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
72          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
73           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
74            HEFF(I,J,3,bi,bj)=HEFF1(I,J)            HEFF(I,J,3,bi,bj)=HEFF1(I,J,bi,bj)
75           ENDDO           ENDDO
76          ENDDO          ENDDO
77         ENDDO         ENDDO
78        ENDDO        ENDDO
79    
80  #endif ALLOW_SEAICE  #endif /* ALLOW_SEAICE */
81    
82        RETURN        RETURN
83        END        END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22