/[MITgcm]/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/kpp_calc.F
ViewVC logotype

Diff of /MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/kpp_calc.F

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

revision 1.4 by atn, Thu May 1 08:09:29 2014 UTC revision 1.5 by atn, Thu May 1 21:30:48 2014 UTC
# Line 639  C     the bottom of the mixing layer. Line 639  C     the bottom of the mixing layer.
639             ENDDO             ENDDO
640          ENDDO          ENDDO
641  #else  #else
642    Catn if decide to include in non-local transport, need to rethink
643    C    how to do. For now, set to zero.
644          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
645           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
646            DO k=1,Nr  C          DO k=1,Nr
647             IF(worka(i,j).LT.rF(k) .AND. work2(i,j).GE.rF(k)) THEN  C           IF(worka(i,j).LT.rF(k) .AND. work2(i,j).GE.rF(k)) THEN
648              KPPplumefrac(i,j,bi,bj) = SPplumek(i,j,k,bi,bj)  Catn:this is wrong KPPplumefrac(i,j,bi,bj) = SPplumek(i,j,k,bi,bj)
649             ENDIF              KPPplumefrac(i,j,bi,bj) = 0. _d 0
650            ENDDO  C           ENDIF
651    C          ENDDO
652           ENDDO           ENDDO
653          ENDDO          ENDDO
654  #endif /* ndef SALT_PLUME_VOLUME */  #endif /* ndef SALT_PLUME_VOLUME */
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22