Index: /trunk/LMDZ.TITAN/deftank/callphys.def.allmuparams
===================================================================
--- /trunk/LMDZ.TITAN/deftank/callphys.def.allmuparams	(revision 4046)
+++ /trunk/LMDZ.TITAN/deftank/callphys.def.allmuparams	(revision 4047)
@@ -104,4 +104,6 @@
 # If yes, number of ices ? (must be compatible with traceur.def AND microphysical model)
 nices      = 1
+# Apply condensation heating rate?
+latent_heat = .true.
 # Use new optics for clouds ?
 opt4clouds = .true.
@@ -137,7 +139,7 @@
 rho_aer                = 600.
 # Enable/disable Haze production process 
-haze_production        = T
+haze_production        = .true.
 # Enable/disable Haze coagulation process 
-haze_coagulation       = T
+haze_coagulation       = .true.
 # Coagulation interactions, a combination of:
 #    0 - no interactions (same as haze_coagulation == F)
@@ -148,7 +150,7 @@
 haze_coag_interactions = 7
 # Enable/disable Haze sedimentation process 
-haze_sedimentation     = T
+haze_sedimentation     = .true.
 # Disable Fiadero correction for sedimentation process
-no_fiadero             = T
+no_fiadero             = .true.
 # Fiadero correction minimum ratio threshold
 fiadero_min_ratio      = 0.1
@@ -156,18 +158,18 @@
 fiadero_max_ratio      = 10.
 # Force settling velocity to M0
-wsed_m0                = T
+wsed_m0                = .true.
 # Force settling velocity to M3
-wsed_m3                = F
+wsed_m3                = .false.
 # Enable/disable clouds sedimentation process
 # (automatically set to F if clouds microphysics is not enabled)
-clouds_sedimentation   = T
+clouds_sedimentation   = .true.
 # Enable/disable clouds nucleation and condensation processes
 # (automatically set to F if clouds microphysics is not enabled)
-clouds_nuc_cond        = T
+clouds_nuc_cond        = .true.
 # Condensible species configuration file
 # (not needed if clouds microphysics is not enabled)
 specie_cfg             = mp2m_species.cfg
 # Enable/disable spherical mode transfert probability
-transfert_probability = T
+transfert_probability = .true.
 # Spherical mode transfert probability look-up tables file
 # (optional if 'transfert_probability' is False) 
@@ -175,5 +177,5 @@
 # Electric charging coagulation correction
 # If set to .false. then no correction is assumed
-electric_charging     = T 
+electric_charging     = .true.
 # File for the electric charging correction factor.
 # (optional if 'electric_charging' is False) 
Index: /trunk/LMDZ.TITAN/libf/phytitan/callkeys_mod.F90
===================================================================
--- /trunk/LMDZ.TITAN/libf/phytitan/callkeys_mod.F90	(revision 4046)
+++ /trunk/LMDZ.TITAN/libf/phytitan/callkeys_mod.F90	(revision 4047)
@@ -21,4 +21,6 @@
       logical,save :: callchim, callmufi, callclouds
 !$OMP THREADPRIVATE(callchim,callmufi,callclouds)
+      logical,save :: latent_heat
+!$OMP THREADPRIVATE(latent_heat)
       logical,save :: global1d
 !$OMP THREADPRIVATE(global1d)
Index: /trunk/LMDZ.TITAN/libf/phytitan/cond_muphy.F90
===================================================================
--- /trunk/LMDZ.TITAN/libf/phytitan/cond_muphy.F90	(revision 4046)
+++ /trunk/LMDZ.TITAN/libf/phytitan/cond_muphy.F90	(revision 4047)
@@ -69,5 +69,5 @@
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 do iq = 1, size(ices_indx)
-    dtlc(:,:) = dtlc(:,:) + (dqmuficond(:,:,iq) * Lc(:,:,iq))
+    dtlc(:,:) = dtlc(:,:) - (dqmuficond(:,:,iq) * Lc(:,:,iq))
 enddo
 
@@ -75,6 +75,6 @@
 ! Condensation heating rate :
 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! If ice formation  : dqmuficond > 0 --> dtlc > 0
-! Else vaporisation : dqmuficond < 0 --> dtlc < 0
+! If ice formation  : dqmuficond < 0 --> dtlc > 0
+! Else vaporisation : dqmuficond > 0 --> dtlc < 0
 dtlc(:,:) = dtlc(:,:) / cpp ! [K.s-1]
 
Index: /trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90
===================================================================
--- /trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90	(revision 4046)
+++ /trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90	(revision 4047)
@@ -490,4 +490,9 @@
      endif
 
+     write(*,*) "Apply condensation heating rate?"
+     latent_heat=.false. ! default value
+     call getin_p("latent_heat",latent_heat)
+     write(*,*)" latent_heat = ",latent_heat
+
      write(*,*) "Disable the coupling of microphysics within rad. transf. ?"
      write(*,*) "If disabled we will assume a planetwide vert. profile of extinction ..."
Index: /trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90
===================================================================
--- /trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90	(revision 4046)
+++ /trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90	(revision 4047)
@@ -196,5 +196,5 @@
       real,intent(in) :: pv(ngrid,nlayer)      ! Meridional wind component (m/s).
       real,intent(in) :: pt(ngrid,nlayer)      ! Temperature (K).
-      real,intent(in) :: pq(ngrid,nlayer,nq)   ! Tracers (kg/kg_of_air).
+      real,intent(in) :: pq(ngrid,nlayer,nq)   ! Tracers (X/kg_of_air). (X is kg or nothing (number of particles))
       real,intent(in) :: flxw(ngrid,nlayer)    ! Vertical mass flux (ks/s) at lower boundary of layer
 
@@ -431,4 +431,5 @@
 
         ! Allocate saved arrays (except for 1D model, where this has already been done)
+        ! For mesoscale it is done in the interface
 #ifndef MESOSCALE
         if (ngrid>1) call phys_state_var_init(nq)
@@ -1177,5 +1178,7 @@
                enddo
                call cond_muphy(ngrid,nlayer,pt,dmuficond,zdtlc)
-               !pdt(:,:) = pdt(:,:) + zdtlc(:,:)
+               if (latent_heat) then
+                  pdt(:,:) = pdt(:,:) + zdtlc(:,:)
+               endif
             endif
          endif ! callmufi