我從一個nc文件中讀取了數(shù)組形式的變量?，F(xiàn)在我已經(jīng)用所有這些變量的第一個元素進行了計算并創(chuàng)建了一個新變量。我想為初始數(shù)組中的每個元素重復相同的一組計算，而不更改我考慮單點的計算代碼。我已經(jīng)嘗試過zip，nditer但在這兩種情況下，變量中的 if 語句a都將更改為.any()or .all()。我也不能這樣做，因為我希望 if 語句只考慮一個點而不是整個數(shù)組。T = AD06_ALL_OMNI.variables['A_TEMP'][:][0] REL_HUM = AD06_ALL_OMNI.variables['HUMIDITY'][:][0]AIR_PRES = AD06_ALL_OMNI.variables['A_PRES'][:][0]a = T-29.65#masking of values so that division by 0 is avoidedif a!=0.0: exponent1 = math.exp(17.67*T-0.16/a) q = REL_HUM*exponent1/(26.3*AIR_PRES) deltaq = 0.98*qs-qprint (deltaq)我需要a對每個點進行計算，以便從 T,和對應的點deltaq中找出相同的點。所有變量的大小相同（一維數(shù)組）。請幫忙！REL_HUMAIR_PRES
1 回答
波斯汪
TA貢獻1811條經(jīng)驗 獲得超4個贊
for var in range(0, length(AD06_ALL_OMNI.variables['A_TEMP'][:])-1):
T = AD06_ALL_OMNI.variables['A_TEMP'][:][var]
REL_HUM = AD06_ALL_OMNI.variables['HUMIDITY'][:][var]
AIR_PRES = AD06_ALL_OMNI.variables['A_PRES'][:][var]
a = T-29.65
#masking of values so that division by 0 is avoided
count = 0
for element in a:
if element!=0.0:
exponent1 = math.exp(17.67*T[count] -0.16/element)
q = REL_HUM[count]*exponent1/(26.3*AIR_PRES[count] )
deltaq = 0.98*qs-q
print (deltaq)
count = count + 1
假設所有數(shù)組的長度相同（如果沒有相等長度的氣壓、氣溫和濕度，這沒有意義），您可以使用循環(huán)遍歷 a 的所有值，檢查每個值都為 0，并計算并打印每個點的 deltaq。我希望這有幫助。
添加回答
舉報
0/150
提交
取消