登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
 - program the skill
- real obs(160,17),model(160,17),Ps(17)
- integer n1(17),n2(17),n0(17)
- parameter (n=160,f1=5,f2=2)
- open(10,file='d:\ssx\two\160data.grb',form='binary')
- do it=1,17
- do j=1,160
- read(10)obs(j,it)
- enddo
- do j=1,160
- read(10)model(j,it)
- enddo
- enddo
- do it=1,17
- do j=1,160
- if((abs(obs(j,it))>=50).and.(abs(model(j,it))>=50)) then
- n1(it)=n1(it)+1
- endif
- if((abs(obs(j,it))>=20).and.(abs(model(j,it))>=20)) then
- if((abs(obs(j,it))<=50).and.(abs(model(j,it))<=50)) then
- n2(it)=n2(it)+1
- endif
- endif
- if(obs(j,it)*model(j,it)>0) then
- n0(it)=n0(it)+1
- else if(obs(j,it)*model(j,it)<0) then
- if(abs(obs(j,it))<=15.and.abs(model(j,it))<=15) then
- n0(it)=n0(it)+1
- endif
- endif
- enddo;enddo
- open(11,file='d:\ssx\ps.txt')
- do it=1,17
- Ps(it)=100*(n0(it)+n1(it)*f1+n2(it)*f2)/(n+n1(it)*f1+n2(it)*f2)
- write(11,*)'年份 预测评分Ps'
- write(11,"(i4,f8.2)") it+1983,Ps(it)
- enddo
- end program
|