tag:blogger.com,1999:blog-4834004831115556445.post5276141359860454961..comments2024-03-13T02:05:52.877-04:00Comments on GrADS-aholic!: Script: cone.gs; Draws "Cone of Uncertainty" for projected storm path using the MFHILO functionUnknownnoreply@blogger.comBlogger16125tag:blogger.com,1999:blog-4834004831115556445.post-27590900548748489102022-09-29T00:13:42.837-04:002022-09-29T00:13:42.837-04:00Hi KC, may I get your code? I got the same error Hi KC, may I get your code? I got the same error Berita Cuacahttps://www.blogger.com/profile/08575345742441976582noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-8690956567027879292016-03-18T17:29:14.911-04:002016-03-18T17:29:14.911-04:00I am attempting to use the mflhilo functions but c...I am attempting to use the mflhilo functions but cannot get it to work. I realize I need to enable the user defined extensions but am failing at getting this setup. I have grads v2.1.a3. Any help would be appreciatedAnonymousnoreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-1418157266903786062015-05-04T02:26:38.106-04:002015-05-04T02:26:38.106-04:00Ok.thanks,i will use the set rgb function in grads...Ok.thanks,i will use the set rgb function in grads 2.1 above. Maybe i can modify your cone.gs to show model date in every point the cone plotKChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-81692844881180157162015-05-04T02:26:20.417-04:002015-05-04T02:26:20.417-04:00Ok.thanks,i will use the set rgb function in grads...Ok.thanks,i will use the set rgb function in grads 2.1 above. Maybe i can modify your cone.gs to show model date in every point the cone plotKChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-25811086445694678072015-04-30T07:37:03.525-04:002015-04-30T07:37:03.525-04:00Sorry I didn't get back to your earlier, I am ...Sorry I didn't get back to your earlier, I am glad you solved the problem and that it works okay! Unfortunately, one of the deficiencies of GrADS is that it does not do transparency unless you have version 2.1 or newer. If you have 2.1 or newer try following the last example on this page. I know it's for shapefiles, but it may work for the polygons as well. <br /><br />http://www.iges.org/grads/gadoc/shapefiles.html<br />Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-59217603274382193122015-04-27T13:36:07.277-04:002015-04-27T13:36:07.277-04:00Hi there,
I managed to solve the problems and had...Hi there,<br /><br />I managed to solve the problems and had successfully plot the cone graphic. Just the -tint and the -maxtime options. But do you know how can i control the transperacy of filled colour?KChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-22244510953847131892015-04-26T23:24:24.338-04:002015-04-26T23:24:24.338-04:00Hi there,
Did you manage to resolve the errors in...Hi there,<br /><br />Did you manage to resolve the errors in my scripts?. Your opinion to solve this issue is highly appreciated.KChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-10788452800369619882015-04-22T22:43:15.205-04:002015-04-22T22:43:15.205-04:00The output :
Tambahan KC :
4.49
101.67
6
Is just...The output : <br />Tambahan KC :<br />4.49<br />101.67<br />6<br /><br />Is just additional lines that i print to check values..<br />KChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-49832555906043208162015-04-22T22:31:36.728-04:002015-04-22T22:31:36.728-04:00Here is my script :
*'sdfopen 'file
'...Here is my script :<br /><br />*'sdfopen 'file<br />'open jma_ssm.ctl'<br />'set display color white'<br />'clear'<br />'set gxout contour'<br />'set mpdset hires'<br />'set t 1'<br />'set lon 92 120'<br />'set lat -6 9'<br />'d PRMSLmsl/100'<br />*'rgbset'<br /><br />'basemap L 74 1 M'<br />'basemap O 45 1 M'<br /><br />'printim indera.png'<br /><br />'clear'<br />'set t 5'<br />'set lon 92 120'<br />'set lat -6 9'<br />*'set lat 5 30'<br />*'set lon -80 1'<br />*'set lev 850'<br />'set gxout contour'<br />'set clevs 10000000'<br />'set ccolor 0'<br />'d PRMSLmsl/100'<br />'set map 1'<br />'draw map'<br />*'cone PRMSLmsl -e 50 -tunit day -dunit km -e_reduct 8 -end 3 -tint 4 -fcol 2'<br />'cone prmslmsl/100 -e 50 -mnlt 4 -maxlat 7 -mnln 94 -mxln 114 -tunit50 -tunit day -dunit km -fcol 2'<br />*-mnlt 10 -maxlat 25 -mnln -70 -mxln -20 -tunit<br />'printim ucone.png -b indera.png -t 0'<br />KChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-60095797142342323552015-04-22T22:30:23.919-04:002015-04-22T22:30:23.919-04:00I modified the script but still received the follo...I modified the script but still received the following error message :<br />ga-> cone_Indera.gs<br />mfhilo CL method N: 1 maxmin: -1 radinf: 1000.0 [nm] cintinf: 300.0 latc: 5.50 lonc: 104.00<br />L 4.49 101.67 M: 948 D: 151.86<br />Tambahan KC :<br />4.49<br />101.67<br />6<br />Data Request Warning: Request is completely outside file limits<br />mfhilo CL method N: 0 maxmin: -1 radinf: 1000.0 [nm] cintinf: 300.0 latc: 5.50 lonc: 104.00<br />Tambahan KC :<br />4.49 grid<br />101.67 contents<br />6<br />Query Error: Syntax is QUERY W2XY Lon Lat<br />Non-numeric args to numeric operation<br /> Error occurred on line 341<br /> In file /opt/opengrads/Contents/Resources/Scripts/cone.gs<br /><br />KChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-57115766685179060332015-04-22T22:27:42.658-04:002015-04-22T22:27:42.658-04:00ok..then let me first look at your comments and tr...ok..then let me first look at your comments and try to resolve the issues...<br />KChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-35187304523100635402015-04-22T10:41:23.936-04:002015-04-22T10:41:23.936-04:00I'll have to look more intently at the script ...I'll have to look more intently at the script to give you a definitive answer, but at first glance is appears that there is an issue when the script calls the "mflhilo" function, which is the function use to find the location of the minimum SLP. The outside of file limits request makes me think you are setting your boundaries wrong, I most commonly do this when I try to set a negative longitude when the data is gridded from 0-360, not -180-180. <br /><br />Also it looks like there is a weird call in this : <br />'cone prmslmsl/100 -e 50 -tunit day -dunit km ucone.png-fcol 2'<br /><br />I don't think the ucone.png should be in there (though honestly, it's been a while since I dug around this script so its entirely possible I put something looking for an image in there...).<br /><br />I will glance at the script later and maybe I can give you a better answer. But at first look, those are my thoughts.<br />Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-53465005967344995182015-04-22T00:15:53.421-04:002015-04-22T00:15:53.421-04:00By the way, below are the error messages i got :
...By the way, below are the error messages i got :<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br /> ******nlo exceeds nmin=2000<br />mfhilo CL method N: 1 maxmin: -1 radinf: 1000.0 [nm] cintinf: 300.0 latc: 1.65 lonc: 106.28<br />L -0.59 96.96 M: 1010.51 D: 937.54<br />Data Request Warning: Request is completely outside file limits<br />mfhilo CL method N: 0 maxmin: -1 radinf: 1000.0 [nm] cintinf: 300.0 latc: 1.65 lonc: 106.28<br />Query Error: Syntax is QUERY W2XY Lon Lat<br />Non-numeric args to numeric operation<br /> Error occurred on line 337<br /> In file /opt/opengrads/Contents/Resources/Scripts/cone.gs<br /><br /><br />Do you have any idea?. Your prompt comment on these is highly appreciated.<br /><br />Thanks,<br />KChangKChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-62851966124524684282015-04-22T00:14:44.489-04:002015-04-22T00:14:44.489-04:00Hi there ,
Great GrADS script to plot Cone Graphic...Hi there ,<br />Great GrADS script to plot Cone Graphic. I try to use that script to plot cone graphic for my imaginary typhoon using my storm surge model. Below are the best track data i'd created :<br /><br />66666 1501 16 0002 1501 0 6 INDERA 20150101<br />15051500 002 2 65 1131 1000 030 00000 0000 60070 0070<br /><br />15051512 002 3 47 1102 998 055 00000 0000 60070 0070<br /><br />15051600 002 4 51 1075 990 085 00000 0000 60100 0080 15051612 002 4 46 1039 970 091 00000 0000 60100 0080<br /><br />15051700 002 5 43 1015 948 128 00000 0000 60120 0090<br />15051712 002 5 54 0973 930 141 00000 0000 60130 0100<br />15051800 002 5 59 0954 950 125 00000 0000 60150 0120<br /><br />15051812 002 4 58 0951 998 078 00000 0000 60150 0120<br /><br />It uses the JMA RSMC Best Track text data format.<br /><br />So i use my control file :<br />dset ^jma_ssm.grb<br />index ^jma_ssm.grb.idx<br />undef 9.999E+20<br />title jma_ssm.grb<br />* produced by grib2ctl v0.9.13<br />* command line options: -verf jma_ssm.grb<br />dtype grib 255<br />ydef 901 linear -6.000000 0.017<br />xdef 1681 linear 92.000000 0.017000<br />tdef 15 linear 00Z15may2015 6hr<br />zdef 1 linear 1 1<br />vars 4<br />DSLMsfc 0 82,1,0 ** surface Deviation of sea level from mean [m]<br />PRMSLmsl 0 2,102,0 ** mean-sea level Pressure reduced to MSL [Pa]<br />UGRD10m 0 33,105,10 ** 10 m above ground u wind [m/s]<br />VGRD10m 0 34,105,10 ** 10 m above ground v wind [m/s]<br />ENDVARS<br /><br />Here are the result when i run the "q file" in grads :<br />ga-> q file<br />File 1 : jma_ssm.grb<br /> Descriptor: jma_ssm.ctl<br /> Binary: jma_ssm.grb<br /> Type = Gridded<br /> Xsize = 1681 Ysize = 901 Zsize = 1 Tsize = 15 Esize = 1<br /> Number of Variables = 4<br /> dslmsfc 0 82 ** surface Deviation of sea level from mean [m]<br /> prmslmsl 0 2 ** mean-sea level Pressure reduced to MSL [Pa]<br /> ugrd10m 0 33 ** 10 m above ground u wind [m/s]<br /> vgrd10m 0 34 ** 10 m above ground v wind [m/s]<br /><br />I use script file cone_indera.gs as below :<br />*'sdfopen 'file<br />'open jma_ssm.ctl'<br />'set display color white'<br />'clear'<br />'set gxout contour'<br />'set mpdset hires'<br />*'set lon 92 120'<br />*'set lat -6 9'<br />'d dslmsfc'<br />*'rgbset'<br /><br />'basemap L 74 1 M'<br />'basemap O 45 1 M'<br /><br />'printim indera.png'<br /><br />'clear'<br />'set t 3'<br />*'set lon 92 120'<br />*'set lat -6 9'<br />*'set lat 5 30'<br />*'set lon -80 1'<br />*'set lev 850'<br />'set gxout contour'<br />*'set clevs 10000000'<br />'set ccolor 0'<br />'d PRMSLmsl/100'<br />'set map 1'<br />'draw map'<br /><br /><br />*'cone PRMSLmsl -e 50 -tunit day -dunit km -e_reduct 8 -end 3 -tint 4 -fcol 2'<br />'cone prmslmsl/100 -e 50 -tunit day -dunit km ucone.png-fcol 2'<br />*'printim ucone.png -b indera.png -t 0'<br /><br />Do you have any idea?. Your prompt comment on these is highly appreciated.<br /><br />Thanks,<br />KChangKChttps://www.blogger.com/profile/12199822834309858214noreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-31547164216039219132014-04-14T07:36:27.041-04:002014-04-14T07:36:27.041-04:00Are you talking about the minlat/maxlat values for...Are you talking about the minlat/maxlat values for this script? Or are you talking about finding the lat/lon coordinates of your low pressure center. <br /><br />If you are talking about the values in this script, the values are your values to choose. What they do is put a "bounding" box around the region you are interested in, so that you avoid jumping from one low pressure area to another. If you used this script over an larger domain, there is a greater chance that you will get crazy results from jumping around to different low pressure areas. But if you restrict the region you are focused on, you will simply track the same LPA. So choose these values as you need them.<br /><br />If you are looking for a way to find the coordinates of you minimum pressure, you can just dig around through this script. The function used is the MFHILO function, and that finds your minimum value in a given domain. I believe you could also nest a "minloc" function in x and y to find the coordinates as well, though as I recall this was less intuitive than using the mfhillo method.<br /><br />Hope this answers your question!<br />Good luck!Anonymousnoreply@blogger.comtag:blogger.com,1999:blog-4834004831115556445.post-4113561037254708652014-04-13T09:42:34.109-04:002014-04-13T09:42:34.109-04:00Good day. I've been using GrADS for quite some...Good day. I've been using GrADS for quite some time now, yet I'm still new to the program. I'm having difficulties on maxlt/mnlat etc. variables. I've been tracking LPAs and tropical cyclones based on the Philippines using GrADS. I don't know the values for the maxlt/mnlat etc. to be used for my location. Hope you would help me with that. Thank you for your further response. God Bless!Anonymousnoreply@blogger.com