/[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.7 by atn, Sat May 3 06:10:54 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 404  c--------------------------------------- Line 404  c---------------------------------------
404          temp1(i,j,1) = saltPlumeFlux(i,j,bi,bj)          temp1(i,j,1) = saltPlumeFlux(i,j,bi,bj)
405          temp2(i,j,1) = 0. _d 0          temp2(i,j,1) = 0. _d 0
406          DO k=2,Nr          DO k=2,Nr
407           temp1(i,j,k) = 0. _ 0           temp1(i,j,k) = 0. _d 0
408           temp2(i,j,k) = 0. _ 0           temp2(i,j,k) = 0. _d 0
409          ENDDO          ENDDO
410  #else /* def SALT_PLUME_VOLUME */  #else /* def SALT_PLUME_VOLUME */
411          DO k=1,Nr          DO k=1,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.7
  ViewVC Help
Powered by ViewVC 1.1.22