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

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

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

revision 1.9 by heimbach, Tue Sep 18 20:30:59 2001 UTC revision 1.10 by heimbach, Sat Jul 13 03:12:30 2002 UTC
# Line 385  c Line 385  c
385    
386           end do           end do
387        end do        end do
388          
389    cph(
390    cph  without this store, there's a recomputation error for
391    cph  rib in adbldepth (probably partial recomputation problem)    
392    CADJ store Rib = comlev1_kpp
393    CADJ &   , key = ikey, shape = (/ (sNx+2*OLx)*(sNy+2*OLy),Nr /)
394    cph)
395    
396        do kl = 2, Nr        do kl = 2, Nr
397           do i = 1, imt           do i = 1, imt
398              if (kbl(i).eq.kmtj(i) .and. Rib(i,kl).gt.Ricr) kbl(i) = kl              if (kbl(i).eq.kmtj(i) .and. Rib(i,kl).gt.Ricr) kbl(i) = kl
# Line 432  c--   ensure bfsfc is never 0 Line 439  c--   ensure bfsfc is never 0
439           bfsfc(i) = sign(eins,bfsfc(i))*max(phepsi,abs(bfsfc(i)))           bfsfc(i) = sign(eins,bfsfc(i))*max(phepsi,abs(bfsfc(i)))
440        end do        end do
441    
442  CADJ store bfsfc = comlev1_kpp  cph(
443    cph  added stable to store list to avoid extensive recomp.
444    CADJ store bfsfc, stable = comlev1_kpp
445  CADJ &   , key = ikey, shape = (/ (sNx+2*OLx)*(sNy+2*OLy) /)  CADJ &   , key = ikey, shape = (/ (sNx+2*OLx)*(sNy+2*OLy) /)
446    cph)
447    
448  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
449  c check hbl limits for hekman or hmonob  c check hbl limits for hekman or hmonob
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22