twin cycloneをプロット

twin cycloneをプロット

~/wkyoshimura2_3/draw4.gs
draw4.gs

*tmp1=0
*tmp=read(data.dat)
*while(tmp1=0)
  tmp=read(data.dat)
  tmp1=sublin(tmp,1)
  tmp2=sublin(tmp,2)
  nlon=subwrd(tmp2,1)
  nlat=subwrd(tmp2,2)
  slon=subwrd(tmp2,3)
  slat=subwrd(tmp2,4)
  say tmp2
  sym = 3
  siz = 0.15
  'q w2xy 'nlon' 'nlat''
  x1 = subwrd(result,3)
  y1 = subwrd(result,6)
  'draw mark 'sym' 'x1' 'y1' 'siz
  'q w2xy 'slon' 'slat''
  x2 = subwrd(result,3)
  y2 = subwrd(result,6)
  'draw mark 'sym' 'x2' 'y2' 'siz
  'draw line 'x1' 'y1' 'x2' 'y2''