/[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.32 by jmc, Sun May 15 03:04:56 2005 UTC revision 1.33 by heimbach, Sat Feb 25 16:29:55 2006 UTC
# Line 182  c     vRef   (nx,ny)       - Reference m Line 182  c     vRef   (nx,ny)       - Reference m
182        _KPP_RL vddiff( ibot:itop    , jbot:jtop    , 0:Nrp1, mdiff )        _KPP_RL vddiff( ibot:itop    , jbot:jtop    , 0:Nrp1, mdiff )
183        _KPP_RL ghat  ( ibot:itop    , jbot:jtop    , Nr            )        _KPP_RL ghat  ( ibot:itop    , jbot:jtop    , Nr            )
184        _KPP_RL hbl   ( ibot:itop    , jbot:jtop                    )        _KPP_RL hbl   ( ibot:itop    , jbot:jtop                    )
185    cph(
186          _KPP_RL TTALPHA( ibot:itop, jbot:jtop, Nrp1 )
187          _KPP_RL SSBETA ( ibot:itop, jbot:jtop, Nrp1 )
188    cph)
189  #ifdef KPP_ESTIMATE_UREF  #ifdef KPP_ESTIMATE_UREF
190        _KPP_RL z0    ( ibot:itop    , jbot:jtop                    )        _KPP_RL z0    ( ibot:itop    , jbot:jtop                    )
191        _KPP_RL zRef  ( ibot:itop    , jbot:jtop                    )        _KPP_RL zRef  ( ibot:itop    , jbot:jtop                    )
# Line 252  c--------------------------------------- Line 256  c---------------------------------------
256        CALL STATEKPP(        CALL STATEKPP(
257       I       ikppkey, bi, bj, myThid       I       ikppkey, bi, bj, myThid
258       O     , work2, dbloc, Ritop       O     , work2, dbloc, Ritop
259       O     , vddiff(ibot,jbot,1,1), vddiff(ibot,jbot,1,2)       O     , TTALPHA, SSBETA
260       &     )       &     )
261        CALL TIMER_STOP ('STATEKPP      [KPP_CALC]', myThid)        CALL TIMER_STOP ('STATEKPP      [KPP_CALC]', myThid)
262    
# Line 292  c     horizontally smooth density relate Line 296  c     horizontally smooth density relate
296       U        Ritop (ibot,jbot,k)  )       U        Ritop (ibot,jbot,k)  )
297           CALL KPP_SMOOTH_HORIZ (           CALL KPP_SMOOTH_HORIZ (
298       I        k, bi, bj,       I        k, bi, bj,
299       U        vddiff(ibot,jbot,k,1) )       U        TTALPHA(ibot,jbot,k) )
300           CALL KPP_SMOOTH_HORIZ (           CALL KPP_SMOOTH_HORIZ (
301       I        k, bi, bj,       I        k, bi, bj,
302       U        vddiff(ibot,jbot,k,2) )       U        SSBETA(ibot,jbot,k) )
303        ENDDO        ENDDO
304  #endif /* KPP_SMOOTH_DENS */  #endif /* KPP_SMOOTH_DENS */
305    
# Line 329  CADJ store work2              = comlev1_ Line 333  CADJ store work2              = comlev1_
333  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE
334  CADJ store dbloc, Ritop, ghat = comlev1_kpp, key = ikppkey  CADJ store dbloc, Ritop, ghat = comlev1_kpp, key = ikppkey
335  CADJ store vddiff             = comlev1_kpp, key = ikppkey  CADJ store vddiff             = comlev1_kpp, key = ikppkey
336    CADJ store TTALPHA, SSBETA    = comlev1_kpp, key = ikppkey
337  #endif  #endif
338  cph)  cph)
339    
# Line 379  cph) Line 384  cph)
384          endif          endif
385    
386          bo(I,J) = - gravity *          bo(I,J) = - gravity *
387       &       ( vddiff(I,J,1,1) * (surfaceForcingT(i,j,bi,bj)+       &       ( TTALPHA(I,J,1) * (surfaceForcingT(i,j,bi,bj)+
388       &       surfaceForcingTice(i,j,bi,bj)) +       &       surfaceForcingTice(i,j,bi,bj)) +
389       &       vddiff(I,J,1,2) * surfaceForcingS(i,j,bi,bj) )       &       SSBETA(I,J,1) * surfaceForcingS(i,j,bi,bj) )
390       &       / work2(I,J)       &       / work2(I,J)
391    
392          bosol(I,J) = gravity * vddiff(I,J,1,1) * Qsw(i,j,bi,bj) *          bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) *
393       &       recip_Cp*recip_rhoConst       &       recip_Cp*recip_rhoConst
394       &       / work2(I,J)       &       / work2(I,J)
395    
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22