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

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

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

revision 1.21 by dimitri, Fri Apr 9 06:08:08 2004 UTC revision 1.22 by mlosch, Wed Apr 28 12:00:53 2004 UTC
# Line 94  C--   Initialize grid info Line 94  C--   Initialize grid info
94           DO I=1-OLx,sNx+OLx           DO I=1-OLx,sNx+OLx
95            CSTICE(i,j,bi,bj) =cos(yC(I,J,bi,bj)*deg2rad)            CSTICE(i,j,bi,bj) =cos(yC(I,J,bi,bj)*deg2rad)
96            CSUICE(i,j,bi,bj) =cos(yG(I,J,bi,bj)*deg2rad)            CSUICE(i,j,bi,bj) =cos(yG(I,J,bi,bj)*deg2rad)
97    CML(
98    C     inverses of CSTICE and CSUICE. Let's hope we are never
99    C     at the poles
100              IF ( CSTICE(I,J,bi,bj) .ne. 0. _d 0 ) THEN
101               RECIP_CSTICE(I,J,bi,bj) = 1./CSTICE(I,J,bi,bj)
102              ELSE  
103               RECIP_CSTICE(I,J,bi,bj) =0. _d 0
104              ENDIF
105              IF ( CSUICE(I,J,bi,bj) .ne. 0. _d 0 ) THEN
106               RECIP_CSUICE(I,J,bi,bj) = 1./CSUICE(I,J,bi,bj)
107              ELSE  
108               RECIP_CSUICE(I,J,bi,bj) =0. _d 0
109              ENDIF
110    CML)
111            SINEICE(i,j,bi,bj)=sin(yC(I,J,bi,bj)*deg2rad)            SINEICE(i,j,bi,bj)=sin(yC(I,J,bi,bj)*deg2rad)
112            TNGTICE(i,j,bi,bj)=SINEICE(i,j,bi,bj)/CSTICE(i,j,bi,bj)            TNGTICE(i,j,bi,bj)=SINEICE(i,j,bi,bj)*RECIP_CSTICE(i,j,bi,bj)
113            SINEICE(i,j,bi,bj)=sin(yG(I,J,bi,bj)*deg2rad)            SINEICE(i,j,bi,bj)=sin(yG(I,J,bi,bj)*deg2rad)
114            TNGICE(i,j,bi,bj) =SINEICE(i,j,bi,bj)/CSUICE(i,j,bi,bj)            TNGICE(i,j,bi,bj) =SINEICE(i,j,bi,bj)*RECIP_CSUICE(i,j,bi,bj)
115            DXTICE(i,j,bi,bj)=dxF(i,j,bi,bj)/CSTICE(i,j,bi,bj)            DXTICE(i,j,bi,bj)=dxF(i,j,bi,bj)*RECIP_CSTICE(i,j,bi,bj)
116            DXUICE(i,j,bi,bj)=dxV(i,j,bi,bj)/CSUICE(i,j,bi,bj)            DXUICE(i,j,bi,bj)=dxV(i,j,bi,bj)*RECIP_CSUICE(i,j,bi,bj)
117            DYTICE(i,j,bi,bj)=dyF(i,j,bi,bj)            DYTICE(i,j,bi,bj)=dyF(i,j,bi,bj)
118            DYUICE(i,j,bi,bj)=dyU(i,j,bi,bj)            DYUICE(i,j,bi,bj)=dyU(i,j,bi,bj)
119           ENDDO           ENDDO
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22