/[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.6 by dimitri, Wed Apr 30 07:04:08 2003 UTC revision 1.7 by dimitri, Fri May 23 20:19:16 2003 UTC
# Line 20  C     === Routine arguments === Line 20  C     === Routine arguments ===
20  C     myThid - Thread no. that called this routine.  C     myThid - Thread no. that called this routine.
21        _RL HEFF       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,3,nSx,nSy)        _RL HEFF       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,3,nSx,nSy)
22        _RL HEFFM      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,  nSx,nSy)        _RL HEFFM      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,  nSx,nSy)
23        _RL DIFFA      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL DIFFA      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,  nSx,nSy)
24        _RL DELTT        _RL DELTT
25        INTEGER myThid        INTEGER myThid
26  CEndOfInterface  CEndOfInterface
# Line 46  C     i,j,bi,bj - Loop counters Line 46  C     i,j,bi,bj - Loop counters
46  C NOW DO DIFFUSION WITH NUIT CONVERSION  C NOW DO DIFFUSION WITH NUIT CONVERSION
47          DO j=1,sNy          DO j=1,sNy
48           DO i=1,sNx           DO i=1,sNx
49            DELTXX1=DELTT*DIFFA(I,J)            DELTXX1=DELTT*DIFFA(I,J,bi,bj)
50            DELTYY1=DELTT*DIFFA(I,J)            DELTYY1=DELTT*DIFFA(I,J,bi,bj)
51            DELTXX=DELTXX1/(DXTICE(I,J,bi,bj)*DXTICE(I,J,bi,bj)            DELTXX=DELTXX1/(DXTICE(I,J,bi,bj)*DXTICE(I,J,bi,bj)
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)
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22