Changeset 2811 for LMDZ5/trunk


Ignore:
Timestamp:
Mar 6, 2017, 2:49:30 PM (7 years ago)
Author:
jyg
Message:

In interp_case_vertical and interp2_case_vertical,
dtrad_mod_cas were not computed (MPL)

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r2716 r2811  
    27832783         hq_mod_cas(l)= hq_prof_cas(k2) - frac*(hq_prof_cas(k2)-hq_prof_cas(k1))
    27842784         vq_mod_cas(l)= vq_prof_cas(k2) - frac*(vq_prof_cas(k2)-vq_prof_cas(k1))
     2785         dtrad_mod_cas(l)= dtrad_prof_cas(k2) - frac*(dtrad_prof_cas(k2)-dtrad_prof_cas(k1))
    27852786     
    27862787         else !play>plev_prof_cas(1)
    … …  
    28092810         hq_mod_cas(l)= frac1*hq_prof_cas(k1) - frac2*hq_prof_cas(k2)
    28102811         vq_mod_cas(l)= frac1*vq_prof_cas(k1) - frac2*vq_prof_cas(k2)
     2812         dtrad_mod_cas(l)= frac1*dtrad_prof_cas(k1) - frac2*dtrad_prof_cas(k2)
    28112813
    28122814         endif ! play.le.plev_prof_cas(1)
    … …  
    28372839         hq_mod_cas(l)= hq_prof_cas(nlev_cas)*fact                            !jyg
    28382840         vq_mod_cas(l)= vq_prof_cas(nlev_cas)*fact                            !jyg
     2841         dtrad_mod_cas(l)= dtrad_prof_cas(nlev_cas)*fact                      !jyg
    28392842 
    28402843        endif ! play
    … …  
    51625165         hq_mod_cas(l)= hq_prof_cas(k2) - frac*(hq_prof_cas(k2)-hq_prof_cas(k1))
    51635166         vq_mod_cas(l)= vq_prof_cas(k2) - frac*(vq_prof_cas(k2)-vq_prof_cas(k1))
     5167         dtrad_mod_cas(l)= dtrad_prof_cas(k2) - frac*(dtrad_prof_cas(k2)-dtrad_prof_cas(k1))
    51645168     
    51655169         else !play>plev_prof_cas(1)
    … …  
    51985202         hq_mod_cas(l)= frac1*hq_prof_cas(k1) - frac2*hq_prof_cas(k2)
    51995203         vq_mod_cas(l)= frac1*vq_prof_cas(k1) - frac2*vq_prof_cas(k2)
     5204         dtrad_mod_cas(l)= frac1*dtrad_prof_cas(k1) - frac2*dtrad_prof_cas(k2)
    52005205
    52015206         endif ! play.le.plev_prof_cas(1)
    … …  
    52345239         hq_mod_cas(l)= hq_prof_cas(nlev_cas)*fact                     !jyg
    52355240         vq_mod_cas(l)= vq_prof_cas(nlev_cas)*fact                     !jyg
     5241         dtrad_mod_cas(l)= dtrad_prof_cas(nlev_cas)*fact               !jyg
    52365242 
    52375243        endif ! play
Note: See TracChangeset for help on using the changeset viewer.