--- trunk/dyn3d/leapfrog.f 2018/03/07 13:46:18 262 +++ trunk/dyn3d/leapfrog.f 2018/03/07 14:41:46 263 @@ -206,9 +206,10 @@ q(:, :, :, 1)) ENDIF - IF (MOD(itau + 1, iecri * day_step) == 0) THEN + IF (MOD(itau + 1, iecri) == 0) THEN CALL geopot(teta, pk, pks, phis, phi) - CALL writehist(vcov, ucov, teta, pk, phi, q, masse, ps, itau) + CALL writehist(vcov, ucov, teta, pk, phi, q, masse, ps, & + itau_w = itau_dyn + itau + 1) END IF end do time_integration