initialize the type hmap_knot, from given parameters as arguments
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | knot_kl(1:2) | |||
| real(kind=wp), | intent(in) | :: | knot_R0 | |||
| real(kind=wp), | intent(in) | :: | knot_delta |
self
FUNCTION hmap_knot_init_params(knot_kl,knot_R0,knot_delta) RESULT(sf) ! MODULES IMPLICIT NONE !----------------------------------------------------------------------------------------------------------------------------------- ! INPUT VARIABLES INTEGER,INTENT(IN) :: knot_kl(1:2) !parameters of the (k,l)-torus REAL(wp),INTENT(IN) :: knot_R0, knot_delta !major radius and shift !----------------------------------------------------------------------------------------------------------------------------------- ! OUTPUT VARIABLES TYPE(t_hmap_knot) :: sf !! self !=================================================================================================================================== SWRITE(UNIT_stdOut,'(4X,A)')'INIT HMAP :: KNOT ON A (k,l)-TORUS ...' sf%k=REAL(knot_kl(1), wp) sf%l=REAL(knot_kl(2), wp) sf%R0=knot_R0 sf%delta=knot_delta IF (.NOT.((sf%R0 - ABS(sf%delta)) > 0.0_wp)) THEN CALL abort(__STAMP__, & "hmap_knot init: condition R0 - |delta| > 0 not fulfilled!") END IF sf%initialized=.TRUE. SWRITE(UNIT_stdOut,'(4X,A)')'...DONE.' IF(.NOT.test_called) CALL hmap_knot_test(sf) END FUNCTION hmap_knot_init_params