V3FIT
getrz.f
1  SUBROUTINE getrz(rmc,rms,zmc,zms,r0c,z0c,rhoc,rhos,m,mrho_in)
2  USE vname1
3 !-----------------------------------------------
4 ! D u m m y A r g u m e n t s
5 !-----------------------------------------------
6  INTEGER m, mrho_in
7  REAL(rprec), DIMENSION(0:mrho-1) :: rhoc, rhos
8  REAL(rprec) :: rmc, rms, zmc, zms, r0c, z0c
9 !-----------------------------------------------
10 ! L o c a l V a r i a b l e s
11 !-----------------------------------------------
12  INTEGER :: mrho1
13 !-----------------------------------------------
14 
15  mrho1 = mrho_in-1
16  rhos(0) = 0 !constraint on poloidal angle
17 
18  IF (m .eq. 0) THEN
19  rmc = r0c + t2m(m)*rhoc(1)
20  zmc = z0c + t2m(m)*rhos(1)
21  rms = 0
22  zms = 0
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))
28  ELSE !Can change highest m constraints here...
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
33  ENDIF
34 
35  END SUBROUTINE getrz