Index: trunk/LMDZ.VENUS/libf/phyvenus/iophy.F90
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/iophy.F90	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/iophy.F90	(revision 1545)
@@ -29,6 +29,6 @@
   USE ioipsl
   implicit none
-    real,dimension(klon),intent(in) :: rlon
-    real,dimension(klon),intent(in) :: rlat
+    real,dimension(klon),intent(in) :: rlon ! longitudes, in degrees
+    real,dimension(klon),intent(in) :: rlat ! latitudes, in degrees
 
     REAL,dimension(klon_glo)        :: rlat_glo
Index: trunk/LMDZ.VENUS/libf/phyvenus/nirco2abs.F
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/nirco2abs.F	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/nirco2abs.F	(revision 1545)
@@ -3,5 +3,5 @@
 
        use dimphy
-       use geometry_mod, only: longitude, latitude
+       use geometry_mod, only: longitude_deg, latitude_deg
        use chemparam_mod, only: i_co2, i_o
 c       use compo_hedin83_mod2 
@@ -236,6 +236,7 @@
             zday_int = (n-1)/float(nstep)
 
-            CALL zenang(0.,zday_int,RDAY/nstep,latitude,longitude,
-     s             mu0_int,fract_int)
+            CALL zenang(0.,zday_int,RDAY/nstep,
+     &                  latitude_deg,longitude_deg,
+     &                  mu0_int,fract_int)
 
             do ig=1,nlon
Index: trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90	(revision 1545)
@@ -13,5 +13,5 @@
       USE iostart
       USE infotrac
-      use geometry_mod, only: longitude, latitude
+      use geometry_mod, only: longitude_deg, latitude_deg
       USE time_phylmdz_mod, only: itau_phy, raz_date
 
@@ -33,4 +33,5 @@
 integer :: i,isoil
 CHARACTER(len=2) :: str2
+REAL :: lon_startphy(klon), lat_startphy(klon)
 
 ! les variables globales lues dans le fichier restart
@@ -64,17 +65,33 @@
       ENDIF
 
-! read latitudes
-call get_field("latitude",latitude,found)
-      IF (.not.found) THEN
-         PRINT*, 'phyetat0: Le champ <latitude> est absent'
-         CALL abort
-      ENDIF
-
-! read longitudes
-call get_field("longitude",longitude,found)
-      IF (.not.found) THEN
-         PRINT*, 'phyetat0: Le champ <longitude> est absent'
-         CALL abort
-      ENDIF
+! read latitudes and make a sanity check (because already known from dyn)
+call get_field("latitude",lat_startphy,found)
+IF (.not.found) THEN
+  PRINT*, 'phyetat0: Le champ <latitude> est absent'
+  CALL abort
+ENDIF
+DO i=1,klon
+  IF (ABS(lat_startphy(i)-latitude_deg(i))>=0.01) THEN
+    WRITE(*,*) "phyetat0: Warning! Latitude discrepancy wrt startphy file:",&
+               " i=",i," lat_startphy(i)=",lat_startphy(i),&
+               " latitude_deg(i)=",latitude_deg(i)
+    CALL abort
+  ENDIF
+ENDDO
+
+! read longitudes and make a sanity check (because already known from dyn)
+call get_field("longitude",lon_startphy,found)
+IF (.not.found) THEN
+  PRINT*, 'phyetat0: Le champ <longitude> est absent'
+  CALL abort
+ENDIF
+DO i=1,klon
+  IF (ABS(lon_startphy(i)-longitude_deg(i))>=0.01) THEN
+    WRITE(*,*) "phyetat0: Warning! Longitude discrepancy wrt startphy file:",&
+               " i=",i," lon_startphy(i)=",lon_startphy(i),&
+               " longitude_deg(i)=",longitude_deg(i)
+    CALL abort
+  ENDIF
+ENDDO
 
 ! read in other variables here ...
@@ -344,5 +361,5 @@
 
 ! do some more initializations
