Innkommende stråling.

[Amatør1]

Det meste er mulig, men hva mener du konkret med "månedsbasert datasett"? Er det et slags middel for måned eller noe?

[la1goa]

Jepp. Gjennomsnittsverdier for periodene man velger. Jeg skjønner at jeg startet på et håpløst prosjekt. Minuttverdiene jeg tok ut i første omgang var egentlig for kontroll, men jeg ser nyttigheten av dem også. Det som jeg også var ute etter var f.eks. forskjellen på 60-80N og 60-80S. Derfor startet jeg med dette.

Jeg kommer til å sette opp ei side der vi skal ta for oss breddegrad for breddegrad. Noen utvalgte lengdegrader avhengig av lengden på breddegraden slik at vi får et "gridsystem" på 104 punkter. Der hadde jeg planer om å legge ut temp, olr, innkommende stråling, skydekke, fuktighet m.m. pr punkt. igjennom tidene. Disse dataene er som regel gjennomsnitt pr. måned. Ja, jeg vet at dette er et huge prosjekt, men med litt hjelp er det ikke umulig.

Forresten var det ganske artige grafer du fikk til der. Mye morro man kan få ut av en liten snutt, som i aller høyeste grad kan bli et effektivt hjelpemiddel i klimaforskning.

[Amatør1]

Så hvis du velger en periode på 3 måneder, vil du ha gjennomsnitt for hele perioden eller for 3x1 måned? Er det kalendermåneder, hva om perioden starter f.eks. den 10. i en måned?

Med utgangspunkt i Julian Date er det litt klønete å lage snitt for kalendermåneder.

Det var litt stunt-programmering dette, men med AA+ er det mulig å gjøre ting fort om man har litt bakgrunn fra før. Tilfeldigvis er det å beregne solhøyden nesten samme problem som det jeg hadde da jeg programmerte et motorisert teleskop til å peke på vilkårlig valgte objekter på himmelen. Det var forresten et interessant problem, siden himmelen roterer mens du forflytter deg, så det finnes ingen deterministisk løsning, man må iterere seg fram.

Her er den essensielle snutten fra bruken av AA+ (dvs. Meeus-metodene)

Kode [Velg] Utvid

   
   // Ecliptical coordinates for the Sun at Julian Date JD
   CAAEllipticalPlanetaryDetails details = CAAElliptical::Calculate(JD,CAAElliptical::SUN);

   // Parallax compensation for observer's location
   CAA2DCoordinate delta;
   double Alpha     = details.ApparentGeocentricRA;
   double Delta     = details.ApparentGeocentricDeclination;
   double Distance  = details.ApparentGeocentricDistance;
   delta = CAAParallax::Equatorial2TopocentricDelta(Alpha,Delta,Distance,Longitude,Latitude,Height,JD);
   double ra_deg = (details.ApparentGeocentricRA + delta.X )*360.0/24.0;
   double de_deg =  details.ApparentGeocentricDeclination + delta.Y;

   // Apparent Greenwich Sidereal Time
   double AST = CAASidereal::ApparentGreenwichSiderealTime(JD);

   // observer longitude as hour angle
   double LongtitudeAsHourAngle = CAACoordinateTransformation::DegreesToHours(Longitude);

   // Local hour angle of object
   double Alpha2 = CAACoordinateTransformation::DegreesToHours(ra_deg);
   double LocalHourAngle = AST - LongtitudeAsHourAngle - Alpha2;
   
   // Compute Alt/Az coordinates
   CAA2DCoordinate Horizontal = CAACoordinateTransformation::Equatorial2Horizontal(LocalHourAngle, de_deg, Latitude);
   double az_deg  = Horizontal.X;
   double alt_deg = Horizontal.Y;

   if(alt_deg < 0.0) alt_deg = 0.0;

   // Solar altitude in radians
   double alt_rad = CAACoordinateTransformation::DegreesToRadians(alt_deg);


Skal prøve å få lagt ut noe imorgen, men det kan jo være en ide å legge inn middel som du foreslo.

[la1goa]

Jeg tenkte at en dataserie fra f.eks. 1970 til d.d. var månedsbasert, da det er slike serier jeg som regel laster ned. Så det hadde vært en flott funksjon.
Så gjennomsnitt måned for måned er det jeg er ute etter.

