source: codes/icosagcm/devel/src/kernels_unst/theta.k90

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

devel : DYSL for compute_theta
File size: 1.0 KB
RevLine 
[614]1   !--------------------------------------------------------------------------
2   !---------------------------- theta ----------------------------------
3   IF(caldyn_eta==eta_mass) THEN ! Compute mass
[917]4      ! compute mass_col from rhodz
[614]5      !$OMP DO SCHEDULE(STATIC)
6      DO ij=1,primal_num
7         mass_col(ij)=0.
8         DO l = 1,llm
9            mass_col(ij) = mass_col(ij) + rhodz(l,ij)
10         END DO
11      END DO
12      !$OMP END DO
13      !$OMP DO SCHEDULE(STATIC)
14      DO ij = 1, primal_num
[658]15         !DIR$ SIMD
[614]16         DO l = 1, llm
[878]17            rhodz(l,ij) = MASS_DAK(l,ij) + mass_col(ij)*MASS_DBK(l,ij)
[614]18         END DO
19      END DO
20      !$OMP END DO
21   END IF
22   DO iq=1,nqdyn
23      !$OMP DO SCHEDULE(STATIC)
24      DO ij = 1, primal_num
[658]25         !DIR$ SIMD
[614]26         DO l = 1, llm
27            theta(l,ij,iq) = theta_rhodz(l,ij,iq)/rhodz(l,ij)
28         END DO
29      END DO
30      !$OMP END DO
31   END DO
32   !---------------------------- theta ----------------------------------
33   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.