35 SUBROUTINE trpaus(PTROP,TTROP,ZTROP,UTROP,VTROP,SHTROP)
41 use vrbls3d, only: pint, t, zint, uh, vh
44 use ctlblk_mod
, only: jsta, jend, spval, im, jm, lm, &
53 real,
PARAMETER :: crtlap=0.002e0, pstart=5.0e4
57 REAL,
dimension(IM,JM),
intent(out) :: ptrop,ttrop,ztrop,utrop, &
59 REAL tlapse(lm),dz2(lm),delt2(lm),tlapse2(lm)
62 real pm,delt,dz,rsqdif
94 loopl:
DO l=llmh-1,2,-1
96 delt = t(i,j,l-1)-t(i,j,l)
97 dz = d50*(zint(i,j,l-1)-zint(i,j,l+1))
100 IF ((tlapse(l)<crtlap).AND.(pm<pstart))
THEN
101 IF (l == 2 .AND. tlapse(l) < crtlap) goto 15
108 dz2(ll) = (2./3.)*(zint(i,j,ll-2)-zint(i,j,l+1))
109 IF ((dz2(ll) > 2000.) .AND. &
110 (dz2(ll+1) > 2000.)) go to 15
111 delt2(ll) = t(i,j,ll-2)-t(i,j,l)
112 tlapse2(ll) = -delt2(ll)/dz2(ll)
114 IF (tlapse2(ll) > crtlap)
THEN
123 15 ptrop(i,j) = d50*(pint(i,j,l)+pint(i,j,l+1))
124 ttrop(i,j) = t(i,j,l)
125 ztrop(i,j)= 0.5*(zint(i,j,l)+zint(i,j,l+1))
127 utrop(i,j) = uh(i,j,l)
128 vtrop(i,j) = vh(i,j,l)
129 dz = zint(i,j,l)-zint(i,j,l+1)
130 rsqdif = sqrt(((uh(i,j,l-1)-uh(i,j,l+1))*0.5)**2 &
131 & +((vh(i,j,l-1)-vh(i,j,l+1))*0.5)**2)
132 shtrop(i,j) = rsqdif/dz
139 dz = d50*(zint(i,j,2)-zint(i,j,3))
140 ptrop(i,j) = d50*(pint(i,j,2)+pint(i,j,3))
141 ttrop(i,j) = t(i,j,2)
142 ztrop(i,j)= d50*(zint(i,j,2)+zint(i,j,3))
143 utrop(i,j) = uh(i,j,2)
144 vtrop(i,j) = vh(i,j,2)
145 rsqdif = sqrt(((uh(i,j,1)-uh(i,j,3))*0.5)**2 &
146 & +((vh(i,j,1)-vh(i,j,3))*0.5)**2)
147 shtrop(i,j) = rsqdif/dz