-call init_iophy_new(latitude,longitude)
+call init_iophy_new(latitude_deg,longitude_deg)
 
 end subroutine phyetat0
Index: trunk/LMDZ.VENUS/libf/phyvenus/phyredem.F90
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/phyredem.F90	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/phyredem.F90	(revision 1545)
@@ -12,5 +12,5 @@
                           put_var,put_field
       USE infotrac
-      use geometry_mod, only: longitude, latitude
+      use geometry_mod, only: longitude_deg, latitude_deg
       USE time_phylmdz_mod, only: day_end, annee_ref, itau_phy, raz_date
 
@@ -62,8 +62,8 @@
 
       CALL put_field("longitude", &
-                     "Longitudes de la grille physique",longitude)
+                     "Longitudes de la grille physique",longitude_deg)
      
       CALL put_field("latitude", &
-                     "Latitudes de la grille physique",latitude)
+                     "Latitudes de la grille physique",latitude_deg)
 
 ! variables
Index: trunk/LMDZ.VENUS/libf/phyvenus/physiq.F
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/physiq.F	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/physiq.F	(revision 1545)
@@ -58,5 +58,7 @@
       USE infotrac
       use dimphy
-      USE geometry_mod, only: longitude, latitude, cell_area, dx, dy
+      USE geometry_mod,only: longitude, latitude, ! in radians
+     &                       longitude_deg,latitude_deg, ! in degrees
+     &                       cell_area,dx,dy
       USE mod_phys_lmdz_para, only : is_parallel,jj_nb
       USE phys_state_var_mod ! Variables sauvegardees de la physique
@@ -665,5 +667,5 @@
       DO k = 1, klev
       DO i = 1, klon      
-         ilat=(latitude(i)/5.625) + 17.
+         ilat=(latitude_deg(i)/5.625) + 17.
          delta_temp(i,k)=mat_dtemp(INT(ilat),k)
       ENDDO
@@ -832,7 +834,8 @@
       IF (cycle_diurne) THEN
         zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
-        CALL zenang(zlongi,gmtime,zdtime,latitude,longitude,rmu0,fract)
+        CALL zenang(zlongi,gmtime,zdtime,latitude_deg,longitude_deg,
+     &              rmu0,fract)
       ELSE
-        call mucorr(klon,zlongi,latitude,rmu0,fract)
+        call mucorr(klon,zlongi,latitude_deg,rmu0,fract)
       ENDIF
       
@@ -859,5 +862,5 @@
      I                   debut,lafin,nqmax,
      I                   nlon,nlev,dtime,paprs,
-     I                   latitude,longitude,
+     I                   latitude_deg,longitude_deg,
      O                   tr_seri)
 
@@ -869,6 +872,6 @@
            DO k = 1, klev
            DO i = 1, klon      
-             ilat=(latitude(i)/5.625) + 17.
-!             PRINT*,INT(ilat),latitude(i),mat_dtemp(INT(ilat),k)
+             ilat=(latitude_deg(i)/5.625) + 17.
+!             PRINT*,INT(ilat),latitude_deg(i),mat_dtemp(INT(ilat),k)
              delta_temp(i,k)=mat_dtemp(INT(ilat),k)
            ENDDO
@@ -888,6 +891,6 @@
      I             nqmax,
      I             klon,
-     I             latitude,
-     I             longitude,
+     I             latitude_deg,
+     I             longitude_deg,
      I             nlev,
      I             dtime,
@@ -902,6 +905,6 @@
      I             nqmax,
      I             klon,
-     I             latitude,
-     I             longitude,
+     I             latitude_deg,
+     I             longitude_deg,
      I             nlev,
      I             dtime,
@@ -1032,5 +1035,5 @@
      $            paprs,pplay,ppk,radsol,falbe,
      e            solsw, sollw, sollwdown, fder,
-     e            longitude, latitude, dx, dy,   
+     e            longitude_deg, latitude_deg, dx, dy,   
      e            debut, lafin,
      s            d_t_vdf,d_u_vdf,d_v_vdf,d_ts,
