Changeset 435


Ignore:
Timestamp:
06/15/16 22:59:19 (8 years ago)
Author:
ymipsl
Message:

Update of physic dcmip2016

YM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/trunk/src/physics_dcmip2016.f90

    r434 r435  
    225225          DO l=1,llm 
    226226           ll=llm+1-l 
    227            rho(l) = pmid(ij,ll)/(287*Temp(ij,l)) 
     227           rho(l) = pmid(ij,ll)/(287*Tfi(ij,ll)*(1+0.608*qfi(ij,ll,1))) 
    228228           z(l)=lastz 
    229229           lastz=lastz+ (p(ij,l)-p(ij,l+1)) /g / rho(l) 
Note: See TracChangeset for help on using the changeset viewer.