grads输出dat数据用matlab读取，矩阵边界为999的问题

chuyingjia

数据是ncepde 2.5*2.5再分析数据，显示计算整层水汽通量，输出.dat数据，然后用matlab读取进行后续计算时发现读出的【23,13,840】数据时间维度是对的，但每个时间维度的第一行和最后一行，第一列和最后一列都是-99900000这样。不知道是哪里的问题，是默认把首尾两个格点数据去掉了吗？按说不应该啊。下面是我的gs、matlab以及读取数据截图'reinit''sdfopen uwnd.mon.mean.nc'
'sdfopen vwnd.mon.mean.nc'
'sdfopen shum.mon.mean.nc'
'sdfopen pres.mon.mean.nc'
*'gxprint shuiqi.png white'
*************Figure 1*********
'set vpage 0 11 0 8.5'
'set parea 0.4 5.5  4.6 8.5'
'set grads off''set grid off'
'set map 1 1 5''set gxout fwrite'
'set fwritezcsandu.dat'
*'set gxout shaded'
'set lat 20 50''set lon 90 145'
'set t 1 840''define uq=shum.3*uwnd.1/9.8'
'define vq=shum.3*vwnd.2/9.8'
'define sandu=hdivg(uq,vq)'
'define zcsandu=vint(pres.4(lev=0),sandu,300)'
'd zcsandu'
'disable fwrite'
'reinit'
+++++++++++++++++++++++++++++++matlab程序+++++++++++++
clear all;
clc;
mx=23;%line
my=13; %Columns
mt=840; %year
fid1=fopen('J:\zcsandu.dat','rb');

sand=fread(fid1,'float32');
sandu=reshape(sand,mx,my,mt);
sandu1=sandu(2:mx-1,2:my-1,:);
fclose(fid1);
++++++++++++++++++++++++++++数据截图++++++++++

chuyingjia

grads刚出了一下图，计算出来边界就是空的。

chuyingjia

经大神指点，找到原因了，在原帖贴上示意图，希望能帮到有同样疑惑的小伙伴

伽蓝鸟

嗯，grads是用中央差的方法算涡度散度的，所以边界上的点没法定义。你要是需要的话，可以把边界扩大一点。

张一鸣

帮助到啦，感谢 {:5_235:} {:5_235:}