Ignore:
Timestamp:
06/18/19 14:17:55 (5 years ago)
Author:
dubos
Message:

devel : added variable Cp to compute_temperature

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/devel/src/kernels_hex/compute_temperature.k90

    r915 r916  
    5252         END DO 
    5353      END IF 
     54   CASE(thermo_variable_Cp) 
     55      DO l = ll_begin, ll_end 
     56         !DIR$ SIMD 
     57         DO ij=ij_begin, ij_end 
     58            p_ik = pmid(ij,l) 
     59            theta_ik = temp(ij,l) 
     60            qv = q(ij,l,1) ! water vapor mixing ratio = mv/md 
     61            Cp_ik = nu*( theta_ik + Rd*log(p_ik/preff) ) 
     62            temp_ik = Treff* (Cp_ik/cpp)**(1./nu) 
     63            temp(ij,l) = temp_ik 
     64         END DO 
     65      END DO 
    5466   CASE(thermo_moist) 
    5567      DO l = ll_begin, ll_end 
Note: See TracChangeset for help on using the changeset viewer.