Index: LMDZ5/trunk/libf/dyn3d/calfis.F
===================================================================
--- LMDZ5/trunk/libf/dyn3d/calfis.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3d/calfis.F	(revision 1615)
@@ -434,10 +434,8 @@
 c
       if (planet_type=="earth") then
-#ifdef CPP_EARTH
 cIM calcul PV a teta=350, 380, 405K
       CALL PVtheta(ngridmx,llm,pucov,pvcov,pteta,
      $           ztfi,zplay,zplev,
      $           ntetaSTD,rtetaSTD,PVteta)
-#endif
       endif
 c
@@ -450,6 +448,4 @@
 
 
-      if (planet_type=="earth") then
-#ifdef CPP_EARTH
 
 !      write(lunout,*) 'PHYSIQUE AVEC NSPLIT_PHYS=',nsplit_phys
@@ -460,5 +456,8 @@
       zdqfic(:,:,:)=0.
 
-      do isplit=1,nsplit_phys
+      if (planet_type=="earth") then
+#ifdef CPP_PHYS
+
+       do isplit=1,nsplit_phys
 
          jH_cur_split=jH_cur+(isplit-1) * dtvr / (daysec *nsplit_phys)
@@ -503,5 +502,10 @@
          zdqfic(:,:,:)=zdqfic(:,:,:)+zdqfi(:,:,:)
 
-      enddo
+       enddo ! of do isplit=1,nsplit_phys
+
+#endif
+! of #ifdef CPP_PHYS
+      endif ! of if (planet_type=="earth")
+
       zdufi(:,:)=zdufic(:,:)/nsplit_phys
       zdvfi(:,:)=zdvfic(:,:)/nsplit_phys
@@ -509,6 +513,4 @@
       zdqfi(:,:,:)=zdqfic(:,:,:)/nsplit_phys
 
-#endif
-      endif !of if (planet_type=="earth")
 
 500   CONTINUE
Index: LMDZ5/trunk/libf/dyn3d/ce0l.F90
===================================================================
--- LMDZ5/trunk/libf/dyn3d/ce0l.F90	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3d/ce0l.F90	(revision 1615)
@@ -28,4 +28,5 @@
   IMPLICIT NONE
 #ifndef CPP_EARTH
+#include "iniprint.h"
   WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
 #else
Index: LMDZ5/trunk/libf/dyn3d/gcm.F
===================================================================
--- LMDZ5/trunk/libf/dyn3d/gcm.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3d/gcm.F	(revision 1615)
@@ -21,6 +21,5 @@
 ! A nettoyer. On ne veut qu'une ou deux routines d'interface 
 ! dynamique -> physique pour l'initialisation
-! Ehouarn: for now these only apply to Earth:
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
       USE dimphy
       USE comgeomphy
@@ -180,11 +179,8 @@
 ! A nettoyer. On ne veut qu'une ou deux routines d'interface 
 ! dynamique -> physique pour l'initialisation
-! Ehouarn : temporarily (?) keep this only for Earth
-      if (planet_type.eq."earth") then
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
       CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
       call InitComgeomphy
 #endif
-      endif
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 c-----------------------------------------------------------------------
@@ -435,14 +431,11 @@
          WRITE(lunout,*)
      .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
-! Earth:
-         if (planet_type.eq."earth") then
-#ifdef CPP_EARTH
+! Physics:
+#ifdef CPP_PHYS
          CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys/nsplit_phys ,
      ,                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp     )
 #endif
-         endif ! of if (planet_type.eq."earth")
          call_iniphys=.false.
       ENDIF ! of IF (call_iniphys.and.(iflag_phys.eq.1))
-!#endif
 
 c  numero de stockage pour les fichiers de redemarrage:
@@ -466,5 +459,5 @@
 #endif
 
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
 ! Create start file (startphy.nc) and boundary conditions (limit.nc) 
 ! for the Earth verstion
Index: LMDZ5/trunk/libf/dyn3d/leapfrog.F
===================================================================
--- LMDZ5/trunk/libf/dyn3d/leapfrog.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3d/leapfrog.F	(revision 1615)
@@ -403,6 +403,8 @@
 ! Ehouarn: be careful, diagedyn is Earth-specific (includes ../phylmd/..)!
            IF (planet_type.eq."earth") THEN
