Changeset 3524 for LMDZ6/trunk/libf/phylmd/dyn1d/lmdz1d.F90

Timestamp:

May 27, 2019, 2:48:00 PM (5 years ago)

Author:

fhourdin

Message:

Correction pour le cas fire, salle en attendant le format commun ... et Godot.

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/lmdz1d.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/lmdz1d.F90

@@ -149 +149 @@
 !                                                            (cf read_tsurf1d.F)
 
+real wwww
 !vertical advection computation
 !       real d_t_z(llm), d_q_z(llm)
@@ -1162 +1163 @@
       endif
 !
+       if (forcing_fire) THEN
+
+!let ww=if ( alt le 1100 ) then alt*-0.00001 else 0
+!let wt=if ( alt le 1100 ) then min( -3.75e-5 , -7.5e-8*alt)  else 0
+!let wq=if ( alt le 1100 ) then max( 1.5e-8 , 3e-11*alt)  else 0
+           d_t_adv=0.
+           d_q_adv=0.
+           teta=temp*(pzero/play)**rkappa
+           d_t_adv=0.
+           d_q_adv=0.
+           do l=2,llm-1
+              if (zlay(l)<=1100) then
+                  wwww=-0.00001*zlay(l)
+                  d_t_adv(l)=-wwww*(teta(l)-teta(l+1))/(zlay(l)-zlay(l+1)) /(pzero/play(l))**rkappa
+                  d_q_adv(l,1:2)=-wwww*(q(l,1:2)-q(l+1,1:2))/(zlay(l)-zlay(l+1))
+                  d_t_adv(l)=d_t_adv(l)+min(-3.75e-5 , -7.5e-8*zlay(l))
+                  d_q_adv(l,1)=d_q_adv(l,1)+max( 1.5e-8 , 3e-11*zlay(l))
+              endif
+           enddo
+
+        endif
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !         call  writefield_phy('dv_age' ,dv_age,llm)