/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_calc_rhs.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_calc_rhs.F

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

revision 1.17 by jmc, Sun Jan 26 21:08:36 2003 UTC revision 1.22 by dimitri, Thu Sep 25 03:01:59 2003 UTC
# Line 62  C  diffKh               :: horizontal di Line 62  C  diffKh               :: horizontal di
62  C  diffK4               :: bi-harmonic diffusion coefficient  C  diffK4               :: bi-harmonic diffusion coefficient
63  C  KappaRT              :: 3-D array for vertical diffusion coefficient  C  KappaRT              :: 3-D array for vertical diffusion coefficient
64  C  Tracer               :: tracer field  C  Tracer               :: tracer field
65  C  tracerIdentity       :: identifier for the tracer (required only for KPP)  C  tracerIdentity       :: identifier for the tracer (required for KPP and GM)
66  C  advectionScheme      :: advection scheme to use  C  advectionScheme      :: advection scheme to use
67  C  calcAdvection        :: =False if Advec terms computed with multiDim scheme  C  calcAdvection        :: =False if Advec terms computed with multiDim scheme
68  C  myThid               :: thread number  C  myThid               :: thread number
# Line 131  C--   Make local copy of tracer array Line 131  C--   Make local copy of tracer array
131    
132  C--   Unless we have already calculated the advection terms we initialize  C--   Unless we have already calculated the advection terms we initialize
133  C     the tendency to zero.  C     the tendency to zero.
134        IF (calcAdvection) THEN  C     <== now done earlier at the beginning of thermodynamics.
135         DO j=1-Oly,sNy+Oly  c     IF (calcAdvection) THEN
136          DO i=1-Olx,sNx+Olx  c      DO j=1-Oly,sNy+Oly
137           gTracer(i,j,k,bi,bj)=0. _d 0  c       DO i=1-Olx,sNx+Olx
138          ENDDO  c        gTracer(i,j,k,bi,bj)=0. _d 0
139         ENDDO  c       ENDDO
140        ENDIF  c      ENDDO
141    c     ENDIF
142    
143  C--   Pre-calculate del^2 T if bi-harmonic coefficient is non-zero  C--   Pre-calculate del^2 T if bi-harmonic coefficient is non-zero
144        IF (diffK4 .NE. 0.) THEN        IF (diffK4 .NE. 0.) THEN
# Line 402  C *note* should update KPP_TRANSPORT_T t Line 403  C *note* should update KPP_TRANSPORT_T t
403       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
404       I     KappaRT,       I     KappaRT,
405       U     df )       U     df )
406    #ifdef ALLOW_PTRACERS
407           ELSEIF (tracerIdentity .GE. GAD_TR1) THEN
408            CALL KPP_TRANSPORT_PTR(
409         I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
410         I     tracerIdentity-GAD_TR1+1,KappaRT,
411         U     df )
412    #endif
413         ELSE         ELSE
414            PRINT*,'invalid tracer indentity: ', tracerIdentity
415          STOP 'GAD_CALC_RHS: Ooops'          STOP 'GAD_CALC_RHS: Ooops'
416         ENDIF         ENDIF
417         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22