forum meteorologie, starea vremii, climatologie

Partners

Portal meteorologie   Galerie meteo   Meteo chat   Model WRF-RoMetEx   RSS  
Acum este Vin Mar 29, 2024 7:49 am

Ora este UTC + 2 [ DST ]




Scrie un subiect nou Răspunde la subiect  [ 19 mesaje ]  Du-te la pagina 1, 2  Următorul
Autor Mesaj
 Subiectul mesajului: Modelul WRF
MesajScris: Dum Oct 02, 2011 8:06 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
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 :lol: 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.htm

Pana 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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 04, 2011 8:20 pm 
Neconectat
Ploaie
Ploaie

Membru din: Joi Sep 15, 2011 4:55 am
Mesaje: 580
Localitate: Darwin-Australia
La ce rezolutie WRF ai rulat simularea pentru Katrina?


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 04, 2011 10:50 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
La rezolutie slaba, folosind nucleul WRF-ARW - Advanced Weather Research. In momentul actual, lucrez la implementarea nucleului NMM - Nonhydrostatic Mesoscale Model, care are o rezolutie de procesare a datelor mult mai buna.

Modelul WRF nu face singur prognozele, adica sa-i dau ca input datele de azi si sa faca prognoza pe maine si poimaine, ci trebuie sa-i dau practic prognoza GFS in format GRIB 1 sau 2 cu datele de azi, maine si poimaine. Modelul doar reexamineaza si recalculeaza datele GFS. Cu toate astea, parametrii de output sunt foarte numerosi, mult mai multi decat are GFS-ul sau ECMWF-ul.
Avand in vedere ca maine plec in Timisoara, daca o sa am acces la un internet destul de bun, o sa lucrez si acolo, dar daca nu, abia vineri voi putea sa ma reapuc de lucru.


Later Edit : Inca nu am ajuns la partea de setare a rezolutiei, setarile fiind foarte complicate.

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mie Oct 05, 2011 12:11 am 
Neconectat
Ploaie
Ploaie

Membru din: Joi Sep 15, 2011 4:55 am
Mesaje: 580
Localitate: Darwin-Australia
Ca sa merite efortul trebuie sa rulezi modelul WRF la o rezolutie foarte ridicata, cel putin 2 km, si daca nu ai un super computer, chestia iti ia si cateva zile, cand prognoza nu-ti mai foloseste la nimic. Iar la rezolutie slaba, modelele regionale dau rezultate f f bune...


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mie Oct 05, 2011 9:51 am 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Deoarece nu este un model propriu-zis, ci doar recalculeaza datele GFS, o rulare de 12 ore imi ia aproximativ 15 minute pe un procesor AMD Athlon II-x2 dual-core cu 3 GB de rami.

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 11, 2011 10:55 am 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Am reusit sa setez modelul WRF sa ruleze pentru Romania, la o rezolutie de 10 km. in momentul actual, fac o rulare de azi de la ora 00 pana in data de 13 oct. la ora 00. O sa postez rezultatele imediat ce se termina rularea :D 8) .

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 11, 2011 11:40 am 
Neconectat
Ploaie
Ploaie

Membru din: Joi Sep 15, 2011 4:55 am
Mesaje: 580
Localitate: Darwin-Australia
Daca merge bine.......prevad o poarta deschisa pentru 'Private sector".


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 11, 2011 1:40 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Rularea a fost un succes, timp total : o ora si jumatate, rezolutia : 10 km, suprafata : jumatatea de SV a Romaniei si o buna parte din Serbia.

Am produs harti de temperatura, precipitatii orare si precipitatii cumulate. Hartile sunt la un interval de o ora, valabile de azi, ora 00z pana in 13 oct, ora 00z si sunt puse in fisiere .pdf :

Temperatura la 2m : http://rometex.devip.net/WRF/teste/init11102011-00z/temp_2m.pdf
Precipitatii totale cumulate : http://rometex.devip.net/WRF/teste/init11102011-00z/precip_cumulative.pdf
Precipitatii totale orare : http://rometex.devip.net/WRF/teste/init11102011-00z/precip_orare.pdf

:D :lol:

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 11, 2011 6:49 pm 
Neconectat
Weather Guru
Weather Guru
Avatar utilizator

