1 SUBROUTINE getrz(rmc,rms,zmc,zms,r0c,z0c,rhoc,rhos,m,mrho_in)
7 REAL(rprec),
DIMENSION(0:mrho-1) :: rhoc, rhos
8 REAL(rprec) :: rmc, rms, zmc, zms, r0c, z0c
19 rmc = r0c + t2m(m)*rhoc(1)
20 zmc = z0c + t2m(m)*rhos(1)
23 ELSE IF (m .lt. mrho1)
THEN
24 rmc = (t1m(m)*rhoc(m-1) + t2m(m)*rhoc(m+1))
25 zms = (t1m(m)*rhoc(m-1) - t2m(m)*rhoc(m+1))
26 rms = (t1m(m)*rhos(m-1) + t2m(m)*rhos(m+1))
27 zmc =-(t1m(m)*rhos(m-1) - t2m(m)*rhos(m+1))
29 rmc = t1m(m)*rhoc(m-1) * hb_parameter
30 zms = rmc * hb_parameter
31 rms = t1m(m)*rhos(m-1) * hb_parameter
32 zmc =-rms * hb_parameter