/[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.16 by jmc, Tue Dec 10 03:02:00 2002 UTC revision 1.17 by jmc, Sun Jan 26 21:08:36 2003 UTC
# Line 427  C--   Divergence of fluxes Line 427  C--   Divergence of fluxes
427         ENDDO         ENDDO
428        ENDDO        ENDDO
429    
430    #ifdef NONLIN_FRSURF
431    C-- account for 3.D divergence of the flow in rStar coordinate:
432          IF (calcAdvection .AND. select_rStar.GT.0) THEN
433           DO j=1-Oly,sNy+Oly-1
434            DO i=1-Olx,sNx+Olx-1
435             gTracer(i,j,k,bi,bj) = gTracer(i,j,k,bi,bj)
436         &     - (rStarExpC(i,j,bi,bj) - 1. _d 0)/deltaTfreesurf
437         &       *tracer(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
438            ENDDO
439           ENDDO
440          ENDIF
441          IF (calcAdvection .AND. select_rStar.LT.0) THEN
442           DO j=1-Oly,sNy+Oly-1
443            DO i=1-Olx,sNx+Olx-1
444             gTracer(i,j,k,bi,bj) = gTracer(i,j,k,bi,bj)
445         &     - rStarDhCDt(i,j,bi,bj)
446         &       *tracer(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
447            ENDDO
448           ENDDO
449          ENDIF
450    #endif /* NONLIN_FRSURF */
451          
452    
453        RETURN        RETURN
454        END        END
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22