caline3遠端計算服務

fname=sys.argv[1]
NS,NR,nms,hgts,lon,lat,lonp,latp,TYPs,VPHs,EMFs,WIDs=rd_kmlLL(fname)

NM=8
parm={'U':1.0,'CLAS':6,'MIXH':100,'AMB':1.0}
for var in 'U CLAS MIXH AMB'.split():
  exec(var+'=[parm["'+var+'"] for i in range(NM)]')
BRG=[45*i for i in range(NM)]

#receptors
recp,zr=nms[NS:],hgts[NS:]
xr,yr=pnyc(lonp, latp, inverse=False)
xr=[int(i+Xcent) for i in xr]
yr=[int(i+Ycent) for i in yr]

#links
NL,(lnks,X1,Y1,X2,Y2,TYP,VPH,EMF,H,W)=0,([] for i in range(10))
for l in range(NS):
  llt=[i for i in lat[l,:] if i !=0]
  lln=[i for i in lon[l,:] if i !=0]
  x,y=pnyc(lln, llt, inverse=False)
  x=[int(i+Xcent) for i in x]
  y=[int(i+Ycent) for i in y]
  X1+=x[:-1]
  Y1+=y[:-1]
  X2+=x[1:]
  Y2+=y[1:]
  ms=len(llt)-1
  lnks+=[nms[l]+'_'+str(i) for i in range(ms)]
  TYP+=[TYPs[l] for i in range(ms)]
  VPH+=[VPHs[l] for i in range(ms)]
  EMF+=[EMFs[l] for i in range(ms)]
  H+=[hgts[l] for i in range(ms)]
  W+=[WIDs[l] for i in range(ms)]
  NL+=ms