/[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.21 by jmc, Tue Oct 28 22:23:25 2003 UTC revision 1.22 by heimbach, Thu Nov 6 22:12:09 2003 UTC
# Line 372  cph) Line 372  cph)
372             ustar(i,j) = SQRT( tempVar2 )             ustar(i,j) = SQRT( tempVar2 )
373          endif          endif
374    
375          if ( .NOT. useSEAICE )  #ifdef ALLOW_SEAICE
376       &       bo(I,J) = - gravity *          if ( .NOT. useSEAICE ) then
377    #endif
378               bo(I,J) = - gravity *
379       &       ( vddiff(I,J,1,1) * SurfaceTendencyT(i,j,bi,bj) +       &       ( vddiff(I,J,1,1) * SurfaceTendencyT(i,j,bi,bj) +
380       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)
381       &       ) *       &       ) *
382       &       drF(1) / work2(I,J)       &       drF(1) / work2(I,J)
383    
384  #ifdef ALLOW_SEAICE  #ifdef ALLOW_SEAICE
385          if ( useSEAICE )          else
386       &       bo(I,J) = - gravity *             bo(I,J) = - gravity *
387       &       ( vddiff(I,J,1,1) * (SurfaceTendencyT(i,j,bi,bj)+       &       ( vddiff(I,J,1,1) * (SurfaceTendencyT(i,j,bi,bj)+
388       &                            SurfaceTendencyTice(i,j,bi,bj)) +       &                            SurfaceTendencyTice(i,j,bi,bj)) +
389       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)
390       &       ) *       &       ) *
391       &       drF(1) / work2(I,J)       &       drF(1) / work2(I,J)
392            endif
393  #endif /* ALLOW_SEAICE */  #endif /* ALLOW_SEAICE */
394    
395          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) *
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22