) else: mindtime, maxdtime, i_mindepth, i_maxdepth, model, sites = get_neracoos_ctl( inputfilename) #get input from input file model = 'aanderaa' sdtime_n = date2num(mindtime) - date2num(dt.datetime( 1858, 11, 17, 0, 0)) #get number type of start time edtime_n = date2num(maxdtime) - date2num(dt.datetime( 1858, 11, 17, 0, 0)) #get number type of end time for index_site in range(len(sites)): period_str, temp_list, depth_list = [], [], [] url = 'http://neracoos.org:8080/opendap/' + sites[ index_site] + '/' + sites[index_site] + '.' + model + '.historical.nc?' id_s, id_e0, id_max_url, maxtime, mintime = get_id_s_id_e_id_max_url( url, sdtime_n, edtime_n) if mintime == '': histvsreal = '1' #"histvsreal" can help us judge if this site has historical data. url = 'http://neracoos.org:8080/opendap/' + sites[ index_site] + '/' + sites[ index_site] + '.' + model + '.realtime.nc?' id_s, id_e0, id_max_url, maxtime, mintime = get_id_s_id_e_id_max_url( url, sdtime_n, edtime_n) print 'realtime from ' + str( num2date(date2num(dt.datetime(1858, 11, 17, 0, 0)) + mintime) ) + 'to' + str( num2date(date2num(dt.datetime(1858, 11, 17, 0, 0)) + maxtime)) else: histvsreal = '' print 'historical from ' + str(
from get_neracoos_ctl import get_neracoos_ctl_py inputfilename='./get_neracoos_ctl.txt' if inputfilename[-2:]=='py': mindtime,maxdtime,i_mindepth,i_maxdepth,model,sites=get_neracoos_ctl_py() else: mindtime,maxdtime,i_mindepth,i_maxdepth,model,sites=get_neracoos_ctl(inputfilename) #get input from input file model='met' #hard code sdtime_n=date2num(mindtime)-date2num(dt.datetime(1858, 11, 17, 0, 0)) #get number type of start time edtime_n=date2num(maxdtime)-date2num(dt.datetime(1858, 11, 17, 0, 0)) #get number type of end time for index_site in range(len(sites)): url='http://neracoos.org:8080/opendap/'+sites[index_site]+'/'+sites[index_site]+'.'+model+'.historical.nc?' id_s,id_e0,id_max_url,maxtime,mintime=get_id_s_id_e_id_max_url(url,sdtime_n,edtime_n) # 'maxtime',the max time in this url, "id_s",the index of start time we want if mintime=='': histvsreal='1' #"histvsreal" can help us judge if this site has historical data. url='http://neracoos.org:8080/opendap/'+sites[index_site]+'/'+sites[index_site]+'.'+model+'.realtime.nc?' id_s,id_e0,id_max_url,maxtime,mintime=get_id_s_id_e_id_max_url(url,sdtime_n,edtime_n) # 'maxtime',the max time in this url, "id_s",the index of start time we want print 'realtime from '+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+mintime))+'to'+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+maxtime)) else: print 'historical from '+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+mintime))+'to'+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+maxtime)) histvsreal='' if id_e0<>'': (period_str,wind_all)=get_neracoos_wind_data(url,id_s,id_e0,id_max_url) #get data from web neracoos df = DataFrame(np.array(wind_all),index=period_str,columns=['wind speed(m/s)','direction(degree)']) else: print "According to your input, there is no data here" if histvsreal<>'1': if maxtime<edtime_n: #make sure if we need a realtime data
import numpy as np from pandas import * pydir='../' sys.path.append(pydir) from neracoos_def import get_neracoos_ctl,get_id_s_id_e_id_max_url,get_neracoos_current_data inputfilename='./get_neracoos_ctl.txt' mindtime,maxdtime,i_mindepth,i_maxdepth,model,sites=get_neracoos_ctl(inputfilename) #get input from input file model='aanderaa' sdtime_n=date2num(mindtime)-date2num(dt.datetime(1858, 11, 17, 0, 0)) #get number type of start time edtime_n=date2num(maxdtime)-date2num(dt.datetime(1858, 11, 17, 0, 0)) #get number type of end time for index_site in range(len(sites)): period_str,temp_list,depth_list=[],[],[] url='http://neracoos.org:8080/opendap/'+sites[index_site]+'/'+sites[index_site]+'.'+model+'.historical.nc?' id_s,id_e0,id_max_url,maxtime,mintime=get_id_s_id_e_id_max_url(url,sdtime_n,edtime_n) if mintime=='': histvsreal='1' #"histvsreal" can help us judge if this site has historical data. url='http://neracoos.org:8080/opendap/'+sites[index_site]+'/'+sites[index_site]+'.'+model+'.realtime.'+depths[index_site]+'m.nc?' id_s,id_e0,id_max_url,maxtime,mintime=get_id_s_id_e_id_max_url(url,sdtime_n,edtime_n) print 'realtime from '+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+mintime))+'to'+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+maxtime)) else: histvsreal='' print 'historical from '+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+mintime))+'to'+str(num2date(date2num(dt.datetime(1858, 11, 17, 0, 0))+maxtime)) if id_e0<>'': (period_str,current_all)=get_neracoos_current_data(url,id_s,id_e0,id_max_url) #get data from web neracoos df = DataFrame(np.array(current_all),index=period_str,columns=['current','direction','u','v' ]) else: print "According to your input, there is no data here" if histvsreal<>'1': if maxtime<edtime_n: #make sure if we need a realtime data
print 'No right site , Please check website and your input data, ' for index_site in range(len(site_d)): if model == 'sbe16': url = 'http://neracoos.org:8080/opendap/' + site_d[ index_site] + '/' + site_d[ index_site] + '.' + model + '.historical.nc?' else: url = 'http://neracoos.org:8080/opendap/' + site_d[ index_site] + '/' + site_d[ index_site] + '.' + model + '.historical.' + depths[ index_site] + 'm.nc?' id_s, id_e0, id_max_url, maxtime, mintime = get_id_s_id_e_id_max_url( url, sdtime_n, edtime_n ) # 'maxtime',the max time in this url, "id_s",the index of start time we want if mintime == '': histvsreal = '1' #"histvsreal" can help us judge if this site has historical data. url = 'http://neracoos.org:8080/opendap/' + site_d[ index_site] + '/' + site_d[ index_site] + '.' + model + '.realtime.' + depths[ index_site] + 'm.nc?' id_s, id_e0, id_max_url, maxtime, mintime = get_id_s_id_e_id_max_url( url, sdtime_n, edtime_n ) # 'maxtime',the max time in this url, "id_s",the index of start time we want #get id of start time and end time,max time id of this url, print 'realtime from ' + str( num2date(date2num(dt.datetime(1858, 11, 17, 0, 0)) + mintime) ) + 'to' + str( num2date(date2num(dt.datetime(1858, 11, 17, 0, 0)) + maxtime))