--- trunk/phylmd/physiq.f	2014/12/18 17:30:24	118
+++ trunk/phylmd/physiq.f	2015/02/13 18:22:38	129
@@ -38,6 +38,7 @@
     USE dimphy, ONLY: klon
     USE dimsoil, ONLY: nsoilmx
     use drag_noro_m, only: drag_noro
+    use dynetat0_m, only: day_ref, annee_ref
     USE fcttre, ONLY: foeew, qsatl, qsats, thermcep
     use fisrtilp_m, only: fisrtilp
     USE hgardfou_m, ONLY: hgardfou
@@ -57,7 +58,7 @@
     use readsulfate_preind_m, only: readsulfate_preind
     use sugwd_m, only: sugwd
     USE suphec_m, ONLY: ra, rcpd, retv, rg, rlvtt, romega, rsigma, rtt
-    USE temps, ONLY: annee_ref, day_ref, itau_phy
+    USE temps, ONLY: itau_phy
     use unit_nml_m, only: unit_nml
     USE ymds2ju_m, ONLY: ymds2ju
     USE yoethf_m, ONLY: r2es, rvtmp2
@@ -66,7 +67,7 @@
     logical, intent(in):: lafin ! dernier passage
 
     REAL, intent(in):: rdayvrai
-    ! (elapsed time since January 1st 0h of the starting year, in days)
+    ! elapsed time since January 1st 0h of the starting year, in days
 
     REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour
     REAL, intent(in):: dtphys ! pas d'integration pour la physique (seconde)
@@ -223,8 +224,8 @@
     ! Variables propres a la physique
 
     INTEGER, save:: radpas
-    ! (Radiative transfer computations are made every "radpas" call to
-    ! "physiq".)
+    ! Radiative transfer computations are made every "radpas" call to
+    ! "physiq".
 
     REAL radsol(klon)
     SAVE radsol ! bilan radiatif au sol calcule par code radiatif
@@ -649,7 +650,7 @@
 
        CALL printflag(radpas, ok_journe, ok_instan, ok_region)
 
-       IF (dtphys * REAL(radpas) > 21600. .AND. cycle_diurne) THEN 
+       IF (dtphys * radpas > 21600. .AND. cycle_diurne) THEN 
           print *, "Au minimum 4 appels par jour si cycle diurne"
           call abort_gcm('physiq', &
                "Nombre d'appels au rayonnement insuffisant", 1)
@@ -680,7 +681,7 @@
        ! Initialisation des sorties
 
        call ini_histins(dtphys, ok_instan, nid_ins)
-       CALL ymds2ju(annee_ref, 1, int(day_ref), 0., date0)
+       CALL ymds2ju(annee_ref, 1, day_ref, 0., date0)
        ! Positionner date0 pour initialisation de ORCHIDEE
        print *, 'physiq date0: ', date0
     ENDIF test_firstcal
@@ -775,7 +776,7 @@
 
     CALL orbite(REAL(julien), longi, dist)
     IF (cycle_diurne) THEN
-       CALL zenang(longi, time, dtphys * REAL(radpas), mu0, fract)
+       CALL zenang(longi, time, dtphys * radpas, mu0, fract)
     ELSE
        mu0 = -999.999
     ENDIF