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

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

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

revision 1.18 by mlosch, Thu Mar 2 16:09:47 2006 UTC revision 1.19 by mlosch, Fri Mar 3 07:52:01 2006 UTC
# Line 18  C     === Global variables === Line 18  C     === Global variables ===
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19  #include "PARAMS.h"  #include "PARAMS.h"
20  #include "GRID.h"  #include "GRID.h"
21    #include "DYNVARS.h"
22  #include "FFIELDS.h"  #include "FFIELDS.h"
23  #include "SEAICE.h"  #include "SEAICE.h"
24  CML#include "SEAICE_GRID.h"  CML#include "SEAICE_GRID.h"
# Line 66  C--   FIRST SET UP BASIC CONSTANTS Line 67  C--   FIRST SET UP BASIC CONSTANTS
67        PSTAR=SEAICE_strength        PSTAR=SEAICE_strength
68    
69  C--   introduce turning angle (default is zero)  C--   introduce turning angle (default is zero)
70        SINWIN=SIN(SEAICE_airTurnAngle)        SINWIN=SIN(SEAICE_airTurnAngle*deg2rad)
71        COSWIN=COS(SEAICE_airTurnAngle)        COSWIN=COS(SEAICE_airTurnAngle*deg2rad)
72        SINWAT=SIN(SEAICE_waterTurnAngle)        SINWAT=SIN(SEAICE_waterTurnAngle*deg2rad)
73        COSWAT=COS(SEAICE_waterTurnAngle)        COSWAT=COS(SEAICE_waterTurnAngle*deg2rad)
74    
75    C--   Compute proxy for geostrophic velocity,
76          DO bj=myByLo(myThid),myByHi(myThid)
77           DO bi=myBxLo(myThid),myBxHi(myThid)
78            DO j=0,sNy+1
79             DO i=0,sNx+1
80              GWATX(I,J,bi,bj)=HALF*(uVel(i,j,KGEO(I,J,bi,bj),bi,bj)
81         &                         +uVel(i,j-1,KGEO(I,J,bi,bj),bi,bj))
82              GWATY(I,J,bi,bj)=HALF*(vVel(i,j,KGEO(I,J,bi,bj),bi,bj)
83         &                         +vVel(i-1,j,KGEO(I,J,bi,bj),bi,bj))
84    #ifdef SEAICE_DEBUG
85    c          write(*,'(2i4,2i2,f7.1,7f12.3)')
86    c     &     ,i,j,bi,bj,UVM(I,J,bi,bj)
87    c     &     ,GWATX(I,J,bi,bj),GWATY(I,J,bi,bj)
88    c     &     ,uVel(i+1,j,3,bi,bj),uVel(i+1,j+1,3,bi,bj)
89    c     &     ,vVel(i,j+1,3,bi,bj),vVel(i+1,j+1,3,bi,bj)
90    #endif
91             ENDDO
92            ENDDO
93           ENDDO
94          ENDDO
95    
96  C--   NOW SET UP MASS PER UNIT AREA AND CORIOLIS TERM  C--   NOW SET UP MASS PER UNIT AREA AND CORIOLIS TERM
97        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22