-import math
-#Bezeichnungen:
-#r = relative Luftfeuchte
-#T = Temperatur in C
-#TK = Temperatur in Kelvin (TK = T + 273.15)
-#TD = Taupunkttemperatur in C
-#DD = Dampfdruck in hPa
-#SDD = Saettigungsdampfdruck in hPa
-
-#Parameter:
-#a = 7.5, b = 237.3 fuer T >= 0
-#a = 7.6, b = 240.7 fuer T < 0 ueber Wasser (Taupunkt)
-#a = 9.5, b = 265.5 fuer T < 0 ueber Eis (Frostpunkt)
-
-Rs = 8314.3 #J/(kmol*K) (universelle Gaskonstante)
-mw = 18.016 #kg/kmol (Molekulargewicht des Wasserdampfes)
-#AF = absolute Feuchte in g Wasserdampf pro m3 Luft
-
-#Formeln:
-
-# SDD(T) = 6.1078 * 10^((a*T)/(b+T))
-# DD(r,T) = r/100 * SDD(T)
-# r(T,TD) = 100 * SDD(TD) / SDD(T)
-# TD(r,T) = b*v/(a-v) mit v(r,T) = log10(DD(r,T)/6.1078)
-# AF(r,TK) = 10^5 * mw/R* * DD(r,T)/TK; AF(TD,TK) = 10^5 * mw/R* * SDD(TD)/TK
-
-
-a=7.5
-b=237.3
-
-#a=r/100.0*SDD(T)/(T+273.15)*mw/Rs*10000
-#r=a*100/(SDD(T)/(T+273.15)*mw/Rs*10000)
-
-def SDD(T):
- return 6.1078 * 10**((a*T)/(b+T))
-
-#def DD(r,T):
-# return r/100 * SDD(T)
-
-#def _r(T,TD):
-# return 100 * SDD(TD) / SDD(T)
-
-#def _v(r,T):
-# return math.log10(DD(r,T)/6.1078)
-
-#def TD(r,T):
-# return b*v/(a-_v(r,T))
-
-
-#def AF(r,T):
-# return 10000 * mw/Rs * DD(r,T)/(T+273.15) #; AF(TD,TK) = 10^5 * mw/R* * SDD(TD)/TK
-
-def AF1(r,T):
- return r/100.0*SDD(T)/(T+273.15)*mw/Rs*10000
-
-def RF1(a,T):
- return a*100/(SDD(T)/(T+273.15)*mw/Rs*10000)
-#print(0.5*SDD(20)/(20+273.15)*mw/Rs*10000)
-d=[]
-for i in range(11):
- s=""
- r=[]
- for j in range(6):
- r.append((i*10)-RF1(AF1(i*10,j*10),j*10+0.5))
- s=s+"%0.2f\t" % (r[len(r)-1])
- print(s)
- d.append(r)
-
-for i in range(11):
- s=""
- for j in range(5):
- s=s+"%0.4f\t" % (d[i][j]-d[i][j+1])
- print(s)
-
+import math\r
+#Bezeichnungen:\r
+#r = relative Luftfeuchte\r
+#T = Temperatur in C\r
+#TK = Temperatur in Kelvin (TK = T + 273.15)\r
+#TD = Taupunkttemperatur in C\r
+#DD = Dampfdruck in hPa\r
+#SDD = Saettigungsdampfdruck in hPa\r
+\r
+#Parameter:\r
+#a = 7.5, b = 237.3 fuer T >= 0\r
+#a = 7.6, b = 240.7 fuer T < 0 ueber Wasser (Taupunkt)\r
+#a = 9.5, b = 265.5 fuer T < 0 ueber Eis (Frostpunkt)\r
+\r
+Rs = 8314.3 #J/(kmol*K) (universelle Gaskonstante)\r
+mw = 18.016 #kg/kmol (Molekulargewicht des Wasserdampfes)\r
+#AF = absolute Feuchte in g Wasserdampf pro m3 Luft\r
+\r
+#Formeln:\r
+\r
+# SDD(T) = 6.1078 * 10^((a*T)/(b+T))\r
+# DD(r,T) = r/100 * SDD(T)\r
+# r(T,TD) = 100 * SDD(TD) / SDD(T)\r
+# TD(r,T) = b*v/(a-v) mit v(r,T) = log10(DD(r,T)/6.1078)\r
+# AF(r,TK) = 10^5 * mw/R* * DD(r,T)/TK; AF(TD,TK) = 10^5 * mw/R* * SDD(TD)/TK\r
+\r
+\r
+a=7.5\r
+b=237.3\r
+\r
+#a=r/100.0*SDD(T)/(T+273.15)*mw/Rs*10000\r
+#r=a*100/(SDD(T)/(T+273.15)*mw/Rs*10000)\r
+\r
+def SDD(T):\r
+ return 6.1078 * 10**((a*T)/(b+T))\r
+\r
+#def DD(r,T):\r
+# return r/100 * SDD(T)\r
+\r
+#def _r(T,TD):\r
+# return 100 * SDD(TD) / SDD(T)\r
+\r
+#def _v(r,T):\r
+# return math.log10(DD(r,T)/6.1078)\r
+ \r
+#def TD(r,T):\r
+# return b*v/(a-_v(r,T)) \r
+\r
+ \r
+#def AF(r,T):\r
+# return 10000 * mw/Rs * DD(r,T)/(T+273.15) #; AF(TD,TK) = 10^5 * mw/R* * SDD(TD)/TK\r
+\r
+def AF1(r,T):\r
+ return r/100.0*SDD(T)/(T+273.15)*mw/Rs*10000\r
+\r
+def RF1(a,T):\r
+ return a*100/(SDD(T)/(T+273.15)*mw/Rs*10000) \r
+#print(0.5*SDD(20)/(20+273.15)*mw/Rs*10000)\r
+d=[]\r
+for i in range(11):\r
+ s=""\r
+ r=[]\r
+ for j in range(6):\r
+ r.append((i*10)-RF1(AF1(i*10,j*10),j*10+0.5))\r
+ s=s+"%0.2f\t" % (r[len(r)-1])\r
+ print(s)\r
+ d.append(r)\r
+ \r
+for i in range(11):\r
+ s=""\r
+ for j in range(5):\r
+ s=s+"%0.4f\t" % (d[i][j]-d[i][j+1])\r
+ print(s)\r
+ \r
tm3d / owSlave2.git / blobdiff