r.sungrass - Online w chmurze

PROGRAM:

IMIĘ

r.słońce - Natężenie promieniowania słonecznego i model napromieniowania.
Oblicza mapy rastrowe bezpośredniego (wiązek), rozproszonego i odbitego promieniowania słonecznego dla danego dnia,
szerokość geograficzna, powierzchnia i warunki atmosferyczne. Parametry słoneczne (np. wschód, zachód słońca
czasy, deklinacja, natężenie promieniowania pozaziemskiego, długość dnia) są zapisywane na mapie
plik historii. Alternatywnie można określić czas lokalny w celu obliczenia padania światła słonecznego
mapy rastrowe kąta i/lub natężenia napromienienia. Efekt cieniowania topografii jest opcjonalny
rejestrowy.

SŁOWA KLUCZOWE

raster, energia słoneczna, energia słońca, cień

STRESZCZENIE

r.słońce
r.słońce --help
r.słońce [-pm] podniesienie=ciąg [aspekt=ciąg] [wartość_aspektu=unosić się] [nachylenie=ciąg]
[wartość_nachylenia=unosić się] [lewo=ciąg] [wartość_powiązania=unosić się] [albedo=ciąg]
[wartość_albedo=unosić się] [lat=ciąg] [długie=ciąg] [współczynnik_bh=ciąg]
[współczynnik_dh=ciąg] [nazwa_bazowa_horyzontu=nazwa podstawowa] [krok_horyzontu=unosić się] [incydent=ciąg]
[promień_rad=ciąg] [różnica_rad=ciąg] [refl_rad=ciąg] [glob_rad=ciąg]
[czas_słoneczny=ciąg] dzień=liczba całkowita [krok=unosić się] [deklinacja=unosić się] [czas=unosić się]
[odległość_kroku=unosić się] [nprzegrody=liczba całkowita] [czas_cywilny=unosić się] [--przepisać]
[--pomoc] [--gadatliwy] [--cichy] [--ui]

Flagi:
-p
Nie uwzględniaj efektu cienia terenu

-m
Użyj wersji programu o małej ilości pamięci

--przepisać
Zezwalaj plikom wyjściowym na zastępowanie istniejących plików

--help
Wydrukuj podsumowanie wykorzystania

--gadatliwy
Pełne wyjście modułu

--cichy
Cichy moduł wyjściowy

--UI
Wymuś uruchomienie okna GUI

Parametry:
podniesienie=ciąg [wymagany]
Nazwa wejściowej mapy rastrowej wysokości [metry]

aspekt=ciąg
Nazwa wejściowej mapy aspektowej (układ terenu lub azymut panelu słonecznego) [dziesiętny
stopni]

wartość_aspektu=unosić się
Pojedyncza wartość orientacji (aspektu), 270 to południe
Zaniedbanie: 270

nachylenie=ciąg
Nazwa mapy rastrowej nachylenia wejściowego (nachylenie terenu lub nachylenie panelu słonecznego) [dziesiętne
stopni]

wartość_nachylenia=unosić się
Pojedyncza wartość nachylenia (nachylenie)
Zaniedbanie: 0.0

lewo=ciąg
Nazwa mapy rastrowej wejściowego współczynnika zmętnienia atmosfery Linkego [-]

wartość_powiązania=unosić się
Pojedyncza wartość współczynnika zmętnienia atmosfery Linkego [-]
Zaniedbanie: 3.0

albedo=ciąg
Nazwa mapy rastrowej wejściowego współczynnika albedo gruntu [-]

wartość_albedo=unosić się
Pojedyncza wartość współczynnika albedo gruntu [-]
Zaniedbanie: 0.2

lat=ciąg
Nazwa wejściowej mapy rastrowej zawierającej szerokości geograficzne [stopnie dziesiętne]

długie=ciąg
Nazwa wejściowej mapy rastrowej zawierającej długości geograficzne [stopnie dziesiętne]

współczynnik_bh=ciąg
Nazwa wejściowej mapy rastrowej współczynnika promieniowania wiązki nieba rzeczywistego (gruba chmura) [0-1]

współczynnik_dh=ciąg
Nazwa wejściowej mapy rastrowej współczynnika promieniowania rozproszonego (zamglenia) rzeczywistego nieba [0-1]

nazwa_bazowa_horyzontu=nazwa podstawowa
Nazwa bazowa mapy wejściowej informacji o horyzoncie