+#ifdef CPP_EARTH
             CALL diagedyn(ztit,2,1,1,dtphys
      &    , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
+#endif
            ENDIF
          ENDIF ! of IF (ip_ebil_dyn.ge.1 )
Index: LMDZ5/trunk/libf/dyn3dpar/bands.F90
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/bands.F90	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/bands.F90	(revision 1615)
@@ -93,6 +93,6 @@
    SUBROUTINE  Set_Bands 
      USE parallel
-#ifdef CPP_EARTH
-! Ehouarn: what follows is only related to // physics; for now only for Earth 
+#ifdef CPP_PHYS
+! Ehouarn: what follows is only related to // physics
      USE mod_phys_lmdz_para, ONLY : jj_para_begin,jj_para_end
 #endif
@@ -106,6 +106,5 @@
       enddo
           
-#ifdef CPP_EARTH
-! Ehouarn: what follows is only related to // physics; for now only for Earth          
+#ifdef CPP_PHYS
       do i=0,MPI_Size-1
         jj_Nb_physic(i)=jj_para_end(i)-jj_para_begin(i)+1
@@ -332,6 +331,6 @@
     subroutine AdjustBands_physic
       use times
-#ifdef CPP_EARTH
-! Ehouarn: what follows is only related to // physics; for now only for Earth 
+#ifdef CPP_PHYS
+! Ehouarn: what follows is only related to // physics
       USE mod_phys_lmdz_para, only : klon_mpi_para_nb
 #endif
@@ -359,6 +358,5 @@
       medium=medium/mpi_size      
       NbTot=0
-#ifdef CPP_EARTH
-! Ehouarn: what follows is only related to // physics; for now only for Earth 
+#ifdef CPP_PHYS
       do i=0,mpi_size-1
         Inc(i)=nint(klon_mpi_para_nb(i)*(medium-value(i))/value(i))
Index: LMDZ5/trunk/libf/dyn3dpar/calfis_p.F
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/calfis_p.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/calfis_p.F	(revision 1615)
@@ -27,9 +27,6 @@
      $                  pdqfi,
      $                  pdpsfi)
-#ifdef CPP_EARTH
-! Ehouarn: For now, calfis_p needs Earth physics
-c
-c    Auteur :  P. Le Van, F. Hourdin 
-c   .........
+#ifdef CPP_PHYS
+! If using physics
       USE dimphy
       USE mod_phys_lmdz_para, mpi_root_xx=>mpi_root 
@@ -146,5 +143,4 @@
       REAL clesphy0( longcles )
 
-#ifdef CPP_EARTH
 c    Local variables :
 c    -----------------
