時間內插取代初始小時濃度

kuang@DEVP /nas1/Data/javascripts/D3js/earthFcst03/public/data/weather/current
$ cat intp_json.py
#!/opt/anaconda3/envs/py37/bin/python
import numpy as np
import json
import sys, os
import subprocess
import datetime

#the json template
grds=subprocess.check_output('pwd',shell=True).decode('utf8').strip('\n').split('/')[5].replace('earthFcst','')
pwd='/nas1/Data/javascripts/D3js/earthFcst'+grds+'/public/data/weather/current/'
bdate=datetime.datetime.strptime(sys.argv[1],"%Y-%m-%d")
dirs=[(bdate+datetime.timedelta(days=i)).strftime("../%Y/%m/%d/") for i in [-1,0,0]]
hhs =[(bdate+datetime.timedelta(hours=i)).strftime("%H00") for i in [-3,3,0]]

for sp in ['no2','ozone','ozone8','pm1','pm10','pm25','so2','vocs']:
  arr=[]
  skip=0
  for id in [0,1]:
    fname=pwd+dirs[id]+hhs[id]+'-'+sp+'-surface-level-fcst-'+grds+'.json'
    if not os.path.isfile(fname):
      print('file '+fname+' not found!')
      skip=1
      break
    with open(fname,'r') as f:
      gfs=json.load(f)
    arr.append(np.array(gfs[0]['data']))
  if skip==1:continue
  fnameO=pwd+dirs[2]+hhs[2]+'-'+sp+'-surface-level-fcst-'+grds+'.json'
  with open(fnameO,'r') as f:
    gfs=json.load(f)
  gfs[0]['data']=list((arr[0]+arr[1])/2.)
  with open(fnameO,'w') as f:
    json.dump(gfs,f)