krok_horyzontu=unosić się
Wielkość kroku kąta dla horyzontu wielokierunkowego [stopnie]

incydent=ciąg
Wyjściowa mapa rastrowa kąta padania (tylko tryb 1)

promień_rad=ciąg
Natężenie promieniowania wiązki wyjściowej [Wm-2] (tryb 1) lub rastrowa mapa napromieniowania [Wh.m-2.day-1] (tryb
2)

różnica_rad=ciąg
Wyjście rozproszonego natężenia napromieniowania [Wm-2] (tryb 1) lub rastrowa mapa napromieniowania [Wh.m-2.day-1]
(tryb 2)

refl_rad=ciąg
Masa wyjściowa odbitego natężenia napromieniowania [Wm-2] (tryb 1) lub rastrowa mapa napromieniowania
[Wh.m-2.day-1] (tryb 2)

glob_rad=ciąg
Wyjściowe globalne (całkowite) natężenie napromieniowania/napromieniowanie [Wm-2] (tryb 1) lub
mapa rastrowa napromieniowania/napromieniowania [Wh.m-2.day-1] (tryb 2)

czas_słoneczny=ciąg
Wyjściowa mapa rastrowa czasu nasłonecznienia [h] (tylko tryb 2)

dzień=liczba całkowita [wymagany]
Nr dnia w roku (1-365)
Opcje: 1-365

krok=unosić się
Krok czasowy przy obliczaniu całodobowych sum promieniowania [godziny dziesiętne]
Zaniedbanie: 0.5

deklinacja=unosić się
Wartość deklinacji (nadpisująca wartość obliczoną wewnętrznie) [radany]

czas=unosić się
Czas lokalny (słoneczny) (do ustawienia tylko dla trybu 1) [godziny dziesiętne]
Opcje: 0-24

odległość_kroku=unosić się
Współczynnik kroku odległości próbkowania (0.5–1.5)
Zaniedbanie: 1.0

nprzegrody=liczba całkowita
Przeczytaj pliki wejściowe w tej liczbie fragmentów
Zaniedbanie: 1

czas_cywilny=unosić się
Wartość cywilnej strefy czasowej, jeśli nie, czas będzie lokalnym czasem słonecznym

OPIS

r.słońce oblicza mapy rastrowe promieniowania słonecznego (bezpośredniego), rozproszonego i odbitego od podłoża
dla danego dnia, szerokości geograficznej, powierzchni i warunków atmosferycznych. Parametry słoneczne (np. czas
wschodu i zachodu słońca, deklinacja, natężenie promieniowania pozaziemskiego, długość światła dziennego).
przechowywane w plikach historii wynikowych map. Alternatywnie może być czas lokalny
przeznaczony do obliczania map rastrowych kąta padania światła słonecznego i/lub natężenia napromienienia. Cienie
Domyślnie uwzględniany jest efekt topografii. Można to zrobić wewnętrznie za pomocą
Obliczenie efektu zacienienia bezpośrednio z cyfrowego modelu elewacji lub poprzez
określenie map rastrowych wysokości horyzontu, co jest znacznie szybsze. Te rastrowe horyzonty
mapy można obliczyć za pomocą r.horizon.

W przypadku współrzędnych szerokości i długości geograficznej wymagane jest, aby mapa wysokości była podana w metrach. The
zasady są następujące:

· współrzędne szerokości/długości: wysokość w metrach;

· Inne współrzędne: wysokość w tej samej jednostce, co współrzędne wschód-północ.
Geometria modelu słońca opiera się na pracach Krcho (1990), później poprawiona przez
Jenco (1992). Równania opisujące położenie Słońca - Ziemi oraz interakcję
Promieniowanie słoneczne z atmosferą pierwotnie opierało się na wzorach zaproponowanych przez
Kitlera i Miklera (1986). Komponent ten został znacznie zaktualizowany przez wyniki i
sugestie grupy roboczej koordynowanej przez Scharmera i Greifa (2000) (niniejszy algorytm
może zostać zastąpiony przez bibliotekę algorytmów SOLPOS zawartą w GRASS w poleceniu r.sunmask).
Model oblicza wszystkie trzy składowe globalnego promieniowania (wiązkę, rozproszenie i odbicie).
dla warunków czystego nieba, czyli bez uwzględnienia wymiaru przestrzennego i czasowego
odmiana chmur. Zasięg i rozdzielczość przestrzenna modelowanego obszaru, a także
integracji w czasie, ograniczone są jedynie zasobami pamięci i przechowywania danych. The
model zbudowany jest z myślą o potrzebach użytkowników z różnych dziedzin nauki (hydrologia, klimatologia,
ekologia i nauki o środowisku, fotowoltaika, inżynieria itp.) dla kontynentalnych,
regionalną aż po skalę krajobrazową.

