# Procedures for G2S
#
PROCEDURE g2s
PARAMETERS
wg = workspace
RESULT ws
LOCAL g2s res
ws = wg
g2s=fields() ; g2s.lptin=wg.ntc
g2s.Xin=wg.x ; g2s.Yin=wg.y ; g2s.Ein=wg.e
g2s.qmax=inquire("g2s> qmax ")
g2s.npt =inquire("g2s> number of points ")
g2s.ic =inquire("g2s> window function code ")
g2s.rho =inquire("g2s> number density ")
module/load "g2s.so" ;
res=module:execute("g2s", g2s)
ws.ntc=res.lptout
ws.x=res.Xout ; ws.y=res.Yout ; ws.e=res.Eout
ws.xlabel=res.xcaptout ; ws.ylabel=res.ycaptout
ENDPROCEDURE