Jeg vet at julien day er problematisk å kjøre den funksjonen på, men ved å bruke vanlig dato konvertert til julian day er det lett: f.eks.

Fra  1970, 1, 1, 00:00   Til   1970,2,1, 00:00   minus et minutt
Fra  1970, 2, 1, 00:00   Til   1970,3,1, 00:00   minus et minutt
Fra  1970, 3, 1, 00:00   Til   1970,4,1, 00:00   minus et minutt
..
..

Legg merke til at det er kun månden som endrer seg. En trenger  ikke en gang ta hensyn til skuddårs problematikken.

Jeg husker fra teleskoptiden. Det var artige saker det. Computer styring og greier. Det var da jeg bestemte meg for å kjøpe boka. Den var litt mer avansert enn jeg forventet så den lå bare i skuffen. Men så plukket jeg den frem i sommer og begynte med å kalkulere en del, men har aldri kommet så langt som deg.

[Amatør1]

Da har jeg laget en versjon med støtte for månedsbasert graf + lagring til fil

Du kan laste ned Windows-versjonen sitt setup program her: WattSun_setup_V1.0-01.exe


Programmet ser nå slik ut
[attachimg=1 width=600]

Følgende ting gjør du
1. Velg From og To dato for tidsperioden du vil se på
2. Velg tidsskritt i minutter mellom hver utregning
3. Trykk på Run Graph for å få en graf med samtlige tidssskritt
4. Trykk på Run Monthly for å få en graf for tidsperioden med én verdi (middelverdien) for hele måneden. Denne blir plottet midt i måneden (ca. 15de i hver måned). For denne funksjonen blir dagen alltid satt til 1 for både From of To, uansett hva du velger.
5. Trykk på Save Watt.txt for å lagre siste graf i en tekstfil. Denne havner på installsjonsfolderen med mindre du justerer Start In-feltet i Windows Startmenyen (men hvis du lager desktop shortcut og starter programmet med den, havner watt.txt på desktopen).
6. Trykk på Clipboard for å kopiere siste graf til system clipboard. Du kan så bruke Paste (evt. "Lim Inn" på norsk) i et hvilket som helst program som støtter bilder.

Tekstfilen ser slik ut

Kode [Velg] Utvid

Date JD Watt at (long,lat)=(10.000,60.000)
2000-08-16T12:00:00 2451772.91668 361.761912679
2000-09-16T00:00:00 2451803.41668 262.768108304
2000-10-16T15:00:00 2451834.04168 142.391351007
2000-11-16T05:00:00 2451864.66668 42.7573264257
2000-12-16T17:00:00 2451895.16668 6.68734575475
2001-01-16T17:00:00 2451926.16668 21.3518872777
2001-02-15T05:00:00 2451955.66668 108.859647329
2001-03-16T14:00:00 2451985.04168 228.168038825
2001-04-16T00:00:00 2452015.41668 336.350744639
2001-05-16T12:00:00 2452045.91668 411.96734628
2001-06-16T00:00:00 2452076.41668 453.592781149
2001-07-16T12:00:00 2452106.91668 427.690690133
2001-08-16T12:00:00 2452137.91668 362.42981745
2001-09-16T00:00:00 2452168.41668 263.660339573
2001-10-16T15:00:00 2452199.04168 143.306328055
2001-11-16T05:00:00 2452229.66668 43.2880530595
2001-12-16T17:00:00 2452260.16668 6.78907933404
2002-01-16T17:00:00 2452291.16668 20.9252161899
2002-02-15T05:00:00 2452320.66668 107.943320983
2002-03-16T14:00:00 2452350.04168 227.200362039
2002-04-16T00:00:00 2452380.41668 335.610149588
2002-05-16T12:00:00 2452410.91668 411.490917701
2002-06-16T00:00:00 2452441.41668 453.524585743
2002-07-16T12:00:00 2452471.91668 428.098883445
2002-08-16T12:00:00 2452502.91668 363.074951752
2002-09-16T00:00:00 2452533.41668 264.555116137
2002-10-16T15:00:00 2452564.04168 144.234435366
2002-11-16T05:00:00 2452594.66668 43.83401141
2002-12-16T17:00:00 2452625.16668 6.89523164398
2003-01-16T17:00:00 2452656.16668 20.5016997097
2003-02-15T05:00:00 2452685.66668 107.051250276
2003-03-16T14:00:00 2452715.04168 226.266382123
2003-04-16T00:00:00 2452745.41668 334.89148533
2003-05-16T12:00:00 2452775.91668 411.001762427
2003-06-16T00:00:00 2452806.41668 453.409862858

