File indexing completed on 2018-03-02 18:43:12 UTC
view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
e10b4e6008 Patr*0001 #include "SALT_PLUME_OPTIONS.h"
0002
0003
0004
0005
0006 SUBROUTINE SALT_PLUME_TENDENCY_APPLY_T(
73b1dccda0 Jean*0007 U gT_arr,
0008 I iMin,iMax,jMin,jMax, k, bi, bj,
0009 I myTime, myIter, myThid )
e10b4e6008 Patr*0010
0011
0012
0013
0014
0015
0016
0017
73b1dccda0 Jean*0018
e10b4e6008 Patr*0019 IMPLICIT NONE
0020 #include "SIZE.h"
0021 #include "EEPARAMS.h"
0022 #include "PARAMS.h"
73b1dccda0 Jean*0023 #include "GRID.h"
0024
e10b4e6008 Patr*0025 #include "SALT_PLUME.h"
0026
73b1dccda0 Jean*0027
0028
0029
0030
0031
0032
0033
0034
0035
0036 _RL gT_arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
0037 INTEGER iMin, iMax, jMin, jMax
0038 INTEGER k, bi, bj
0039 _RL myTime
0040 INTEGER myIter
0041 INTEGER myThid
e10b4e6008 Patr*0042
0043
0044 #ifdef ALLOW_SALT_PLUME
0045 #ifdef SALT_PLUME_VOLUME
0046
0047
73b1dccda0 Jean*0048 INTEGER i, j
e10b4e6008 Patr*0049
0050 _RL plumeEtend(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
0051
0052
0053 DO j=jMin,jMax
0054 DO i=iMin,iMax
73b1dccda0 Jean*0055 IF ( SaltPlumeDepth(i,j,bi,bj) .GT. ABS(rF(k)) ) THEN
0056 plumeEtend(i,j) = SPforcingT(i,j,k,bi,bj)
0057 gT_arr(i,j) = gT_arr(i,j) + plumeEtend(I,J)
0058 & *recip_drF(k)*mass2rUnit*_recip_hFacC(i,j,k,bi,bj)
e10b4e6008 Patr*0059 & /HeatCapacity_Cp
0060 ELSE
0061 plumeEtend(i,j) = 0. _d 0
0062 ENDIF
0063 ENDDO
0064 ENDDO
0065
0066 #ifdef ALLOW_DIAGNOSTICS
0067 IF ( useDiagnostics ) THEN
0068 CALL DIAGNOSTICS_FILL (
73b1dccda0 Jean*0069 & plumeEtend,'oceEPtnd',k,1,2,bi,bj,myThid )
e10b4e6008 Patr*0070 ENDIF
0071 #endif /* ALLOW_DIAGNOSTICS */
73b1dccda0 Jean*0072
e10b4e6008 Patr*0073 #endif /* SALT_PLUME_VOLUME */
0074 #endif /* ALLOW_SALT_PLUME */
73b1dccda0 Jean*0075
e10b4e6008 Patr*0076 RETURN
0077 END