Model uwzględnia efekt cienia lokalnej topografii, chyba że zostanie wyłączony
dotychczasowy -p flag. r.słońce działa w dwóch trybach: W pierwszym trybie dokonuje obliczeń dla zadanego lokalnego
czas i kąt padania światła słonecznego [stopnie] oraz wartości natężenia promieniowania słonecznego [Wm-2]. W sekundę
w trybie dobowe sumy promieniowania słonecznego [Wh.m-2.day-1] obliczane są w obrębie zadanego dnia. przez A
skryptów, te dwa tryby mogą być używane oddzielnie lub w kombinacji w celu zapewnienia szacunków
na dowolny żądany przedział czasu. Model uwzględnia zasłanianie nieba przez lokalną rzeźbę terenu
cechy. W powstałych plikach historii map zapisywanych jest kilka parametrów oświetlenia słonecznego, które
można wyświetlić za pomocą polecenia r.info.

Mapa rastrowa kąta padania słońca incydent jest obliczana z określeniem mapy rastrowej wysokości
podniesienie, aspekt mapy rastrowej aspekt, mapa rastrowa stromości zboczy nachylenie, dany dzień dzień
i czasu lokalnego czas. Nie ma potrzeby definiowania szerokości geograficznej dla lokalizacji o znanych i
zdefiniowany rzut/układ współrzędnych (sprawdź to komendą g.proj). Jeśli masz
niezdefiniowany rzut, układ (x, y) itp., dla którego można jawnie zdefiniować szerokość geograficzną
duże obszary według wejściowej mapy rastrowej łacina z interpolowanymi wartościami szerokości geograficznej. Wszystkie rastry wejściowe
mapy muszą być mapami rastrowymi zmiennoprzecinkowymi (FCELL). Dane zerowe na mapach są wyłączone z
obliczeń (a także przyspieszenie obliczeń), więc każda wyjściowa mapa rastrowa będzie zawierać
puste dane w komórkach zgodnie ze wszystkimi wejściowymi mapami rastrowymi. Użytkownik może użyć polecenia r.null
utwórz/zresetuj plik zerowy dla wejściowych map rastrowych.
Określony dzień dzień to numer dnia w roku ogólnym, w którym dniem jest 1 stycznia
nr 1 i 31 grudnia to 365. Czas czas musi być czasem lokalnym (słonecznym) (tj. NIE czasem strefowym,
np. GMT, CET) w systemie dziesiętnym, np. 7.5 (= 7h 30m), 16.1 = 4h 6m PM.

Słoneczna deklinacja parametr jest opcją zastąpienia wartości obliczonej przez
Wewnętrzna rutyna na dzień roku. Wartość szerokości geograficznej można ustawić jako
stała dla całego obliczanego obszaru lub opcjonalnie siatka reprezentująca przestrzennie
rozłożone wartości na dużym obszarze. Szerokość geograficzna musi być również podana w postaci dziesiętnej
układ z wartościami dodatnimi dla półkuli północnej i ujemnymi dla półkuli południowej. W
podobna zasada współczynnika zmętnienia Linkego (lewo, lin ) i albedo zmielone (albedo, biały)
może być ustawiony.

Oprócz promieniowania czystego nieba użytkownik może obliczyć promieniowanie nieba rzeczywistego (wiązka, rozproszone)
za pomocą współczynnik_bh i współczynnik_dh wejściowe mapy rastrowe określające ułamek odpowiedniego
promieniowanie czystego nieba zmniejszone przez czynniki atmosferyczne (np. zachmurzenie). Wartość jest
pomiędzy 0-1. Zwykle współczynniki te można uzyskać z długoterminowej prognozy meteorologicznej
pomiary dostarczone w postaci map rastrowych z rozkładem przestrzennym tych współczynników
oddzielnie dla promieniowania wiązkowego i rozproszonego (patrz Suri i Hofierka, 2004, rozdział 3.2).

