Index: LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h
===================================================================
--- LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h	(revision 2989)
+++ LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h	(revision 2990)
@@ -5085,4 +5085,5 @@
        implicit none
  
+#include "YOMCST.h"
 #include "dimensions.h"
 
@@ -5160,4 +5161,5 @@
          t_mod_cas(l)= t_prof_cas(k2) - frac*(t_prof_cas(k2)-t_prof_cas(k1))
          theta_mod_cas(l)= th_prof_cas(k2) - frac*(th_prof_cas(k2)-th_prof_cas(k1))
+         t_mod_cas(l)= theta_mod_cas(l)*(play(l)/100000.)**(RD/RCPD)
          thv_mod_cas(l)= thv_prof_cas(k2) - frac*(thv_prof_cas(k2)-thv_prof_cas(k1))
          thl_mod_cas(l)= thl_prof_cas(k2) - frac*(thl_prof_cas(k2)-thl_prof_cas(k1))
@@ -5197,4 +5199,5 @@
          t_mod_cas(l)= frac1*t_prof_cas(k1) - frac2*t_prof_cas(k2)
          theta_mod_cas(l)= frac1*th_prof_cas(k1) - frac2*th_prof_cas(k2)
+         t_mod_cas(l)= theta_mod_cas(l)*(play(l)/100000.)**(RD/RCPD)
          thv_mod_cas(l)= frac1*thv_prof_cas(k1) - frac2*thv_prof_cas(k2)
          thl_mod_cas(l)= frac1*thl_prof_cas(k1) - frac2*thl_prof_cas(k2)
@@ -5245,4 +5248,5 @@
          vg_mod_cas(l)= vg_prof_cas(nlev_cas)*fact                     !jyg
          w_mod_cas(l)= 0.0                                             !jyg
+         omega_mod_cas(l)= 0.0                                         !jyg
          du_mod_cas(l)= du_prof_cas(nlev_cas)*fact 
          hu_mod_cas(l)= hu_prof_cas(nlev_cas)*fact                     !jyg
@@ -5266,9 +5270,4 @@
        enddo ! l
 
-!       do l = 1,llm
-!       print *,'t_mod_cas(l),q_mod_cas(l),ht_mod_cas(l),hq_mod_cas(l) ',
-!     $        l,t_mod_cas(l),q_mod_cas(l),ht_mod_cas(l),hq_mod_cas(l)
-!       enddo
- 
           return
           end