Changeset 588

Timestamp:

10/16/17 15:23:49 (7 years ago)

Author:

dubos

Message:

devel : output time- averaged tracer fluxes

File:

• codes/icosagcm/devel/src/diagnostics/observable.f90 (modified) (2 diffs)

codes/icosagcm/devel/src/diagnostics/observable.f90

@@ -49 +49 @@
     USE theta2theta_rhodz_mod
     USE omega_mod
+    USE diagflux_mod
     LOGICAL, INTENT(IN) :: init
     INTEGER :: l
@@ -160 +161 @@
     END IF
 
+    IF(.NOT. init) THEN
+       IF(diagflux_on) THEN
+          CALL flux_centered_lonlat(1./(itau_out*dt) , f_qfluxt, f_qfluxt_lon, f_qfluxt_lat)
+          CALL output_field("qflux_lon",f_qfluxt_lon)
+          CALL output_field("qflux_lat",f_qfluxt_lat)
+          CALL zero_qfluxt  ! restart accumulating fluxes
+       END IF
+    END IF
   END SUBROUTINE write_output_fields_basic