/[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.22 by heimbach, Thu Nov 6 22:12:09 2003 UTC revision 1.23 by dimitri, Thu Nov 13 06:35:14 2003 UTC
# Line 100  c     to OLx=1, OLy=1. Line 100  c     to OLx=1, OLy=1.
100  #include "FFIELDS.h"  #include "FFIELDS.h"
101  #include "GRID.h"  #include "GRID.h"
102    
 #ifdef ALLOW_SEAICE  
 #include "SEAICE_EXTERNAL.h"  
 #endif /* ALLOW_SEAICE */  
   
103  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
104  #include "tamc.h"  #include "tamc.h"
105  #include "tamc_keys.h"  #include "tamc_keys.h"
# Line 372  cph) Line 368  cph)
368             ustar(i,j) = SQRT( tempVar2 )             ustar(i,j) = SQRT( tempVar2 )
369          endif          endif
370    
371  #ifdef ALLOW_SEAICE          bo(I,J) = - gravity *
         if ( .NOT. useSEAICE ) then  
 #endif  
            bo(I,J) = - gravity *  
      &       ( vddiff(I,J,1,1) * SurfaceTendencyT(i,j,bi,bj) +  
      &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)  
      &       ) *  
      &       drF(1) / work2(I,J)  
   
 #ifdef ALLOW_SEAICE  
         else  
            bo(I,J) = - gravity *  
372       &       ( vddiff(I,J,1,1) * (SurfaceTendencyT(i,j,bi,bj)+       &       ( vddiff(I,J,1,1) * (SurfaceTendencyT(i,j,bi,bj)+
373       &                            SurfaceTendencyTice(i,j,bi,bj)) +       &       SurfaceTendencyTice(i,j,bi,bj)) +
374       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)       &       vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj) ) *
      &       ) *  
375       &       drF(1) / work2(I,J)       &       drF(1) / work2(I,J)
         endif  
 #endif /* ALLOW_SEAICE */  
376    
377          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) *
378       &       recip_Cp*recip_rhoConst*recip_dRf(1) *       &       recip_Cp*recip_rhoConst*recip_dRf(1) *
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22