Hi all, I need a suggestion, i got geomagnetic and geographic coordinates for lats and lons.
How can I use map_set in geomagnetic coordinates system
Geographic I used map_set correct, I got a right map, but geomagntical, stations(cities) are located in a wrong place:
map_set, 72.5,-45,/orthographic,/isotropic, limit=[55,-75,90,-15],/label,latlab=0,lonlab=0,/continents,/ grid
;geomagnetic coordinates
;lons=[-42.68,-38.56,-37.56,-36.66,-40.25,-37.56,-38.72,-41. 87,-39.46,-41.11,-32.05,-28.35]
;
;lats=[65.22,66.92,69.52,70.98,72.18,73.56,74.83,75.97,78.57 ,80.35,82.71,84.46]
;geographic coordinates
lons=[314.56,310.32,308.27,307.10,309.28,306.43,306.47,307.8 7,303.85,302.82,294.90,290.77]
lats=[61.16,62.00,64.17,65.42,67.02,67.93,69.25,70.68,72.78, 74.57,76.02,77.47]
!p.charsize=0.8
cgplot,lons,lats,psym=6,/overplot
cities=['NAQ','FHB','GHB','SKT','STF','ATU','GDH','UMQ','UPN ','KUV','SVS','THL']
citycolor=['lawn green','blue','royal blue','lime green','green','crimson','chartreuse','red','tomato', $
'orange red','green','navy']
for j=0,n_elements(cities)-1 do cgtext,lons(j),lats(j),cities(j),color=citycolor(j)
Any kind of suggestion is welcome
Many thanks
Zolile
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
