Difference between revisions of "Stationsseite Lebensministerium"

From hms.sternhell.at
Jump to: navigation, search
(Kalibration 28./29. August)
(AT_LebMin_2 - Lightrouter, Mark 2.3l)
Line 22: Line 22:
 
ab 8. August 2011
 
ab 8. August 2011
  
15. August 2015 - Ausfall. Vermutlich Speicher voll oder UMTS Kontingent zu erneuern!
+
== Kontrolle 22. Jan. 2016 ==
 +
 
 +
Die Daten von 2015 sind komplett und auch alles seit Anfang 2016.
 +
 
 +
Der LightRouter ist gelaufen und hat aufgezeichnet. Nur der UMTS-WLAN Router war auf "Standby". Router aus, Daten kopiert,
 +
2014 und 2015 auf dem Router gelöscht damit er nicht voll wird.
 +
 
 +
    UMTS-Router Einschalten; Warten auf WLAN; Router Starten;
 +
    Heartbeat; alles scheint OK.
 +
 
 +
Der "Ausfall" vom 15. August war nur ein Verbindungsverlust. Daten wurden weiter gewonnen.
 +
 
 +
--[[User:GuentherWuchterl|GW]] ([[User talk:GuentherWuchterl|talk]]) 22:51, 23 February 2016 (CET)
 +
 
 +
'''15. August 2015''' - Ausfall. Vermutlich Speicher voll oder UMTS Kontingent zu erneuern!
  
 
==Kalibration 28./29. August==
 
==Kalibration 28./29. August==

Revision as of 22:51, 23 February 2016

Metadata

Koordinaten:
48° 12' 33,35" N
16° 22' 54,32" O
209 m Höhe

Sensor:
Mark 2.3L, S/N ....028
Liefert Daten seit dem 9. Aug. 2011. ; validiert --GuentherWuchterl 19:35, 30. Aug. 2011 (UTC)
--Nero2401 19:36, 30. Aug. 2011 (UTC)

Messdaten

Nahezu-Echtzeit-Daten finden sich auf http://www.lightdata.at.vu


AT_LebMin_2 - Lightrouter, Mark 2.3l

Station: AT_LebMin Sensor ID: 918028 ab 8. August 2011

Kontrolle 22. Jan. 2016

Die Daten von 2015 sind komplett und auch alles seit Anfang 2016.

Der LightRouter ist gelaufen und hat aufgezeichnet. Nur der UMTS-WLAN Router war auf "Standby". Router aus, Daten kopiert, 2014 und 2015 auf dem Router gelöscht damit er nicht voll wird.

    UMTS-Router Einschalten; Warten auf WLAN; Router Starten;
    Heartbeat; alles scheint OK.

Der "Ausfall" vom 15. August war nur ein Verbindungsverlust. Daten wurden weiter gewonnen.

--GW (talk) 22:51, 23 February 2016 (CET)

15. August 2015 - Ausfall. Vermutlich Speicher voll oder UMTS Kontingent zu erneuern!

Kalibration 28./29. August

Kalibration erstellt am 1.9.2011

Kalibration Lebensministerium neue Station.

