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

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

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

revision 1.12 by mlosch, Wed Apr 18 18:06:52 2007 UTC revision 1.13 by mlosch, Fri Apr 20 18:29:58 2007 UTC
# Line 41  C     i,j,bi,bj - Loop counters Line 41  C     i,j,bi,bj - Loop counters
41        _RL  fuIceLoc, fvIceLoc, FX, FY        _RL  fuIceLoc, fvIceLoc, FX, FY
42        _RL  areaW, areaS        _RL  areaW, areaS
43    
44          _RL e11         (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
45          _RL e22         (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
46          _RL e12         (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
47        _RL press       (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)        _RL press       (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
48        _RL etaPlusZeta (1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL etaPlusZeta (1-Olx:sNx+Olx,1-Oly:sNy+Oly)
49        _RL zetaMinusEta(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL zetaMinusEta(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
# Line 82  C     use an intergral over ice and ocea Line 85  C     use an intergral over ice and ocea
85  C     surface stresses on ocean following Hibler and Bryan (1987, JPO)  C     surface stresses on ocean following Hibler and Bryan (1987, JPO)
86  C      C    
87  C     recompute viscosities from updated ice velocities  C     recompute viscosities from updated ice velocities
88         CALL SEAICE_CALC_VISCOSITIES(        CALL SEAICE_CALC_STRAINRATES(
89       I      uIce(1-Olx,1-Oly,1,1,1), vIce(1-Olx,1-Oly,1,1,1),       I     uIce(1-Olx,1-Oly,1,1,1), vIce(1-Olx,1-Oly,1,1,1),
90       I      zMin, zMax, hEffM, press0,       O     e11, e22, e12,
91       O      eta, zeta, press,       I     myThid )
92  #ifdef SEAICE_ALLOW_EVP  
93       O      seaice_div, seaice_tension, seaice_shear,        CALL SEAICE_CALC_VISCOSITIES(
94  #endif /* SEAICE_ALLOW_EVP */       I     e11, e22, e12, zMin, zMax, hEffM, press0,
95       I      myThid )       O     eta, zeta, press,
96         I     myThid )
97  C     re-compute internal stresses with updated ice velocities  C     re-compute internal stresses with updated ice velocities
98         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
99          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22