Public Member Functions
subroutine	init_corfield_trace ()
	Input trace_field parameters initialisation. Read input namelist, allocate seed_phi and seed_the. More...

subroutine	close_corfield_trace ()
	Deallocate trace_field parameters. More...

subroutine	trace_field
	Master fieldline tracing subroutine. More...

subroutine	pfss_line_tracer_grid
	Given a Magcor structure, trace magnetic fieldlines from all seed point of a regular grid. Curvilinear abscissa increases with constant ds. When the upper/bottom limit is overpassed, the integration stops just before. The last point is interpolated from previous point. More...

subroutine	pfss_line_tracer_list
	Given a Magcor structure, trace magnetic fieldlines from a list of seed points. Curvilinear abscissa increases with constant ds. When the upper/bottom limit is overpassed, the integration stops just before. The last point is interpolated from previous point. More...

Public Attributes
type(map)	ssmagmap

type(map)	ssfssmap

type(map)	sunlinemap

type(t_fieldline), dimension(:), allocatable	fieldlines

integer, parameter	nmax_integration = 999

real(kind=dp)	eps

real(kind=dp)	ds

character(3)	seed_pos = 'top'

character(3)	seed_unit = 'deg'

real(kind=dp), dimension(:), allocatable	seed_phi

real(kind=dp), dimension(:), allocatable	seed_the

Member Function/Subroutine Documentation

subroutine corfield_trace::close_corfield_trace ( )

Deallocate trace_field parameters.

Author

seed_phi
seed_the
fieldlines

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine corfield_trace::init_corfield_trace ( )

Input trace_field parameters initialisation. Read input namelist, allocate seed_phi and seed_the.

Author

Todo:: Allocatable arrays cannot be passed on a namelist, even if allocated when

Here is the caller graph for this function:

subroutine corfield_trace::pfss_line_tracer_grid ( )

Given a Magcor structure, trace magnetic fieldlines from all seed point of a regular grid.
Curvilinear abscissa increases with constant ds.
When the upper/bottom limit is overpassed, the integration stops just before.
The last point is interpolated from previous point.

Author

Output files with Carrington map option (magmaptype = CRMAP)

Magmapsrc_CRXXXX[_RXXX]_line_NNNNN.dat: fieldline number NNNNN (RXXX is ADAPT realization number). This format is used.
Magmapsrc_CRXXXX[_RXXX]_line_key.dat: list of all fieldlines seed points
Magmapsrc_crmap_XXXX[_RXXX].dat: magnetogram

Output files with Synchronic map option (magmaptype = SMAP)

Magmapsrc_YYYYMMDDHHMMSS[_RXXX]_line_NNNNN.dat: fieldline number NNNNN (RXXX is ADAPT realization number). This format is used.
Magmapsrc_YYYYMMDDHHMMSS[_RXXX]_line_key.dat: list of all fieldlines seed points
Magmapsrc_smap_YYYYMMDDHHMMSS[_RXXX].dat: magnetogram

Last point interpolation on source surface or sun surface: Linear interpolation using polintdouble with 2 points

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine corfield_trace::pfss_line_tracer_list ( )

Given a Magcor structure, trace magnetic fieldlines from a list of seed points. Curvilinear abscissa increases with constant ds.
When the upper/bottom limit is overpassed, the integration stops just before.
The last point is interpolated from previous point.

Author

Output files with Carrington map option (magmaptype = CRMAP)

Magmapsrc_CRXXXX[_RXXX]_linelist_NNNNN.dat: fieldline number NNNNN (RXXX is ADAPT realization number). This format is used.
Magmapsrc_CRXXXX[_RXXX]_linelist_key.dat: list of all fieldlines seed points
Magmapsrc_crmap_linelist_XXXX[_RXXX].dat: magnetogram

Output files with Synchronic map option (magmaptype = SMAP)

Magmapsrc_YYYYMMDDHHMMSS[_RXXX]_linelist_NNNNN.dat: fieldline number NNNNN (RXXX is ADAPT realization number). This format is used.
Magmapsrc_YYYYMMDDHHMMSS[_RXXX]_linelist_key.dat: list of all fieldlines seed points
Magmapsrc_smap_linelist_YYYYMMDDHHMMSS[_RXXX].dat: magnetogram

Last point interpolation on source surface or sun surface: Linear interpolation using polintdouble with 2 points

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine corfield_trace::trace_field ( )

Master fieldline tracing subroutine.

Author

Todo:: I've separated the two types of tracing in two almost identical routines (but should, in principle, merge all into a main subroutine with calls to simpler embedded functions). To do...

Here is the call graph for this function: