You can try MeteoInfo/MeteoInfoLab to plot the meteorological data. MeteoInfo is free and open source (http://www.meteothinker.com/).

Yaqiang Wang, thank you for the alternative solution. I will try to use your suggestion. Thank you very much. Best Regards.

I use McIDAS to plot various meteorological data (grib, satellite, ..).

http://www.ssec.wisc.edu/mcidas/software/v/

It is open source and available on Linux and Windows platforms. Various tutorials on how to use it. Regularly updated and mantained.

Best,

Stefano

Dear Stefano Costa, thank you very much for your suggestion.

I will try this software after I used MeteoInfoLab. Best Regards.

Dear Yaqiang Wang,

Now I try to use MeteoInfoLab since I found much software related to the netCDF data format (http://www.unidata.ucar.edu/software/netcdf/software.html#ANDX), however, many of them did not work in Window OS (Windows 10).

Now, I used the data from http://www.esrl.noaa.gov/psd/data/composites/day/ for wind vector and the result can be seen in the attached file (Fig. (a)). As a comparison in Fig. (b), the graph was plotted using LDI commercial software for geopotential height and wind vector as well as SO2 (Kim et al., 2016).

I need more time to familiar with this software, but can you tell me in the short comment that this software can visualize the data u-wind and v-wind as a wind vector, and can I also added geopotential height in the same graphs?

Thank you very much for your help.

Best Regards,

Indra Chandra

An example MeteoInfoLab script:

hfn = 'C:/Temp/compday.5XTNsoLwus.nc'

ufn = 'C:/Temp/compday.5YkRep4hMt.nc'

vfn = 'C:/Temp/vcompday.5YkRep4hMt.nc'

hf = addfile(hfn)

uf = addfile(ufn)

vf = addfile(vfn)

hgt = hf['hgt'][0,[20,50],[100,150]]

u = uf['uwnd'][0,[20,50],[100,150]]

v = vf['vwnd'][0,[20,50],[100,150]]

#Plot

axesm()

lworld = shaperead('D:/Temp/map/country1.shp')

geoshow(lworld, edgecolor='k')

lhgt = contourfm(hgt, 20)

lwind = quiverm(u, v, color='b')

quiverkey(lwind, 0.85, 0.18, 10, bbox={'edge':True, 'fill':True})

colorbar(lhgt, extendrect=False)

title('Geopotential height and wind vector at 850hPa')

#savefig('C:/Temp/uv_test.png')

Dear Yaqiang Wang,

I would like to thank you for the script. Fig. (c) was visualized by MeteoInfoLab software which uses your script. The geopotential height result (Fig. (c) was similar with NCEP/NCAR Reanalysis (Fig. (b)) and Kim et al. (2016) paper (Fig. (d)). However, regarding the wind vector, it was a little bit different as shown in Fig. (a), (c), and (d). I will ask the raw data to the authors of this paper, but, how about the difference (wind speed (length of the arrow) and wind direction (arrowhead)) between Fig. (a) and (c)? 

Thank you very much.

Best Regards,

Indra Chandra

Please make sure you have downloaded correct uwind/vwind data (time, level).

Dear Yaqiang Wang,

Yes, I will confirm the raw data from the NCEP/NCAR Reanalysis and the author of this paper. Thank you very much.

Best Regards,

Indra Chandra

I suggest two different packages:

* OpenGrADS (opegrads.org). Just load the files you've downloaded and try one of the examples. You've also lots of scripts to run (grads recipes). Is perhaps the most used to plot model outputs

* NCL (https://www.ncl.ucar.edu) The most powerful language to plot any output from a model. It needs some time to learn but it's very well documented and you can ask to a very active and helpful community

Luis Cana, thank you for your suggestion.

Best Regards, Indra Chandra.

Hi, I see you are (or were) working with ERA-Interim. I am too, and I need to convert the special lat/long coordinates of one dataset back into the regular lat/long (-180/180 for long, -90/90 for lat, Reference Meridian=Greenwich). So far I have not been able to find a proper way to make this conversion. Are you perhaps able to give me any suggestion?

Thanks,

Francesco

Francesco, you have several examples here if you're working with NCL:

http://www.ncl.ucar.edu/Applications/era40.shtml

It worths a try

Luis

http://www.ncl.ucar.edu/Applications/era40.shtml

My dataset is in NETCDF format though.

Netcdf  is ok for working with NCL