/[MITgcm]/MITgcm_contrib/dcarroll/iceberg/code/shelfice_thermodynamics.F
ViewVC logotype

Diff of /MITgcm_contrib/dcarroll/iceberg/code/shelfice_thermodynamics.F

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

revision 1.2 by dcarroll, Thu Aug 29 19:47:27 2019 UTC revision 1.3 by dcarroll, Tue Sep 10 05:33:18 2019 UTC
# Line 85  C     iceFrontWidth    :: the width of t Line 85  C     iceFrontWidth    :: the width of t
85        _RL tLoc        _RL tLoc
86        _RL sLoc        _RL sLoc
87        _RL pLoc        _RL pLoc
88          _RL uLoc
89  #ifndef SHI_USTAR_WETPOINT        _RL vLoc
90        _RL uLoc(1-olx:snx+olx,1-oly:sny+oly)        _RL wLoc
91        _RL vLoc(1-olx:snx+olx,1-oly:sny+oly)        _RL speedLoc
92  #endif  
93        _RL velSq(1-olx:snx+olx,1-oly:sny+oly)  C#ifndef SHI_USTAR_WETPOINT
94    C      _RL uLoc(1-olx:snx+olx,1-oly:sny+oly)
95    C      _RL vLoc(1-olx:snx+olx,1-oly:sny+oly)
96    C#endif
97    C      _RL velSq(1-olx:snx+olx,1-oly:sny+oly)
98                
99        _RL freshWaterFlux        _RL freshWaterFlux
100                
# Line 134  C     (Holland and Jenkins, 1999, eq.21) Line 138  C     (Holland and Jenkins, 1999, eq.21)
138          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
139            DO J = 1-OLy,sNy+OLy            DO J = 1-OLy,sNy+OLy
140              DO I = 1-OLx,sNx+OLx              DO I = 1-OLx,sNx+OLx
141                shelfIceHeatFlux      (I,J,bi,bj) = 0. _d 0               shelfIceHeatFlux      (I,J,bi,bj) = 0. _d 0
142                shelfIceFreshWaterFlux(I,J,bi,bj) = 0. _d 0               shelfIceFreshWaterFlux(I,J,bi,bj) = 0. _d 0
143                SHIICFHeatFlux      (I,J,bi,bj) = 0. _d 0               SHIICFHeatFlux      (I,J,bi,bj) = 0. _d 0
144                SHIICFFreshWaterFlux(I,J,bi,bj) = 0. _d 0               SHIICFFreshWaterFlux(I,J,bi,bj) = 0. _d 0
145                shelficeForcingT      (I,J,bi,bj) = 0. _d 0               shelficeForcingT      (I,J,bi,bj) = 0. _d 0
146                shelficeForcingS      (I,J,bi,bj) = 0. _d 0               shelficeForcingS      (I,J,bi,bj) = 0. _d 0
147                shelficeForcingTR     (I,J,bi,bj) = 0. _d 0               shelficeForcingTR     (I,J,bi,bj) = 0. _d 0
148  #ifndef ALLOW_shiTransCoeff_3d  
               shiTransCoeffS(I,J,bi,bj) = 5.05 _d -3 *  
      &              shiTransCoeffT(I,J,bi,bj)  
 #endif  
149                DO K = 1, NR                DO K = 1, NR
150                   iceFrontHeatFlux(I,J,K,bi,bj)       = 0. _d 0
151                   iceFrontFreshWaterFlux(I,J,K,bi,bj) = 0. _d 0
152                   iceFrontForcingT(I,J,K,bi,bj)       = 0. _d 0
153                   iceFrontForcingS(I,J,K,bi,bj)       = 0. _d 0
154                   iceFrontForcingTR(I,J,K,bi,bj)      = 0. _d 0
155                  ENDDO /* K */
156                  
157                ENDDO /* I */
158              ENDDO /* J */
159    
160    C 3-D velocity-dependent turbulent transfer coefficients
161    
162              DO J = 1-OLy,sNy+OLy
163                DO I = 1-OLx,sNx+OLx
164                  DO K = 1, NR
165    
166    C make local copies of velocity
167                   uLoc = 0.5*(uVel(I,J,K,bi,bj)+uVel(I+1,J,K,bi,bj))
168                   vLoc = 0.5*(vVel(I,J,K,bi,bj)+vVel(I,J+1,K,bi,bj))
169                   wLoc = 0.5*(wVel(I,J,K,bi,bj)+wVel(I,J,K+1,bi,bj))
170                   speedLoc = SQRT(uLoc*uLoc + vLoc*vLoc + wLoc*wLoc)
171    
172  #ifdef ALLOW_shiTransCoeff_3d  #ifdef ALLOW_shiTransCoeff_3d
173                  shiTransCoeffS(I,J,K,bi,bj) = 5.05 _d -3 *                 shiTransCoeffT(I,J,K,bi,bj) =  1.1 _d -3 *
174       &                shiTransCoeffT(I,J,K,bi,bj)       &          abs(speedLoc) + 1. _d -4
175    
176                   shiTransCoeffS(I,J,K,bi,bj) = 3.1 _d -5 *
177         &          abs(speedLoc) + 5.05 _d -7
178  #endif  #endif
                 iceFrontHeatFlux(I,J,K,bi,bj)       = 0. _d 0  
                 iceFrontFreshWaterFlux(I,J,K,bi,bj) = 0. _d 0  
                 iceFrontForcingT(I,J,K,bi,bj)       = 0. _d 0  
                 iceFrontForcingS(I,J,K,bi,bj)       = 0. _d 0  
                 iceFrontForcingTR(I,J,K,bi,bj)      = 0. _d 0  
179                ENDDO /* K */                ENDDO /* K */
                 
180              ENDDO /* I */              ENDDO /* I */
181            ENDDO /* J */            ENDDO /* J */
182    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22