[gradsusr] lon lat to km

Tue Feb 5 10:54:16 EST 2013

Wendi,
Use gxout contour and use a single color for the lines.  Turn off axis
labels.  The contours you draw will become your new grid.  You can change
the appearance of grid lines using set cstyle and set cthick commands.
You can change the tick mark interval by setting the contour interval.

Jeff

On Tue, Feb 5, 2013 at 6:17 AM, wendi harjupa <wendiharjupa at gmail.com>wrote:

>
> Dear Mr. Jeff
>
> Thank you very much for your explanation,
> I tried that script, but the result (plot) not in KM but still in lon and
> lat
> Here I attached the result (try.gif).
>
> I am sorry, I do not understand about ((lon-100) and ((lat-0) in script
> below,
> because when I change the number, for example ((lat-0) become ((lat-1),
> the result was same.
>
> 'define xdist = ((lon-100)/0.004492369)*(dx/2)'
> 'define ydist = ((lat-0)/0.004522556)*(dy/2)'
>
> I tried to display xdist, but the result was very strange as I attached in
> this email (xdist.gs).
>
> I am very sorry I asked you many time. But please help me to get more
> closer to solve this problem.
>
> Thank you very much
> Best regards
>
>
> On Tue, Feb 5, 2013 at 2:30 AM, Jeff Duda <jeffduda319 at gmail.com> wrote:
>
>> Don't change these four statements from what I have suggested:
>> 'define dx = *6.371e3**cos(-1*(3.14159/180))*cdiff(lon,x)*(3.14159/180)'
>> 'define dy = *6.371e3**cdiff(lat,y)*(3.14159/180))'
>> 'define xdist = ((lon-100)/0.004492369)*(dx/2)'
>> 'define ydist = ((lat-0)/0.004522556)*(dy/2)'
>>
>> In the definitions of dx and dy, the grids lon and lat, and dimension
>> expressions x and y must be entered exactly as that.  Do not replace any of
>> those values with anything else.  The rest is fine.
>>
>> Jeff
>>
>>
>> On Mon, Feb 4, 2013 at 1:49 AM, wendi harjupa <wendiharjupa at gmail.com>wrote:
>>
>>> Dear Mr Jeff
>>>
>>> Thank you very much for your guidance
>>>
>>> Here the script
>>>
>>> 'reinit'
>>> 'open try2.ctl'
>>> 'set grads off'
>>> 'set mproj latlon'
>>> 'set mpdset hires'
>>> 'set gxout shaded'
>>> 'set grid on'
>>> 'set lon 100  101'
>>> 'set lat -1  0'
>>>  'define dx = *6.371e3*
>>> *cos(-1*(3.14159/180))*cdiff(100,0)*(3.14159/180)'
>>> 'define dy = *6.371e3**cdiff(-1,0)*(3.14159/180))'
>>> 'define xdist = ((lon-100)/0.004492369)*(dx/2)'
>>> 'define ydist = ((lat-0)/0.004522556)*(dy/2)'
>>> 'set z 5'
>>> 'set t 1'
>>> 'd ref'
>>>
>>> The script is similar with the script that you suggested.
>>> But, I tried to make changes in many points
>>>
>>> for example :
>>>
>>>  'define dx = *6.371e3*
>>> *cos(-1*(3.14159/180))*cdiff(100,0)*(3.14159/180)'
>>> 'define dy = *6.371e3**cdiff(-1,0)*(3.14159/180))'
>>> 'define xdist = ((100 - 101)/0.004492369)*(dx/2)'
>>> 'define ydist = ((-1-0)/0.004522556)*(dy/2)'
>>>
>>> or I tried to make plot
>>>
>>> 'd xdist'
>>>
>>> But, always I got error message.
>>>
>>> Please teach me more.
>>>
>>> Best regards,
>>> wendi
>>>
>>>
>>> On Mon, Feb 4, 2013 at 3:25 AM, Jeff Duda <jeffduda319 at gmail.com> wrote:
>>>
>>>> What does your script look like now?
>>>>
>>>> Jeff
>>>>
>>>>
>>>> On Sun, Feb 3, 2013 at 6:55 AM, wendi harjupa <wendiharjupa at gmail.com>wrote:
>>>>
>>>>> Dear Mr.  Jeff
>>>>>
>>>>> Thank you very much for your guidance as far..
>>>>> But, I am very sorry for asking you many times..
>>>>> since it is very difficult for me..
>>>>>
>>>>> I tried to change some parts of script, but still this did not run..
>>>>> I got error message as below :
>>>>>
>>>>> Error from CDIFF:  Specified dimension non varying
>>>>> Operation Error:  Error from cdiff function
>>>>>   Error ocurred at column 31
>>>>> DEFINE error:  Invalid expression.
>>>>> *** glibc detected *** grads: corrupted double-linked list: 0x089e2918
>>>>> ***
>>>>> Would you mind to explain more??
>>>>>
>>>>> Thank you very much
>>>>>
>>>>> Best regards,
>>>>> wendi
>>>>>
>>>>> On Sun, Feb 3, 2013 at 1:46 AM, Jeff Duda <jeffduda319 at gmail.com>wrote:
>>>>>
>>>>>> Use this code:
>>>>>>
>>>>>> 'reinit'
>>>>>> 'open try2.ctl'
>>>>>> 'set grads off'
>>>>>> 'set mproj latlon'
>>>>>> 'set mpdset hires'
>>>>>> 'set gxout shaded'
>>>>>> 'set grid on'
>>>>>> 'set lon 100  101'
>>>>>> 'set lat -1  0'
>>>>>> 'define dx = *6.371e3*
>>>>>> *cos(-1*(3.14159/180))*cdiff(100,0)*(3.14159/180)'
>>>>>> 'define dy = *6.371e3**cdiff(-1,0)*(3.14159/180))'
>>>>>> 'define xdist = ((lon-100)/0.004492369)*(dx/2)'
>>>>>> 'define ydist = ((lat-0)/0.004522556)*(dy/2)'
>>>>>> 'set z 5'
>>>>>> 'set t 1'
>>>>>> 'd ref'
>>>>>>
>>>>>> The square brackets were just to show different groupings.  When it
>>>>>> comes time to code it, you must use only parentheses.
>>>>>>
>>>>>> The above should work.
>>>>>>
>>>>>> Jeff
>>>>>>
>>>>>> On Sat, Feb 2, 2013 at 1:41 AM, wendi harjupa <wendiharjupa at gmail.com
>>>>>> > wrote:
>>>>>>
>>>>>>> Dear Mr. Jeff
>>>>>>>
>>>>>>> Thank you very much for your explanation,
>>>>>>> but I am very sorry, I still do not get it,
>>>>>>>
>>>>>>> I made new script
>>>>>>>
>>>>>>> 'reinit'
>>>>>>> 'open try2.ctl'
>>>>>>> 'set grads off'
>>>>>>> 'set mproj latlon'
>>>>>>> 'set mpdset hires'
>>>>>>> 'set gxout shaded'
>>>>>>> 'set grid on'
>>>>>>> 'set lon 100  101'
>>>>>>> 'set lat -1  0'
>>>>>>> 'define dx = *6.371e3*
>>>>>>> *cos(-1*(3.14159/180))*cdiff(100,0)*(3.14159/180)'
>>>>>>> 'define dy = *6.371e3**cdiff(-1,0)*(3.14159/180))'
>>>>>>> 'define xdist = ((lon-100)/0.004492369)*(dx/2]'
>>>>>>> 'define ydist = ((lat-0)/0.004522556)*dy/2]'
>>>>>>> 'set z 5'
>>>>>>> 'set t 1'
>>>>>>> 'd ref'
>>>>>>>
>>>>>>> I checked grads document and I found that 6.371e3 is radius of the
>>>>>>> earth (6371 to 6374 km)*.*
>>>>>>> *
>>>>>>> *
>>>>>>> *I tried to change "*6.371e3" to "0e10" which is mean 0 to 10 km)
>>>>>>> But still the program did not run normally.
>>>>>>>
>>>>>>> I changed the value of lon and lat in script with 100 and -1.
>>>>>>> but still I got error.
>>>>>>>
>>>>>>> Please help me to solve it.
>>>>>>>
>>>>>>> Thank you very much
>>>>>>> Best regards,
>>>>>>> wendi
>>>>>>>
>>>>>>> On Sat, Feb 2, 2013 at 2:21 AM, Jeff Duda <jeffduda319 at gmail.com>wrote:
>>>>>>>
>>>>>>>> Wendi,
>>>>>>>> You need to run those two define commands at the top of my previous
>>>>>>>> email before attempting to define the xdist and ydist fields.  Also, since
>>>>>>>> I can see your control file information, I can see that dlon is 0.004492369
>>>>>>>> and dlat is  0.004522556.  So substitute those two values into your
>>>>>>>> expressions for xdist and ydist also.
>>>>>>>>
>>>>>>>> Jeff
>>>>>>>>
>>>>>>>>
>>>>>>>> On Fri, Feb 1, 2013 at 4:36 AM, wendi harjupa <
>>>>>>>> wendiharjupa at gmail.com> wrote:
>>>>>>>>
>>>>>>>>> Dear Mr. Jeff
>>>>>>>>>
>>>>>>>>> I have tried to use the scipt which you suggested,
>>>>>>>>> But, still I could not do it. Would you help me more...
>>>>>>>>>
>>>>>>>>> here I send my gs file :
>>>>>>>>>
>>>>>>>>>  'reinit'
>>>>>>>>> 'open try2.ctl'
>>>>>>>>> 'set grads off'
>>>>>>>>> 'set mproj latlon'
>>>>>>>>> 'set mpdset hires'
>>>>>>>>> 'set gxout shaded'
>>>>>>>>> 'set grid on'
>>>>>>>>> 'set lon 100  101'
>>>>>>>>> 'set lat 0  1'
>>>>>>>>> 'define xdist = [(lon 99.7)/dlon]*[dx/2]'
>>>>>>>>> 'define ydist = [(lat -1.02)/dlat]*[dy/2]'
>>>>>>>>> 'set z 5'
>>>>>>>>> 'set t 1'
>>>>>>>>> 'd ref'
>>>>>>>>>
>>>>>>>>> Below is my ctl file :
>>>>>>>>>
>>>>>>>>> DSET 1013.dat
>>>>>>>>> TITLE  RADAR
>>>>>>>>> OPTIONS  TEMPLATE LITTLE_ENDIAN
>>>>>>>>> UNDEF  -1000
>>>>>>>>> XDEF    321 LINEAR    99.688120996    0.004492369
>>>>>>>>> YDEF    321 LINEAR    -1.083708891    0.004522556
>>>>>>>>> ZDEF    39  LINEAR   1.0 0.5
>>>>>>>>> TDEF    15 LINEAR 13:00Z10apr2004  4mn
>>>>>>>>> VARS    2
>>>>>>>>> ref 39  99  Reflectivity(dBZ)
>>>>>>>>> vel 39  99  Doppler Velocity(m/s)
>>>>>>>>> ENDVARS
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Would you check my gs file, is there something that I have to
>>>>>>>>> revise??
>>>>>>>>>
>>>>>>>>> Thank you very much.
>>>>>>>>>
>>>>>>>>> Best regards,
>>>>>>>>> wendi
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Fri, Feb 1, 2013 at 10:34 AM, Jeff Duda <jeffduda319 at gmail.com>wrote:
>>>>>>>>>
>>>>>>>>>> These commands will help.
>>>>>>>>>>
>>>>>>>>>> 'define dx =
>>>>>>>>>> 6.371e3*cos(lat*(3.14159/180))*cdiff(lon,x)*(3.14159/180)'
>>>>>>>>>> 'define dy = 6.371e3*cdiff(lat,y)*(3.14159/180))'
>>>>>>>>>>
>>>>>>>>>> They define grid increments in terms of spatial distance rather
>>>>>>>>>> than geographic distance.  If you're trying to change the axes and tick
>>>>>>>>>> marks, you'll have to define an origin and define your x and y distances
>>>>>>>>>> based on that.  For instance, if you set your lower left domain corner
>>>>>>>>>> (latitude 0 and longitude 100), then you can define your coordinate
>>>>>>>>>> distances using
>>>>>>>>>>
>>>>>>>>>> 'define xdist = [(lon-100)/dlon]*[dx/2]'
>>>>>>>>>> 'define ydist = [(lat-0)/dlat]*[dy/2]'
>>>>>>>>>>
>>>>>>>>>> I think this should work, but I haven't tested it myself, so I
>>>>>>>>>> don't know.  dlon and dlat are the step values from the XDEF and YDEF lines
>>>>>>>>>> of the control file, respectively.
>>>>>>>>>>
>>>>>>>>>> Jeff Duda
>>>>>>>>>>
>>>>>>>>>> On Thu, Jan 31, 2013 at 7:02 AM, wendi harjupa <
>>>>>>>>>> wendiharjupa at gmail.com> wrote:
>>>>>>>>>>
>>>>>>>>>>> Dear Sir/Madam
>>>>>>>>>>>
>>>>>>>>>>> I need help
>>>>>>>>>>> I have problem with grads command,
>>>>>>>>>>> I am not well understand with command gr2XY,
>>>>>>>>>>>
>>>>>>>>>>> Actually I want to convert the lon 100  to 101, and lat 0 to 1
>>>>>>>>>>> to distance (km)
>>>>>>>>>>> let think that lon 100 is 0 km, and lat 0 is 0 km,
>>>>>>>>>>>
>>>>>>>>>>> what command I can use to convert lon and lat to km?
>>>>>>>>>>>
>>>>>>>>>>> Thank you very much for your help
>>>>>>>>>>>
>>>>>>>>>>> best regards,
>>>>>>>>>>> --
>>>>>>>>>>> -------------
>>>>>>>>>>> Wendi Harjupa. ST
>>>>>>>>>>> Shimane University Graduate School of Engineering Sciences
>>>>>>>>>>> Department of Electronic Control Systems
>>>>>>>>>>> Remote Sensing Laboratory
>>>>>>>>>>> s119447
>>>>>>>>>>> Cellphone : +81-080-4268-6676
>>>>>>>>>>> email : wendi at rslab.riko.shimane-u.ac.jp
>>>>>>>>>>>
>>>>>>>>>>> _______________________________________________
>>>>>>>>>>> gradsusr mailing list
>>>>>>>>>>> gradsusr at gradsusr.org
>>>>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> --
>>>>>>>>>> Jeff Duda
>>>>>>>>>> Graduate research assistant
>>>>>>>>>> University of Oklahoma School of Meteorology
>>>>>>>>>> Center for Analysis and Prediction of Storms
>>>>>>>>>>
>>>>>>>>>> _______________________________________________
>>>>>>>>>> gradsusr mailing list
>>>>>>>>>> gradsusr at gradsusr.org
>>>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> --
>>>>>>>>> -------------
>>>>>>>>> Wendi Harjupa. ST
>>>>>>>>> Shimane University Graduate School of Engineering Sciences
>>>>>>>>> Department of Electronic Control Systems
>>>>>>>>> Remote Sensing Laboratory
>>>>>>>>> s119447
>>>>>>>>> Cellphone : +81-080-4268-6676
>>>>>>>>> email : wendi at rslab.riko.shimane-u.ac.jp
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> gradsusr mailing list
>>>>>>>>> gradsusr at gradsusr.org
>>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> --
>>>>>>>> Jeff Duda
>>>>>>>> Graduate research assistant
>>>>>>>> University of Oklahoma School of Meteorology
>>>>>>>> Center for Analysis and Prediction of Storms
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> gradsusr mailing list
>>>>>>>> gradsusr at gradsusr.org
>>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> --
>>>>>>> -------------
>>>>>>> Wendi Harjupa. ST
>>>>>>> Shimane University Graduate School of Engineering Sciences
>>>>>>> Department of Electronic Control Systems
>>>>>>> Remote Sensing Laboratory
>>>>>>> s119447
>>>>>>> Cellphone : +81-080-4268-6676
>>>>>>> email : wendi at rslab.riko.shimane-u.ac.jp
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> gradsusr mailing list
>>>>>>> gradsusr at gradsusr.org
>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>> --
>>>>>> Jeff Duda
>>>>>> Graduate research assistant
>>>>>> University of Oklahoma School of Meteorology
>>>>>> Center for Analysis and Prediction of Storms
>>>>>>
>>>>>> _______________________________________________
>>>>>> gradsusr mailing list
>>>>>> gradsusr at gradsusr.org
>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> -------------
>>>>> Wendi Harjupa. ST
>>>>> Shimane University Graduate School of Engineering Sciences
>>>>> Department of Electronic Control Systems
>>>>> Remote Sensing Laboratory
>>>>> s119447
>>>>> Cellphone : +81-080-4268-6676
>>>>> email : wendi at rslab.riko.shimane-u.ac.jp
>>>>>
>>>>> _______________________________________________
>>>>> gradsusr mailing list
>>>>> gradsusr at gradsusr.org
>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> Jeff Duda
>>>> Graduate research assistant
>>>> University of Oklahoma School of Meteorology
>>>> Center for Analysis and Prediction of Storms
>>>>
>>>> _______________________________________________
>>>> gradsusr mailing list
>>>> gradsusr at gradsusr.org
>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>
>>>>
>>>
>>>
>>> --
>>> -------------
>>> Wendi Harjupa. ST
>>> Shimane University Graduate School of Engineering Sciences
>>> Department of Electronic Control Systems
>>> Remote Sensing Laboratory
>>> s119447
>>> Cellphone : +81-080-4268-6676
>>> email : wendi at rslab.riko.shimane-u.ac.jp
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>>
>>
>>
>> --
>> Jeff Duda
>> Graduate research assistant
>> University of Oklahoma School of Meteorology
>> Center for Analysis and Prediction of Storms
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> -------------
> Wendi Harjupa. ST
> Shimane University Graduate School of Engineering Sciences
> Department of Electronic Control Systems
> Remote Sensing Laboratory
> s119447
> Cellphone : +81-080-4268-6676
> email : wendi at rslab.riko.shimane-u.ac.jp
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>

-- 
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130205/118cc869/attachment-0003.html 