ConnectTool
 All Classes Namespaces Files Functions Variables Pages
Variables
db_corfield_map2D_smap Namespace Reference

Variables

string help_message
 
string map2D_type = 'NSO'
 Step 1.1. More...
 
tuple date_min = datetime.datetime(2019,1,1)
 
list date_max = interface_read_map2D.info[map2D_type]
 
 tmp = date_min
 
tuple corfield_exec = os.path.join(module_corona.path_corfield_pfss,'test_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_test = os.path.join('dir.test','COMPARE_DB_CORFIELD')
 
tuple fig_size = (25,15)
 
int fig_resolution = 100
 
list map2D_types = ['NSO_DATE','ADAPT_DATE']
 
list map2D_numsub = [1,4]
 
 delta = date_max-date_min
 
 ndays = delta.days
 
list dates = [ date_min + i*datetime.timedelta(days=1) 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...
 
tuple corfield_output_name = ('-lat_%s')
 
tuple corfield_output_file = os.path.join('dir.output',map2D_type.lower() + corfield_output_name + '_txt_tmp.dat')
 
tuple map2D_path = interface_read_map2D.date2map(map2D_type,date)
 
string map2D_time = 'smap'
 
tuple simu_path = model_corfield_pfss.create_simu_folder()
 Step 1.2. More...
 
tuple sub = plt.subplot(230+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 img = plot_common.map2D(sub,map2D_type,map2D_corfield,0)
 
tuple eq = plot_common.equator(sub,'black')
 
tuple map2D_python = interface_read_map2D.read(map2D_type,map2D_path)
 Step 2. More...
 
tuple map2D_diff = module_map2D.rel_diff(map2D_python,map2D_corfield)
 DIFFERENCE. More...
 

Variable Documentation

string db_corfield_map2D_smap.axis_color = 'k'
string db_corfield_map2D_smap.axis_name_x = 'Carrington Longitude'
string db_corfield_map2D_smap.axis_name_y = 'Latitude'
int db_corfield_map2D_smap.axis_size = 10
tuple db_corfield_map2D_smap.connect_tool_path = os.getcwd()
tuple db_corfield_map2D_smap.corfield_exec = os.path.join(module_corona.path_corfield_pfss,'test_read_map')
tuple db_corfield_map2D_smap.corfield_output_file = os.path.join('dir.output',map2D_type.lower() + corfield_output_name + '_txt_tmp.dat')
tuple db_corfield_map2D_smap.corfield_output_name = ('-lat_%s')
list db_corfield_map2D_smap.date_max = interface_read_map2D.info[map2D_type]
tuple db_corfield_map2D_smap.date_min = datetime.datetime(2019,1,1)
list db_corfield_map2D_smap.dates = [ date_min + i*datetime.timedelta(days=1) for i in range(0,ndays+1) ]
db_corfield_map2D_smap.delta = date_max-date_min
int db_corfield_map2D_smap.dphi = 40
tuple db_corfield_map2D_smap.eq = plot_common.equator(sub,'black')
tuple db_corfield_map2D_smap.fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)
int db_corfield_map2D_smap.fig_resolution = 100
tuple db_corfield_map2D_smap.fig_size = (25,15)
tuple db_corfield_map2D_smap.folder_test = os.path.join('dir.test','COMPARE_DB_CORFIELD')
string db_corfield_map2D_smap.help_message
Initial value:
1 = """
2 Usage :\n
3 ./db_corfield_map2D_smap.py
4 [-h] : brief help
5 [d1 d2] : from date d1 to d2 with format YYYYMMDD
6 """
int db_corfield_map2D_smap.i = -1

Step 1.

Loop for all map2D_types

tuple db_corfield_map2D_smap.ifile = open(os.path.join(simu_path,corfield_output_file))
tuple db_corfield_map2D_smap.img = plot_common.map2D(sub,map2D_type,map2D_corfield,0)
tuple db_corfield_map2D_smap.map2D_corfield = ascii_map2D.read(ifile)
tuple db_corfield_map2D_smap.map2D_diff = module_map2D.rel_diff(map2D_python,map2D_corfield)

DIFFERENCE.

list db_corfield_map2D_smap.map2D_numsub = [1,4]
tuple db_corfield_map2D_smap.map2D_path = interface_read_map2D.date2map(map2D_type,date)
tuple db_corfield_map2D_smap.map2D_python = interface_read_map2D.read(map2D_type,map2D_path)

Step 2.

Python reading.

string db_corfield_map2D_smap.map2D_time = 'smap'
string db_corfield_map2D_smap.map2D_type = 'NSO'

Step 1.1.

Initialisation Corfield map2D outputs filename

list db_corfield_map2D_smap.map2D_types = ['NSO_DATE','ADAPT_DATE']
db_corfield_map2D_smap.ndays = delta.days
tuple db_corfield_map2D_smap.simu_path = model_corfield_pfss.create_simu_folder()

Step 1.2.

Corfield reading/writing. Corfield simulation name

tuple db_corfield_map2D_smap.sub = plt.subplot(230+map2D_numsub[i])

Step 1.3.

Plot map2D.

db_corfield_map2D_smap.tmp = date_min
list db_corfield_map2D_smap.xtick_label = [ str(int(xtick_loc[i])) for i in range(0,len(xtick_loc)) ]
list db_corfield_map2D_smap.xtick_loc = [ i*dphi for i in range(0,int(360./dphi)+1) ]
list db_corfield_map2D_smap.ytick_label = [ str(int(ytick_loc[i])) for i in range(0,len(ytick_loc)) ]
list db_corfield_map2D_smap.ytick_loc = [ -90,-80,-60,-40,-20,0,20,40,60,80,90 ]