登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
数据文件为 文件1. 40 个站的 站点信息、经度、纬度 文件2 两列数据 1列为最高气温,1列为最低气温
Fortran, Ctl, gs 文件如下, grads中提示的错误信息如下图所示,不胜感激~
Fortran 程序 PARAMETER(num=40,ks=2) c num=stn,ks=nt it=2(1=eof;2=reof) DIMENSIONxlon(num),xlat(num),vis(num,ks) CHARACTERSTA(num)*8,FL*35,F5*5,F2*2,FL1*35 FL='c:\abc\40stn.txt' OPEN (11,FILE=FL,STATUS='OLD') Fl='win-ave.txt';FL1='win-ave.grd' OPEN (1,FILE=FL,status='old') open(3,file=fl1,form='binary') DO i=1,num READ(11,*) ist,xlat(i),xlon(i) WRITE(sta(i),'(I8.8)')ist read(1,*)(vis(i,j),j=1,ks) end do close(11) close(1) tim=0.0 nlev=1 nflag=1 do i=1,num write(3)sta(i),xlat(i),xlon(i),tim,nlev,nflag,(vis(i,k),k=1,ks) enddo nlev=0 write(3)sta(num),xlat(num),xlon(num),tim,nlev,nflag close(3) STOP 999 END Ctl 为 dset C:\abc\win-ave.grd dtype station stnmap C:\abc\zj.map undef 3276.6 *undef 1.0 title average tdef 1 linear 01Jan1979 1yr vars 2 V1 0 99 eigenvalue 1 V2 0 99 eigenvalue 2 endvars Gs 为 'REINIT' 'enable print C:\abc\win-ave.gmf' i=1 while(i<=2) 'set mpdset cnworld' 'set map 1' 'OPEN C:\abc\win-ave.ctl';a0='10 7';a1='1.57' 'open C:\abc\prc5.ctl' 'set parea 1. 10.5 1.0 8.' 'set lon 108 119';'set lat 20 27' 'set grid off';'set grads off' 'set dignum 0' 'SET DIGSIZ 0.1' 'set xlopts 1 5 0' 'set ylopts 1 5 0' 'set clab forced' *'set clevs -28.0 -24.0 -20.0 -16.0 -12.0-8.0 -4.0 4.0 8.0 12.0 16.0 20.0 24.0 28.0' 'set ccolor 1' 'set clopts 1 5 0.14';'set cthick 5';'setcint 1' 'Dsmth9(MASKOUT(OACRES(prc.2(T=1),v'i',20,16,12,8,4,2),prc.2(T=1)-0.5))' *'xylab2 2 1 0.16' if(i=1);fs='a';endif if(i=2);fs='b';endif 'set string 1 c 6' 'set strsiz 0.35';'draw string 1.7 7.5`1'fs 'print' 'reset' 'c' 'close 2' 'close 1' i=i+1 endwhile 'disable print'
|