Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_big_step_utilities_em.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_big_step_utilities_em.F	(revision 1728)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_big_step_utilities_em.F	(revision 1729)
@@ -14,5 +14,4 @@
 
    USE module_domain, ONLY : domain
-   USE module_configure
    USE module_model_constants
    USE module_state_description
@@ -4476,6 +4475,7 @@
       p_phy(i,k,j) = p(i,k,j) + pb(i,k,j)
       pi_phy(i,k,j) = (p_phy(i,k,j)/p1000mb)**rcp
+      !! TAKE INTO ACCOUNT cp=f(T) on Venus
       IF (planet.eq. "venus" ) THEN 
-      !call tpot2t(1,th_phy(i,k,j), t_phy(i,k,j), pi_phy(i,k,j))
+        !call tpot2t(1,th_phy(i,k,j), t_phy(i,k,j), pi_phy(i,k,j))
         t_phy(i,k,j)= (th_phy(i,k,j)**nu - nu*(TT00**nu)*log((p1000mb/p_phy(i,k,j))**rcp))**(1/nu)
       ELSE
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_big_step_utilities_em.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_big_step_utilities_em.F	(revision 1728)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_big_step_utilities_em.F	(revision 1729)
@@ -18,4 +18,5 @@
    USE module_configure
    USE module_wrf_error
+
 CONTAINS
 
@@ -568,5 +569,5 @@
       ENDDO
       ENDDO
-      !print*,'muu(i,j)',muu(1,:)
+
    ELSE IF (name .EQ. 'v')THEN
 
@@ -587,5 +588,5 @@
       ENDDO
       ENDDO
-      !print*,'muv(i,j)',muv(1,:)
+
    ELSE IF (name .EQ. 'w')THEN
       itf=MIN(ite,ide-1)
@@ -620,6 +621,5 @@
       ENDDO
       ENDDO
-      !print*,'field(i,k,j)',field(1,:,1)
-      !print*,'(mu(i,j)+mub(i,j))',(mu(1,:)+mub(1,:))   
+
    ENDIF
 
@@ -1052,8 +1052,7 @@
 
   IF (non_hydrostatic) THEN
-      print*,'n_moist',n_moist
-      print*,'PARAM_FIRST_SCALAR',PARAM_FIRST_SCALAR
+
       IF (n_moist >= PARAM_FIRST_SCALAR ) THEN  
-        print*,'p(i,k,j)_before',p(1,:,1)
+
         DO j=jts,jtf
         DO k=kts,ktf
@@ -1076,5 +1075,5 @@
         ENDDO
         ENDDO
-        print*,'p(i,k,j)_after',p(1,:,1)
+
       ELSE
 
@@ -4076,5 +4075,5 @@
    INTEGER :: i_start, i_end, j_start, j_end, k_start, k_end, i_startu, j_startv
    INTEGER :: i, j, k
-   REAL    :: w1, w2, z0, z1, z2,temp
+   REAL    :: w1, w2, z0, z1, z2
 
 !-----------------------------------------------------------------------
@@ -4100,4 +4099,5 @@
 
 !  compute thermodynamics and velocities at pressure points
+
     do j = j_start,j_end
     do k = k_start, k_end
@@ -4107,4 +4107,5 @@
       p_phy(i,k,j) = p(i,k,j) + pb(i,k,j)
       pi_phy(i,k,j) = (p_phy(i,k,j)/p1000mb)**rcp
+      !! TAKE INTO ACCOUNT cp=f(T) on Venus
       IF (planet .eq. "venus" ) THEN
         t_phy(i,k,j)= (th_phy(i,k,j)**nu - nu*(TT00**nu)*log((p1000mb/p_phy(i,k,j))**rcp))**(1/nu)
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F	(revision 1728)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F	(revision 1729)
@@ -160,5 +160,5 @@
   ijds = min(ids, jds)
   ijde = max(ide, jde)
-  !num_moist = 3
+
   num_3d_m        = num_moist
   num_3d_c        = num_chem