GrADS：出现 logic errors 8 in gxshad 的原因

黎大页子

最近想对一块区域做时间的平均，用Fortran处理数据之后写成dat文件，放进GrADS里面画图，set  gxout contour时没有问题，可以出图，但是set gxout shaded总是出现‘logic errors 8 in gxshad’的提示错误，并且只能填色部分，不能正常显示。





后来发现是fortran处理数据时出现了问题，一些缺测的点在处理时被记录成了NaN（Fortran中对无效数据的记录），所以把这些点挑出来，改成32766就解决了。

可以另外写一个txt文件，看看是否数据中混有NaN，据说GrADS有命令可以查看数据中有没有NaN，但是我不会用，也没有搜到可用的正确语句。

小问题卡了一下午，可见严谨在编程中的重要。

这是之前出现错误的部分

1. do x=1,xnum
   
2. do y=1,ynum
   
3.   do i=1,12 !From Jan to Dec
   
4.     count=0
   
5.     do t=i,tnum-12,12
   
6.           if(yy(x,y,t)/=32766)then
   
7.           a(x,y,i)=a(x,y,i)+yy(x,y,t)
   
8.           count=count+1
   
9.           end if
   
10.         end do
    
11.         a(x,y,i)=a(x,y,i)/count !当全部缺测的站点，此处的count=0，所以a（x,y,i）会被记录成NaN，而不是32766
    
12.         end if
    
13.   end do
    
14. end do
    
15. end do
    

do x=1,xnum<br /> do y=1,ynum<br />   do i=1,12 !From Jan to Dec<br />     count=0<br />     do t=i,tnum-12,12<br />           if(yy(x,y,t)/=32766)then<br />           a(x,y,i)=a(x,y,i)+yy(x,y,t)<br />           count=count+1 <br />           end if<br />         end do<br />         a(x,y,i)=a(x,y,i)/count !当全部缺测的站点，此处的count=0，所以a（x,y,i）会被记录成NaN，而不是32766<br />         end if<br />   end do<br /> end do<br /> end do<br />

修改后的部分

1. do x=1,xnum
   
2. do y=1,ynum
   
3.   do i=1,12 !From Jan to Dec
   
4.     count=0
   
5.     do t=i,tnum-12,12
   
6.           if(yy(x,y,t)/=32766)then
   
7.           a(x,y,i)=a(x,y,i)+yy(x,y,t)
   
8.           count=count+1
   
9.           end if
   
10.         end do
    
11.         if(count==0)then  !加上一个判断，避免出现NaN的记录
    
12.     a(x,y,i)=32766
    
13.     else
    
14.     a(x,y,i)=a(x,y,i)/count
    
15.         end if
    
16.   end do
    
17. end do
    
18. end do
    

do x=1,xnum<br /> do y=1,ynum<br />   do i=1,12 !From Jan to Dec<br />     count=0<br />     do t=i,tnum-12,12<br />           if(yy(x,y,t)/=32766)then<br />           a(x,y,i)=a(x,y,i)+yy(x,y,t)<br />           count=count+1 <br />           end if<br />         end do<br />         if(count==0)then  !加上一个判断，避免出现NaN的记录<br />     a(x,y,i)=32766<br />     else<br />     a(x,y,i)=a(x,y,i)/count<br />         end if<br />   end do<br /> end do<br /> end do<br />

一个不小心

我没用FORTRAN，直接用GRADS处理NC资料也遇到过这种情况。难道是GRADS不能自动识别缺测值？

mofangbao

谢谢页子提醒，我好像也遇到过...

黎大页子

mofangbao 发表于 2012-5-12 20:58

                               
登录/注册后可看大图

谢谢页子提醒，我好像也遇到过...

嘿嘿，不客气，为什么我自己发的帖的代码部分不会像你编辑过的帖子一样呢？

mofangbao

黎大页子 发表于 2012-5-12 21:37

                               
登录/注册后可看大图

嘿嘿，不客气，为什么我自己发的帖的代码部分不会像你编辑过的帖子一样呢？

在你的代码两头加上
[code=fortran]
你的代码
[/fortran]
我先给你加上，你点一下编辑帖子就能看到样子了

Really

果然解决了，谢谢~~~

xiexieixiexjie

1217342256

好贴，拜谢

逸儿~~

好贴~

小小一只杨

我想请问楼主，如果是wrfoutput的数据，用GRADS画图时出现这个问题要到哪修改？