ps评分程序

无敌三脚猫

program ps
implicit none
integer,parameter::f1=2,f2=1,m=17,n=3*m     !!旬为5和2，m为总站数
integer i,j
real for(n),shi(n),ave(m)
real a,b,c,d,aa,da,db,bb,pps,bridge
open(1,file='F:\ch\emd\m\20088\1.1.txt')     !!!for为各站预报值
read(1,'()')
read(1,*)(for(i),i=1,n)
do i=1,m
print*,"ok"
bridge=for(3*i)
for(i)=bridge
end do
close(1)
open(2,file='F:\ch\combine\result\20088.txt')
read(2,'()')          !!!shi为各站实况值
do i=1,m
bridge=shi(3*i)
shi(i)=bridge
end do   
read(2,*)(shi(i),i=1,n)
close(2)
open(3,file='F:\ch\combine\result\a08.txt')     !!!ave为各站多年平均值
read(3,*)(ave(i),i=1,m)
close(3)   
aa=0
da=0
db=0
do i=1,m
  a=for(i)-ave(i)
  b=shi(i)-ave(i)              
  c=a*b
  da=abs(a/ave(i))
  db=abs(b/ave(i))
  if(c.GT.0)then
  aa=aa+1
  else if((0.15>da).AND.(0.15>db))then
  aa=aa+1
  else
  aa=aa+0
  end if
end do
bb=0
do i=1,m
  a=for(i)-ave(i)
  b=shi(i)-ave(i)              
  c=a*b
  da=abs(a/ave(i))
  db=abs(b/ave(i))
  if ((c>0).AND.(((da>0.2).AND.(da<0.5)).AND.((db>0.2).AND.(db<0.5))))then
  bb=bb+1
  else if((c>0).AND.((da>0.5).AND.(db>0.5)))then
  bb=bb+2
  else
  bb=bb+0
  endif
end do
pps=(aa+bb)/(m+bb)
print*,pps
end

随缘

哎呀！无偿献血啊！风格高尚的很嘛！

deeli

这个是那个降水的评分那个吗？

无敌三脚猫

deeli 发表于 2013-4-1 21:01

                               
登录/注册后可看大图

这个是那个降水的评分那个吗？

是的。气温的分成7级。

筑梦成真

这个对么？f1,f2系数都没有用到啊。。。

无敌三脚猫

筑梦成真 发表于 2013-5-4 16:26

                               
登录/注册后可看大图

这个对么？f1,f2系数都没有用到啊。。。

F1，F2就是bb=bb+1，bb=bb+2了。不好意思。写的太乱。

aa22556688

楼主无私，灰常给力。。。{:eb515:}{:eb515:}

故乡的云

新的评分标准不是变了吗？

无敌三脚猫

故乡的云 发表于 2014-7-13 12:34
新的评分标准不是变了吗？

可以自己调整吧？

故乡的云