rain_grads.txt

气象单要素资料:如雨量气压 温度、温度露点差、 相对湿度 等转为grads数据

!  ***********************************************
!    单站要素转换为GRADS数据格式
!    重庆：王中   2007年4月11日
!    [e-mail] 
!  ***********************************************
!    数据说明：
!    StaN  站点数
!	 T_M   数据列数
!	 Sta   站号
!	 Lat   纬度  Lon  经度
!	 DZ_Data   要素
!  ***********************************************
     program DZData
    integer,parameter ::StaN=9,T_M=1
     character*8 Sta(StaN)
	 real Lat(StaN)
	 real Lon(StaN)
	 real DZ_Data(StaN,T_M)
!     character*9 flnm(32)

!------------------------------------------------------------------------
      do jj=1,StaN
      do ii=1,T_M
         DZ_Data(StaN,T_M)=-9999.9
      enddo
      enddo
!    打开原始数据文件
      open(10,file='d:\资料处理\DZ_Data1.dat', form='formatted')
!    读取数据
      do kk=1,StaN
       read(10,'(a8,f7.2,2f6.2)') sta(kk),lon(kk),lat(kk),DZ_Data(KK,1)
	   ! write(*,*)sta(kk),lon(kk),lat(kk),DZ_Data(KK,1)
	   ! write(*,*)DZ_Data(KK,1)  ,DZ_Data(KK,3)
	   
	  enddo
	  close(10)
	  write(*,*)" ----- 读取数据----> OK !"
!    写入数据文件中
      open(11,file='d:\资料处理\DZ_Data0.dat', form='binary')
	  open(12,file='d:\资料处理\DZ_Data0.txt')

	   TIM=0.0
	   NLEV=1
	   NFLAG=1
	   DO JJ=1,T_M
	   DO KK=1,StaN
	     WRITE(11) sta(kk),lon(kk),lat(kk),TIM,NLEV,NFLAG
	     WRITE(11) DZ_Data(KK,JJ)
	
	     WRITE(12,*) sta(kk),lon(kk),lat(kk),TIM,NLEV,NFLAG
	     WRITE(12,*) DZ_Data(KK,JJ)

	   ENDDO
	   ENDDO
	   NLEV=0
	   sta(kk)="    57000"
	   WRITE(11) sta(kk),lon(kk),lat(kk),TIM,NLEV,NFLAG
	   WRITE(12,*) sta(kk),lon(kk),lat(kk),TIM,NLEV,NFLAG
	   WRITE(*,*) sta(kk),lon(kk),lat(kk),TIM,NLEV,NFLAG

	   CLOSE(11)
	   CLOSE(12)
   	  write(*,*)" ----- 写入数据----> OK !"

!------------------------------------------------------------------
          end