Log:

 In [4]: x=lies_Verzeichnis('/home/lightmeter1/Desktop/LIGHT/atlebmin2011-08-30/selected/', Datenschritt=3, Format="Lightrouter");
 In [5]: x.Observer.elev=209; x.Observer.lat="48:12:33.35"; x.Observer.long="16:22:54.32"; # LEBmin
 In [6]: x.plot_NatLight(); x.plot_Mond_GR_Lux_h_extinction()
 In [63]: ssh1=(x.Sunheight_deg_from_JD(x.JD)>30)|((x.Sunheight_deg_from_JD(x.JD)<0)&(x.Sunheight_deg_from_JD(x.JD)>-18))
 In [64]: sjd6=(x.JD>mx.DateTime.DateTime(2011,8,28,14,7,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,28,16,43,0).jdn)
 In [65]: sjd5=(x.JD>mx.DateTime.DateTime(2011,8,28,18,11,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,28,18,30,0).jdn)
 In [66]: sjd4=(x.JD>mx.DateTime.DateTime(2011,8,28,7,57,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,28,9,21,0).jdn)
 In [68]: sjd2=(x.JD>mx.DateTime.DateTime(2011,8,29,5,19,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,29,5,59,0).jdn)
 In [69]: sjd1=(x.JD>mx.DateTime.DateTime(2011,8,29,7,13,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,29,7,46,0).jdn)
 In [70]: sjd7=(x.JD>mx.DateTime.DateTime(2011,8,29,8,07,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,29,8,30,0).jdn)
 In [71]: sjd8=(x.JD>mx.DateTime.DateTime(2011,8,29,10,39,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,29,11,10,0).jdn)
 In [86]: sjd3=(x.JD>mx.DateTime.DateTime(2011,8,29,2,40,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,8,29,3,20,0).jdn)
 In [87]: SELECT=(sjd1|sjd2|sjd3|sjd4|sjd5|sjd6|sjd7|sjd8)&ssh1
 In [88]: x.fit_em1_to_natLight(JD_select=SELECT, Lux_range=[0.001,300000]);
 Fit: a,b,c,x0,d,: 2.501e+05,3.681e-03,3.637e-06,-1.094e-02,4.075e-03, (res2/(N-n-1))^1/2 = 0.078 ,N=5138,n= 5
 In [89]: a,b,c,x0,d= 2.501e+05,3.681e-03,3.637e-06,-1.094e-02,4.075e-03
 In [90]: execfile("draw_all_plots.py")
 In [91]: x.plot_Sunlight_Watt_per_square_meter_from_JD_DM(surface_pressure=991.0, T_dew_C=15.0, factor=110.9, c='c', marker=',' , linestyle=':')
 In [92]: legend(prop=matplotlib.font_manager.FontProperties(size='small'),loc='center left')
 In [93]: x.set_simple_time_axis_format()

a , b , c , x0 , d = 2.501e+05 , 3.681e-03 , 3.637e-06 , -1.094e-02 , 4.075e-03

--Nero2401 22:01, 1. Sep. 2011 (UTC)

22.August Abenddämmerung

In [128]: SJD1=(X.JD > mx.DateTime.DateTime(2011,8,22,14,0,0).jdn) & (X.JD < mx.DateTime.DateTime  (2011,8,24,19,10,0).jdn)
In [129]: SSH=(X.Sunheight_deg_from_JD(X.JD)>30)|((X.Sunheight_deg_from_JD(X.JD)<-4)&(X.Sunheight_deg_from_JD(X.JD)>-18))
In [130]: SJD=(SJD1)&SSH
In [131]: X.fit_em1_to_natLight(JD_select=SJD, Lux_range=[0.01,300000],type='em1c0T')
Fit: a,b,c,x0,d,: 2.630e+05,2.307e-03,5.728e-06,-5.280e-02,5.248e-03, (res2/(N-n-1))^1/2 = 0.0529 ,N=5900,n= 5
Out[131]: (array([  2.63024701e+05,   2.30669604e-03,   5.72755715e-06, -5.27981786e-02,   5.24831315e-03]), 1)
In [132]: clf()
In [133]: a,b,c,x0,d=2.630e+05,2.307e-03,5.728e-06,-5.280e-02,5.248e-03
In [134]: X.plot_em1c0T(a,b,c,x0,d,marker=',',color='g')
In [135]: X.plot_em1T(a,b,c,d,marker=',',color='k')
In [136]: X.plot_em1c0T(a,b,c,x0,d,marker=',',color='r', JD_select=SJD)
In [137]: X.plot_NatLight(); X.plot_Mond_GR_Lux_h_extinction()
In [138]: X.plot_em1T(a,b,c,d,marker=',',color='k')

AT_Lebensministerium_1, Mark 2.3

Sensor ID: 908065

... bis 2011-08-07

Bei Auswertung zu Beachten

In der Nacht von 26. -> 27. März '11 erfolgt die Sommerzeitumstellung. Davor: timezone=1, danach timezone=2.


Kalibration Ende April '11

Zur Anwendbarkeit dieser Kalibration im März. Kalibration by GuentherWuchterl

Kalibration April

Log:

a,b,c,x0,d = 1.206e+05,3.493e-02,1.393e-07,-1.164e-02,5.056e-03# Lebensmini cal April 2011, 
a,b,c,x0,d = 1.206e+05,3.493e-02,1.2898e-09,-1.164e-02,5.056e-03# Lebensmini cal April 2011, f=108
 ==> beste und mondschwach , Mittagsvariation minimiert
 SJD1=(X.JD > mx.DateTime.DateTime(2011,4,22,0,0).jdn) & (X.JD < mx.DateTime.DateTime(2011,4,24,0,00).jdn)
 SSH=(X.Sunheight_deg_from_JD(X.JD)>35)|((X.Sunheight_deg_from_JD(X.JD)<0)&(X.Sunheight_deg_from_JD(X.JD)>-18))
 SJD=(SJD1)&SSH
 X.fit_em1_to_natLight(JD_select=SJD, Lux_range=[0.01,300000],type='em1c0T',time_correction_seconds=0)
 a,b,c,x0,d,: 1.206e+05,3.493e-02,1.393e-07,-1.164e-02,5.056e-03, (res2/(N-n-1))^1/2 = 0.107 ,N=67186,n= 5
 X.plot_Sunlight_Watt_per_square_meter_from_JD_DM(surface_pressure=1013.0,T_dew_C=4.0,factor=110,c='c',marker=',',linestyle=':') 

Wetterquelle: Weather underground 1030 Wien http://www.wunderground.com/weatherstation/WXDailyHistory.asp?ID=IWIENWIE4&month=4&day=23&year=2011

 a,b,c,x0,d = 1.206e+05,3.493e-02,1.393e-07,-1.164e-02,5.056e-03# Lebensmini cal April 2011, 
 a,b,c,x0,d = 1.206e+05,3.493e-02,1.2898e-09,-1.164e-02,5.056e-03# Lebensmini cal April 2011, f=108
 Halbschaltung (bei den perfekten Kal-Beindungen): 1.-0.024/0.026      Out[53]: 0.0769

Log:

import os
os.chdir('/home/ke/Desktop/GW/Licht/python')
from Luxmeterdaten import *
import mx.DateTime
X=Luxmeterdaten('/home/ke/Desktop/GW/Licht/WienKartierung/Monitor    Netz/L0xx_Lebensministerium/AT_VIENNA_1_2011_03_17_15_06.csv' ,timezone=1, Datenschritt=1)
X.Observer.long=" 16:22:54.12"; X.Observer.lat=" 48:12:33.72"; X.Observer.elev=209 # Lebensministerium
X.Berechne_Sonne(); 
X.Berechne_Mond();
a,b,c,x0,d = 1.206e+05,3.493e-02,1.393e-07,-1.164e-02,5.056e-03# Lebensmini cal April 2011, (res2/(N-n-1))^1/2 = 0.107 ,N=67186,n= 5
a,b,c,x0,d = 1.206e+05,3.493e-02,1.2898e-09,-1.164e-02,5.056e-03# Lebensmini W/m2 cal April 2011, f=108
Schön und wenig Monschein
SJD1=(X.JD > mx.DateTime.DateTime(2011,4,20,3,30).jdn) & (X.JD < mx.DateTime.DateTime(2011,4,24,0,00).jdn)
niedrigste Mittagsvariation
SJD1=(X.JD > mx.DateTime.DateTime(2011,4,22,0,0).jdn) & (X.JD < mx.DateTime.DateTime(2011,4,24,0,00).jdn)
SSH=(X.Sunheight_deg_from_JD(X.JD)>35)|((X.Sunheight_deg_from_JD(X.JD)<0)&(X.Sunheight_deg_from_JD(X.JD)>-18))
SJD=(SJD1)&SSH
FIT=X.fit_em1_to_natLight(JD_select=SJD, Lux_range=[0.01,300000],type='em1c0T',time_correction_seconds=0)
figure(3); clf()
X.plot_NatLight(); X.plot_Mond_GR_Lux_h_extinction()
X.plot_em1c0T(a,b,c,x0,d,marker=',',color='g')
X.plot_em1T(a,b,c,d,marker=',',color='k')
X.plot_em1T(a,b,c,d,marker=',',color='r', JD_select=SJD)
X.plot_Sunlight_Watt_per_square_meter_from_JD_DM(surface_pressure=1013.0,T_dew_C=4.0,factor=110,c='c',marker=',',linestyle=':') 
Weather underground 1030 Wien http://www.wunderground.com/weatherstation/WXDailyHistory.asp? ID=IWIENWIE4&month=4&day=23&year=2011
title("Lightmeter Wien-Zentrum, Lebensministerium, cal. April 2011"); 
suptitle(X.Datum[X.JD.argmin()]+" to "+X.Datum[X.JD.argmax()])
xlabel("Time [UTC]"); ylabel("Illumination / [Lux]");
legend(prop=matplotlib.font_manager.FontProperties(size='small'),loc='upper right')
lebenstit="Wien Zentrum, Lebensministerium "+X.Datum[X.JD.argmin()]+" to "+X.Datum[X.JD.argmax()]
figure(); 
X.plot_dndLux_all_skies(X.em1T(X.Strom,X.Temperatur,a,b,c,d),bins=50,lgLux_range=[-5.5,-2],unit="Watt/(m*m)")
figure(); 
X.plot_dndLux_all_skies(X.em1T(X.Strom,X.Temperatur,a,b,c,d),bins=200,title=lebenstit)
--GuentherWuchterl 15:16, 25. Jun. 2011 (UTC)

