ruby-dcl 1.6.2

data/demo/grph2/uspack/uspk03.rb

@@ -0,0 +1,38 @@
+#
+# uspk03.rb
+#
+# $Id: uspk03.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 30
+y = NArray.sfloat(n)
+
+#---------------------------- DATA DEFINITION --------------------------
+r = 0.2
+a = 3.7
+for i in 0..n-1
+  r = a*r*(1.0-r)
+  y[i] = r
+end
+
+#----------------------------- GRAPH -----------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+rundef = DCL::glrget('RUNDEF')
+DCL::gropn iws
+
+DCL::grfrm
+DCL::grswnd(0.0, 1.0, rundef, rundef)
+
+DCL::usgrph(rundef, y)
+
+DCL::grcls
+
+#-----------------------------------------------------------------------
+

data/demo/grph2/uspack/uspk04.rb

@@ -0,0 +1,64 @@
+#
+# uspk04.rb
+#
+# $Id: uspk04.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 200
+m = 5
+x = NArray.sfloat(n)
+y0 = NArray.sfloat(n)
+y1 = NArray.sfloat(n)
+y2 = NArray.sfloat(n)
+a = NArray.sfloat(m)
+
+#-----------------------------------------------------------------------
+dt = 1.0/(n-1)
+for j in 0..m-1
+  jj = j*2+1.0
+  a[j] = (-1.0)**(j+1.0) *2.0/(jj*PI)
+end
+
+x = x.indgen * dt
+
+for i in 0..n-1
+  t = dt*i*2*PI
+  y2[i] = 0.0
+  for j in 0..m-1
+    jj = j*2+1
+    yy = a[j]*cos(jj*t)
+    y2[i] = y2[i] + yy
+  end
+  y1[i] = a[0]*cos(t)
+
+  if (t < PI/2.0 || t >= PI*3.0/2.0)
+    y0[i] = -0.5
+  else
+    y0[i] = 0.5
+  end
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::usspnt(x, y1)
+DCL::usspnt(x, y2)
+
+DCL::uuslni(5)
+DCL::usgrph(x, y0)
+
+DCL::uulinz(x, y1, 3, 1)
+DCL::uulinz(x, y2, 2, 2)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk05.rb

@@ -0,0 +1,68 @@
+#
+# uspk05.rb
+#
+# $Id: uspk05.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 200
+m = 5
+x = NArray.sfloat(n)
+y0 = NArray.sfloat(n)
+y1 = NArray.sfloat(n)
+y2 = NArray.sfloat(n)
+a = NArray.sfloat(m)
+
+#-----------------------------------------------------------------------
+dt = 1.0/(n-1)
+for j in 0..m-1
+  jj = j*2+1
+  a[j] = (-1.0)**(j+1) *2.0/(jj*PI)
+end
+
+x = x.indgen * dt
+for i in 0..n-1
+  t = dt*i*2*PI
+  y2[i] = 0.0
+  for j in 0..m-1
+    jj = j*2+1
+    yy = a[j]*cos(jj*t)
+    y2[i] = y2[i] + yy
+  end
+  y1[i] = a[0]*cos(t)
+
+  if (t < PI/2.0 || t >= PI*3.0/2.0)
+    y0[i] = -0.5
+  else
+    y0[i] = 0.5
+  end
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::usspnt(x, y0)
+DCL::usspnt(x, y1)
+DCL::usspnt(x, y2)
+
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uulinz(x, y0, 1, 5)
+DCL::uulinz(x, y1, 3, 1)
+DCL::uulinz(x, y2, 2, 2)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk06.rb

@@ -0,0 +1,43 @@
+#
+# uspk06.rb
+#
+# $Id: uspk06.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 200
+t = NArray.sfloat(n)
+a = NArray.sfloat(n)
+x, y, z, s, r, b, dt = [0.0, 1.0, 1.0, 10.0, 26.0, 2.6, 0.01]
+
+#-----------------------------------------------------------------------
+t = t.indgen(0.0, 1000)
+for i in 0..n-1
+  for j in 0..7
+    dx = -s*x + s*y
+    dy = -x*z + r*x - y
+    dz =  x*y - b*z
+    x = x + dx*dt
+    y = y + dy*dt
+    z = z + dz*dt
+  end
+  a[i] = y + 20.0
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::ussttl('TIME', 'SEC', 'HEAT FLUX', 'W/m|2"')
+DCL::usgrph(t, a)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk07.rb

@@ -0,0 +1,62 @@
+#
+# uspk07.rb
+#
+# $Id: uspk07.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 200
+m = 10
+x = NArray.sfloat(n)
+y = NArray.sfloat(n)
+a = NArray.sfloat(m)
+
+#-----------------------------------------------------------------------
+dt = 1.0/(n-1)
+for j in 0..m-1
+  jj = j*2+1
+  a[j] = (-1.0)**(j+1) *2.0/(jj*PI)
+end
+
+x = x.indgen * dt
+for i in 0..n-1
+  t = dt*i*2*PI
+  yy = 0.0
+  for j in 0..m-1
+    jj = j*2+1
+    yd = a[j]*cos(jj*t)
+    yy = yy + yd
+  end
+  y[i] = yy*2.0
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+#      --- X AXIS ---
+DCL::uzrset('UYUSER', 0.0)
+DCL::uscset('CXSIDE', 'U')
+
+#      --- Y AXIS ---
+DCL::uziset('IROTLYL', 1)
+DCL::uziset('ICENTYL', 0)
+DCL::uscset('CYSIDE', 'L')
+
+#      --- etc. ---
+DCL::uziset('INNER', -1)
+DCL::uscset('CBLKT', '[]')
+
+DCL::ussttl('TIME', 'SEC', 'VOLTAGE', 'mV')
+DCL::usgrph(x, y)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk08.rb

