Last change
on this file since 1010 was
917,
checked in by dubos, 5 years ago
|
devel : DYSL for compute_theta
|
File size:
805 bytes
|
Rev | Line | |
---|
[917] | 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.