Kalibration Ende März '11

Die Nächte sind leider nicht sehr sauber.

Log:

 x=lies_Verzeichnis('/home/lightmeter1/Desktop/LIGHT/atlebmin2011-09-04/', Datenschritt=30, timezone=1);
 a,b,c,x0,d= 1.296e+05,3.508e-02,2.289e-07,-2.910e-03,6.276e-03
 x.Observer.elev=209; x.Observer.lat="48:12:33.35"; x.Observer.long="16:22:54.32"; # LEBmin
 execfile("_hello.py")
 Calculating Lunar Ephemeris; please wait
 Calculating Solar Ephemeris; please wait
 ssh1=(x.Sunheight_deg_from_JD(x.JD)>30)|((x.Sunheight_deg_from_JD(x.JD)<0)&(x.Sunheight_deg_from_JD(x.JD)>-18))
 sjd1=(x.JD>mx.DateTime.DateTime(2011,4,3,0,0,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,4,3,4,48,0).jdn)
 sjd2=(x.JD>mx.DateTime.DateTime(2011,4,3,17,40,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,4,3,18,10,0).jdn)
 sjd3=(x.JD>mx.DateTime.DateTime(2011,4,2,15,25,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,4,2,18,25,0).jdn)
 sjd4=(x.JD>mx.DateTime.DateTime(2011,4,4,3,5,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,4,4,4,10,0).jdn)
 sjd5=(x.JD>mx.DateTime.DateTime(2011,4,3,9,0,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,4,3,14,48,0).jdn)
 sjd7=(x.JD>mx.DateTime.DateTime(2011,3,30,5,35,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,3,30,14,0,0).jdn)
 sjd6=(x.JD>mx.DateTime.DateTime(2011,4,10,5,15,0).jdn) & (x.JD<mx.DateTime.DateTime(2011,4,10,6,15,0).jdn)
 SELECT=(sjd1|sjd2|sjd3|sjd4|sjd5|sjd6|sjd7)&ssh1
 x.fit_em1_to_natLight(JD_select=SELECT, Lux_range=[0.00001,300000]);
 Fit: a,b,c,x0,d,: 1.214e+05,3.952e-02,1.248e-07,-1.833e-02,5.393e-03, (res2/(N-n-1))^1/2 = 0.173 ,N=1802,n= 5
 a,b,c,x0,d= 1.214e+05,3.952e-02,1.248e-07,-1.833e-02,5.393e-03
 x.plot_Sunlight_Watt_per_square_meter_from_JD_DM(surface_pressure=997.0, T_dew_C=3.0, factor=106, c='c', marker=',' , linestyle=':')
 legend(prop=matplotlib.font_manager.FontProperties(size='small'),loc='center left')
 x.set_simple_time_axis_format()

a , b , c , x0 , d = 1.214e+05 , 3.952e-02 , 1.248e-07 , -1.833e-02 , 5.393e-03 #hier c_Lux! f=106

--Nero2401 12:47, 4. Sep. 2011 (UTC)