/[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.6 by atn, Fri May 2 05:46:01 2014 UTC
# Line 156  c     temp?  (nx,ny,Nr)    - 3d working Line 156  c     temp?  (nx,ny,Nr)    - 3d working
156  c     ustar  (nx,ny)       - surface friction velocity                  (m/s)  c     ustar  (nx,ny)       - surface friction velocity                  (m/s)
157  c     bo     (nx,ny)       - surface turbulent buoyancy forcing     (m^2/s^3)  c     bo     (nx,ny)       - surface turbulent buoyancy forcing     (m^2/s^3)
158  c     bosol  (nx,ny)       - surface radiative buoyancy forcing     (m^2/s^3)  c     bosol  (nx,ny)       - surface radiative buoyancy forcing     (m^2/s^3)
159  c     boplume(nx,ny,Nr)    - surface haline buoyancy forcing        (m^2/s^3)  c     boplume(nx,ny,Nrp1)  - surface haline buoyancy forcing        (m^2/s^3)
160  c     shsq   (nx,ny,Nr)    - local velocity shear squared  c     shsq   (nx,ny,Nr)    - local velocity shear squared
161  c                            at interfaces for ri_iwmix             (m^2/s^2)  c                            at interfaces for ri_iwmix             (m^2/s^2)
162  c     dVsq   (nx,ny,Nr)    - velocity shear re surface squared  c     dVsq   (nx,ny,Nr)    - velocity shear re surface squared
# Line 185  c     vRef   (nx,ny)       - Reference m Line 185  c     vRef   (nx,ny)       - Reference m
185  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
186        _RL temp1   ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr          )        _RL temp1   ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr          )
187        _RL temp2   ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr          )        _RL temp2   ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr          )
188        _RL boplume ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr          )        _RL boplume ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, 0:Nr        )
189  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
190        _RL shsq  ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr            )        _RL shsq  ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr            )
191        _RL dVsq  ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr            )        _RL dVsq  ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr            )
# 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.6
  ViewVC Help
Powered by ViewVC 1.1.22