Changeset 595 for codes/icosagcm/devel/src/time
- Timestamp:
- 10/19/17 10:32:52 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/time/timeloop_gcm.f90
r594 r595 304 304 305 305 IF(MOD(it,itau_adv)==0) THEN 306 IF(diagflux_on) CALL diagflux_energy(adv_over_out, f_phis,f_rhodz,f_theta_rhodz,f_u, f_geopot,f_theta, f_hfluxt)307 306 CALL advect_tracer(f_hfluxt,f_wfluxt,f_u, f_q,f_rhodz, & ! update q and rhodz after RK step 308 307 adv_over_out, f_masst,f_qmasst,f_massfluxt, f_qfluxt) ! accumulate mass and fluxes if diagflux_on 309 308 fluxt_zero=.TRUE. ! restart accumulation of hfluxt and qfluxt at next call to explicit_scheme / HEVI_scheme 309 ! At this point advect_tracer has obtained the halos of u and rhodz, 310 ! needed for correct computation of kinetic energy 311 IF(diagflux_on) CALL diagflux_energy(adv_over_out, f_phis,f_rhodz,f_theta_rhodz,f_u, f_geopot,f_theta, f_hfluxt) 312 310 313 IF (check_rhodz) THEN ! check that rhodz is consistent with ps 311 314 DO ind=1,ndomain
Note: See TracChangeset
for help on using the changeset viewer.