Mapy rastrowe nasłonecznienia lub natężenia napromieniowania promień_rad, różnica_rad, refl_rad są obliczane
na dany dzień dzień, szerokość łacina, elewacja podniesienie, nachylenie nachylenie i aspekt aspekt
mapy rastrowe. Dla wygody raster wyjściowy podano jako glob_rad wyświetli sumę
trzech składników promieniowania. Program wykorzystuje współczynnik zmętnienia atmosfery Linkego i
współczynnik albedo przyziemnego. Domyślną, pojedynczą wartością współczynnika Linkego jest lin= 3.0 i jest blisko
średnioroczna dla obszarów wiejskich i miejskich. Czynnik Linkego jest absolutnie jasny
atmosfera jest lin=1.0. Zobacz uwagi poniżej, aby dowiedzieć się więcej na temat tego czynnika. Częstość występowania
kąt słoneczny to kąt między horyzontem a wektorem wiązki słonecznej.

Mapy nasłonecznienia na dany dzień obliczane są poprzez całkowanie odpowiednich danych
natężenie napromienienia pomiędzy wschodem a zachodem słońca w danym dniu. Użytkownik może ustawić dokładniejszy lub
zgrubniejszy krok czasowy używany do obliczeń promieniowania całodobowego za pomocą krok opcja.
domyślna wartość krok wynosi 0.5 godziny. Większe kroki (np. 1.0-2.0) mogą przyspieszyć obliczenia
ale dają mniej wiarygodne (i bardziej postrzępione) wyniki. Gdy słońce przechodzi przez ok. 15°
nieba w ciągu godziny, ustawienie domyślne krok pół godziny spowoduje powstanie kroków co 7.5°
dane. Dla stosunkowo płynnego wyjścia ze słońcem umieszczonym dla każdego stopnia ruchu w
niebo powinieneś ustawić krok do 4 minut lub mniej. krok=0.05 odpowiada każdemu 3
minuty. Oczywiście ustawienie bardzo małego kroku czasowego proporcjonalnie zwiększa
czas działania modułu.

Jednostkami wyjściowymi są Wh na metr kwadratowy na dany dzień [Wh/(m*m)/dzień]. Częstość występowania
mapy kąta i napromieniowania/napromieniowania są obliczane z uwzględnieniem wpływu cienia na relief
domyślnie. Możliwe jest również ich obliczenie bez tego wpływu za pomocą
flaga płaska (-p). Na obszarach górskich może to prowadzić do bardzo różnych wyników! Użytkownik
Należy mieć świadomość, że uwzględnienie efektu cieniowania reliefu może spowolnić
szybkość obliczeń, zwłaszcza gdy wysokość słońca jest niska.

Biorąc pod uwagę efekt cieniowania, można uwzględnić szybkość i precyzję obliczeń
kontrolowany przez odległość_kroku parametr określający gęstość próbkowania, przy której
widoczność komórki siatki jest obliczana w kierunku ścieżki przepływu światła słonecznego. To
definiuje również metodę obliczania wysokości przeszkody. Wybierając A
odległość_kroku mniejsza niż 1.0 (tj. punkty pobierania próbek zostaną obliczone w odległość_kroku *
odległość wielkości komórki), r.słońce pobiera wysokość z najbliższego punktu siatki. Wartości powyżej 1.0
użyje maksymalnej wartości wysokości znalezionej w najbliższych 4 otaczających punktach siatki. The
domyślna wartość odległość_kroku=1.0 powinno dać rozsądne wyniki w większości przypadków (np
DEM). The odległość_kroku wartość określa współczynnik mnożenia odległości próbkowania.
Ta podstawowa odległość próbkowania jest równa średniej arytmetycznej obu rozmiarów komórek. The
rozsądne wartości mieszczą się w przedziale 0.5-1.5. Wartości poniżej 0.5 będą się zmniejszać i wartości
powyżej 1.0 zwiększy prędkość obliczeń. Wartości większe niż 2.0 mogą dawać szacunki
z mniejszą dokładnością w przypadku silnie rozciętego reliefu. Zapisywane są całkowicie zacienione obszary
mapy wyjściowe mają wartości zerowe. Obszary z danymi NULL są uważane za nie stanowiące bariery
efekt cieniowania.

