/[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.39 by jmc, Mon Apr 30 13:49:40 2007 UTC revision 1.40 by mlosch, Tue May 1 04:09:25 2007 UTC
# Line 251  c--------------------------------------- Line 251  c---------------------------------------
251    
252        CALL TIMER_START('STATEKPP      [KPP_CALC]', myThid)        CALL TIMER_START('STATEKPP      [KPP_CALC]', myThid)
253        CALL STATEKPP(        CALL STATEKPP(
254       I       ikppkey, bi, bj, myThid       O     work2, dbloc, Ritop,
255       O     , work2, dbloc, Ritop       O     TTALPHA, SSBETA,
256       O     , TTALPHA, SSBETA       I     ikppkey, bi, bj, myThid )
      &     )  
257        CALL TIMER_STOP ('STATEKPP      [KPP_CALC]', myThid)        CALL TIMER_STOP ('STATEKPP      [KPP_CALC]', myThid)
258    
259        DO k = 1, Nr        DO k = 1, Nr
# Line 274  c     levels therefore k+1 mask must be Line 273  c     levels therefore k+1 mask must be
273        DO k = 1, Nr-1        DO k = 1, Nr-1
274           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
275       I        k+1, bi, bj,       I        k+1, bi, bj,
276       U        ghat (1-OLx,1-OLy,k) )       U        ghat (1-OLx,1-OLy,k),
277         I        myThid )
278        ENDDO        ENDDO
279    
280  #endif /* KPP_SMOOTH_DBLOC */  #endif /* KPP_SMOOTH_DBLOC */
# Line 283  c     levels therefore k+1 mask must be Line 283  c     levels therefore k+1 mask must be
283  c     horizontally smooth density related quantities with 121 filters  c     horizontally smooth density related quantities with 121 filters
284        CALL SMOOTH_HORIZ (        CALL SMOOTH_HORIZ (
285       I     1, bi, bj,       I     1, bi, bj,
286       U     work2 )       U     work2,
287         I     myThid )
288        DO k = 1, Nr        DO k = 1, Nr
289           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
290       I        k+1, bi, bj,       I        k+1, bi, bj,
291       U        dbloc (1-OLx,1-OLy,k) )       U        dbloc (1-OLx,1-OLy,k),
292         I        myThid )
293           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
294       I        k, bi, bj,       I        k, bi, bj,
295       U        Ritop (1-OLx,1-OLy,k)  )       U        Ritop (1-OLx,1-OLy,k),
296         I        myThid )
297           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
298       I        k, bi, bj,       I        k, bi, bj,
299       U        TTALPHA(1-OLx,1-OLy,k) )       U        TTALPHA(1-OLx,1-OLy,k),
300         I        myThid )
301           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
302       I        k, bi, bj,       I        k, bi, bj,
303       U        SSBETA(1-OLx,1-OLy,k) )       U        SSBETA(1-OLx,1-OLy,k),
304         I        myThid )
305        ENDDO        ENDDO
306  #endif /* KPP_SMOOTH_DENS */  #endif /* KPP_SMOOTH_DENS */
307    
# Line 711  c     matching diffusivities at bottom o Line 716  c     matching diffusivities at bottom o
716        END DO        END DO
717        CALL TIMER_START('KPPMIX [KPP_CALC]', myThid)        CALL TIMER_START('KPPMIX [KPP_CALC]', myThid)
718        CALL KPPMIX (        CALL KPPMIX (
719       I       mytime, mythid       I       work1, shsq, dVsq, ustar
720       I     , work1, shsq, dVsq, ustar       I     , maskC(1-Olx,1-Oly,1,bi,bj)
721       I     , bo, bosol, dbloc, Ritop, work2       I     , bo, bosol, dbloc, Ritop, work2
722       I     , KPPdiffKzS(1-Olx,1-Oly,1,bi,bj)       I     , KPPdiffKzS(1-Olx,1-Oly,1,bi,bj)
723       I     , KPPdiffKzT(1-Olx,1-Oly,1,bi,bj)       I     , KPPdiffKzT(1-Olx,1-Oly,1,bi,bj)
724       I     , ikppkey       I     , ikppkey
725       O     , vddiff       O     , vddiff
726       U     , ghat       U     , ghat
727       O     , hbl )       O     , hbl
728         I     , mytime, mythid )
729        CALL TIMER_STOP ('KPPMIX [KPP_CALC]', myThid)        CALL TIMER_STOP ('KPPMIX [KPP_CALC]', myThid)
730    
731  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
# Line 753  c     horizontal smoothing of vertical v Line 759  c     horizontal smoothing of vertical v
759        DO k = 1, Nr        DO k = 1, Nr
760           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
761       I        k, bi, bj,       I        k, bi, bj,
762       U        KPPviscAz(1-OLx,1-OLy,k,bi,bj) )       U        KPPviscAz(1-OLx,1-OLy,k,bi,bj),
763         I        myThid )
764        END DO        END DO
765        _EXCH_XYZ_R8(KPPviscAz  , myThid )        _EXCH_XYZ_R8(KPPviscAz  , myThid )
766  #endif /* KPP_SMOOTH_VISC */  #endif /* KPP_SMOOTH_VISC */
# Line 763  c     horizontal smoothing of vertical d Line 770  c     horizontal smoothing of vertical d
770        DO k = 1, Nr        DO k = 1, Nr
771           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
772       I        k, bi, bj,       I        k, bi, bj,
773       U        KPPdiffKzS(1-OLx,1-OLy,k,bi,bj) )       U        KPPdiffKzS(1-OLx,1-OLy,k,bi,bj),
774         I        myThid )
775           CALL SMOOTH_HORIZ (           CALL SMOOTH_HORIZ (
776       I        k, bi, bj,       I        k, bi, bj,
777       U        KPPdiffKzT(1-OLx,1-OLy,k,bi,bj) )       U        KPPdiffKzT(1-OLx,1-OLy,k,bi,bj),
778         I        myThid )
779        END DO        END DO
780        _EXCH_XYZ_R8(KPPdiffKzS , myThid )        _EXCH_XYZ_R8(KPPdiffKzS , myThid )
781        _EXCH_XYZ_R8(KPPdiffKzT , myThid )        _EXCH_XYZ_R8(KPPdiffKzT , myThid )
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22