Changeset 594 for codes/icosagcm/devel/src/time

Timestamp:

10/18/17 18:46:11 (7 years ago)

Author:

dubos

Message:

devel : diagnose energy fluxes

File:

• codes/icosagcm/devel/src/time/timeloop_gcm.f90 (modified) (2 diffs)

codes/icosagcm/devel/src/time/timeloop_gcm.f90

@@ -273 +273 @@
           CALL HEVI_scheme(it, fluxt_zero)
        END SELECT
-       
+
        IF (MOD(it,itau_dissip)==0) THEN
           
@@ -304 +304 @@
        
        IF(MOD(it,itau_adv)==0) THEN
+          IF(diagflux_on) CALL diagflux_energy(adv_over_out, f_phis,f_rhodz,f_theta_rhodz,f_u, f_geopot,f_theta, f_hfluxt)
           CALL advect_tracer(f_hfluxt,f_wfluxt,f_u, f_q,f_rhodz, & ! update q and rhodz after RK step
                adv_over_out, f_masst,f_qmasst,f_massfluxt, f_qfluxt)  ! accumulate mass and fluxes if diagflux_on
-          fluxt_zero=.TRUE.
-          ! FIXME : check that rhodz is consistent with ps
-          IF (check_rhodz) THEN
+          fluxt_zero=.TRUE. ! restart accumulation of hfluxt and qfluxt at next call to explicit_scheme / HEVI_scheme
+          IF (check_rhodz) THEN ! check that rhodz is consistent with ps
              DO ind=1,ndomain
                 IF (.NOT. assigned_domain(ind)) CYCLE