Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/dyn1d/lmdz1d.F90

    r3442 r3524  
    149149!                                                            (cf read_tsurf1d.F)
    150150
     151real wwww
    151152!vertical advection computation
    152153!       real d_t_z(llm), d_q_z(llm)
     
    11621163      endif
    11631164!
     1165       if (forcing_fire) THEN
     1166
     1167!let ww=if ( alt le 1100 ) then alt*-0.00001 else 0
     1168!let wt=if ( alt le 1100 ) then min( -3.75e-5 , -7.5e-8*alt)  else 0
     1169!let wq=if ( alt le 1100 ) then max( 1.5e-8 , 3e-11*alt)  else 0
     1170           d_t_adv=0.
     1171           d_q_adv=0.
     1172           teta=temp*(pzero/play)**rkappa
     1173           d_t_adv=0.
     1174           d_q_adv=0.
     1175           do l=2,llm-1
     1176              if (zlay(l)<=1100) then
     1177                  wwww=-0.00001*zlay(l)
     1178                  d_t_adv(l)=-wwww*(teta(l)-teta(l+1))/(zlay(l)-zlay(l+1)) /(pzero/play(l))**rkappa
     1179                  d_q_adv(l,1:2)=-wwww*(q(l,1:2)-q(l+1,1:2))/(zlay(l)-zlay(l+1))
     1180                  d_t_adv(l)=d_t_adv(l)+min(-3.75e-5 , -7.5e-8*zlay(l))
     1181                  d_q_adv(l,1)=d_q_adv(l,1)+max( 1.5e-8 , 3e-11*zlay(l))
     1182              endif
     1183           enddo
     1184
     1185        endif
     1186
    11641187!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    11651188!         call  writefield_phy('dv_age' ,dv_age,llm)
Note: See TracChangeset for help on using the changeset viewer.