Generowane są pliki historii map zawierające następujące parametry używane w
obliczenia:
- Stała słoneczna 1367 Wm-2
- Natężenie promieniowania pozaziemskiego w płaszczyźnie prostopadłej do wiązki światła słonecznego [Wm-2]
- Dzień w roku
- Deklinacja [radany]
- Godzina dziesiętna (tylko alternatywa 1)
- Wschód i zachód słońca (min.-maks.) w płaszczyźnie poziomej
- Długości światła dziennego
- Szerokość geograficzna (min.-maks.)
- Współczynnik zmętnienia Linkego (min.-maks.)
- Albedo naziemne (min.-maks.)

Użytkownik może użyć ładnego skryptu powłoki ze zmiennym dniem do obliczenia promieniowania przez pewien czas
odstępie czasu w ciągu roku (np. okres wegetacyjny lub zimowy). Wysokość, aspekt i nachylenie
wartości wejściowych nie należy reklasyfikować do bardziej ogólnych kategorii. Może to prowadzić do
nieprawidłowe wyniki.

OPCJE

Obecnie istnieją dwa tryby r.sun. W pierwszym trybie oblicza padanie światła słonecznego
mapy rastrowe kąta i natężenia promieniowania słonecznego z wykorzystaniem ustawionego czasu lokalnego. W drugim trybie codziennie
sumy nasłonecznienia [Wh.m-2.day-1] obliczane są dla określonego dnia.

UWAGI

Energia słoneczna jest ważnym parametrem wejściowym w różnych modelach dotyczących energii
przemysł, krajobraz, roślinność, ewapotranspiracja, topnienie śniegu lub teledetekcja. Słoneczny
Mapy kąta padania promieni mogą być skutecznie wykorzystywane w radiometrii i topografii
poprawki w terenie górzystym i pagórkowatym, gdzie powinny być bardzo dokładne badania
wykonane.

Model promieniowania słonecznego przy czystym niebie zastosowany w r.sun opiera się na podjętych pracach
na rzecz opracowania Europejskiego Atlasu Promieniowania Słonecznego (Scharmer i Greif 2000, Page i in.
2001, Rigollier 2001). Model czystego nieba szacuje globalne promieniowanie na podstawie sumy
jego wiązki, elementy rozproszone i odbite. Główna różnica między promieniowaniem słonecznym
modeli dla powierzchni pochyłych w Europie jest obróbka elementu rozproszonego. w
Klimat europejski Składnik ten jest często największym źródłem błędu szacunków. Nabierający
pod uwagę istniejące modele i ich ograniczenia w europejskim promieniowaniu słonecznym
Zespół Atlas wybrał model Muneera (1990), ponieważ ma on solidne podstawy teoretyczne i tym samym
większy potencjał do późniejszej poprawy.

Szczegóły podstawowych równań używanych w tym programie można znaleźć w odnośniku
literatura cytowana poniżej lub książka opublikowana przez Neteler i Mitasova: Open Source GIS: A GRASS
Podejście GIS (opublikowane w Kluwer Academic Publishers w 2002 r.).

Średnie miesięczne wartości współczynnika zmętnienia Linkego dla klimatu łagodnego w woj
półkula północna (zobacz literaturę dotyczącą swojego obszaru badań):

Miesiąc styczeń luty mar kwi maj cze lipiec sie wrzesień paź lis grudzień roczny

góry 1.5 1.6 1.8 1.9 2.0 2.3 2.3 2.3 2.1 1.8 1.6 1.5 1.90

wieś 2.1 2.2 2.5 2.9 3.2 3.4 3.5 3.3 2.9 2.6 2.3 2.2 2.75

miasto 3.1 3.2 3.5 4.0 4.2 4.3 4.4 4.3 4.0 3.6 3.3 3.1 3.75

przemysłowy 4.1 4.3 4.7 5.3 5.5 5.7 5.8 5.7 5.3 4.9 4.5 4.2 5.00

Planowane ulepszenia obejmują wykorzystanie algorytmu SOLPOS do geometrii Słońca
obliczenia i wewnętrzne obliczenia aspektu i nachylenia.

