提取单个省份数据后画图无法显示，求各位帮忙看下~（已解决）

sunshinezhpb

我有一个全国的0.25格点的dat数据（中国区域外全是缺测），我只想要西南地区的，于是就把这个dat转成txt文件在arcgis里裁剪出西南几个省份的数据，把其余部分全部设为了缺测，这样得出一个txt文件再转成dat文件。但是我用处理后的画图，有显示变量值的范围，但在图形上完全不显示，请问这是怎么回事呢？
情况如下图：

dat转txt的程序：
parameter(m=283,n=163,t=20)
dimension  a1(m,n,t)
real lon(m),lat(n)
integer i,j

open(11,file='H:\ccsm4\t-1mon.dat',form='binary')
read(11) (((a1(i,j,k),j=1,n),i=1,m),k=1,t)
close(11)

do i = 1,283
   lon(i) = 69.75 + (i-1)*0.25
enddo
print *,lon(280)
do i = 1,163
   lat(i) = 14.75 + (i-1)*0.25
enddo

open(13,file='C:\Users\Administrator\Desktop\xishu\area\arcgis.txt')
do i=1,m
  do j=1,n
      write(13,*) a1(i,j,1),lon(i),lat(j)
  enddo
enddo
end

我的ctl：
dset C:\Users\Administrator\Desktop\xishu\area\xinan1.dat                                                        
title CN05.1 0.25x0.25 grided daily obs data provided by 2400 weather stations                     
undef -9999.                                                                                       
options little_endian                                                                              
XDEF  283 LINEAR 69.75 0.25                                                                        
YDEF  163 LINEAR 14.75 0.25                                                                        
ZDEF 1 levels 1000                                                                                 
tdef 1 linear  01jan1986 1yr                                                                  
vars 1                                                                                               
tm 0 99 CN05.1 month, mean temperature,c                                                        
endvars   

sunshinezhpb

已解决，fortran程序问题，grads最内层循环应该是x，这点老是搞混。

月影青枫

请问下你这个cn05.1数据是在哪儿下载的呀，可以方便告知下么？

handangjh

月影青枫 发表于 2016-10-14 17:44
请问下你这个cn05.1数据是在哪儿下载的呀，可以方便告知下么？

您好。我有同样的问题，想用一下CN05.1数据集，但不知道在哪儿下载

白夜桃夭

同求CN05.1数据集