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

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

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

revision 1.5 by dimitri, Tue Feb 18 05:33:55 2003 UTC revision 1.6 by dimitri, Wed Apr 30 07:04:08 2003 UTC
# Line 6  CStartOfInterface Line 6  CStartOfInterface
6        SUBROUTINE seaice_model( myTime, myIter, myThid )        SUBROUTINE seaice_model( myTime, myIter, myThid )
7  C     /==========================================================\  C     /==========================================================\
8  C     | SUBROUTINE SEAICE_MODEL                                  |  C     | SUBROUTINE SEAICE_MODEL                                  |
9  C     | o Time stepping of sea ice model.                        |  C     | o Time stepping of a dynamic thermodynamic sea ice model |
10    C     |  Dynamics solver: Zhang/Hibler, JGR, 102, 8691-8702, 1997|
11    C     |    and          Zhang/Rothrock, JGR, 105, 3325-3338, 2000|
12    C     |  Thermodynamics:        Hibler, MWR, 108, 1943-1973, 1980|
13    C     |  Rheology:              Hibler, JPO,   9,  815- 846, 1979|
14    C     |  Written by Jinlun Zhang (PSC/UW), Feb-2001              |
15    C     |  Incorporated into MITgcm by Dimitris Menemenlis (JPL)   |
16  C     |==========================================================|  C     |==========================================================|
17  C     \==========================================================/  C     \==========================================================/
18        IMPLICIT NONE        IMPLICIT NONE
# Line 67  C--   Third level model velocity is used Line 73  C--   Third level model velocity is used
73         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
74          DO j=0,sNy+1          DO j=0,sNy+1
75           DO i=0,sNx+1           DO i=0,sNx+1
76            GWATX(I,J,bi,bj)=HALF*(uVel(i+1,j,3,bi,bj)            GWATX(I,J,bi,bj)=HALF*(uVel(i,j,3,bi,bj)
77       &                         +uVel(i+1,j+1,3,bi,bj))       &                         +uVel(i,j-1,3,bi,bj))
78            GWATY(I,J,bi,bj)=HALF*(vVel(i,j+1,3,bi,bj)            GWATY(I,J,bi,bj)=HALF*(vVel(i,j,3,bi,bj)
79       &                         +vVel(i+1,j+1,3,bi,bj))       &                         +vVel(i-1,j,3,bi,bj))
80  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
81  c          write(*,'(2i4,2i2,f7.1,7f12.3)')  c          write(*,'(2i4,2i2,f7.1,7f12.3)')
82  c     &     ,i,j,bi,bj,UVM(I,J,bi,bj)  c     &     ,i,j,bi,bj,UVM(I,J,bi,bj)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22