initialize the type hmap_cyl, reading from parameterfile and then call init_params
self
FUNCTION hmap_cyl_init() RESULT(sf) ! MODULES USE MODgvec_Globals, ONLY:TWOPI USE MODgvec_ReadInTools, ONLY: GETREAL IMPLICIT NONE !----------------------------------------------------------------------------------------------------------------------------------- ! OUTPUT VARIABLES TYPE(t_hmap_cyl) :: sf !! self !----------------------------------------------------------------------------------------------------------------------------------- ! LOCAL VARIABLES REAL(wp) :: cyl_len !=================================================================================================================================== SWRITE(UNIT_stdOut,'(4X,A)')'INIT HMAP :: CYLINDER WITH X1:=x, X2:=z, zeta := -2*pi*(y/cyl_len) , GET PARAMETERS:' cyl_len=GETREAL("hmap_cyl_len",TWOPI) sf= hmap_cyl_init_params(cyl_len) END FUNCTION hmap_cyl_init