source: codes/icosagcm/devel/src/kernels_hex/compute_theta.k90 @ 935

Last change on this file since 935 was 917, checked in by dubos, 5 years ago

devel : DYSL for compute_theta
File size: 805 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- compute_theta ----------------------------------
3   IF(caldyn_eta==eta_mass) THEN ! Compute mass
4      DO l = ll_begin, ll_end
5         !DIR$ SIMD
6         DO ij=ij_begin_ext, ij_end_ext
7            m = MASS_DAK(ij,l)+(mass_col(ij)*g+ptop)*MASS_DBK(ij,l)
8            rhodz(ij,l) = m/g
9         END DO
10      END DO
11   END IF
12   DO iq=1,nqdyn
13      DO l = ll_begin, ll_end
14         !DIR$ SIMD
15         DO ij=ij_begin_ext, ij_end_ext
16            theta(ij,l,iq) = theta_rhodz(ij,l,iq)/rhodz(ij,l)
17         END DO
18      END DO
19   END DO
20   !---------------------------- compute_theta ----------------------------------
21   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.