Solar czas
Domyślnie r.sun oblicza czasy jako prawdziwy czas słoneczny, przy czym południe słoneczne jest zawsze dokładne
12:XNUMX wszędzie w bieżącym regionie. W zależności od tego, gdzie znajduje się strefa zainteresowania
zlokalizowane w odpowiedniej strefie czasowej, w niektórych przypadkach może to powodować różnice sięgające nawet godziny
przypadków (takich jak zachodnia Hiszpania) jest jeszcze więcej. Poza tym offset zmienia się w ciągu roku
zgodnie z równaniem czasu.

Aby przezwyciężyć ten problem, użytkownik może skorzystać z tej opcji czas_cywilny= in
r.sun, aby używać czasu rzeczywistego (zegar ścienny). Na przykład dla Europy Środkowej
przesunięcie strefy czasowej wynosi +1, +2, gdy obowiązuje czas letni.

Ekstrakcja of cień mapy
Mapę cieni można wyodrębnić z mapy kąta padania światła słonecznego (incidout). Obszary
z wartościami zerowymi są cieniowane. To nie zadziała, jeśli -p flaga była używana.

Duży mapy i na zewnątrz of pamięć problemy
W przypadku dużej liczby lub kolumn i wierszy, r.słońce może zużywać znaczną ilość pamięci.
Chociaż wyjściowe mapy rastrowe nie podlegają partycjonowaniu, wejściowe mapy rastrowe korzystają z
nprzegrody parametr. W przypadku błędu braku pamięci (BŁĄD: G_malloc: brak pamięci),
dotychczasowy nprzegrody parametru można użyć do uruchomienia obliczeń segmentowych, które zużywają mniej
pamięć podczas obliczeń. Ilość pamięci wg r.słońce szacuje się w następujący sposób:
# bez wejściowego podziału mapy rastrowej:
# wymagania dotyczące pamięci: 4 bajty na komórkę rastrową
# wiersze, kolumny: wiersze i kolumny bieżącego regionu (dowiedz się za pomocą g.region)
# IR: liczba wejściowych map rastrowych bez map horyzontalnych
# OR: liczba wyjściowych map rastrowych
bajty_pamięci = wiersze*kolumny*(IR*4 + kroki_horyzontu + LUB*4)
# z wejściowym podziałem mapy rastrowej:
bajty_pamięci = wiersze*kols*((IR*4+kroki_horyzontu)/npartycje + LUB*4)

PRZYKŁADY

Przykład z Karoliny Północnej (biorąc pod uwagę również rzucane cienie):
g.region raster=elewacja -p
# oblicz kąty horyzontu (aby przyspieszyć późniejsze obliczenia r.słońce)
r.horyzont wzniesienie=stopień wzniesienia=30 strefa buforowa=200 nazwa podstawowa=horkąt \
maksymalna odległość=5000
# nachylenie + aspekt
r.slope.aspect elewacja=elewacja aspekt=aspekt.dem nachylenie=nachylenie.dem
# oblicz globalne promieniowanie dla dnia 180 o godzinie 2:XNUMX, korzystając z wyjścia r.horizon
r.słońce wzniesienie=wzniesienie nazwa_bazy_horyzontu=horkąt stopień_horyzontu=30 \
aspekt=aspekt.dem nachylenie=nachylenie.dem glob_rad=global_rad dzień=180 czas=14
# wynik: globalne (całkowite) natężenie napromieniowania/napromieniowanie [Wm-2] dla danego dnia/godziny
r.univar global_rad

Obliczenie zintegrowanego dziennego napromieniowania dla regionu w Północnej Karolinie dla danego zagadnienia
dzień roku w rozdzielczości 30 m. Tutaj dzień 172 (tj. 21 czerwca w latach nieprzestępnych):
g.region raster=wysokość_ned_30m -str
# biorąc pod uwagę rzucane cienie
r. wysokość słońca=wysokość_ned_30m wartość_linke=2.5 wartość_albedo=0.2 dzień=172 \
promień_rad=b172 diff_rad=d172 \
refl_rad=r172 insol_time=it172
d.mon wx0
# pokaż mapę rastrową napromieniowania [Wh.m-2.day-1]
d.rast.leg b172
# pokaż mapę rastrową czasu nasłonecznienia [h]
d.rast.leg it172
Możemy obliczyć dzień roku na podstawie określonej daty w powłoce Pythona:
>>> importuj datę i godzinę
>>> datetime.datetime(2014, 6, 21).timetuple().tm_yday
172

Korzystaj z r.sungrass online, korzystając z usług onworks.net