Index: trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F	(revision 4176)
+++ trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F	(revision 4185)
@@ -210,5 +210,4 @@
       end do
 
-
 !     update mass mixing ratio tracers
 
@@ -216,5 +215,4 @@
          do ig = 1,klon
             do i = 1,nqmx
-!               iq = niq(i) 
                zq(ig,l,i) = max(1.e-30, tr_seri(ig,l,i))
             end do
@@ -232,5 +230,4 @@
          do ig = 1,klon
             do i = 1,nqmx-nmicro
-c               iq = niq(i) 
                mmean(ig,l) = mmean(ig,l) + zq(ig,l,i)/M_tr(i)
             end do
@@ -248,31 +245,24 @@
 
       do l = 1,klev
-          do ig = 1,klon
-
-            ntot = pplay(ig,l)/(RKBOL*zt(ig,l))*1.e-6  ! in #/cm3 
+         do ig = 1,klon
+
+            ntot = pplay(ig,l)/(RKBOL*zt(ig,l))*1.e-6       ! in #/cm3 
             rho(ig,l) = (ntot * mmean(ig,l))/RNAVO*1.e3     ! in kg/m3
-
-c            write(*,*),'Air density: ',ig, l, rho(0,l)            
-
-!!  WARNING -> Cp here below doesn't depend on T (cpdet)
 
             do i = 1,nbq
                iq = niq(i)
+               if (iq == i_co2) then
+                 cpnew(ig,l) = cpnew(ig,l) + zq(ig,l,iq)*cpdet(zt(ig,l))
+               else
+                 cpnew(ig,l) = cpnew(ig,l) + zq(ig,l,iq)*cpi(i)
+               end if
                ni(i) = ntot*zq(ig,l,iq)*mmean(ig,l)/M_tr(iq)
-!! On a une super formule pour calculer cp_co2 sur Venus
-               if (iq == i_co2) then
-                  cpnew(ig,l) = cpnew(ig,l) + ni(i)*cpdet(zt(ig,l))
-               else
-                  cpnew(ig,l) = cpnew(ig,l) + ni(i)*cpi(i)
-               end if
                akknew(ig,l) = akknew(ig,l) + ni(i)*aki(i)
             end do 
  
-            cpnew(ig,l) = cpnew(ig,l)/ntot
             akknew(ig,l)= akknew(ig,l)/ntot
 
-
-          end do 
-       end do
+         end do 
+      end do
 
       return