SWDatabase
 All Classes Namespaces Files Functions Variables Enumerations Enumerator Macros Pages
Functions | Variables
test_map2D_nso_gong Namespace Reference

Functions

def get_nso_gong_jpg
 Find in NSO website the jpg corresponding to the date. More...
 

Variables

string help_message
 
string map2D_type = 'NSO'
 
string mag_table = 'OBS_MAPS_NSO'
 
 date_max = module_timespace.today_mid
 
tuple date_min = module_timespace.today_mid-datetime.timedelta(days=2)
 
 tmp = date_min
 
tuple folder_test = os.path.join('dir.test','TEST')
 
list files_test = []
 
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 fig_size = (12,6)
 
int fig_resolution = 100
 
 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 time_in = module_timespace.Time()
 
tuple fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)
 
tuple mag_db = database_map2D.time2map(mag_table,time_in)
 NSO DB. More...
 
tuple mag_data = interface_read_map2D.read(mag_table,os.path.join(mag_db['dirname'],mag_db['filename']))
 
tuple min_value = np.nanmin(mag_data.data)
 
tuple max_value = np.nanmax(mag_data.data)
 
tuple mean_value = np.nanmean(mag_data.data)
 
tuple sub = plt.subplot(222)
 
tuple img = plot_common.map2D(sub,mag_table,mag_data,mag_db['vmin'],mag_db['vmax'],1)
 
tuple eq = plot_common.equator(sub,'black')
 
tuple remote_path = get_nso_gong_jpg(date_in)
 NSO REF Download and plot reference image. More...
 
tuple imgplot = sub.imshow(img)
 
tuple filepath = os.path.join(folder_test,'test_' + date_in.strftime('%Y%m%dT%H%M') + '.png')
 
tuple command = ('eog ' + '%s '*len(files_test))
 

Function Documentation

def test_map2D_nso_gong.get_nso_gong_jpg (   date_in)

Find in NSO website the jpg corresponding to the date.

Returns
Complete url of jpg file.

Variable Documentation

string test_map2D_nso_gong.axis_color = 'k'
string test_map2D_nso_gong.axis_name_x = 'Carrington Longitude'
string test_map2D_nso_gong.axis_name_y = 'Latitude'
int test_map2D_nso_gong.axis_size = 10
tuple test_map2D_nso_gong.command = ('eog ' + '%s '*len(files_test))
list test_map2D_nso_gong.date_max = module_timespace.today_mid
tuple test_map2D_nso_gong.date_min = module_timespace.today_mid-datetime.timedelta(days=2)
list test_map2D_nso_gong.dates = [ date_min + i*datetime.timedelta(days=1) for i in range(0,ndays+1) ]
test_map2D_nso_gong.delta = date_max-date_min
int test_map2D_nso_gong.dphi = 40
tuple test_map2D_nso_gong.eq = plot_common.equator(sub,'black')
tuple test_map2D_nso_gong.fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)
int test_map2D_nso_gong.fig_resolution = 100
tuple test_map2D_nso_gong.fig_size = (12,6)
tuple test_map2D_nso_gong.filepath = os.path.join(folder_test,'test_' + date_in.strftime('%Y%m%dT%H%M') + '.png')
list test_map2D_nso_gong.files_test = []
tuple test_map2D_nso_gong.folder_test = os.path.join('dir.test','TEST')
string test_map2D_nso_gong.help_message
Initial value:
1 = """
2 Test NSO/GONG reading routine.
3 For a given date, read the NSO magnetogram and compare with a reference image taken from: https://gong.nso.edu/data/magmap/QR/bq1\n
4 Usage :\n
5 ./test_map2D_nso_read.py
6 [-h] : brief help
7 [d1 d2] : from date d1 to d2 with format YYYYMMDD
8 """
tuple test_map2D_nso_gong.img = plot_common.map2D(sub,mag_table,mag_data,mag_db['vmin'],mag_db['vmax'],1)
tuple test_map2D_nso_gong.imgplot = sub.imshow(img)
tuple test_map2D_nso_gong.mag_data = interface_read_map2D.read(mag_table,os.path.join(mag_db['dirname'],mag_db['filename']))
tuple test_map2D_nso_gong.mag_db = database_map2D.time2map(mag_table,time_in)

NSO DB.

string test_map2D_nso_gong.mag_table = 'OBS_MAPS_NSO'
string test_map2D_nso_gong.map2D_type = 'NSO'
tuple test_map2D_nso_gong.max_value = np.nanmax(mag_data.data)
tuple test_map2D_nso_gong.mean_value = np.nanmean(mag_data.data)
tuple test_map2D_nso_gong.min_value = np.nanmin(mag_data.data)
test_map2D_nso_gong.ndays = delta.days
tuple test_map2D_nso_gong.remote_path = get_nso_gong_jpg(date_in)

NSO REF Download and plot reference image.

tuple test_map2D_nso_gong.sub = plt.subplot(222)
tuple test_map2D_nso_gong.time_in = module_timespace.Time()
test_map2D_nso_gong.tmp = date_min
list test_map2D_nso_gong.xtick_label = [ str(int(xtick_loc[i])) for i in range(0,len(xtick_loc)) ]
list test_map2D_nso_gong.xtick_loc = [ i*dphi for i in range(0,int(360./dphi)+1) ]
list test_map2D_nso_gong.ytick_label = [ str(int(ytick_loc[i])) for i in range(0,len(ytick_loc)) ]
list test_map2D_nso_gong.ytick_loc = [ -90,-80,-60,-40,-20,0,20,40,60,80,90 ]