@@ -489,6 +485,5 @@
 
 
-      IF (is_sequential) THEN
-c
+      IF (is_sequential.and.(planet_type=="earth")) THEN
 cIM calcul PV a teta=350, 380, 405K
         CALL PVtheta(ngridmx,llm,pucov,pvcov,pteta,
@@ -627,7 +622,4 @@
 c$OMP BARRIER
       
-      if (planet_type=="earth") then
-#ifdef CPP_EARTH
-
 !$OMP MASTER
 !      write(lunout,*) 'PHYSIQUE AVEC NSPLIT_PHYS=',nsplit_phys
@@ -639,4 +631,6 @@
       zdqfic_omp(:,:,:)=0.
 
+      if (planet_type=="earth") then
+#ifdef CPP_PHYS
       do isplit=1,nsplit_phys
 
@@ -687,11 +681,12 @@
       enddo
 
+#endif
+! of #ifdef CPP_PHYS
+      endif !of if (planet_type=="earth")
+
       zdufi_omp(:,:)=zdufic_omp(:,:)/nsplit_phys
       zdvfi_omp(:,:)=zdvfic_omp(:,:)/nsplit_phys
       zdtfi_omp(:,:)=zdtfic_omp(:,:)/nsplit_phys
       zdqfi_omp(:,:,:)=zdqfic_omp(:,:,:)/nsplit_phys
-
-#endif
-      endif !of if (planet_type=="earth")
 c$OMP BARRIER
 
@@ -1110,10 +1105,4 @@
       firstcal = .FALSE.
 
-#else
-      write(lunout,*)
-     & "calfis_p: for now can only work with parallel physics"
-      stop
-#endif
-! of #ifdef CPP_EARTH
       RETURN
       END
Index: LMDZ5/trunk/libf/dyn3dpar/ce0l.F90
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/ce0l.F90	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/ce0l.F90	(revision 1615)
@@ -31,4 +31,5 @@
   IMPLICIT NONE
 #ifndef CPP_EARTH
+#include "iniprint.h"
   WRITE(lunout,*)'limit_netcdf: Earth-specific routine, needs Earth physics'
 #else
Index: LMDZ5/trunk/libf/dyn3dpar/gcm.F
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/gcm.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/gcm.F	(revision 1615)
@@ -20,6 +20,5 @@
       USE control_mod
 
-! Ehouarn: for now these only apply to Earth:
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
       USE mod_grid_phy_lmdz
       USE mod_phys_lmdz_para, ONLY : klon_mpi_para_nb
@@ -187,13 +186,10 @@
       call ini_getparam("out.def")
       call Read_Distrib
-! Ehouarn : temporarily (?) keep this only for Earth
-      if (planet_type.eq."earth") then
-#ifdef CPP_EARTH
+
+#ifdef CPP_PHYS
         CALL Init_Phys_lmdz(iim,jjp1,llm,mpi_size,distrib_phys)
 #endif
-      endif ! of if (planet_type.eq."earth")
       CALL set_bands
-#ifdef CPP_EARTH
-! Ehouarn: For now only Earth physics is parallel
+#ifdef CPP_PHYS
       CALL Init_interface_dyn_phys
 #endif
@@ -207,12 +203,9 @@
 c$OMP END PARALLEL
 
-! Ehouarn : temporarily (?) keep this only for Earth
-      if (planet_type.eq."earth") then
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
 c$OMP PARALLEL
       call InitComgeomphy
 c$OMP END PARALLEL 
 #endif
-      endif ! of if (planet_type.eq."earth")
 
 c-----------------------------------------------------------------------
@@ -451,11 +444,9 @@
          WRITE(lunout,*)
      .       'GCM: WARNING!!! vitesse verticale nulle dans la physique'
-! Earth:
-         if (planet_type.eq."earth") then
-#ifdef CPP_EARTH
+! Physics:
+#ifdef CPP_PHYS
          CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys/nsplit_phys ,
      ,                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp     )
 #endif
-         endif ! of if (planet_type.eq."earth")
          call_iniphys=.false.
       ENDIF ! of IF (call_iniphys.and.(iflag_phys.eq.1))
@@ -490,4 +481,12 @@
  301  FORMAT('1'/,15x,'run du ', i2,'/',i2,'/',i4)
  302  FORMAT('1'/,15x,'    au ', i2,'/',i2,'/',i4)
+#endif
+
+#ifdef CPP_PHYS
+! Create start file (startphy.nc) and boundary conditions (limit.nc) 
+! for the Earth verstion
+       if (iflag_phys>=100) then
+          call iniaqua(ngridmx,latfi,lonfi,iflag_phys)
+       endif
 #endif
 
Index: LMDZ5/trunk/libf/dyn3dpar/gr_dyn_fi_p.F
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/gr_dyn_fi_p.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/gr_dyn_fi_p.F	(revision 1615)
@@ -3,7 +3,6 @@
 !
       SUBROUTINE gr_dyn_fi_p(nfield,im,jm,ngrid,pdyn,pfi)
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
 ! Interface with parallel physics,
-! for now this routine only works with Earth physics
       USE mod_interface_dyn_phys
       USE dimphy
@@ -40,10 +39,6 @@
       ENDDO
 c$OMP END DO NOWAIT
-#else
-      write(lunout,*) "gr_fi_dyn_p : This routine should not be called",
-     &   "without parallelized physics"
-      stop
 #endif
-! of #ifdef CPP_EARTH
+! of #ifdef CPP_PHYS
       RETURN
       END
