10 |
|
|
11 |
USE dimens_m |
USE dimens_m |
12 |
USE dimphy |
USE dimphy |
13 |
|
use nr_util, only: pi |
14 |
USE suphec_m |
USE suphec_m |
15 |
USE yoegwd |
USE yoegwd |
16 |
|
|
65 |
!old zcons2=g**2/cpd |
!old zcons2=g**2/cpd |
66 |
zcons2 = rg**2/rcpd |
zcons2 = rg**2/rcpd |
67 |
!old zcons3=1.5*api |
!old zcons3=1.5*api |
68 |
zcons3 = 1.5*rpi |
zcons3 = 1.5*pi |
69 |
|
|
70 |
! 2. |
! 2. |
71 |
|
|
204 |
zu = pulow(jl) |
zu = pulow(jl) |
205 |
END IF |
END IF |
206 |
zphi = atan(pvlow(jl)/zu) |
zphi = atan(pvlow(jl)/zu) |
207 |
ppsi(jl, klev+1) = ptheta(jl)*rpi/180. - zphi |
ppsi(jl, klev+1) = ptheta(jl)*pi/180. - zphi |
208 |
zb(jl) = 1. - 0.18*pgamma(jl) - 0.04*pgamma(jl)**2 |
zb(jl) = 1. - 0.18*pgamma(jl) - 0.04*pgamma(jl)**2 |
209 |
zc(jl) = 0.48*pgamma(jl) + 0.3*pgamma(jl)**2 |
zc(jl) = 0.48*pgamma(jl) + 0.3*pgamma(jl)**2 |
210 |
pd1(jl) = zb(jl) - (zb(jl)-zc(jl))*(sin(ppsi(jl, klev+1))**2) |
pd1(jl) = zb(jl) - (zb(jl)-zc(jl))*(sin(ppsi(jl, klev+1))**2) |
338 |
zrhop = prho(jl, jk+1) |
zrhop = prho(jl, jk+1) |
339 |
znup(jl) = znup(jl) + (zdelp/rg)*((zstabp/zrhop+zstabm/zrhom)/2.)/ & |
znup(jl) = znup(jl) + (zdelp/rg)*((zstabp/zrhop+zstabm/zrhom)/2.)/ & |
340 |
zwind |
zwind |
341 |
IF ((znum(jl)<=rpi/2.) .AND. (znup(jl)>rpi/2.) .AND. (kkcrith( & |
IF ((znum(jl)<=pi/2.) .AND. (znup(jl)>pi/2.) .AND. (kkcrith( & |
342 |
jl)==klev)) kkcrith(jl) = jk |
jl)==klev)) kkcrith(jl) = jk |
343 |
|
|
344 |
END IF |
END IF |
363 |
zu = pum1(jl, jk) |
zu = pum1(jl, jk) |
364 |
END IF |
END IF |
365 |
zphi = atan(pvm1(jl, jk)/zu) |
zphi = atan(pvm1(jl, jk)/zu) |
366 |
ppsi(jl, jk) = ptheta(jl)*rpi/180. - zphi |
ppsi(jl, jk) = ptheta(jl)*pi/180. - zphi |
367 |
END IF |
END IF |
368 |
end DO |
end DO |
369 |
end DO |
end DO |