|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
# W% T4 }/ E+ O- clear,clc
- nt=735;
- ny=73; %2.5*2.5格点的
- nx=144; %2.5*2.5格点的
- f=netcdf('air.mon.mean.nc','nowrite');
- tt=f{'air'}(:);
- close(f)
- fid=fopen('a5.grd','w');
- for it=1:nt
- for j=1:ny
- for i=1:nx
- count= fwrite(fid,tt(it,j,i),'float');%根据例子这里用原始的数据
- end
- end
- end
- fclose(fid)
- a5.grd 对应的ctl文件是
- dset F:a5.grd
- title Monthly mean air temperature NCEP Reanalysis
- options sequential yrev
- undef -9.96921e+36
- xdef 145 linear 0 2.5
- ydef 73 linear -90 2.5
- zdef 1 linear 0 0
- tdef 735 linear 00Z01JAN1948 1mo
- vars 1
- a 0 -999 Monthly Mean Air Temperature
- endvars% n8 @, K" o% _, H: \$ Y" `
1 |. X e9 r0 A) l
% n6 m$ F( z7 i/ l
% k Q- _" w" L3 z F x6 K1 B1 e+ ^: Q: @: D
|
|