initialize the type hmap_knot, reading from parameter file and then call init_params
self
FUNCTION hmap_knot_init() RESULT(sf) ! MODULES USE MODgvec_ReadInTools, ONLY: GETINTARRAY, GETREAL IMPLICIT NONE !----------------------------------------------------------------------------------------------------------------------------------- ! OUTPUT VARIABLES TYPE(t_hmap_knot) :: sf !! self !----------------------------------------------------------------------------------------------------------------------------------- ! LOCAL VARIABLES INTEGER :: knot_kl(1:2) !parameters of the (k,l)-torus REAL(wp) :: knot_R0, knot_delta !major radius and shift !=================================================================================================================================== SWRITE(UNIT_stdOut,'(4X,A)')'INIT HMAP :: KNOT ON A (k,l)-TORUS, GET PARAMETERS:' knot_kl=GETINTARRAY("hmap_knot_kl",2,proposal=(/2,3/)) knot_R0=GETREAL("hmap_knot_major_radius",1.0_wp) knot_delta=GETREAL("hmap_knot_delta_shift",0.4_wp) sf=hmap_knot_init_params(knot_kl,knot_R0,knot_delta) END FUNCTION hmap_knot_init