ELTE Meteorológiai Tanszék => Oktatás => Magics++ oktatóanyag => Fortran oktatóanyag
Az első példában egy európai térképkivágatot rajzolunk.
program coast_1
c Magics megnyitása és a kimenő postscript fájl beállítása: A4-es lap
call popen
call psetc ('ps_device', 'ps_a4')
call psetc ('ps_file_name', 'coast_1a.ps')
c A térkép beállítása:
call psetr ('subpage_lower_left_latitude', 30.)
call psetr ('subpage_lower_left_longitude', -15.)
call psetr ('subpage_upper_right_latitude', 70.)
call psetr ('subpage_upper_right_longitude', 40.)
c A partvonalak és a rácshálózat beállításai:
call psetc ('map_coastline_colour', 'tan')
call psetc ('map_grid_colour', 'tan')
call psetc ('map_grid_line_style', 'dot' )
c Partvonal kirajzolása:
call pcoast
c Magics bezárása:
call pclose
end
A polársztereografikus térkép paraméterei:
'subpage_lower_left_latitude': a térkép bal alsó szélességi koordinátája'subpage_lower_left_longitude': a térkép bal alsó hosszúsági koordinátája'subpage_upper_right_latitude': a térkép jobb felső szélességi koordinátája'subpage_upper_right_longitude': a térkép jobb felső hosszúsági koordinátája'map_coastline_colour': a partvonalak színe'map_grid_colour': a rácshálózat színe'map_grid_line_style': a rácshálózat vonal típusa (SOLID/DASH/DOT/CHAIN_DASH/CHAIN_DOT)Alapértelmezésben hengervetületi térképet kapunk.
Ha ezen változtatni szeretnénk, a
call psetr ('SUBPAGE_MAP_PROJECTION',...) paranccsal tudunk
CYLINDRICAL/POLAR_STEREOGRAPHIC/MERCATOR/SATELLITE/LAMBERT/AITOFF/OCEAN_SECTION/NONE
vetületeket kirajzolni. A polár-sztereografikus vetülettel a 2. példában foglalkozunk.
Az eredmény:
További beállítási lehetőségek, paraméterek: az ECMWF honlapján.
Második példánkban a térkép elkészítéséhez polár-sztereografikus vetületet alkalmazunk.
program coast_2
call popen
call psetc ('ps_device', 'ps_a4')
call psetc ('ps_file_name', 'coast_2.ps')
c A polár-sztereografikus térkép:
call psetc ('subpage_map_projection','polar_stereographic')
call psetc ('subpage_map_area_definition','centre')
call psetr ('subpage_map_centre_longitude',-95.)
call psetr ('subpage_map_centre_latitude',40.)
call psetr ('subpage_map_vertical_longitude',-95.)
call psetr ('subpage_map_scale', 27.e6 )
c A partvonalak és rácshálózat:
call psetc ('map_coastline_colour','mustard')
call psetc ('map_grid_colour','grey')
call pseti ('map_coastline_thickness',6)
call psetc ('map_grid_line_style','dot')
call psetr ('map_grid_latitude_increment',5.0)
call psetr ('map_grid_longitude_increment',5.0)
call pcoast
call psetc ('text_line_1', 'Eszak-Amerika polar-sztereografikus ve
xtuleten')
call psetc ('text_line_2', 'Ez a kovetkezo sor.')
call pseti ('text_line_count', 2)
c A szövegek kiírása:
call ptext
call pclose
end
A polársztereografikus térkép paraméterei:
'subpage_map_area_definition': a polár vetület fajtája (lehet még CORNERS is)'subpage_map_centre_longitude': a középpont földrajzi hosszúsága'subpage_map_centre_latitude': a középpont földrajzi szélessége'subpage_map_vertical_longitude': a vetület vertikális középpontja'subpage_map_scale': a vetület mérete (minél nagyobb, annál nagyobb területet láthatunk)'map_grid_latitude_increment': hosszúsági körök közötti intervallum'map_grid_longitude_increment': szélességi körök közötti intervallumA PCOAST további beállítási lehetőségei, paraméterek: az
ECMWF honlapján.
A szövegek paraméterei:
'text_line_1': a kiírandó szöveg (n=1,10)'text_line_count': hány sorba írja ki a szövegetA szöveg további beállítási lehetőségei, paraméterek: az ECMWF honlapján.
Az eredmény:
program coast_3
call popen
call psetc ('ps_device', 'ps_a4')
call psetc ('ps_file_name', 'coast_3.ps')
c 2 térkép egy A4- es oldalon (portrait formátumban)
c 1. globális térkép színezéssel
c 2. Északi-félteke polársztereografikus vetületben
call psetr ('super_page_y_length', 29.7)
call psetr ('super_page_x_length', 21.)
call psetr ('page_y_length', 14.5)
call psetr ('page_x_length', 21.)
c Az 1. térkép beállításai:
call psetc ('subpage_map_projection','cylindrical')
c Partvonalak tulajdonságai:
call psetc ('map_coastline_land_shade','on')
call psetc ('map_coastline_land_shade_colour','cream')
call psetc ('map_coastline_sea_shade','on')
call psetc ('map_coastline_sea_shade_colour','blue')
call psetc ('map_grid_colour', 'grey')
call psetc ('map_coastline_colour', 'evergreen')
call pcoast
c Szöveg kiírása:
call psetc ('text_line_1', 'Hengervetületi')
call ptext
call pnew("PAGE")
c A 2. térkép beállításai (amiket az elsőhöz képest megváltoztatunk):
call psetc ('subpage_map_projection','polar_stereographic')
call psetc ('map_coastline_land_shade_colour','cream')
call psetc ('map_grid_colour','red')
call pseti ('map_label_latitude_frequency',2)
call psetr ('map_label_height', 0.6)
call pcoast
c Szöveg kiírása:
call psetc ('text_line_1', 'Polár-sztereografikus vetület')
call ptext
call pclose
end
Az új beállítások magyarázatai:
'super_page_y_length','super_page_x_length': a fő oldal hosszúsága és szélessége cm-ben'page_y_length','page_x_length': a kisebb "oldalak" (amire a térképeket rajzoljuk) hosszúsága és szélessége cm-ben'map_coastline_land_shade': kontinensek beszínezése'map_coastline_land_shade_colour':kontinensek színének beállítása'map_coastline_sea_shade': tengerek beszínezése'map_coastline_sea_shade_colour':tengerek színének beállítása'map_grid_colour': a rácsvonalak színe'map_label_latitude_frequency': minden n. szélességi kört feliratoz'map_label_longitude_frequency': minden n. hosszúsági kört feliratoz'map_label_height': a feliratok vastagsága cm-ben (az alapértelmezett 0.2 cm)Az eredmény:
Készítette: Haszpra Tímea, V. éves meteorológus hallgató, 2009/2010 I. félév
Az elektronikus oktatóanyag a Korszerű meteorológiai adatformátumok című
speciális kollégium keretén belül készült.
Oktató: dr. Barcza Zoltán