天诊实习报告一.docVIP

实习报告一

实习要求

利用所提供的观测资料、模拟资料和相应的资料读取程序，编写气候极值检查的程序，并对观测资料进行质量控制试验，并根据各个站点剔除资料量画出台站资料质量分布特征图。

二．实习内容

Read.f:

Programmain

implicitnone

integerim,km,i,k

parameter(im=3197,km=239)

realobs(im,km)

reallat(im),lon(im)

open(1,file=f:\zhiliang\observation.grd,form=binary)

open(2,file=f:\zhiliang\lat_all.txt)

open(3,file=f:\zhiliang\lon_all.txt)

open(4,file=f:\zhiliang\tcd.txt)

dok=1,km

doi=1,im

read(1)obs(i,k)

enddo

enddo

doi=1,im

read(2,*)lat(i)

read(3,*)lon(i)

enddo

calltichu(obs,lon,lat)

close(1)

close(2)

close(3)

close(4)

end

subroutinetichu(a,b,c)

implicitnone

integerim,km,i,k

parameter(im=3197,km=239)

integermiss(im)

reala(im,km)

realb(im),c(im)

doi=1,im

miss(i)=0

dok=1,km

if(a(i,k)(-96.0))then

miss(i)=miss(i)+1

elseif((c(i)45.0).and.(a(i,k)50.0))then

miss(i)=miss(i)+1

elseif((c(i)=45.0).and.(a(i,k)60.0))then

miss(i)=miss(i)+1

endif

enddo

write(4,*)c(i),b(i),miss(i)

enddo

end

Parameter(im=3197,km=1)

realwint(im,km)

reallat(im,km),lon(im,km),owint(im,km)

realolat(im),olon(im),a,b,c

character*8zh(im)

integeri_qzk

STATION.for:

open(9,file=f:\zhilaing\rs160.grd,form=binary)

open(1,file=f:\zhiliang\tcd.txt)

dok=1,km

doisk=1,im

read(1,*,end=1000)a,b,c

olat(isk)=a

olon(isk)=b

owint(isk,k)=c

enddo

enddo

close(1)

1000continue

doi=1,im

zh(i)=char(i)

enddo

do101k=1,km

ii=0

doi=1,im

print*,i

if(owint(i,k).le.90000)then

ii=ii+1

lat(ii,k)=olat(i)

lon(ii,k)=olon(i)

wint(ii,k)=owint(i,k)

endif

enddo

print*,ii,k

72format(f5.2,f7.2)

tim=0.0

do100i=1,im

nlev=1

nflag=1

write(9)zh(i),lat(i,k),lon(i,k),tim,nlev,nflag

write(9)wint(i,k)

100continue

nlev=0

write(9)zh(i-1),lat(i-1,1),lon