请帮忙查看程序，如何用GRads画规则圆和不规则圆

雾由心生

不好意思，还没掌握上传图片。
很明显在画内圆中，bear=0与bear=360时应该在同一个点上，我觉得问题出在第二个循环的内循环中变量att并没有随着tlat和tlon实时变化

mofangbao

雾由心生 发表于 2013-5-12 15:42

                               
登录/注册后可看大图

图片在附件中，修改过的程序画出的是图片a，但是应该画出类似于图片b

放弃你现在的脚本把  画这种图很简单，大圆可以用我前两天发帖的那个画椭圆的脚本，画不规则的，只要你知道了坐标，用draw line画出来就行了，甚至可以用draw mark点出来

雾由心生

mofangbao 发表于 2013-5-12 16:02

                               
登录/注册后可看大图

放弃你现在的脚本把  画这种图很简单，大圆可以用我前两天发帖的那个画椭圆的脚本，画不规则的，只要你知 ...

关键是要想求得所需经纬度需要对integ、 fid进行计算并利用其控制循环次数，为了找出问题，我在程序开始添加file1='f:/fog/model/att.txt'，在'set lon ' tlon     'set lat ' tlat之后加rc=write (file1,att' ',append)运行得到att.txt中只有att，没有具体数值；更进一步，我在'set lat ' tlat之后先加'define att=25*qcwsfc'
    'q define '
    att=sublin(result,10)
    att=subwrd (att,2)
再加rc=write运行的处的att.txt中，只有一个固定的值。
很显然程序运行中att没有随经纬度变化而变化，现在就是想法解决这个问题了，正在努力中
谢谢帮忙，请再帮忙调调

mofangbao

雾由心生 发表于 2013-5-12 20:04

                               
登录/注册后可看大图

关键是要想求得所需经纬度需要对integ、 fid进行计算并利用其控制循环次数，为了找出问题，我在程序开始添 ...

怎么和我说不到一路上去。你非得用这个脚本我就帮不了你了，我实在看不下这个脚本。
祝你成功~

雾由心生

mofangbao 发表于 2013-5-12 20:08

                               
登录/注册后可看大图

怎么和我说不到一路上去。你非得用这个脚本我就帮不了你了，我实在看不下这个脚本。
祝你成功~

还是得谢谢你，我再调调吧，实在调出再说吧，谢谢了