Oppdatering: Endret nedlastingslinken øverst, rettet en feil. Versjonen er nå V1.0-01

[la1goa]

Da har jeg laget en versjon med støtte for månedsbasert graf + lagring til fil

Kjempeflott!!!!! 

Det er tydelig at astronomiske algoritmer slår inn da effekten er forskjellig fra år til år på samme tidspunkt. Men det som er meget interresant å se er at Sydpolen mottar ca 561 Watt på det meste mens Nordpolen mottar 526 Watt på det meste et halvår før/senere. dvs. ca 6.5% mer på maxnivåene. Men det må kjøres gjennomsnittlig for å være reelt.

Ellers en liten bugs som forvirret meg med det første. Lat/lon boksene har byttet plass, det tok litt tid før jeg oppdaget det.

Men nå skal det generes serier som settes opp mot OLR, temp, skydekke, m.m. for hvert punkt jeg har definert på kloden. Ca 104 stykker som skal gåes igjennom.

Takker så meget. Håper at det er flere som får nytte av denne.

/Knuta

[Amatør1]

Det er tydelig at astronomiske algoritmer slår inn da effekten er forskjellig fra år til år på samme tidspunkt. Men det som er meget interresant å se er at Sydpolen mottar ca 561 Watt på det meste mens Nordpolen mottar 526 Watt på det meste et halvår før/senere. dvs. ca 6.5% mer på maxnivåene. Men det må kjøres gjennomsnittlig for å være reelt.

Jepp. Jorda befiner seg i Perihel i Januar, dvs. den sydlige sommeren. Så dette er å forvente når man tar hensyn til den elliptiske banen.

Ellers en liten bugs som forvirret meg med det første. Lat/lon boksene har byttet plass, det tok litt tid før jeg oppdaget det.

Den eneste gjenstående feilen jeg ser er at (lat,long) i tittelen på grafen skulle vært (long,lat).

Men nå skal det generes serier som settes opp mot OLR, temp, skydekke, m.m. for hvert punkt jeg har definert på kloden. Ca 104 stykker som skal gåes igjennom.

Takker så meget. Håper at det er flere som får nytte av denne.

/Knuta

Vel bekomme  Ta noen kvalitetssjekker på kjente steder, så vi får luket ut eventuelle feil.

Det skal bli artig å se hva du kan bruke dette til.

[Jostemikk]

Aldeles fabelaktig! Dette forumet besitter sannelig en brukbar samlet kunnskap, for å si det mildt. Hva blir det neste? Albedometer etter snøstatistikker og is på innsjøene?

[Amatør1]

Aldeles fabelaktig! Dette forumet besitter sannelig en brukbar samlet kunnskap, for å si det mildt. Hva blir det neste? Albedometer etter snøstatistikker og is på innsjøene?

Det er kun fantasien som setter grenser, pleier jeg å si når jeg ikke kommer på noe annet     

Linux-versjonen fungerer jamen også:
[attachimg=1 width=600]

[la1goa]

Nå blir det i allefall fart i vellinga

Et eksempler på lon: 0   Lat: 80  der jeg har dataserier for inn og ut.



rød er ståling fra sola fra 0-500watt og blå er utgående 180-225watt. En ser med en gang at utgående stråling øker samtidig sol solen er mer eller mindre konstant. Men i følge gjeldende teorier skulle ikke OLR kunne øke. Men det blir artig etterhvert å kunne plukke punkter på jorden, sette dem sammen med diverse serier. Og legge puslespillet bit for bit, og dermed avlive hysteriet en gang for alle. Noen forslag på hva som kan gjøres?

[Amatør1]

Tidsserier med nettostråling kunne jo gjort seg. Hvis man hadde tilstrekkelig mange, kunne man lage animerte kart (med interpolasjon imellom punktene), der grader av blått representerer ulike verdier på netto utstråling. Ditto grader av rødt for innstråling.

Man kan jo drømme.