/[MITgcm]/MITgcm/pkg/kpp/kpp_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_calc.F

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

revision 1.9 by adcroft, Thu Sep 27 18:08:57 2001 UTC revision 1.10 by jmc, Thu Mar 7 14:27:47 2002 UTC
# Line 302  c     so that the subroutine "bldepth" w Line 302  c     so that the subroutine "bldepth" w
302  c------------------------------------------------------------------------  c------------------------------------------------------------------------
303  c     friction velocity, turbulent and radiative surface buoyancy forcing  c     friction velocity, turbulent and radiative surface buoyancy forcing
304  c     -------------------------------------------------------------------  c     -------------------------------------------------------------------
305  c     taux / rho = SurfaceTendencyU * delZ(1)                     (N/m^2)  c     taux / rho = SurfaceTendencyU * drF(1)                     (N/m^2)
306  c     tauy / rho = SurfaceTendencyV * delZ(1)                     (N/m^2)  c     tauy / rho = SurfaceTendencyV * drF(1)                     (N/m^2)
307  c     ustar = sqrt( sqrt( taux^2 + tauy^2 ) / rho )                 (m/s)  c     ustar = sqrt( sqrt( taux^2 + tauy^2 ) / rho )                (m/s)
308  c     bo    = - g * ( alpha*SurfaceTendencyT +  c     bo    = - g * ( alpha*SurfaceTendencyT +
309  c                     beta *SurfaceTendencyS ) * delZ(1) / rho  (m^2/s^3)  c                     beta *SurfaceTendencyS ) * drF(1) / rho  (m^2/s^3)
310  c     bosol = - g * alpha * Qsw * delZ(1) / rho                 (m^2/s^3)  c     bosol = - g * alpha * Qsw * drF(1) / rho                 (m^2/s^3)
311  c------------------------------------------------------------------------  c------------------------------------------------------------------------
312    
313  c initialize arrays to zero  c initialize arrays to zero
# Line 329  c initialize arrays to zero Line 329  c initialize arrays to zero
329       &   (SurfaceTendencyV(i,j,bi,bj) + SurfaceTendencyV(i,jp1,bi,bj)) *       &   (SurfaceTendencyV(i,j,bi,bj) + SurfaceTendencyV(i,jp1,bi,bj)) *
330       &   (SurfaceTendencyV(i,j,bi,bj) + SurfaceTendencyV(i,jp1,bi,bj))       &   (SurfaceTendencyV(i,j,bi,bj) + SurfaceTendencyV(i,jp1,bi,bj))
331          if ( tempVar1 .lt. (phepsi*phepsi) ) then          if ( tempVar1 .lt. (phepsi*phepsi) ) then
332             ustar(i,j) = SQRT( phepsi * p5 * delZ(1) )             ustar(i,j) = SQRT( phepsi * p5 * drF(1) )
333          else          else
334             tempVar2 =  SQRT( tempVar1 ) * p5 * delZ(1)             tempVar2 =  SQRT( tempVar1 ) * p5 * drF(1)
335             ustar(i,j) = SQRT( tempVar2 )             ustar(i,j) = SQRT( tempVar2 )
336          endif          endif
337          bo(I,J) = - gravity *          bo(I,J) = - gravity *
338       &       ( vddiff(I,J,1,1) * SurfaceTendencyT(i,j,bi,bj) +       &       ( vddiff(I,J,1,1) * SurfaceTendencyT(i,j,bi,bj) +
339       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)
340       &       ) *       &       ) *
341       &       delZ(1) / work2(I,J)       &       drF(1) / work2(I,J)
342          bosol(I,J) = gravity * vddiff(I,J,1,1) * Qsw(i,j,bi,bj) *          bosol(I,J) = gravity * vddiff(I,J,1,1) * Qsw(i,j,bi,bj) *
343       &       recip_Cp*recip_rhoNil*recip_dRf(1) *       &       recip_Cp*recip_rhoNil*recip_dRf(1) *
344       &       delZ(1) / work2(I,J)       &       drF(1) / work2(I,J)
345         END DO         END DO
346        END DO        END DO
347    
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22