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