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_unst/compute_temperature.k90

    r915 r916  
    6060         !$OMP END DO 
    6161      END IF 
     62   CASE(thermo_variable_Cp) 
     63      !$OMP DO SCHEDULE(STATIC) 
     64      DO ij = 1, primal_num 
     65         !DIR$ SIMD 
     66         DO l = 1, llm 
     67            p_ik = pmid(l,ij) 
     68            theta_ik = temp(l,ij) 
     69            qv = q(l,ij,1) ! water vapor mixing ratio = mv/md 
     70            Cp_ik = nu*( theta_ik + Rd*log(p_ik/preff) ) 
     71            temp_ik = Treff* (Cp_ik/cpp)**(1./nu) 
     72            temp(l,ij) = temp_ik 
     73         END DO 
     74      END DO 
     75      !$OMP END DO 
    6276   CASE(thermo_moist) 
    6377      !$OMP DO SCHEDULE(STATIC) 
Note: See TracChangeset for help on using the changeset viewer.