2 ! { dg-require-effective-target pthread }
3 ! { dg-options "-Ofast -ftree-parallelize-loops=4" }
15 REAL, DIMENSION( its
:ite
, kts
:kte
), &
18 REAL, DIMENSION( ims
:ime
, kms
:kme
), &
23 REAL, DIMENSION( its
:ite
, kts
:kte
) :: &
30 diffus(x
,y
) = 8.794e-5*x
**1.81/y
31 diffac(a
,b
,c
,d
,e
) = d
*a
*a
/(xka(c
,d
)*rv
*c
*c
)+1./(e
*diffus(c
,b
))
32 venfac(a
,b
,c
) = (viscos(b
,c
)/diffus(b
,a
))**(.3333333) &
33 /viscos(b
,c
)**(.5)*(den0
/c
)**0.25
39 if(t(i
,k
).lt
.ttp
) then
40 qs(i
,k
) =psat
*(tr
**xa
)*exp(xb
*(1.-tr
))
42 qs0(i
,k
) =psat
*(tr
**xa
)*exp(xb
*(1.-tr
))
45 if(t(i
,k
).ge
.t0c
) then
46 work1(i
,k
) = diffac(xl(i
,k
),p(i
,k
),t(i
,k
),den(i
,k
),qs(i
,k
))
48 work2(i
,k
) = venfac(p(i
,k
),t(i
,k
),den(i
,k
))