NCL 画台风路径图问题

尽头的尽头

指甲钳 发表于 2013-11-14 18:59

                               
登录/注册后可看大图

怎么画啊？求指教！！

一种方法是在你的数据最后敲一个回车，让数据多一行，这样脚本就可以不改，图就出来了，还有就是不改数据，把行数直接用数字写到脚本，都可以，图就出来了

dataiyang

我使用相同的数据，相同的脚本，并没有提示任何错误，但是图中并没有任何台风信息！所以问题应该不是这里的吧！当然这个里面我觉得跟版本可能也有一些关系吧！

dataiyang

我使用相同的数据、相同的脚本，运行时并未出现任何报错，正常出图，但是图中并未绘上任何台风信息！
贴图：
虽然，没有解决问题，慢慢解析是什么原因！
我觉得跟ncl版本可能也有一定关系！
不知道有没有大神找到问题所在啊！

指甲钳

尽头的尽头 发表于 2013-11-14 19:37

                               
登录/注册后可看大图

一种方法是在你的数据最后敲一个回车，让数据多一行，这样脚本就可以不改，图就出来了，还有就是不改数据 ...

谢谢~~我明天试试！

指甲钳

dataiyang 发表于 2013-11-14 20:32

                               
登录/注册后可看大图

我使用相同的数据、相同的脚本，运行时并未出现任何报错，正常出图，但是图中并未绘上任何台风信息！
贴图 ...

我也遇到同样的问题了，不知道你解决了没有？

dataiyang

指甲钳 发表于 2013-11-15 09:15

                               
登录/注册后可看大图

我也遇到同样的问题了，不知道你解决了没有？

问题已解决，还是从脚本本身入手啊！
主要小小修改两三个地方
首先拷贝的别人的脚本这个默认是世界时，但是我们的txt文件中是北京时间，所以致使后面绘图时判断语句使我们需要的信息无法进入下一步；
另外一个就是txt文件读取的时候可能没有留意吧，经度读取的应该是第二例，原来的脚本中可能疏忽了，错误读取第一列，肯定是不对的。这里将图贴上

dataiyang

刚才的贴图啦！！！

指甲钳

dataiyang 发表于 2013-11-15 11:23

                               
登录/注册后可看大图

刚才的贴图啦！！！

我也发现我的经纬度读取出错了。但是那个时间应该怎么改呢？我是新手，也是照着脚本画的，那个世界时的怎么改？
do i = 1, nrow-1
    if (HH(i) .ne. "08" ) then
      continue
    end if
    dumDate(i) = gsn_add_text(wks,plot, MM(i)+DD(i), lon(i)+0.5, lat(i), resTx)
  end do

这一段的意思不是很懂啊

dataiyang

指甲钳 发表于 2013-11-15 11:59

                               
登录/注册后可看大图

我也发现我的经纬度读取出错了。但是那个时间应该怎么改呢？我是新手，也是照着脚本画的，那个世界时的怎 ...

俺要是新手啊，不断摸索中啊！
    if (HH(i) .eq. "08") then
      dumDot(i) = gsn_add_polymarker(wks, plot, lon(i), lat(i), resDot)
    end if
  do i = 1, nrow-1
    if (HH(i) .ne. "08" ) then
      continue
    end if
    dumDate(i) = gsn_add_text(wks,plot, MM(i)+DD(i), lon(i)+0.5, lat(i), resTx)
  end do

这两个地方"00"改为"08"
以下面为例：
do i = 1, nrow-1
    if (HH(i) .ne. "08" ) then
      continue
    end if
    dumDate(i) = gsn_add_text(wks,plot, MM(i)+DD(i), lon(i)+0.5, lat(i), resTx)
  end do
如果时间不是北京时08时时（这里是string类型），那么i 增加1，直到HH(i )=“08”为true时执行dumDate(i) = gsn_add_text(wks,plot, MM(i)+DD(i), lon(i)+0.5, lat(i), resTx)

指甲钳

dataiyang 发表于 2013-11-15 14:00

                               
登录/注册后可看大图

俺要是新手啊，不断摸索中啊！
    if (HH(i) .eq. "08") then
      dumDot(i) = gsn_add_polymarker( ...

这个我也改了，但是还是只出线，没有标出点和时间，你能帮我看看脚本吗？