Membru din: Vin Dec 05, 2008 1:00 pm
Mesaje: 36951
Localitate: Bistrita
CrazyFrog merita felicitari pentru efort :)
Ramane de vazut cum se verifica.
Eu nu le-as mai pune in fisier .pdf, as face un site special.


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mar Oct 11, 2011 7:18 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Merci de aprecieri 8) . Asa mi le scoate programul NCAR Graphics :)) , in format pdf .O sa fac cumva sa automatizez rularile si sa extraga automat imaginile din fisierul pdf. Pana una alta, am pornit o rulare pentru intervalul 17.10.2011, 18z - 18.10.2011, 06z :lol: . O sa incerc sa fac harta cu freezing level si cu stratul de zapada la o rezolutie de 10 km.

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Dum Oct 16, 2011 12:13 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Astazi, daca reusesc, am sa implementez indicele ESI (European Snow Index) http://forum.meteorologie.ro/viewtopic.php?f=9&t=496 in programul de procesare grafica a datelor de la modelul WRF si sa afisez pe harta probabilitatea de a cadea precipitatii solide.

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mie Oct 19, 2011 2:48 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Am reusit sa implementez indicele ESI si am produs ca exemplu o harta a indicelui esi pentru data de 16.10.2011, ora 06z. Pot implementa o groaza de indici de instabilitate si nu numai :D .


Fişiere ataşate:
esi-11.jpg
esi-11.jpg [ 37.98 KiB | Vizualizat de 2392 ori ]

_________________
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
Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Mie Oct 19, 2011 7:19 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Am gasit un mic "bug" in scriptul ce l-am scris pentru plotarea indicelui ESI, dar am rezolvat problema :D
Cum am facut ? Raspunsul :

Am calculat indicele ESI cu formula de calcul, apoi mi-am dat seama ca nu e destul de corect afisarea lui in alte zone decat cele cu precipitatii, asa ca am scris un script ce suprapune zonele cu precipitatii peste 0.1 mm peste indicele ESI
Daca pp < 0.1 mm, pp ia valoarea 0, iar daca pp > 0.1, pp ia valoarea 1. Apoi inmultesc indicele ESI cu valoarea variabilei pp. Daca nu sunt precipitatii, deci pp = 0, oricat ar fi indicele ESI, va lua valoarea 0. Daca sunt precipitatii, indicele ESI, inmultit cu 1 va lua valoarea indicelui ESI original. :lol: .... multa bataie de cap pentru aflarea sansei de ninsoare :)). Eu initial am uitat sa activez algoritmii de conditionare si am inmultit indicele ESI direct cu valoarea PP..Cam asa arata corect harta cu indicele ESI :

L.E. : Trebuie sa ajustez culorile :))


Fişiere ataşate:
esi-11.jpg
esi-11.jpg [ 57.61 KiB | Vizualizat de 2379 ori ]

_________________
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
Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Dum Noi 06, 2011 10:38 am 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Am reusit instalarea nucleului WRF-NMM, cel mai comun utilizat si in scurt timp am sa fac un program care sa ruleze modelul periodic ( de 4 ori pe zi, la actualizarea GFS-ului ). Cel mai probabil, timpul de prognoza va fi de 3 sau 4 zile, la o rezolutie intre 5 si 10 km ( inca nu m-am decis ). Parametrii prognozati vor fi despre temperatura, nebulozitate, precipitatii, indici de instabilitate severa, presiune si vant. Suprafata ( domeniul ) modelului va fi Romania si zonele inconjuratoare. Pe viitor, poate am sa produc si meteograme, pentru majoritatea oraselor.

_________________
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


Sus
 Profil  
 
 Subiectul mesajului: Re: Modelul WRF
MesajScris: Joi Apr 26, 2012 7:01 pm 
Neconectat
Tornada
Tornada
Avatar utilizator

Membru din: Mie Feb 20, 2008 1:00 pm
Mesaje: 2761
Localitate: Resita / Timisoara
Din ce ati observat, prognozele oferite de modelul WRF pe care-l rulez sunt cat de cat exacte ? Momentan folosesc modulul de microfizica si setarile implicite, dar vreau sa-l personalizez pentru tara noastra.

_________________
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


Sus
 Profil  
 
Afişează mesajele de la anteriorul:  Sortează după  
Scrie un subiect nou Răspunde la subiect  [ 19 mesaje ]  Du-te la pagina 1, 2  Următorul

Ora este UTC + 2 [ DST ]


Cine este conectat

Utilizatorii ce navighează pe acest forum: Niciun utilizator înregistrat şi 0 vizitatori


Nu puteţi scrie subiecte noi în acest forum
Nu puteţi răspunde subiectelor din acest forum
Nu puteţi modifica mesajele dumneavoastră în acest forum
Nu puteţi şterge mesajele dumneavoastră în acest forum
Nu puteţi publica fişiere ataşate în acest forum

Căutare după:
Mergi la:  
Termeni şi condiţii
Powered by phpBB® Forum Software © phpBB Group
Translation/Traducere: phpBB România
phpBB SEO