/[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.45 by dimitri, Tue Dec 18 09:05:59 2007 UTC revision 1.46 by atn, Fri Dec 21 02:54:34 2007 UTC
# Line 588  C     the bottom of the mixing layer. Line 588  C     the bottom of the mixing layer.
588  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
589  C     Compute fraction of saltplume (flux) penetrating to  C     Compute fraction of saltplume (flux) penetrating to
590  C     the bottom of the mixing layer.  C     the bottom of the mixing layer.
591        DO j=1-OLy,sNy+OLy        IF ( useSALT_PLUME ) THEN
592           DO i=1-OLx,sNx+OLx          DO j=1-OLy,sNy+OLy
593              worka(i,j) = KPPhbl(i,j,bi,bj)             DO i=1-OLx,sNx+OLx
594           ENDDO                worka(i,j) = KPPhbl(i,j,bi,bj)
595        ENDDO             ENDDO
596        CALL SALT_PLUME_FRAC(          ENDDO
597       I     (sNx+2*OLx)*(sNy+2*OLy), minusone,SaltPlumeDepth,          CALL SALT_PLUME_FRAC(
598       U     worka,       I       (sNx+2*OLx)*(sNy+2*OLy), minusone,SaltPlumeDepth,
599       I     myTime, myIter, myThid )       U       worka,
600        DO j=1-OLy,sNy+OLy       I       myTime, myIter, myThid )
601           DO i=1-OLx,sNx+OLx          DO j=1-OLy,sNy+OLy
602              KPPplumefrac(i,j,bi,bj) = worka(i,j)             DO i=1-OLx,sNx+OLx
603           ENDDO                KPPplumefrac(i,j,bi,bj) = worka(i,j)
604        ENDDO             ENDDO
605            ENDDO
606          ENDIF
607  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
608    
609        ENDIF        ENDIF
Legend:
Removed from v.1.45  
changed lines
  Added in v.1.46
  ViewVC Help
Powered by ViewVC 1.1.22