--- trunk/phylmd/physiq.f	2015/02/13 18:22:38	129
+++ trunk/Sources/phylmd/physiq.f	2015/05/22 23:13:19	138
@@ -4,7 +4,7 @@
 
 contains
 
-  SUBROUTINE physiq(lafin, rdayvrai, time, dtphys, paprs, play, pphi, pphis, &
+  SUBROUTINE physiq(lafin, dayvrai, time, dtphys, paprs, play, pphi, pphis, &
        u, v, t, qx, omega, d_u, d_v, d_t, d_qx)
 
     ! From phylmd/physiq.F, version 1.22 2006/02/20 09:38:28
@@ -25,7 +25,7 @@
          ok_orodr, ok_orolf
     USE clmain_m, ONLY: clmain
     use clouds_gno_m, only: clouds_gno
-    USE comgeomphy, ONLY: airephy, cuphy, cvphy
+    USE comgeomphy, ONLY: airephy
     USE concvl_m, ONLY: concvl
     USE conf_gcm_m, ONLY: offline, raz_date
     USE conf_phys_m, ONLY: conf_phys
@@ -66,8 +66,8 @@
 
     logical, intent(in):: lafin ! dernier passage
 
-    REAL, intent(in):: rdayvrai
-    ! elapsed time since January 1st 0h of the starting year, in days
+    integer, intent(in):: dayvrai
+    ! current day number, based at value 1 on January 1st of annee_ref
 
     REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour
     REAL, intent(in):: dtphys ! pas d'integration pour la physique (seconde)
@@ -741,7 +741,7 @@
 
     ! Incrémenter le compteur de la physique
     itap = itap + 1
-    julien = MOD(NINT(rdayvrai), 360)
+    julien = MOD(dayvrai, 360)
     if (julien == 0) julien = 360
 
     forall (k = 1: llm) zmasse(:, k) = (paprs(:, k) - paprs(:, k + 1)) / rg
@@ -1229,8 +1229,8 @@
     ! Introduce the aerosol direct and first indirect radiative forcings:
     IF (ok_ade .OR. ok_aie) THEN
        ! Get sulfate aerosol distribution :
-       CALL readsulfate(rdayvrai, firstcal, sulfate)
-       CALL readsulfate_preind(rdayvrai, firstcal, sulfate_pi)
+       CALL readsulfate(dayvrai, time, firstcal, sulfate)
+       CALL readsulfate_preind(dayvrai, time, firstcal, sulfate_pi)
 
        CALL aeropt(play, paprs, t_seri, sulfate, rhcl, tau_ae, piz_ae, cg_ae, &
             aerindex)
@@ -1475,11 +1475,11 @@
     ! Si c'est la fin, il faut conserver l'etat de redemarrage
     IF (lafin) THEN
        itau_phy = itau_phy + itap
-       CALL phyredem("restartphy.nc", rlat, rlon, pctsrf, ftsol, ftsoil, &
-            tslab, seaice, fqsurf, qsol, fsnow, falbe, falblw, fevap, &
-            rain_fall, snow_fall, solsw, sollw, dlw, radsol, frugs, &
-            agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &
-            q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)
+       CALL phyredem("restartphy.nc", pctsrf, ftsol, ftsoil, tslab, seaice, &
+            fqsurf, qsol, fsnow, falbe, falblw, fevap, rain_fall, snow_fall, &
+            solsw, sollw, dlw, radsol, frugs, agesno, zmea, zstd, zsig, zgam, &
+            zthe, zpic, zval, t_ancien, q_ancien, rnebcon, ratqs, clwcon, &
+            run_off_lic_0, sig1, w01)
     ENDIF
 
     firstcal = .FALSE.