Changeset 583 for codes/icosagcm/devel/src/transport/advect.F90

Timestamp:

10/13/17 21:21:27 (7 years ago)

Author:

dubos

Message:

devel : accumulate tracer fluxes over time for diagnostics

File:

• codes/icosagcm/devel/src/transport/advect.F90 (modified) (2 diffs)

codes/icosagcm/devel/src/transport/advect.F90

@@ -425 +425 @@
   ! Horizontal transport (S. Dubey, T. Dubos)
   ! Slope-limited van Leer approach with hexagons
-  SUBROUTINE compute_advect_horiz(update_mass,hfluxt,cc,gradq3d, mass,qi) 
+  SUBROUTINE compute_advect_horiz(update_mass,diagflux_on, hfluxt,cc,gradq3d, mass,qi,qfluxt)
     USE trace
     USE omp_para
     IMPLICIT NONE
-    LOGICAL, INTENT(IN)       :: update_mass
+    LOGICAL, INTENT(IN)       :: update_mass, diagflux_on
     REAL(rstd), INTENT(IN)    :: gradq3d(iim*jjm,llm,3) 
     REAL(rstd), INTENT(IN)    :: hfluxt(3*iim*jjm,llm) ! mass flux
@@ -435 +435 @@
     REAL(rstd), INTENT(INOUT) :: mass(iim*jjm,llm)
     REAL(rstd), INTENT(INOUT) :: qi(iim*jjm,llm)
+    REAL(rstd), INTENT(INOUT) :: qfluxt(3*iim*jjm,llm) ! time-integrated tracer flux
 
     REAL(rstd) :: dq,dmass,qe,ed(3), newmass