FII.5. Megjelenítés az R program segítségével

A CHAPEAU kimeneti adatok ún. FA formátumban állnak elő ezek megjelenítése az „R” program, illetve az „Rfa” csomag segítségével lehetséges. Az R program ingyenesen elérhető az interneten és könnyen installálható Linux platformokra. Az Rfa csomag nem publikus, a csomagot az Országos Meteorológiai Szolgálaton (OMSZ) keresztül érhetjük el (szepszo.g@met.hu, boloni.g@met.hu). Az Rfa csomag telepítéséhez kövessük a csomagban elhelyezett README fájlban megadott leírást. Sikeres telepítés után kövessük az alábbi lépéseket:

  1. Indítsuk el az R programot:

    you@yourpc:~> R

  2. Az R programban kiadott parancsokat a továbbiakban „>” jellel vezetjük be.  Töltsük be az Rfa csomagot:

    > library (Rfa)

    Loading required package: geogrid

    Loading required package: mapNew

    Loading required package: maps

    Warning message:

    'Sys.putenv' is deprecated.

    Use 'Sys.setenv' instead.

    See help("Deprecated")

  3. Nyissunk meg a CHAPEAU által előállított két FA fájlt (pl. a 6 előrejelzést):

    > fa_6=FAopen("/yourpath/ICMSHABOF+0006")

  4. Listázzuk ki az egyik FA fájl tartalmát:

    > fa_6$list

      [1] "SPECSURFGEOPOTEN" "SURFPRESSION " "S001WIND.U.PHYS "

      [4] "S002WIND.U.PHYS " "S003WIND.U.PHYS " "S004WIND.U.PHYS "

      [7] "S005WIND.U.PHYS " "S006WIND.U.PHYS " "S007WIND.U.PHYS "

     [10] "S008WIND.U.PHYS " "S009WIND.U.PHYS " "S010WIND.U.PHYS "

     [13] "S011WIND.U.PHYS " "S012WIND.U.PHYS " "S013WIND.U.PHYS "

     [16] "S014WIND.U.PHYS " "S015WIND.U.PHYS " "S016WIND.U.PHYS "

     [19] "S017WIND.U.PHYS " "S018WIND.U.PHYS " "S019WIND.U.PHYS "

     [22] "S020WIND.U.PHYS " "S021WIND.U.PHYS " "S022WIND.U.PHYS "

     [25] "S023WIND.U.PHYS " "S024WIND.U.PHYS " "S025WIND.U.PHYS "

  5. Kódoljuk ki az FA fájl egy mezőjét (pl. a szél u komponensét a 20. és 21. modellszinteken) az "U_lev20" és "U_lev21" változókba):

    > U_lev20=FAdec(fa_6,"S020WIND.U.PHYS")

    > U_lev21=FAdec(fa_6,"S021WIND.U.PHYS")

  6. Rajzoljuk ki az "U_lev20" változót:

    > fcview(U_lev20)

  7. Rajzoljuk ki az "U_lev20" és "U_lev21" változók különbségét:

    > fcview(U_lev20-U_lev21)

  8. Rajzoljuk ki a szelet mint vektormezőt a 20. modellszinten:

    > V_lev20=FAdec(fa_6,"S020WIND.V.PHYS")

    > vview(U_lev20,V_lev20)

  9. Rajzoljuk ki a szél u komponensét a 20. szinten mint skalár mezőt, és a két komponens vektormezejét egyaránt:

    > fcview(U_lev20)

    > vview(add=TRUE,U_lev20,V_lev20)

  10. Irányítsuk fájlba (pl. postscript) megjelenítésünk eredményét (alapesetben egy külön ablakban jelenik meg az ábra):

    > postscript(file="fig.ps")

    > fcview(U_lev20)

    > dev.off()

  11. Tekintsük meg az elmentett képfájlt:

    > system("ghostview fig.ps")

Megjegyezzük, hogy az R szoftver segítségével korlátozottak a lehetőségek vertikális metszetek és idősorok megjelenítésére, mivel az FA fájlokban külön tömbökben találhatók az egyes szintek és időpontok. Fentiek okán (és a WRF-fel való könnyebb közös feldolgozás érdekében, valamint azért is, mert az egyetemi meteorológus oktatásban a WMO GRIB formátumú fájlok feldolgozása szerepel tananyagként) az FA fájlok GRIB formátumra való konverziója javasolt. Ez például a gribeuse konvertáló program segítségével lehetséges, ami az Országos Meteorológiai Szolgálaton (OMSZ) keresztül érhető el (szepszo.g@met.hu, boloni.g@met.hu).