/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_advection.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_advection.F

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

revision 1.21 by edhill, Mon Mar 29 03:33:51 2004 UTC revision 1.22 by jmc, Fri Jun 25 18:19:20 2004 UTC
# Line 398  c       kp1=min(Nr,k+1) Line 398  c       kp1=min(Nr,k+1)
398          if (k.EQ.Nr) kp1Msk=0.          if (k.EQ.Nr) kp1Msk=0.
399    
400  C-- Compute Vertical transport  C-- Compute Vertical transport
401          IF (k.EQ.1) THEN  #ifdef ALLOW_AIM
402    C- a hack to prevent Water-Vapor vert.transport into the stratospheric level Nr
403            IF ( k.EQ.1 .OR.
404         &     (useAIM .AND. tracerIdentity.EQ.GAD_SALINITY .AND. k.EQ.Nr)
405         &              ) THEN
406    #else
407            IF ( k.EQ.1 ) THEN
408    #endif
409    
410  C- Surface interface :  C- Surface interface :
411           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
412            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
413             rTransKp1(i,j) = rTrans(i,j)             rTransKp1(i,j) = kp1Msk*rTrans(i,j)
414             rTrans(i,j) = 0.             rTrans(i,j) = 0.
415             fVerT(i,j,kUp) = 0.             fVerT(i,j,kUp) = 0.
416             af(i,j) = 0.             af(i,j) = 0.
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22