ConnectTool
 All Classes Namespaces Files Functions Variables Pages
Functions | Variables
corfield_read_map2D Namespace Reference

Functions

def pattern_wso_wso_gif
 
def get_wso_wso_gif
 

Variables

string help_message
 
int ncar_min = 2207
 
int ncar_max = 2207
 
 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,'main_read_map')
 
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)
 
list local_paths = [os.path.join(folder_ref,pattern_wso_wso_gif(time.ncar))]
 Step 1. More...
 
tuple remote_path = get_wso_wso_gif(time.ncar)
 
tuple sub = plt.subplot(322)
 Step 2.3. More...
 
tuple img = mpimg.imread(local_paths[0])
 
tuple imgplot = sub.imshow(img)
 
int i = -1
 Step 2. More...
 
string map2D_type = 'WSO'
 Step 2.1. More...
 
tuple corfield_output_name = ('_CR%4d')
 
tuple corfield_output_file = os.path.join('dir.output',map2D_type.lower() + corfield_output_name + '_txt_tmp.dat')
 
tuple map2D_path = interface_read_map2D.ncar2map(map2D_type,time.ncar)
 
string map2D_time = 'crmap'
 
tuple simu_path = model_corfield_pfss.create_simu_folder()
 Step 2.2. More...
 
tuple ifile = open(os.path.join(simu_path,corfield_output_file))
 
tuple map2D_corfield = ascii_map2D.read(ifile)
 
tuple eq = plot_common.equator(sub,'black')
 

Function Documentation

def corfield_read_map2D.get_wso_wso_gif (   ncar)
Find in NSO website the png corresponding to the Carrington rotation.
Return its complete url.

Here is the call graph for this function:

def corfield_read_map2D.pattern_wso_wso_gif (   ncar)
Return the corresponding file pattern for WSO gif picture.

Here is the caller graph for this function:

Variable Documentation

string corfield_read_map2D.axis_color = 'k'
string corfield_read_map2D.axis_name_x = 'Carrington Longitude'
string corfield_read_map2D.axis_name_y = 'Latitude'
int corfield_read_map2D.axis_size = 10
tuple corfield_read_map2D.connect_tool_path = os.getcwd()
tuple corfield_read_map2D.corfield_exec = os.path.join(module_corona.path_corfield_pfss,'main_read_map')
tuple corfield_read_map2D.corfield_output_file = os.path.join('dir.output',map2D_type.lower() + corfield_output_name + '_txt_tmp.dat')
tuple corfield_read_map2D.corfield_output_name = ('_CR%4d')
corfield_read_map2D.delta = time_max.datetime-time_min.datetime
int corfield_read_map2D.dphi = 40
tuple corfield_read_map2D.eq = plot_common.equator(sub,'black')
tuple corfield_read_map2D.fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)
int corfield_read_map2D.fig_resolution = 100
tuple corfield_read_map2D.fig_size = (20,15)
tuple corfield_read_map2D.folder_ref = os.path.join('dir.test','REF_WSO' )
tuple corfield_read_map2D.folder_test = os.path.join('dir.test','TEST_CORFIELD')
string corfield_read_map2D.help_message
Initial value:
1 = """
2 Usage :\n
3 ./corfield_read_map2D.py
4 [-h] : brief help
5 [n1 n2] : from Carrington n1 to n2
6 """
int corfield_read_map2D.i = -1

Step 2.

Loop for all map2D_types

tuple corfield_read_map2D.ifile = open(os.path.join(simu_path,corfield_output_file))
tuple corfield_read_map2D.img = mpimg.imread(local_paths[0])
tuple corfield_read_map2D.imgplot = sub.imshow(img)
list corfield_read_map2D.local_paths = [os.path.join(folder_ref,pattern_wso_wso_gif(time.ncar))]

Step 1.

Download and plot a reference image.

tuple corfield_read_map2D.map2D_corfield = ascii_map2D.read(ifile)
list corfield_read_map2D.map2D_numsub = [1,3,4,6]
tuple corfield_read_map2D.map2D_path = interface_read_map2D.ncar2map(map2D_type,time.ncar)
string corfield_read_map2D.map2D_time = 'crmap'
string corfield_read_map2D.map2D_type = 'WSO'

Step 2.1.

Initialisation Corfield map2D outputs filename

list corfield_read_map2D.map2D_types = ['WSO_CR','NSO_CR','NSO_DATE','ADAPT_DATE']
tuple corfield_read_map2D.ncar_max = 2207
tuple corfield_read_map2D.ncar_min = 2207
corfield_read_map2D.ndays = delta.days
tuple corfield_read_map2D.remote_path = get_wso_wso_gif(time.ncar)
tuple corfield_read_map2D.simu_path = model_corfield_pfss.create_simu_folder()

Step 2.2.

Corfield reading/writing. Corfield simulation name

tuple corfield_read_map2D.sub = plt.subplot(322)

Step 2.3.

Plot map2D.

tuple corfield_read_map2D.time_max = module_timespace.Time()
tuple corfield_read_map2D.time_min = module_timespace.Time()
list corfield_read_map2D.times = [ module_timespace.Time() for i in range(0,ndays+1) ]
corfield_read_map2D.tmp = ncar_min
list corfield_read_map2D.xtick_label = [ str(int(xtick_loc[i])) for i in range(0,len(xtick_loc)) ]
list corfield_read_map2D.xtick_loc = [ i*dphi for i in range(0,int(360./dphi)+1) ]
list corfield_read_map2D.ytick_label = [ str(int(ytick_loc[i])) for i in range(0,len(ytick_loc)) ]
list corfield_read_map2D.ytick_loc = [ -90,-80,-60,-40,-20,0,20,40,60,80,90 ]