Index: LMDZ5/trunk/libf/dyn3dpar/gr_fi_dyn_p.F
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/gr_fi_dyn_p.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/gr_fi_dyn_p.F	(revision 1615)
@@ -3,7 +3,6 @@
 !
       SUBROUTINE gr_fi_dyn_p(nfield,ngrid,im,jm,pfi,pdyn)
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
 ! Interface with parallel physics,
-! for now this routine only works with Earth physics
       USE mod_interface_dyn_phys
       USE dimphy
@@ -52,10 +51,6 @@
       ENDDO
 c$OMP END DO NOWAIT
-#else
-      write(lunout,*) "gr_fi_dyn_p : This routine should not be called",
-     &   "without parallelized physics"
-      stop
 #endif
-! of #ifdef CPP_EARTH
+! of #ifdef CPP_PHYS
       RETURN
       END
Index: LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F	(revision 1615)
@@ -728,6 +728,8 @@
 ! Ehouarn: be careful, diagedyn is Earth-specific (includes ../phylmd/..)!
            IF (planet_type.eq."earth") THEN
+#ifdef CPP_EARTH
             CALL diagedyn(ztit,2,1,1,dtphys
      &    , ucov    , vcov , ps, p ,pk , teta , q(:,:,1), q(:,:,2))
+#endif
            ENDIF
       ENDIF 
Index: LMDZ5/trunk/libf/dyn3dpar/mod_interface_dyn_phys.F90
===================================================================
--- LMDZ5/trunk/libf/dyn3dpar/mod_interface_dyn_phys.F90	(revision 1614)
+++ LMDZ5/trunk/libf/dyn3dpar/mod_interface_dyn_phys.F90	(revision 1615)
@@ -7,7 +7,6 @@
   
   
-#ifdef CPP_EARTH
+#ifdef CPP_PHYS
 ! Interface with parallel physics,
-! for now this routine only works with Earth physics
 CONTAINS
   
@@ -56,4 +55,4 @@
   END SUBROUTINE Init_interface_dyn_phys 
 #endif
-! of #ifdef CPP_EARTH
+! of #ifdef CPP_PHYS
 END MODULE mod_interface_dyn_phys
