Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F	(revision 1387)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F	(revision 1388)
@@ -123,5 +123,7 @@
 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ) :: &
      dz8w,p8w,p,exner,t,t8w,rho,u,v,w,z,th
-REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
+!REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
+!     M_Q2
+REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: &
      M_Q2
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -164,4 +166,5 @@
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: PSFC,TSK
    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(OUT)  :: RTHBLTEN,RUBLTEN,RVBLTEN
+   REAL,DIMENSION(:),ALLOCATABLE,SAVE :: dtsw,dtlw
 #endif
 !-------------------------------------------
@@ -855,4 +858,5 @@
 ! specific treatment for tracers !
 !--------------------------------!
+#ifndef NOPHYS
 #ifdef NEWPHYS
 IF (MARS_MODE .EQ. 0) THEN
@@ -886,4 +890,5 @@
 END SELECT
 #endif
+#endif
 
 !!**********************************************************!!
@@ -1043,4 +1048,5 @@
 ! Tracer at surface !
 !-------------------!
+#ifndef NOPHYS
 SELECT CASE (MARS_MODE)
     CASE(0)  
@@ -1080,4 +1086,5 @@
 
 END SELECT
+#endif
 
 !-------------------!
@@ -1219,4 +1226,5 @@
 ! Tracers !  
 !---------!
+#ifndef NOPHYS
 #ifdef NEWPHYS
 pq(subs,:,:) = q_prof(:,:)  !! traceurs generiques, seuls noms sont specifiques
@@ -1234,4 +1242,5 @@
     STOP    
 END SELECT
+#endif
 #endif
 ENDDO
@@ -1486,4 +1495,22 @@
 ENDIF
 
+#ifdef NOPHYS
+!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!
+IF (firstcall .EQV. .true.) THEN
+  ALLOCATE(dtsw(nlayer))
+  ALLOCATE(dtlw(nlayer))
+  dtsw(1:nlayer) = M_Q2(1,1:nlayer,1)
+  dtlw(1:nlayer) = M_Q2(2,1:nlayer,2)
+ENDIF
+!print *, "SW", dtsw(kps:kpe) ! sw
+!print *, "LW", dtlw(kps:kpe) ! lw
+DO i= 1,ngrid
+  pdt(i,kps:kpe) = dtsw(1:nlayer) + dtlw(1:nlayer) 
+ENDDO
+!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!
+#endif
+
 DO j = jps,jpe
 DO i = ips,ipe
@@ -1548,4 +1575,5 @@
 ! Tracers !  
 !---------!
+#ifndef NOPHYS
 #ifdef NEWPHYS
 SCALAR(i,kps:kpe,j,1)=0.
@@ -1591,4 +1619,5 @@
 END SELECT
 #endif
+#endif
         !!TODO: check if adding the whole tendency once, and set the
         !!TODO: following tendencies to 0 until physics is called again