@@ -0,0 +1,46 @@
+#
+# uspk08.rb
+#
+# $Id: uspk08.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 200
+t = NArray.sfloat(n)
+z = NArray.sfloat(n)
+
+#-----------------------------------------------------------------------
+r = 0.2
+a = 4.0
+z.indgen(2.0, 2.0)
+for i in 0..n-1
+  r = a*r*(1.0-r)
+  z2  = ((i-4).to_f/40.0)**2
+  t[i] = 20.0*exp(-z2) + r*2.0
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+#      --- OMAJINAI ---
+DCL::uslset('LYINV', true)
+DCL::uzlset('LABELXT', true)
+DCL::uscset('CYSPOS', 'B')
+
+DCL::uscset('CXSIDE', 'T')
+DCL::uscset('CYSIDE', 'L')
+
+DCL::ussttl('TEMPERATURE', 'C|'+ DCL::csgi(4) + '"', 'DEPTH', 'm')
+DCL::usgrph(t, z)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk09.rb

@@ -0,0 +1,89 @@
+#
+# uspk09.rb
+#
+# $Id: uspk09.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 200
+m = 5
+x = NArray.sfloat(n)
+y = NArray.sfloat(n)
+yc = NArray.sfloat(n, m)
+a = NArray.sfloat(m)
+
+#-----------------------------------------------------------------------
+dt = 1.0/(n-1)
+for j in 0..m-1
+  jj = j*2+1
+  a[j] = (-1.0)**(j+1) *2.0/(jj*PI)
+end
+
+x = x.indgen * dt
+for i in 0..n-1
+  t = dt*i*2*PI
+  y[i] = 0.0
+  for j in 0..m-1
+    jj = j*2+1
+    yc[i,j] = a[j]*cos(jj*t)
+    y[i]  = y[i] + yc[i,j]
+  end
+end
+
+#--------------------------- 1ST PAGE ----------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+rundef = DCL::glrget('RUNDEF')
+DCL::gropn iws
+
+DCL::grfrm
+DCL::grsvpt(0.2, 0.8, 0.2, 0.6)
+DCL::usspnt(rundef, yc)
+
+DCL::ussttl('X-AXIS', ' ', 'COMPONENTS', ' ')
+DCL::usgrph(x, yc[true,0])
+
+for j in 1..m-1
+  ip = j.divmod(4)[1] + 1
+  DCL::uulinz(x, yc[true, j], 1, ip)
+end
+
+#     --- NEW FIG ---
+DCL::grfig
+DCL::grsvpt(0.2, 0.8, 0.62, 0.82)
+DCL::uzlset('LABELXB', false)
+
+DCL::ussttl('X-AXIS', ' ', 'TOTAL', ' ')
+DCL::usgrph(x, y)
+
+#--------------------------- 2ND PAGE ----------------------------------
+DCL::grfrm
+DCL::grsvpt(0.2, 0.8, 0.2, 0.6)
+
+DCL::uzfact(0.5)
+DCL::uzlset('LABELXB', true)
+
+DCL::usspnt(rundef, yc)
+DCL::ussttl('X-AXIS', ' ', 'COMPONENTS', ' ')
+DCL::usgrph(x, yc[true,0])
+
+for j in 1..m-1
+  ip = j.divmod(4)[1] + 1
+  DCL::uulinz(x, yc[true, j], 1, ip)
+end
+
+#     --- NEW FIG ---
+DCL::grfig
+DCL::grsvpt(0.2, 0.8, 0.62, 0.82)
+DCL::uzlset('LABELXB', false)
+
+DCL::ussttl('X-AXIS', ' ', 'TOTAL', ' ')
+DCL::usgrph(x, y)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk10.rb

@@ -0,0 +1,48 @@
+#
+# uspk10.rb
+#
+# $Id: uspk10.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include NMath
+
+
+n = 400
+x = NArray.sfloat(n)
+y = NArray.sfloat(n)
+
+#-----------------------------------------------------------------------
+dt = 2*PI/(n-1)
+x = x.indgen * dt
+y = y.indgen * 2*dt
+y = sin(y) * 10.0
+
+#---------------------------- 1ST PAGE ---------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+DCL::usgrph(x, y)
+
+#---------------------------- 2ND PAGE ---------------------------------
+rundef = DCL::glrget('RUNDEF')
+
+DCL::grfrm
+DCL::grswnd(0.0, 2.0*PI, rundef, rundef)
+
+DCL::usrset('DXL', 1.25)
+DCL::usrset('DXT', 0.25)
+
+DCL::usrset('XFAC', 0.1)
+DCL::usrset('YFAC', 10.0)
+
+DCL::usgrph(x, y)
+
+#-----------------------------------------------------------------------
+
+DCL::grcls
+