)
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))
