corfield_compute_cube3D Namespace Reference

Functions
def	pattern_wso_coronalfield_gif
def	get_wso_coronalfield_gif

Variables
string	help_message
int	ncar_min = 2209
int	ncar_max = 2209
	tmp = ncar_min
tuple	time_min = module_timespace.Time()
tuple	time_max = module_timespace.Time()
tuple	corfield_exec = os.path.join(module_corona.path_corfield_pfss,'corfield')
tuple	connect_tool_path = os.getcwd()
string	axis_name_y = 'Latitude'
string	axis_name_x = 'Carrington Longitude'
int	axis_size = 10
string	axis_color = 'k'
int	dphi = 40
list	xtick_loc = [ i*dphi for i in range(0,int(360./dphi)+1) ]
list	xtick_label = [ str(int(xtick_loc[i])) for i in range(0,len(xtick_loc)) ]
list	ytick_loc = [ -90,-80,-60,-40,-20,0,20,40,60,80,90 ]
list	ytick_label = [ str(int(ytick_loc[i])) for i in range(0,len(ytick_loc)) ]
tuple	folder_ref = os.path.join('dir.test','REF_WSO' )
tuple	folder_test = os.path.join('dir.test','TEST_CORFIELD')
tuple	fig_size = (20,15)
int	fig_resolution = 100
list	map2D_types = ['WSO_CR','NSO_CR','NSO_DATE','ADAPT_DATE']
list	map2D_numsub = [1,3,4,6]
	delta = time_max.datetime-time_min.datetime
	ndays = delta.days
list	times = [ module_timespace.Time() for i in range(0,ndays+1) ]
tuple	fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)
int	i = -1
	Step 1. More...
string	map2D_type = 'WSO'
	Step 1.1. More...
tuple	corfield_output_name = ('_CR%4d')
tuple	corfield_datacube_name = ('CR%4d')
tuple	map2D_path = interface_read_map2D.ncar2map(map2D_type,time.ncar)
tuple	time_str = ('CR%4d')
string	map2D_time = 'crmap'
string	pattern = map2D_type+'_'
tuple	corfield_datacube_file = os.path.join('dir.output',map2D_type.lower() + '_' + corfield_datacube_name + '_datacube.dat')
tuple	corfield_output_file = os.path.join('dir.output',map2D_type.lower() + corfield_output_name + '_txt_tmp.dat')
	find = False
	Step 1.2. More...
	simu_name = idir
tuple	simu_path = os.path.join(model_corfield_pfss.db_connect_tool_path,simu_name)
tuple	sub = plt.subplot(320+map2D_numsub[i])
	Step 1.3. More...
tuple	ifile = open(os.path.join(simu_path,corfield_output_file))
tuple	map2D_corfield = ascii_map2D.read(ifile)
tuple	cube3D_corfield = interface_read_cube3D.read('corfield',os.path.join(simu_path,corfield_datacube_file))
string	label_hcs = 'HCS'
tuple	hcs = module_plot.plot_hcs(sub,cube3D_corfield,[0.],[label_hcs],'red')
tuple	img = plot_common.map2D(sub,map2D_type,map2D_corfield,0)
tuple	eq = plot_common.equator(sub,'black')
list	local_paths = [os.path.join(folder_ref,pattern_wso_coronalfield_gif(time.ncar))]
	Step 2. More...
tuple	remote_path = get_wso_coronalfield_gif(time.ncar)
tuple	imgplot = sub.imshow(img)

Function Documentation

def corfield_compute_cube3D.get_wso_coronalfield_gif

(

ncar

)

Find in NWO website the gif corresponding to the Carrington rotation.
Return its complete url.

Here is the call graph for this function:

def corfield_compute_cube3D.pattern_wso_coronalfield_gif

(

ncar

)

Return the corresponding file pattern for WSO gif picture.

Here is the caller graph for this function:

Variable Documentation

string corfield_compute_cube3D.axis_color = 'k'

string corfield_compute_cube3D.axis_name_x = 'Carrington Longitude'

string corfield_compute_cube3D.axis_name_y = 'Latitude'

int corfield_compute_cube3D.axis_size = 10

tuple corfield_compute_cube3D.connect_tool_path = os.getcwd()

tuple corfield_compute_cube3D.corfield_datacube_file = os.path.join('dir.output',map2D_type.lower() + '_' + corfield_datacube_name + '_datacube.dat')

tuple corfield_compute_cube3D.corfield_datacube_name = ('CR%4d')

