/[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.28 by jmc, Mon Sep 27 14:45:09 2004 UTC revision 1.29 by jmc, Wed Oct 20 21:54:46 2004 UTC
# Line 386  C *note* should update GMREDI_RTRANSPORT Line 386  C *note* should update GMREDI_RTRANSPORT
386        ENDDO        ENDDO
387    
388  #ifdef ALLOW_KPP  #ifdef ALLOW_KPP
389  C-    Add non local KPP transport term (ghat) to diffusive T flux.  C-    Set non local KPP transport term (ghat):
390        IF (useKPP) THEN        IF ( useKPP .AND. k.GE.2 ) THEN
391         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
392          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
393           df(i,j) = 0. _d 0           df(i,j) = 0. _d 0
394          ENDDO          ENDDO
395         ENDDO         ENDDO
396         IF (tracerIdentity.EQ.GAD_TEMPERATURE) THEN         IF (tracerIdentity.EQ.GAD_TEMPERATURE) THEN
 C *note* should update KPP_TRANSPORT_T to set df  *aja*  
397          CALL KPP_TRANSPORT_T(          CALL KPP_TRANSPORT_T(
398       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
399       I     KappaRT,       O     df )
      U     df )  
400         ELSEIF (tracerIdentity.EQ.GAD_SALINITY) THEN         ELSEIF (tracerIdentity.EQ.GAD_SALINITY) THEN
401          CALL KPP_TRANSPORT_S(          CALL KPP_TRANSPORT_S(
402       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
403       I     KappaRT,       O     df )
      U     df )  
404  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
405         ELSEIF (tracerIdentity .GE. GAD_TR1) THEN         ELSEIF (tracerIdentity .GE. GAD_TR1) THEN
406          CALL KPP_TRANSPORT_PTR(          CALL KPP_TRANSPORT_PTR(
407       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,       I     iMin,iMax,jMin,jMax,bi,bj,k,km1,
408       I     tracerIdentity-GAD_TR1+1,KappaRT,       I     tracerIdentity-GAD_TR1+1,
409       U     df )       O     df )
410  #endif  #endif
411         ELSE         ELSE
412          PRINT*,'invalid tracer indentity: ', tracerIdentity          PRINT*,'invalid tracer indentity: ', tracerIdentity
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.22