Instaland pas cu pas modelul WRF, am obtinut informatii foarte faine despre cum functioneaza si ce poate WRF-ul.
Pe scurt, modelul WRF poate rula in 3 moduri, cele mai folosite, 2 la numar fiind : modul de date reale si modul de date ideale.
Cel mai des se foloseste modul de date reale, datele de initializare obtinandu-se din multe surse gratuite, dar cel mai des folosite sunt datele GRIB1 si GRIB2, ale GFS-ului.
Cand ruleaza in modul real, wrf-ul este compus din 2 mari module, si anume : WRF si WPS. WRF este modulul ce se socupa de partea de simulare, iar WPS-ul este modulul ce se ocupa de obtinerea datelor de initializare, interpolarea lor, convertirea in formate pe care WRF poate sa le foloseasca. Prin anumite comenzi in consola csh sau bash a linux-ului, acest proces de preprocesare al datelor este realizat cu usurinta. Urmeaza partea interesanta, cand se defineste domeniul de simulare - suprafata, nivele verticale si alte informatii, apoi se trece la partea efectiva de simulare. Intai se va defini prin editarea unui fisier datele despre domeniu, intervalul de rulare si alte chestii pe care inca nu le-am deslusit pe deplin, apoi se va face legatura - prin comenzi - intre datele de initializare si modulul de simulare. dupa cateva operatii de conversie a fisierelor se ruleaza simularea. Pe laptopul meu, cu un procesor AMD athlon II x2 dual core si cu 3 GB de memorie RAM, simulez 12 ore in aprox. 15 minute. La final se va genera un fisier de output care de fapt reprezinta o matrice 3D cu date.
Partea cea mai complicata insa, spre surprinderea mea
este generarea hartilor. Exista un program, numit NCAR Graphics sau NCL, care instalat si initializat, va putea produce harti, procesand datele din fisierul de output. Dar asta nu se face prin click-uri si next-uri, ci prin rularea unor scripturi destul de complicate. Un exemplu de script pentru a extrage temperatura pentru toti timpii prognozati din fisierul de output este :
Cod:
; Example script to produce plots for a WRF real-data run,
; with the ARW coordinate dynamics option.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
a = addfile("./wrfout_d01_2000-01-24_12:00:00.nc","r")
; We generate plots, but what kind do we prefer?
type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
wks = gsn_open_wks(type,"plt_Surface1")
; Set some basic resources
res = True
res@MainTitle = "REAL-TIME WRF"
pltres = True
mpres = True
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; What times and how many time steps are in the data set?
times = wrf_user_list_times(a) ; get times in the file
ntimes = dimsizes(times) ; number of times in the file
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
do it = 0,ntimes-1,2 ; TIME LOOP
print("Working on time: " + times(it) )
res@TimeLabel = times(it) ; Set Valid time to use on plots
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need
tc2 = wrf_user_getvar(a,"T2",it) ; T2 in Kelvin
tc2 = tc2-273.16 ; T2 in C
; Plotting options for T
opts = res
opts@cnFillOn = True
opts@ContourParameters = (/ -20., 90., 5./)
opts@gsnSpreadColorEnd = -3 ; End third from the last color in color map
contour_tc = wrf_contour(a,wks,tf2,opts)
delete(opts)
; MAKE PLOTS
plot = wrf_map_overlays(a,wks,(/contour_tc/),pltres,mpres)
end do ; END OF TIME LOOP
end
Exista tone de exemple de scripturi pt producerea de harti, diagrame, etc pe site-ul userilor WRF :
http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/NCL_examples.htmPana acum, am rulat cu succes o simulare a uraganului Katrina, pentru 24 ore, acum ruland una pt 3 zile.
Indata ce am instalat modelul pe un calculator, am sa incep simularile cu date actuale.
O chestie interesanta este ca pot folosi ca date de initializare, datele de prognoza ale GFS-ului, deci voi putea continua prognoza GFS, in stil WRF. De asemenea, se pot folosi seturi de date de initializare.
Mai multe date, puteti gasi pe site-ul userilor WRF :
http://www.mmm.ucar.edu/wrf/users/O sa va tin la curent cu progresul facut.
_________________
Oană Liviu
PhD Student - West University of Timisoara, Faculty of Mathematics and Computer Science
Research Assistant - West University of Timisoara, Faculty of Mathematics and Computer Science
http://wrf.rometex.org, e-mail :
liviu.oana88@gmail.com