stest/fdiff.py

   1 import math
   2 #Bezeichnungen:
   3 #r = relative Luftfeuchte
   4 #T = Temperatur in C
   5 #TK = Temperatur in Kelvin (TK = T + 273.15)
   6 #TD = Taupunkttemperatur in C
   7 #DD = Dampfdruck in hPa
   8 #SDD = Saettigungsdampfdruck in hPa
   9
  10 #Parameter:
  11 #a = 7.5, b = 237.3 fuer T >= 0
  12 #a = 7.6, b = 240.7 fuer T < 0 ueber Wasser (Taupunkt)
  13 #a = 9.5, b = 265.5 fuer T < 0 ueber Eis (Frostpunkt)
  14
  15 Rs = 8314.3 #J/(kmol*K) (universelle Gaskonstante)
  16 mw = 18.016 #kg/kmol (Molekulargewicht des Wasserdampfes)
  17 #AF = absolute Feuchte in g Wasserdampf pro m3 Luft
  18
  19 #Formeln:
  20
  21 #    SDD(T) = 6.1078 * 10^((a*T)/(b+T))
  22 #    DD(r,T) = r/100 * SDD(T)
  23 #    r(T,TD) = 100 * SDD(TD) / SDD(T)
  24 #    TD(r,T) = b*v/(a-v) mit v(r,T) = log10(DD(r,T)/6.1078)
  25 #    AF(r,TK) = 10^5 * mw/R* * DD(r,T)/TK; AF(TD,TK) = 10^5 * mw/R* * SDD(TD)/TK
  26
  27
  28 a=7.5
  29 b=237.3
  30
  31 #a=r/100.0*SDD(T)/(T+273.15)*mw/Rs*10000
  32 #r=a*100/(SDD(T)/(T+273.15)*mw/Rs*10000)
  33
  34 def SDD(T):
  35         return 6.1078 * 10**((a*T)/(b+T))
  36
  37 #def DD(r,T):
  38 #       return r/100 * SDD(T)
  39
  40 #def _r(T,TD):
  41 #       return 100 * SDD(TD) / SDD(T)
  42
  43 #def _v(r,T):
  44 #       return math.log10(DD(r,T)/6.1078)
  45
  46 #def TD(r,T):
  47 #       return b*v/(a-_v(r,T))
  48
  49
  50 #def AF(r,T):
  51 #       return 10000 * mw/Rs * DD(r,T)/(T+273.15) #; AF(TD,TK) = 10^5 * mw/R* * SDD(TD)/TK
  52
  53 def AF1(r,T):
  54         return r/100.0*SDD(T)/(T+273.15)*mw/Rs*10000
  55
  56 def RF1(a,T):
  57         return a*100/(SDD(T)/(T+273.15)*mw/Rs*10000)
  58 #print(0.5*SDD(20)/(20+273.15)*mw/Rs*10000)
  59 d=[]
  60 for i in range(11):
  61         s=""
  62         r=[]
  63         for j in range(6):
  64                 r.append((i*10)-RF1(AF1(i*10,j*10),j*10+0.5))
  65                 s=s+"%0.2f\t" % (r[len(r)-1])
  66         print(s)
  67         d.append(r)
  68
  69 for i in range(11):
  70         s=""
  71         for j in range(5):
  72                 s=s+"%0.4f\t" % (d[i][j]-d[i][j+1])
  73         print(s)
  74