天診實(shí)習(xí)報(bào)告一

實(shí)習(xí)報(bào)告一實(shí)習(xí)要求利用所提供的觀測(cè)資料、模擬資料和相應(yīng)的資料讀取程序，編寫氣候極值檢查的程序，并對(duì)觀測(cè)資料進(jìn)行質(zhì)量控制試驗(yàn)，并根據(jù)各個(gè)站點(diǎn)剔除資料量畫出臺(tái)站資料質(zhì)量分布特征圖。二．實(shí)習(xí)內(nèi)容Read.f:Programmain implicitnoneintegerim,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,kmdoi=1,imread(1)obs(i,k) enddo enddo doi=1,imread(2,*)lat(i) read(3,*)lon(i)enddocalltichu(obs,lon,lat)close(1) close(2)close(3)close(4) endsubroutinetichu(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 endParameter(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(i-1,1),tim,nlev,nflag101continue close(9) end gs：'reinit''!stnmap-irs160.ctl''reinit'*****blue'setrgb16470255''setrgb512428240''setrgb533642248''setrgb554856255''setrgb576070255''setrgb597284255''setrgb618498255''setrgb6396112255''setrgb65108126255''setrgb67120140255''setrgb69132154255'*****red'setrgb48255250205''setrgb50255244165''setrgb52255238125''setrgb54255214101''setrgb5625519077''setrgb5825516653''setrgb6025514229''setrgb622551185''setrgb64255940''setrgb66255700''setrgb68255460''setrgb9525500'********************************************'OPENf:/zhiliang/grid.ctl''OPENf:/zhiliang/rs160.ctl'*************************k=1while(k<=1)'setparea1101.57.5''sett'k'qtime'res=subwrd(result,3)'setlon50150''setlat070''setmpdsethires''setmap116''DEFINEA1=OACRES(g,RS.2)''DEFINEb1=MASKOUT(A1,g-0.5)''DEFINEBB1=(b1)''SETGRIDOFF';'SETGRADSOFF';'setframeon''setxlopts1100.24''setylopts1100.24''setxlint30''setylint10''setgxoutshaded''setclevs-8-6-4-202468'*'setccols16535557596367545660646895''d(smth9(bb1))''cbarn1.005.52.5''gxyatsta.png-x1200-y1200'pulldummy'c'k=k+1endwhile'reinit'; 臺(tái)站資料質(zhì)量分布特征圖：結(jié)果分析根據(jù)上圖可做出如下分析：1.總體上看，除西部部分地區(qū)外，我國(guó)大部分地區(qū)資料質(zhì)量較高，剔除都相對(duì)較少，這和地形以及設(shè)備的先進(jìn)有很大的關(guān)系，。2.中國(guó)西部有部分顯示無色區(qū)域，這個(gè)區(qū)域位于我國(guó)青海西藏地區(qū)，在這個(gè)地區(qū)，站點(diǎn)設(shè)的相對(duì)較少，觀測(cè)條件和設(shè)備較落后，人員也少，所以剔除的資料最多。3.在西藏南部，四川地區(qū)被剔除的資料也較多，可能是由于該地區(qū)的地形