Changeset 935 for codes/icosagcm/devel/src/kernels_unst
- Timestamp:
- 07/03/19 17:15:11 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/kernels_unst/caldyn_fast.k90
r837 r935 19 19 !DIR$ SIMD 20 20 DO l = 1, llm 21 berni(l,ij) = .5*(geopot(l,ij)+geopot(l+1,ij)) 21 Phi_ik = .5*(geopot(l,ij)+geopot(l+1,ij)) 22 berni(l,ij) = Phi_ik 22 23 END DO 23 24 END DO … … 28 29 !DIR$ SIMD 29 30 DO l = 1, llm 30 berni(l,ij)= .5*(geopot(l,ij)+geopot(l+1,ij))31 berni(l,ij) = berni(l,ij)+ pk(l,ij)*(cpp-theta(l,ij,1)) ! Gibbs = Cp.T-Ts = T(Cp-s)31 Phi_ik = .5*(geopot(l,ij)+geopot(l+1,ij)) 32 berni(l,ij) = Phi_ik + pk(l,ij)*(cpp-theta(l,ij,1)) ! Gibbs = Cp.T-Ts = T(Cp-s) 32 33 END DO 33 34 END DO 34 35 !$OMP END DO 35 36 CASE(thermo_variable_Cp) 36 ! thermodynamics with variable Cp37 ! Cp(T) = Cp0 * (T/T0)^nu38 ! => h = Cp(T).T/(nu+1)39 37 !$OMP DO SCHEDULE(STATIC) 40 38 DO ij = 1, primal_num 41 39 !DIR$ SIMD 42 40 DO l = 1, llm 43 berni(l,ij) = .5*(geopot(l,ij)+geopot(l+1,ij)) 44 cp_ik = cpp*(pk(l,ij)/Treff)**nu 45 berni(l,ij) = berni(l,ij) + pk(l,ij)*(cp_ik/(nu+1.)-theta(l,ij,1)) ! Gibbs = h-Ts = T(Cp/(nu+1)-s) 41 Phi_ik = .5*(geopot(l,ij)+geopot(l+1,ij)) 42 ! thermodynamics with variable Cp 43 ! Cp(T) = Cp0 * (T/T0)^nu 44 ! => h = Cp(T).T/(nu+1) 45 cp_ik = cpp*(pk(ij,l)/Treff)**nu 46 berni(l,ij) = Phi_ik + pk(l,ij)*(cp_ik/(nu+1.)-theta(l,ij,1)) ! Gibbs = h-Ts = T(Cp/(nu+1)-s) 46 47 END DO 47 48 END DO
Note: See TracChangeset
for help on using the changeset viewer.