Changeset 3008 for LMDZ6/trunk/libf/phylmd

Timestamp:

Oct 6, 2017, 9:31:00 AM (7 years ago)

Author:

Laurent Fairhead

Message:

Modification in case theta exists in the forcings but not temperature.
We compute temp with theta

MPL

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h (modified) (3 diffs)

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

@@ -5161 +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)
+         if(theta_mod_cas(l).NE.0) 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))
@@ -5199 +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)
+         if(theta_mod_cas(l).NE.0) 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)
@@ -5275 +5275 @@
 
 
+
+