Index: trunk/LMDZ.GENERIC/libf/phystd/dyn1d/iniprint.h
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/dyn1d/iniprint.h	(revision 1539)
+++ 	(revision )
@@ -1,1 +1,0 @@
-link ../../dyn3d/iniprint.h
Index: trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F	(revision 1539)
+++ trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F	(revision 1542)
@@ -11,7 +11,6 @@
 !      use comsaison_h
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat, volcapa
-      USE comgeomfi_h, only: lati, long, area
       use phyredem, only: physdem0,physdem1
-      use comgeomphy, only: initcomgeomphy, airephy
+      use comgeomphy, only: initcomgeomphy, cell_area
       use slab_ice_h, only: noceanmx
       use planete_mod, only: apoastr,periastr,year_day,peri_day,
@@ -159,9 +158,4 @@
       IF (.not.ALLOCATED(sinlon)) ALLOCATE(sinlon(1))
       IF (.not.ALLOCATED(coslon)) ALLOCATE(coslon(1))
-      !! those are defined in comgeomfi_h
-      IF (.not. ALLOCATED(lati)) ALLOCATE(lati(1))
-      IF (.not. ALLOCATED(long)) ALLOCATE(long(1))
-      IF (.not. ALLOCATED(area)) ALLOCATE(area(1))
-
 
       saveprofile=.false.
@@ -326,6 +320,5 @@
 !!! GEOGRAPHICAL INITIALIZATIONS
      !!! AREA. useless in 1D
-      area(1)=1.E+0
-      airephy(1)=area(1) !JL+EM to have access to the area in the diagfi.nc files. area in comgeomfi.h and aire in comgeom.h
+      cell_area(1)=1.E+0 !JL+EM to have access to the area in the diagfi.nc files.
      !!! GEOPOTENTIAL. useless in 1D because control by surface pressure
       phisfi(1)=0.E+0
@@ -511,5 +504,5 @@
       cpp=1.d-7 !JL because we divide by cpp in inifis, there may be a more elegant solution
       CALL inifis(1,llm,nq,day0,daysec,nday,dtphys,
-     .            latitude,longitude,area,rad,g,r,cpp)
+     .            latitude,longitude,cell_area,rad,g,r,cpp)
 
       nsoil=nsoilmx
@@ -673,5 +666,5 @@
          IF (qsurf(i_co2_ice).ge.1.E+0) THEN
             ! if we have some CO2 ice on the surface, change emissivity
-            if (lati(1).ge.0) then ! northern hemisphere
+            if (latitude(1).ge.0) then ! northern hemisphere
               emis(1)=emisice(1)
             else ! southern hemisphere
@@ -834,6 +827,6 @@
 c  It is needed to transfert physics variables to "physiq"...
 
-      call physdem0("startfi.nc",long,lati,nsoilmx,1,llm,nq,
-     &              dtphys,real(day0),time,area,
+      call physdem0("startfi.nc",longitude,latitude,nsoilmx,1,llm,nq,
+     &              dtphys,real(day0),time,cell_area,
      &              albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe)
       call physdem1("startfi.nc",nsoilmx,1,llm,nq,
@@ -843,10 +836,4 @@
      &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
 
-!      call physdem1(1,"startfi.nc",long,lati,nsoilmx,nq,
-!     &              dtphys,float(day0),
-!     &              time,tsurf,tsoil,emis,q2,qsurf,
-!     &              area,albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe,
-!     &              cloudfrac,totcloudfrac,hice,nametrac)
-
 c=======================================================================
 c  BOUCLE TEMPORELLE DU MODELE 1D 
@@ -861,9 +848,9 @@
          call stellarlong(day*1.0,zls)
 !         write(103,*) 'Ls=',zls*180./pi
-!         write(103,*) 'Lat=', lati(1)*180./pi
+!         write(103,*) 'Lat=', latitude(1)*180./pi
 !         write(103,*) 'RunEnd - Atmos. Temp. File'
 !         write(103,*) 'RunEnd - Atmos. Temp. File'
 !         write(104,*) 'Ls=',zls*180./pi
-!         write(104,*) 'Lat=', lati(1)
+!         write(104,*) 'Lat=', latitude(1)
 !         write(104,*) 'RunEnd - Atmos. Temp. File'
         ENDIF
@@ -923,5 +910,5 @@
 
 c       Pour certain test : pas de coriolis a l'equateur
-c       if(lati(1).eq.0.) then
+c       if(latitude(1).eq.0.) then
           DO ilayer=1,nlayer
              du(ilayer)=du(ilayer)+ (gru-u(ilayer))/1.e4