tuple corfield_compute_cube3D.corfield_exec = os.path.join(module_corona.path_corfield_pfss,'corfield')

tuple corfield_compute_cube3D.corfield_output_file = os.path.join('dir.output',map2D_type.lower() + corfield_output_name + '_txt_tmp.dat')

tuple corfield_compute_cube3D.corfield_output_name = ('_CR%4d')

tuple corfield_compute_cube3D.cube3D_corfield = interface_read_cube3D.read('corfield',os.path.join(simu_path,corfield_datacube_file))

corfield_compute_cube3D.delta = time_max.datetime-time_min.datetime

int corfield_compute_cube3D.dphi = 40

tuple corfield_compute_cube3D.eq = plot_common.equator(sub,'black')

tuple corfield_compute_cube3D.fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)

int corfield_compute_cube3D.fig_resolution = 100

tuple corfield_compute_cube3D.fig_size = (20,15)

corfield_compute_cube3D.find = False

Step 1.2.

Corfield reading/writing.

tuple corfield_compute_cube3D.folder_ref = os.path.join('dir.test','REF_WSO' )

tuple corfield_compute_cube3D.folder_test = os.path.join('dir.test','TEST_CORFIELD')

tuple corfield_compute_cube3D.hcs = module_plot.plot_hcs(sub,cube3D_corfield,[0.],[label_hcs],'red')

string corfield_compute_cube3D.help_message

Initial value:

= """
Usage :\n
./corfield_compute_cube3D.py
[-h]         : brief help
[n1 n2]      : from Carrington n1 to n2
"""

int corfield_compute_cube3D.i = -1

Step 1.

Loop for all map2D_types

tuple corfield_compute_cube3D.ifile = open(os.path.join(simu_path,corfield_output_file))

tuple corfield_compute_cube3D.img = plot_common.map2D(sub,map2D_type,map2D_corfield,0)

tuple corfield_compute_cube3D.imgplot = sub.imshow(img)

string corfield_compute_cube3D.label_hcs = 'HCS'

list corfield_compute_cube3D.local_paths = [os.path.join(folder_ref,pattern_wso_coronalfield_gif(time.ncar))]

Step 2.

Download and plot a reference image.

tuple corfield_compute_cube3D.map2D_corfield = ascii_map2D.read(ifile)

list corfield_compute_cube3D.map2D_numsub = [1,3,4,6]

tuple corfield_compute_cube3D.map2D_path = interface_read_map2D.ncar2map(map2D_type,time.ncar)

string corfield_compute_cube3D.map2D_time = 'crmap'

string corfield_compute_cube3D.map2D_type = 'WSO'

Step 1.1.

Initialisation Corfield map2D outputs filename

list corfield_compute_cube3D.map2D_types = ['WSO_CR','NSO_CR','NSO_DATE','ADAPT_DATE']

tuple corfield_compute_cube3D.ncar_max = 2209

tuple corfield_compute_cube3D.ncar_min = 2209

corfield_compute_cube3D.ndays = delta.days

string corfield_compute_cube3D.pattern = map2D_type+'_'

tuple corfield_compute_cube3D.remote_path = get_wso_coronalfield_gif(time.ncar)

corfield_compute_cube3D.simu_name = idir

tuple corfield_compute_cube3D.simu_path = os.path.join(model_corfield_pfss.db_connect_tool_path,simu_name)

tuple corfield_compute_cube3D.sub = plt.subplot(320+map2D_numsub[i])

Step 1.3.

Plot map2D and HCS.

tuple corfield_compute_cube3D.time_max = module_timespace.Time()

tuple corfield_compute_cube3D.time_min = module_timespace.Time()

tuple corfield_compute_cube3D.time_str = ('CR%4d')

list corfield_compute_cube3D.times = [ module_timespace.Time() for i in range(0,ndays+1) ]

corfield_compute_cube3D.tmp = ncar_min

list corfield_compute_cube3D.xtick_label = [ str(int(xtick_loc[i])) for i in range(0,len(xtick_loc)) ]

list corfield_compute_cube3D.xtick_loc = [ i*dphi for i in range(0,int(360./dphi)+1) ]

list corfield_compute_cube3D.ytick_label = [ str(int(ytick_loc[i])) for i in range(0,len(ytick_loc)) ]

list corfield_compute_cube3D.ytick_loc = [ -90,-80,-60,-40,-20,0,20,40,60,80,90 ]