def example():
pc = csmPowerCurve()
drivetrain = csmDriveEfficiency(1)
pc.drivetrain = drivetrain
# initialize input parameters
pc.hubHt = 90.0
pc.ratedPower = 5000.0
pc.maxTipSpd = 80
pc.rotorDiam = 126
pc.maxCp = 0.488
pc.maxTipSpdRatio = 7.525
pc.cutInWS = 3.0
pc.cutOutWS = 25.0
pc.altitude = 0.0
pc.airDensity = 0.0
pc.reg2pt5slope = 0.05
pc.compute(pc.drivetrain, pc.hubHt, pc.ratedPower,pc.maxTipSpd,pc.rotorDiam, \
pc.maxCp, pc.maxTipSpdRatio, pc.cutInWS, pc.cutOutWS, \
pc.altitude, pc.airDensity)
print 'Rated Speed: %9.3f mps' % pc.getRatedWindSpeed()
print 'Rated RPM: %9.3f rpm' % pc.getRatedRotorSpeed()
pcurve = np.array(pc.getPowerCurve())
print 'Power Curve:'
for i in xrange(0,pcurve.shape[1]):
print (pcurve[0][i], pcurve[1][i])
'''# plot
def example2():
pc = csmPowerCurve()
drivetrain = csmDriveEfficiency(1)
pc.drivetrain = drivetrain
# initialize input parameters
pc.hubHt = 90.0
pc.maxTipSpd = 80
pc.rotorDiam = 126
pc.maxCp = 0.488
pc.maxTipSpdRatio = 7.525
pc.cutInWS = 3.0
pc.cutOutWS = 25.0
pc.altitude = 0.0
pc.airDensity = 0.0
pc.reg2pt5slope = 0.05
pc.ratedPower = 750.0
pc.compute(pc.drivetrain, pc.hubHt, pc.ratedPower,pc.maxTipSpd,pc.rotorDiam, \
pc.maxCp, pc.maxTipSpdRatio, pc.cutInWS, pc.cutOutWS, \
pc.altitude, pc.airDensity)
print
print 'RatedSpd %9.3f mps' % pc.getRatedWindSpeed()
print 'RatedRPM %9.3f ' % pc.getRatedRotorSpeed()
print 'Power curve values'
pcurve = pc.getPowerCurve()
print pcurve
pc.ratedPower = 1000.0
pc.compute(pc.drivetrain, pc.hubHt, pc.ratedPower,pc.maxTipSpd,pc.rotorDiam, \
pc.maxCp, pc.maxTipSpdRatio, pc.cutInWS, pc.cutOutWS, \
pc.altitude, pc.airDensity)
print
print 'RatedSpd %9.3f mps' % pc.getRatedWindSpeed()
print 'RatedRPM %9.3f ' % pc.getRatedRotorSpeed()
print 'Power curve values'
pcurve = pc.getPowerCurve()
print pcurve
pc.ratedPower = 2000.0
pc.compute(pc.drivetrain, pc.hubHt, pc.ratedPower,pc.maxTipSpd,pc.rotorDiam, \
pc.maxCp, pc.maxTipSpdRatio, pc.cutInWS, pc.cutOutWS, \
pc.altitude, pc.airDensity)
print
print 'RatedSpd %9.3f mps' % pc.getRatedWindSpeed()
print 'RatedRPM %9.3f ' % pc.getRatedRotorSpeed()
print 'Power curve values'
pcurve = pc.getPowerCurve()
print pcurve