Index: LMDZ5/trunk/libf/phydev/comcstphy.h
===================================================================
--- LMDZ5/trunk/libf/phydev/comcstphy.h	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/comcstphy.h	(revision 1615)
@@ -0,0 +1,2 @@
+      real rradius,rr,rg,rcpp
+      common/comcstphy/rradius,rr,rg,rcpp
Index: LMDZ5/trunk/libf/phydev/comgeomphy.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/comgeomphy.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/comgeomphy.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/comgeomphy.F90
Index: LMDZ5/trunk/libf/phydev/dimphy.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/dimphy.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/dimphy.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/dimphy.F90
Index: LMDZ5/trunk/libf/phydev/iniphysiq.F
===================================================================
--- LMDZ5/trunk/libf/phydev/iniphysiq.F	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/iniphysiq.F	(revision 1615)
@@ -0,0 +1,114 @@
+!
+! $Id: iniphysiq.F 1403 2010-07-01 09:02:53Z fairhead $
+!
+c
+c
+      SUBROUTINE iniphysiq(ngrid,nlayer,
+     $           punjours,
+     $           pdayref,ptimestep,
+     $           plat,plon,parea,pcu,pcv,
+     $           prad,pg,pr,pcpp)
+      USE dimphy
+      USE mod_grid_phy_lmdz
+      USE mod_phys_lmdz_para
+      USE comgeomphy
+
+      IMPLICIT NONE
+c
+c=======================================================================
+c
+c   subject:
+c   --------
+c
+c   Initialisation for the physical parametrisations of the LMD 
+c   martian atmospheric general circulation modele.
+c
+c   author: Frederic Hourdin 15 / 10 /93
+c   -------
+c
+c   arguments:
+c   ----------
+c
+c   input:
+c   ------
+c
+c    ngrid                 Size of the horizontal grid.
+c                          All internal loops are performed on that grid.
+c    nlayer                Number of vertical layers.
+c    pdayref               Day of reference for the simulation
+c    firstcall             True at the first call
+c    lastcall              True at the last call
+c    pday                  Number of days counted from the North. Spring
+c                          equinoxe.
+c
+c=======================================================================
+c
+c-----------------------------------------------------------------------
+c   declarations:
+c   -------------
+ 
+cym#include "dimensions.h"
+cym#include "dimphy.h"
+cym#include "comgeomphy.h"
+#include "comcstphy.h"
+      REAL prad,pg,pr,pcpp,punjours
+ 
+      INTEGER ngrid,nlayer
+      REAL plat(ngrid),plon(ngrid),parea(klon_glo)
+      REAL pcu(klon_glo),pcv(klon_glo)
+      INTEGER pdayref
+      INTEGER :: ibegin,iend,offset
+ 
+      REAL ptimestep
+      CHARACTER (LEN=20) :: modname='iniphysiq'
+      CHARACTER (LEN=80) :: abort_message
+ 
+      IF (nlayer.NE.klev) THEN
+         PRINT*,'STOP in inifis'
+         PRINT*,'Probleme de dimensions :'
+         PRINT*,'nlayer     = ',nlayer
+         PRINT*,'klev   = ',klev
+         abort_message = ''
+         CALL abort_gcm (modname,abort_message,1)
+      ENDIF
+
+      IF (ngrid.NE.klon_glo) THEN
+         PRINT*,'STOP in inifis'
+         PRINT*,'Probleme de dimensions :'
+         PRINT*,'ngrid     = ',ngrid
+         PRINT*,'klon   = ',klon_glo
+         abort_message = ''
+         CALL abort_gcm (modname,abort_message,1)
+      ENDIF
+c$OMP PARALLEL PRIVATE(ibegin,iend) 
+c$OMP+         SHARED(parea,pcu,pcv,plon,plat)
+      
+      offset=klon_mpi_begin-1
+      airephy(1:klon_omp)=parea(offset+klon_omp_begin:
+     &                          offset+klon_omp_end)
+      cuphy(1:klon_omp)=pcu(offset+klon_omp_begin:offset+klon_omp_end)
+      cvphy(1:klon_omp)=pcv(offset+klon_omp_begin:offset+klon_omp_end)
+      rlond(1:klon_omp)=plon(offset+klon_omp_begin:offset+klon_omp_end)
+      rlatd(1:klon_omp)=plat(offset+klon_omp_begin:offset+klon_omp_end)
+
+!     call suphel
+!     prad,pg,pr,pcpp
+      rradius=prad
+      rg=pg
+      rr=pr
+      rcpp=pcpp
+
+      
+
+c$OMP END PARALLEL
+
+      print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
+      print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'
+
+      RETURN
+9999  CONTINUE
+      abort_message ='Cette version demande les fichier rnatur.dat
+     & et surf.def'
+      CALL abort_gcm (modname,abort_message,1)
+
+      END
Index: LMDZ5/trunk/libf/phydev/init_phys_lmdz.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/init_phys_lmdz.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/init_phys_lmdz.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/init_phys_lmdz.F90
Index: LMDZ5/trunk/libf/phydev/iophy.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/iophy.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/iophy.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/iophy.F90
Index: LMDZ5/trunk/libf/phydev/mod_grid_phy_lmdz.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_grid_phy_lmdz.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_grid_phy_lmdz.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_grid_phy_lmdz.F90
Index: LMDZ5/trunk/libf/phydev/mod_phys_lmdz_mpi_data.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_phys_lmdz_mpi_data.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_phys_lmdz_mpi_data.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_phys_lmdz_mpi_data.F90
Index: LMDZ5/trunk/libf/phydev/mod_phys_lmdz_mpi_transfert.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_phys_lmdz_mpi_transfert.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_phys_lmdz_mpi_transfert.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_phys_lmdz_mpi_transfert.F90
Index: LMDZ5/trunk/libf/phydev/mod_phys_lmdz_omp_data.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_phys_lmdz_omp_data.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_phys_lmdz_omp_data.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_phys_lmdz_omp_data.F90
Index: LMDZ5/trunk/libf/phydev/mod_phys_lmdz_omp_transfert.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_phys_lmdz_omp_transfert.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_phys_lmdz_omp_transfert.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_phys_lmdz_omp_transfert.F90
Index: LMDZ5/trunk/libf/phydev/mod_phys_lmdz_para.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_phys_lmdz_para.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_phys_lmdz_para.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_phys_lmdz_para.F90
Index: LMDZ5/trunk/libf/phydev/mod_phys_lmdz_transfert_para.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/mod_phys_lmdz_transfert_para.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/mod_phys_lmdz_transfert_para.F90	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/mod_phys_lmdz_transfert_para.F90
Index: LMDZ5/trunk/libf/phydev/phyaqua.F
===================================================================
--- LMDZ5/trunk/libf/phydev/phyaqua.F	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/phyaqua.F	(revision 1615)
@@ -0,0 +1,26 @@
+! Routines complementaires pour la physique planetaire.
+
+
+      subroutine iniaqua(nlon,latfi,lonfi,iflag_phys)
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!  Creation d'un etat initial et de conditions aux limites 
+!  (resp startphy.nc et limit.nc) pour des configurations idealisees 
+! du modele LMDZ dans sa version terrestre.
+!  iflag_phys est un parametre qui controle
+!  iflag_phys = N  
+!    de 100 a 199 : aqua planetes avec SST forcees
+!                 N-100 determine le type de SSTs
+!    de 200 a 299 : terra planetes avec Ts calcule
+!        
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+
+      integer nlon,iflag_phys
+cIM ajout latfi, lonfi
+      REAL, DIMENSION (nlon) :: lonfi, latfi
+
+
+      return
+      end
+
Index: LMDZ5/trunk/libf/phydev/physiq.F90
===================================================================
--- LMDZ5/trunk/libf/phydev/physiq.F90	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/physiq.F90	(revision 1615)
@@ -0,0 +1,124 @@
+! $Id: physiq.F 1565 2011-08-31 12:53:29Z jghattas $
+!#define IO_DEBUG
+
+      SUBROUTINE physiq (nlon,nlev, &
+     &            debut,lafin,jD_cur, jH_cur,pdtphys, &
+     &            paprs,pplay,pphi,pphis,presnivs,clesphy0, &
+     &            u,v,t,qx, &
+     &            flxmass_w, &
+     &            d_u, d_v, d_t, d_qx, d_ps &
+     &            , dudyn &
+     &            , PVteta)
+
+      USE dimphy
+      USE infotrac
+      USE comgeomphy
+
+      IMPLICIT none
+!======================================================================
+! Objet: Moniteur general de la physique du modele
+!======================================================================
+!
+!  Arguments:
+!
+! nlon----input-I-nombre de points horizontaux
+! nlev----input-I-nombre de couches verticales, doit etre egale a klev
+! debut---input-L-variable logique indiquant le premier passage
+! lafin---input-L-variable logique indiquant le dernier passage
+! jD_cur       -R-jour courant a l'appel de la physique (jour julien)
+! jH_cur       -R-heure courante a l'appel de la physique (jour julien)
+! pdtphys-input-R-pas d'integration pour la physique (seconde)
+! paprs---input-R-pression pour chaque inter-couche (en Pa)
+! pplay---input-R-pression pour le mileu de chaque couche (en Pa)
+! pphi----input-R-geopotentiel de chaque couche (g z) (reference sol)
+! pphis---input-R-geopotentiel du sol
+! presnivs-input_R_pressions approximat. des milieux couches ( en PA)
+! u-------input-R-vitesse dans la direction X (de O a E) en m/s
+! v-------input-R-vitesse Y (de S a N) en m/s
+! t-------input-R-temperature (K)
+! qx------input-R-humidite specifique (kg/kg) et d'autres traceurs
+! d_t_dyn-input-R-tendance dynamique pour "t" (K/s)
+! d_q_dyn-input-R-tendance dynamique pour "q" (kg/kg/s)
+! flxmass_w -input-R- flux de masse verticale
+! d_u-----output-R-tendance physique de "u" (m/s/s)
+! d_v-----output-R-tendance physique de "v" (m/s/s)
+! d_t-----output-R-tendance physique de "t" (K/s)
+! d_qx----output-R-tendance physique de "qx" (kg/kg/s)
+! d_ps----output-R-tendance physique de la pression au sol
+!IM
+! PVteta--output-R-vorticite potentielle a des thetas constantes
+!======================================================================
+#include "dimensions.h"
+#include "comcstphy.h"
+
+      integer jjmp1
+      parameter (jjmp1=jjm+1-1/jjm)
+      integer iip1
+      parameter (iip1=iim+1)
+
+      INTEGER ivap          ! indice de traceurs pour vapeur d'eau
+      PARAMETER (ivap=1)
+      INTEGER iliq          ! indice de traceurs pour eau liquide
+      PARAMETER (iliq=2)
+
+!
+!
+! Variables argument:
+!
+      INTEGER nlon
+      INTEGER nlev
+      REAL, intent(in):: jD_cur, jH_cur
+
+      REAL pdtphys
+      LOGICAL debut, lafin
+      REAL paprs(klon,klev+1)
+      REAL pplay(klon,klev)
+      REAL pphi(klon,klev)
+      REAL pphis(klon)
+      REAL presnivs(klev)
+      REAL znivsig(klev)
+      real pir
+
+      REAL u(klon,klev)
+      REAL v(klon,klev)
+      REAL t(klon,klev),theta(klon,klev)
+      REAL qx(klon,klev,nqtot)
+      REAL flxmass_w(klon,klev)
+      REAL omega(klon,klev) ! vitesse verticale en Pa/s
+      REAL d_u(klon,klev)
+      REAL d_v(klon,klev)
+      REAL d_t(klon,klev)
+      REAL d_qx(klon,klev,nqtot)
+      REAL d_ps(klon)
+      real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
+!IM definition dynamique o_trac dans phys_output_open
+!      type(ctrl_out) :: o_trac(nqtot)
+!FH! REAL PVteta(klon,nbteta)
+      REAL PVteta(klon,1)
+      REAL dudyn(iim+1,jjmp1,klev)
+
+    INTEGER        longcles
+    PARAMETER    ( longcles = 20 )
+
+real temp_newton(klon,klev)
+integer k
+logical, save :: first=.true.
+
+      REAL clesphy0( longcles      )
+
+d_u=0.
+d_v=0.
+d_t=0.
+d_qx=0.
+d_ps=0.
+
+     d_u(:,1)=-u(:,1)/86400.
+     do k=1,klev
+        temp_newton(:,k)=280.+cos(rlatd(:))*40.-pphi(:,k)/rg*6.e-3
+        d_t(:,k)=(temp_newton(:,k)-t(:,k))/1.e5
+     enddo
+
+
+      print*,'COUCOU PHYDEV'
+      return
+      end
Index: LMDZ5/trunk/libf/phydev/tetalevel.F
===================================================================
--- LMDZ5/trunk/libf/phydev/tetalevel.F	(revision 1615)
+++ LMDZ5/trunk/libf/phydev/tetalevel.F	(revision 1615)
@@ -0,0 +1,1 @@
+link ../phylmd/tetalevel.F
Index: LMDZ5/trunk/makegcm
===================================================================
--- LMDZ5/trunk/makegcm	(revision 1614)
+++ LMDZ5/trunk/makegcm	(revision 1615)
@@ -568,6 +568,10 @@
    set phys="L_PHY= LIBPHY="
 else
-   #Default planet type is Earth
-   set cppflags="$cppflags -DCPP_EARTH"
+   #We'll use some physics
+   set cppflags="$cppflags -DCPP_PHYS"
+   if ( `echo $physique | grep -i "^lmd"` != "" ) then
+     #For lmd physics, default planet type is Earth
+     set cppflags="$cppflags -DCPP_EARTH"
+   endif
 endif
 
Index: LMDZ5/trunk/makelmdz_fcm
===================================================================
--- LMDZ5/trunk/makelmdz_fcm	(revision 1614)
+++ LMDZ5/trunk/makelmdz_fcm	(revision 1615)
@@ -190,6 +190,11 @@
 if [[ "$physique" != "nophys" ]]
 then
-   #Default planet type is Earth
+   #We'll use some physics
+   CPP_KEY="$CPP_KEY CPP_PHYS"
+   if [[ "${physique:0:3}" == "lmd" ]]
+   then
+   #For lmd physics, default planet type is Earth
    CPP_KEY="$CPP_KEY CPP_EARTH"
+   fi
 fi