/[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.6.1 by heimbach, Thu Nov 28 13:00:32 2002 UTC revision 1.9.6.2 by dimitri, Thu Feb 13 19:28:38 2003 UTC
# Line 183  c--------------------------------------- Line 183  c---------------------------------------
183    
184  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
185  c combine interior and boundary layer coefficients and nonlocal term  c combine interior and boundary layer coefficients and nonlocal term
186    c !!!NOTE!!! In shallow (2-level) regions and for shallow mixed layers
187    c (< 1 level), diffusivity blmc can become negative.  The max's below
188    c are a hack until this problem is properly diagnosed and fixed.
189  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
   
190        do k = 1, Nr        do k = 1, Nr
191           do i = 1, imt           do i = 1, imt
192              if (k .lt. kbl(i)) then              if (k .lt. kbl(i)) then
193                 do md = 1, mdiff                 diffus(i,k,1) = max ( blmc(i,k,1), viscAr  )
194                    diffus(i,k,md) = blmc(i,k,md)                 diffus(i,k,2) = max ( blmc(i,k,2), diffKrS )
195                 end do                 diffus(i,k,3) = max ( blmc(i,k,3), diffKrT )
196              else              else
197                 ghat(i,k) = 0.                 ghat(i,k) = 0.
198              endif              endif
# Line 661  c     fRi, fcon             function of Line 663  c     fRi, fcon             function of
663        _KPP_RL Rig        _KPP_RL Rig
664        _KPP_RL fRi, fcon        _KPP_RL fRi, fcon
665        _KPP_RL ratio        _KPP_RL ratio
666        integer i, ki, mr        integer i, ki
667        _KPP_RL c1, c0        _KPP_RL c1, c0
668    
669  #ifdef ALLOW_KPP_VERTICALLY_SMOOTH  #ifdef ALLOW_KPP_VERTICALLY_SMOOTH
670          integer mr
671  CADJ INIT kpp_ri_tape_mr = common, 1  CADJ INIT kpp_ri_tape_mr = common, 1
672  #endif  #endif
673    
Legend:
Removed from v.1.9.6.1  
changed lines
  Added in v.1.9.6.2
  ViewVC Help
Powered by ViewVC 1.1.22