grads画台风路径（Grads draw typhoon tracks）.docVIP

grads画台风路径（Grads draw typhoon tracks） Draw the path of Typhoon reinit; * first - from SLP (or 1000hPa height) to find the location of typhoon center, the output to a file RM! Weizhi-w.txt; * is written, it must first remove old data files open 928s.ctl; *WRF prediction results I=1 While (i29); * should use =, but the last position is wrong, not set t I T=subwrd (result, 4); * date Tc=substr (T, 8,5); * partial date string d * slvl; sea level pressure define locx=minloc (min (SLVL, lat=10, lat=22), lon=106, lon=121 ); * search for the lowest pressure longitude define locy=minloc (min (SLVL, lon=106, lon=121), lat=10, lat=22 ); * search for the lowest Latitude Atmospheric Pressure d locx Locxx=subwrd (result, 4) d locy Locyy=subwrd (result, 4) q gr2wlocxxlocyy; *GRADS coordinates into geodetic coordinates X=subwrd (result, 3) Y=subwrd (result, 6) Rc=write (weizhi-w.txt, Xy TC, append); * the latitude and longitude of the location and time stored in the file *say RC I=i+4; * 08 and 20 output only when the data interval of 3 hours Endwhile Close (weizhi-w.txt); * closing output file CP! Weizhi-w.txt weizhi-r.txt; GS * in the same file can not read and write, so make 1 copies reinit; * second - from the data file (you can also edit their own reading position) draw the path of Typhoon open 928s.ctl set mproj latlon set mpdset hires; * set high resolution map set lon 105 118; * coordinate display part missing too lax set lat 15.5 23.5 set Cmax 0; * to show the maximum design, is to map out set xlint 2; * definition of abscissa marker intervals set ylint 2; * define the longitudinal coordinate marking interval set xlopts 14 0.22; * abscissa color width size set ylopts 14 0.20; * ordinate color width size set clopts -1 -1 0.20; * contour color marker size width * use the above definition, basically meet the layout requirements -- about 5*7CM axis label is about 6 character size d slvlis a map with only *; Read_file=read (weizhi-r.txt); * open the file for reading Read_line=