import warnings
warnings.filterwarnings("ignore")
print("start program:", datetime.datetime.now())
# intsance the windfield model
georgiou = Georgiou.Georgiou()
starttime = datetime.datetime.now()
# getting parameter
print("getting parameters")
parameter = parameter.SiteInfo()
begYear = parameter.begYear()
endYear = parameter.endYear()
totalYear = endYear - begYear + 1
radiusInflu = parameter.radiusInflu()
returnPeriod = parameter.returnPeriod()
dictInfo = parameter.allWeatherStaionInfo()
#dictInfo = parameter.allWindFarmInfo()
deltaT1 = 10
deltaT2 = 5
for iKey in dictInfo.keys():
### read typhoon key parameter
inputFileName = r"data_allTyKeyParameter/" + iKey + "KeyParameters.csv"
print("reading data from", inputFileName)
dataset = pd.read_csv(inputFileName, header=None, sep=',')
dataset = np.array(dataset)
m, n = np.shape(dataset)
allTyNum = dataset[:, 0]
allDate = dataset[:, 1]
plt.ylabel(u'风速(m/s)', fontproperties=myfont,fontsize=fontsize)
plt.show()
plt.tight_layout()
fig.savefig(figName)
plt.close()
return None
if __name__ == '__main__':
# get parameter
print("getting parameter")
parameter = parameter.SiteInfo()
begYear = parameter.begYear()
endYear = parameter.endYear()
totalYear = endYear-begYear+1
radiusInflu = parameter.radiusInflu() # influence radius ,unit:KM
returnPeriod = parameter.returnPeriod()
allWeatherStaionInfo = parameter.allWeatherStaionInfo()
allWindFarmInfo = parameter.allWindFarmInfo()
specialWindFarmInfo = parameter.specialWindFarmInfo()
allObsInfo = parameter.allObsInfo()
# plotTemporalSpatial
plotTemporal = plotTemporalSpatial()
#plt0 = plotTemporal.plotSiteInterannualVariability(allWeatherStaionInfo,begYear,endYear)
#plt1 = plotTemporal.plotSiteSeasonalVariability(allWeatherStaionInfo,begYear,endYear)
#plt2 = plotTemporal.plotSiteInterannualVariabilitySubplot(specialWindFarmInfo,begYear,endYear)
#plt3 = plotTemporal.plotSiteSeasonalVariability(specialWindFarmInfo,begYear,endYear)
#plt4 = plotTemporal.plotSiteIntensitySpatialDistribution(specialWindFarmInfo,begYear,endYear,minLat=16,maxLat=30,minLon=106,maxLon=128,figName='SiteIntensitySpatialDistribution.png')
#plt5 = plotTemporal.plotSiteIntensityFrequencyDistribution(specialWindFarmInfo,begYear,endYear,figName='SiteIntensityFrequencyDistribution.png')
# plotVmax
plotVmax = plotVmax()
