Search found 49 matches

by yaqiang
November 1st, 2018, 4:06 am
Forum: MeteoInfo Software
Topic: How to calculate potential vorticity (PV) along the trajectory
Replies: 2
Views: 1463

Re: How to calculate potential vorticity (PV) along the trajectory

You can calculate PV array from the existing variables and than interpolate to the end point positions of the trajectories. There is a MPV example online (http://meteothink.org/examples/meteoinfolab/meteo_analysis/moisture_potential_vorticity.html). Replace pseudo-equivalent potential temperature to...
by yaqiang
October 26th, 2018, 11:31 am
Forum: MeteoInfo Software
Topic: Web domain changed to meteothink.org
Replies: 4
Views: 1633

Re: Web domain changed to meteothink.org

Great, thanks!
by yaqiang
October 19th, 2018, 9:18 pm
Forum: MeteoInfo Software
Topic: Web domain changed to meteothink.org
Replies: 4
Views: 1633

Re: Web domain changed to meteothink.org

Glenn,

Thanks for your quickly response on it! Also TrajStat links need to be updated on the READY web site.

Regards
Yaqiang
by yaqiang
October 19th, 2018, 12:52 pm
Forum: MeteoInfo Software
Topic: Read trajectory data as 2D array
Replies: 0
Views: 1644

Read trajectory data as 2D array

The trajectory endpoint data file can be read as 2D arrays following netCDF data model from MeteoInfo version 1.6. The dimensions are 'trajectory' (trajectory number) and 'obs' (endpoint number of one trajectory). Detail and example can be found here: http://www.meteothink.org/examples/meteoinfolab/...
by yaqiang
October 19th, 2018, 12:19 pm
Forum: MeteoInfo Software
Topic: Web domain changed to meteothink.org
Replies: 4
Views: 1633

Web domain changed to meteothink.org

The previous meteothinker.com domain is no appropriate for MeteoInfo software due to its open source nature, and the previous web hosting was not stable from several weeks ago. So I moved the web site to GoDaddy with a new domain of meteothink.org (http://meteothink.org/). Please revise the web link...
by yaqiang
October 19th, 2018, 10:50 am
Forum: Dispersion Model
Topic: ASCII to cdump file
Replies: 3
Views: 1053

Re: ASCII to cdump file

You may try MeteoInfoLab script to get the max value and location from cdump file, and plot it on map. #Read data array from HYSPLIT concentration output data fn = r'D:\Working\MIScript\JapanNuclear\cdump15' f = addfile_hyconc(fn) data = f['C137'][:,0,:,:] # 3 dimension array - time, lat, lon #Get t...
by yaqiang
August 23rd, 2018, 3:42 am
Forum: Users
Topic: mixing height data
Replies: 1
Views: 950

Re: mixing height data

If you mean to extract PBLH data at a specific location from existing meteorological data (such as GDAS data), you can follow the example here: http://www.meteothinker.com/examples/me ... rp_st.html .
interpn_st.png
by yaqiang
August 8th, 2018, 8:32 pm
Forum: Conversion programs
Topic: Connverting tgs file to shape file
Replies: 1
Views: 1496

Re: Connverting tgs file to shape file

Firstly you should make sure that you have calculated trajectories correctly. Or you may got empty .tgs file so no .shp file can be created. You can learn MeteoInfo and TrajStat from my presentation in 2018 HYSPLIT workshop (http://www.meteothinker.com/docs/introd ... esentation).
by yaqiang
May 31st, 2018, 11:39 pm
Forum: MeteoInfo Software
Topic: TrajStat 1.4.5 was released
Replies: 0
Views: 3739

TrajStat 1.4.5 was released

Add Check point number button in cluster calculation dialog to check the end point number of each trajectory. Add 3D trajectory chart to explore the trajectories in 3D view. Add Save to grid menu item to save PSCF/CWT data in a Sufer ASCII grid data file, which can be loaded and plot as a contour la...
by yaqiang
March 19th, 2018, 9:08 am
Forum: General questions on meteorological data needed for HYSPLIT
Topic: Clarification on meteorological data
Replies: 20
Views: 9643

Re: Clarification on meteorological data

You can use MeteoInfoMap to explore the ARL data file. Or use MeteoInfoLab script to extract PBLH data at a specific location from GDAS data. An example script: lat = 40.0 lon = 120.5 fn = 'D:/Temp/arl/20140301_gdas0p5' f = addfile(fn) pblh = f['PBLH'][:,str(lat),str(lon)] #Plot plot(pblh, '-bo') yl...