Changeset 2990


Ignore:
Timestamp:
Sep 19, 2017, 11:25:13 AM (7 years ago)
Author:
fhourdin
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h

    r2933 r2990  
    50855085       implicit none
    50865086 
     5087#include "YOMCST.h"
    50875088#include "dimensions.h"
    50885089
     
    51605161         t_mod_cas(l)= t_prof_cas(k2) - frac*(t_prof_cas(k2)-t_prof_cas(k1))
    51615162         theta_mod_cas(l)= th_prof_cas(k2) - frac*(th_prof_cas(k2)-th_prof_cas(k1))
     5163         t_mod_cas(l)= theta_mod_cas(l)*(play(l)/100000.)**(RD/RCPD)
    51625164         thv_mod_cas(l)= thv_prof_cas(k2) - frac*(thv_prof_cas(k2)-thv_prof_cas(k1))
    51635165         thl_mod_cas(l)= thl_prof_cas(k2) - frac*(thl_prof_cas(k2)-thl_prof_cas(k1))
     
    51975199         t_mod_cas(l)= frac1*t_prof_cas(k1) - frac2*t_prof_cas(k2)
    51985200         theta_mod_cas(l)= frac1*th_prof_cas(k1) - frac2*th_prof_cas(k2)
     5201         t_mod_cas(l)= theta_mod_cas(l)*(play(l)/100000.)**(RD/RCPD)
    51995202         thv_mod_cas(l)= frac1*thv_prof_cas(k1) - frac2*thv_prof_cas(k2)
    52005203         thl_mod_cas(l)= frac1*thl_prof_cas(k1) - frac2*thl_prof_cas(k2)
     
    52455248         vg_mod_cas(l)= vg_prof_cas(nlev_cas)*fact                     !jyg
    52465249         w_mod_cas(l)= 0.0                                             !jyg
     5250         omega_mod_cas(l)= 0.0                                         !jyg
    52475251         du_mod_cas(l)= du_prof_cas(nlev_cas)*fact
    52485252         hu_mod_cas(l)= hu_prof_cas(nlev_cas)*fact                     !jyg
     
    52665270       enddo ! l
    52675271
    5268 !       do l = 1,llm
    5269 !       print *,'t_mod_cas(l),q_mod_cas(l),ht_mod_cas(l),hq_mod_cas(l) ',
    5270 !     $        l,t_mod_cas(l),q_mod_cas(l),ht_mod_cas(l),hq_mod_cas(l)
    5271 !       enddo
    5272  
    52735272          return
    52745273          end
Note: See TracChangeset for help on using the changeset viewer.