@@ -1508,5 +1511,5 @@
 c A ADAPTER POUR VENUS!!!
 c            CALL lift_noro(klon,klev,dtime,paprs,pplay,
-c     e                   latitude,zmea,zstd,zpic,zgam,zthe,zpic,zval,
+c     e                   latitude_deg,zmea,zstd,zpic,zgam,zthe,zpic,zval,
 c     e                   igwd,idx,itest,
 c     e                   t_seri, u_seri, v_seri,
@@ -1567,5 +1570,6 @@
 c====================================================================
       if (ballons.eq.1) then
-        CALL ballon(30,pdtphys,rjourvrai,gmtime*RDAY,latitude,longitude,
+        CALL ballon(30,pdtphys,rjourvrai,gmtime*RDAY,
+     &              latitude_deg,longitude_deg,
 c    C              t,pplay,u,v,pphi)   ! alt above surface (smoothed for GCM)
      C              t,pplay,u,v,zphi)   ! alt above planet average radius
@@ -1601,5 +1605,5 @@
       CALL aaam_bud (27,klon,klev,rjourvrai,gmtime*RDAY,
      C               ra,rg,romega,
-     C               latitude,longitude,pphis,
+     C               latitude_deg,longitude_deg,pphis,
      C               zustrdr,zustrli,zustrcl,
      C               zvstrdr,zvstrli,zvstrcl,
@@ -1687,6 +1691,6 @@
 c     DO k = 1, klev
 c     DO i = 1, klon
-c       mang(i,k) = RA*cos(latitude(i)*RPI/180.)
-c    .     *(u_seri(i,k)+RA*cos(latitude(i)*RPI/180.)*ROMEGA)
+c       mang(i,k) = RA*cos(latitude(i))
+c    .     *(u_seri(i,k)+RA*cos(latitude(i))*ROMEGA)
 c    .     *cell_area(i)*(paprs(i,k)-paprs(i,k+1))/RG
 c       mangtot=mangtot+mang(i,k)
Index: trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F	(revision 1545)
@@ -27,5 +27,5 @@
 c======================================================================
       use dimphy
-      USE geometry_mod, ONLY: latitude
+      USE geometry_mod, ONLY: latitude_deg
       USE phys_state_var_mod, only: heat,cool,radsol,
      .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
@@ -164,11 +164,11 @@
          lat=1
        else
-         if (abs(latitude(j)).le.50.) then
+         if (abs(latitude_deg(j)).le.50.) then
 	   lat=1
-	 elseif (abs(latitude(j)).le.60.) then
+	 elseif (abs(latitude_deg(j)).le.60.) then
 	   lat=2
-	 elseif (abs(latitude(j)).le.70.) then
+	 elseif (abs(latitude_deg(j)).le.70.) then
 	   lat=3
-	 elseif (abs(latitude(j)).le.80.) then
+	 elseif (abs(latitude_deg(j)).le.80.) then
 	   lat=4
 	 else
Index: trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool	(revision 1543)
+++ trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool	(revision 1545)
@@ -24,5 +24,5 @@
 c======================================================================
       use dimphy
-      USE geometry_mod, ONLY: latitude
+      USE geometry_mod, ONLY: latitude ! in radians
       USE phys_state_var_mod, only: heat,cool,radsol,
      .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
@@ -98,6 +98,6 @@
 	    ztemp = tempCLee(level)*(1-fact)+tempCLee(level+1)*fact
 	    zdt   = dt_epCLee(level)*(1-fact)+dt_epCLee(level+1)*fact
-c           zt_eq(i,k) = ztemp + zdt*(cos(latitude(i)*RPI/180.)-2./RPI)
-	    zt_eq(i,k) = ztemp + zdt*(cos(latitude(i)*RPI/180.)-RPI/4.)
+c           zt_eq(i,k) = ztemp + zdt*(cos(latitude(i))-2./RPI)
+	    zt_eq(i,k) = ztemp + zdt*(cos(latitude(i))-RPI/4.)
 	    
 	  enddo
