Changeset 2990

Timestamp:

Sep 19, 2017, 11:25:13 AM (7 years ago)

Author:

fhourdin

Message:

Modif 1D, lecture format standard. Interpolation en temperature potentielle.

File:

• LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -5085 +5085 @@
        implicit none
  
+#include "YOMCST.h"
 #include "dimensions.h"
 
@@ -5160 +5161 @@
          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 +5199 @@
          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 +5248 @@
          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 +5270 @@
        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