Packs all relevant parameters into a flux_surface_type instance. Note we can't fully populate the surface instance as we don't hold delrho
function get_parameters_as_surf() result(surf)
use geo_utils, only: flux_surface_type
implicit none
type(flux_surface_type) :: surf
surf = flux_surface_type(geotype, &
mMode=mMode, nMode=nMode, &
rmaj=rmaj, r_geo=R_geo, r=rhoc, &
aSurf=aSurf, sHorz=shift, &
sVert=shiftVert, delm=deltam,&
deln=deltan, delmp=deltampri, &
delnp=deltanpri, thm=thetam, &
thn=thetan, q=qinp, shat=shat, &
nt=ntheta,& ! Might not want to pack this here
n_mxh=n_mxh, c_mxh=c_mxh, s_mxh=s_mxh,&
dc_mxh_dr=dc_mxh_dr, ds_mxh_dr=ds_mxh_dr)
end function get_parameters_as_surf