Micaps格式站点数据转换为Grads格式数据.docVIP

将Micaps格式站点数据转换为Grads格式数据 !------------------------------------------------------------------------ ! 本程序将Micaps格式站点数据转换为Grads格式数据 ! ----------------------------------------------------------------------- program GKSTA parameter (inum=1000,kmax=9) ! inum 站点数 kmax 层次数 dimension nyear(kmax),nmon(kmax),nday(kmax),ntim(kmax) dimension lev(kmax),nsta(kmax) dimension slon(inum,kmax),slat(inum,kmax),h(inum,kmax) dimension t(inum,kmax), td(inum,kmax),orit(inum,kmax) dimension rr(inum,kmax), uu(inum,kmax),vv(inum,kmax) dimension state(inum,kmax) dimension pp(kmax) character*8 stid character*80 file1, file2,file3,file4,file5,file6,file7 character*80 file8, file9,file00 !file10,file00 character*8 yymmddhh character*2 mon(12),da(31) character*4 RIQI real YYYY integer mon0,cc,ys,ll0,ll1 character*2 tim0 logical alive character*50 FilePath character*11 FilePath0 character*4 cc0(9) data cc0/1000,925 ,850 ,700 ,500 ,400 ,300 ,200 ,100/ data pp/1000.,925.,850.,700.,500.,400.,300.,200.,100./ data mon/01,02,03,04,05,06,07,08,09,10,11,12/ data da/01,02,03,04,05,06,07,08,09,10, 11,12,13,14,15,16,17,18,19,20, 21,22,23,24,25,26,27,28,29,30,31/ !------------------------------------------------------------------------ FilePath0=d:\micaps write(*,*) 将Micaps格式 站点 数据转换为Grads格式数据 开始。。。 OPEN (3,FILE=D:\HIS_DATA\RIQI.DAT,status=old) READ(3,*) RIQI CLOSE(3) OPEN (3,FILE=D:\HIS_DATA\riqi.DAT) read(3,*) yyyy CLOSE(3) open(2,file=D:\HIS_DATA\DATA\Up//RIQI//.sta,form=binary)!存数据的文件 do mon0=1,12 !month begain DDM=31 if(mon0.eq.4)then DDM=30 ENDIF if(mon0.eq.6)then DDM=30 ENDIF if(mon0.eq.9)then DDM=30 ENDIF if(mon0.eq.11)then DDM=30 ENDIF if(mon0.eq.2)then DDM=28 if(YYYY/4.eq.INT(YYYY/4))then DDM=29 ENDIF ENDIF do ddd=1,DDM !day begain do ti=1,2 !ti