本文整理汇总了Python中Data.dalecData方法的典型用法代码示例。如果您正苦于以下问题:Python Data.dalecData方法的具体用法?Python Data.dalecData怎么用?Python Data.dalecData使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Data
的用法示例。
在下文中一共展示了Data.dalecData方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: ACM
# 需要导入模块: import Data [as 别名]
# 或者: from Data import dalecData [as 别名]
import numpy as np
from scipy import misc
import matplotlib.pyplot as plt
import Data as D
d=D.dalecData()
from matplotlib.mlab import csv2rec
dat=csv2rec('NiwotRidgeModelDat.csv')
def ACM(Cf,i):
q=d.a_3-d.a_4
L=(Cf/111.)
e_0=((d.a_7*L**2)/(L**2+d.a_9))
g_c=((np.abs(float(d.phi_d)))**(d.a_10))/(0.5*float(max(dat.ta_fill[(dat.jd==d.D[i])]-min(dat.ta_fill[(dat.jd==d.D[i])])+d.a_6*float(d.R_tot))))
p=(((d.a_1*d.N*L)/g_c)*np.exp(d.a_8*float(max(dat.ta_fill[(dat.jd==d.D[i])]))))
C_i=(0.5*(d.C_a+q-p+np.sqrt((d.C_a+q-p)**2-4*(d.C_a*q-p*d.a_3))))
delta=(-0.408*np.cos(((360*(float(d.D[i])+10)*np.pi)/(365*180))))
s=(24*np.arccos((-np.tan(d.bigdelta)*np.tan(delta)))/np.pi)
def GPPfunc(Cf):
#L=(Cf/111.)
#e_0=((d.a_7*L**2)/(L**2+d.a_9))
#GPP=((e_0*float(d.I[i])*g_c*(d.C_a-C_i)*(d.a_2*s+d.a_5))/(e_0*float(d.I[i])+g_c*(d.C_a-C_i))) #0.2*float(d.I[i])*(1-np.exp(-0.5*Cf/111.))
#GPP_diff=((2*float(d.I[i])*d.a_7*d.a_9*((111.0*g_c*(d.C_a-C_i))**2)*(d.a_2*s+d.a_5)*Cf)/(((g_c*(d.C_a-C_i)+d.a_7*float(d.I[i]))*Cf**2+d.a_9*(111.0**2)*g_c*(d.C_a-C_i))**2)) #(0.2*0.5*float(d.I[i])/111.)*np.exp(-0.5*Cf/111.)
GPP= (d.a_7*(Cf**2)*float(d.I[i])*g_c*(d.C_a-C_i)*(d.a_2*s+d.a_5))/((d.a_7*float(d.I[i])+g_c*(d.C_a-C_i))*(Cf**2)+d.a_9*(111.**2)*g_c*(d.C_a-C_i))
return GPP
def GPPdifffunc(Cf):