1111111

 VALUE=new((/1,8,53/),string) 
   file_obs="china_sites_20200630.csv"
  ;file_obs="china_sites_20200901.csv"
 row=361
 col=1673
 readvalue=new((/row,col/),string)
 f_h = asciiread("sitelist.csv",-1,"string")
 delim=","
 site=str_get_field(f_h,1,delim)
 allvalue=asciiread(file_obs,-1,"string")
 allvalue@_Fillvalue="-999"
  do u=1,col
    readvalue(0:360,u-1)=str_get_field(allvalue,u,delim)
    readvalue@_Fillvalue="-999"
  end do
 Col=0
 Row=0
   do j=241,346,15
  ;do j=1,107,15
      do k=3,col-1 
           if (any(site .eq. readvalue(0,k)  )) then
             VALUE(0,Row,Col)=readvalue(j,k)
             Col=Col+1
            end if
       end do
       Row=Row+1
       Col=0
  end do
VALUE!0="day"
VALUE!1="hour"
VALUE!2="site"
reValue=VALUE(site|:,day|:,hour|:)
;write_table("AQI8_27.csv","a",[/reValue(0,0,:),reValue(1,0,:),reValue(2,0,:),reValue(3,0,:),reValue(4,0,:),reValue(5,0,:),reValue(6,0,:),reValue(7,0,:),reValue(8,0,:),reValue(9,0,:),reValue(10,0,:),reValue(11,0,:),reValue(12,0,:),reValue(13,0,:),reValue(14,0,:),reValue(15,0,:),reValue(16,0,:),reValue(17,0,:),reValue(18,0,:),reValue(19,0,:),reValue(20,0,:),reValue(21,0,:),reValue(22,0,:),reValue(23,0,:),reValue(24,0,:),reValue(25,0,:),reValue(26,0,:)/],"%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s")
;write_table("AQI8_26.csv","a",[/reValue(27,0,:),reValue(28,0,:),reValue(29,0,:),reValue(30,0,:),reValue(31,0,:),reValue(32,0,:),reValue(33,0,:),reValue(34,0,:),reValue(35,0,:),reValue(36,0,:),reValue(37,0,:),reValue(38,0,:),reValue(39,0,:),reValue(40,0,:),reValue(41,0,:),reValue(42,0,:),reValue(43,0,:),reValue(44,0,:),reValue(45,0,:),reValue(46,0,:),reValue(47,0,:),reValue(48,0,:),reValue(49,0,:),reValue(50,0,:),reValue(51,0,:),reValue(52,0,:)/],"%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s")
write_table("AQI8_27.csv","w",[/reValue(0,0,16:23),reValue(1,0,16:23),reValue(2,0,16:23),reValue(3,0,16:23),reValue(4,0,16:23),reValue(5,0,16:23),reValue(6,0,16:23),reValue(7,0,16:23),reValue(8,0,16:23),reValue(9,0,16:23),reValue(10,0,16:23),reValue(11,0,16:23),reValue(12,0,16:23),reValue(13,0,16:23),reValue(14,0,16:23),reValue(15,0,16:23),reValue(16,0,16:23),reValue(17,0,16:23),reValue(18,0,16:23),reValue(19,0,16:23),reValue(20,0,16:23),reValue(21,0,16:23),reValue(22,0,16:23),reValue(23,0,16:23),reValue(24,0,16:23),reValue(25,0,16:23),reValue(26,0,16:23)/],"%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s")
write_table("AQI8_26.csv","w",[/reValue(27,0,16:23),reValue(28,0,16:23),reValue(29,0,16:23),reValue(30,0,16:23),reValue(31,0,16:23),reValue(32,0,16:23),reValue(33,0,16:23),reValue(34,0,16:23),reValue(35,0,16:23),reValue(36,0,16:23),reValue(37,0,16:23),reValue(38,0,16:23),reValue(39,0,16:23),reValue(40,0,16:23),reValue(41,0,16:23),reValue(42,0,16:23),reValue(43,0,16:23),reValue(44,0,16:23),reValue(45,0,16:23),reValue(46,0,16:23),reValue(47,0,16:23),reValue(48,0,16:23),reValue(49,0,16:23),reValue(50,0,16:23),reValue(51,0,16:23),reValue(52,0,16:23)/],"%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s")
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值