ruby-dcl 1.6.2

data/demo/grph2/uspack/uspk11.rb

@@ -0,0 +1,71 @@
+#
+# uspk11.rb
+#
+# $Id: uspk11.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 181
+N1 = 26
+N2 = 7
+nc = 4
+t = NArray.sfloat(n)
+y = NArray.sfloat(n)
+ux1 = NArray.sfloat(N1)
+
+#-----------------------------------------------------------------------
+ch = ['JAN', 'FEB', 'MAR', 'APL', 'JUN', 'JULY', '   ']
+ux2 = NArray[0.0, 31.0, 59.0, 90.0, 120.0, 151.0, 181.0]
+
+r = 0.2
+a = 4.0
+y[0] = 120.0
+t.indgen(0.5, 1.0)
+for i in 1..n-1
+  r = a*r*(1.0-r)
+  y[i] = y[i-1] + (r-0.46)*2
+end
+
+ux1.indgen(1.0, 7.0)
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+
+rundef = DCL::glrget('RUNDEF')
+DCL::uslset('LYINV', true)
+
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(0.0, 181.0, rundef, rundef)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+
+DCL::usspnt(rundef, y)
+DCL::uspfit
+
+DCL::grstrf
+
+#------------------------------- Y-AXIS --------------------------------
+DCL::usyaxs('L')
+DCL::usyaxs('R')
+DCL::uysttl('L', 'YEN/DOLLAR', 0.0)
+
+#------------------------------- X-AXIS --------------------------------
+DCL::uzlset('LBTWN', true)
+DCL::uxaxlb('B', ux1, ux2, ch, nc)
+DCL::uxaxlb('T', ux1, ux2, ch, nc)
+DCL::uxsttl('T', 'EXCHANGE RATE', 0.0)
+
+#-------------------------------- LINE ---------------------------------
+DCL::sgplu(t, y)
+
+DCL::grcls
+

data/demo/grph2/uspack/uspk12.rb

@@ -0,0 +1,51 @@
+#
+# uspk12.rb
+#
+# $Id: uspk12.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nt = 51
+nz = 21
+zmin = 20.0
+zmax = 50.0
+tmax = 5.0
+dz = (zmax-zmin)/(nz-1)
+dt = tmax/(nt-1)
+u = NArray.sfloat(nt, nz)
+
+#-----------------------------------------------------------------------
+for j in 0..nz-1
+  z  = dz*j
+  uz = exp(-0.2*z)*(z**0.5)
+  for i in 0..nt-1
+    t  = dt*i - 2*exp(-0.1*z)
+    u[i,j] = uz*sin(3.0*t)
+  end
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(0.0, tmax, zmin, zmax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::ussttl('TIME', 'YEAR', 'HEIGHT', 'km')
+
+DCL::usdaxs
+
+DCL::udcntr(u)
+
+DCL::grcls
+

data/demo/grph2/uupack/uupk01.rb

@@ -0,0 +1,50 @@
+#
+# uupk01.rb
+#
+# $Id: uupk01.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
+
+y[n/2-1] = 999.0
+DCL::gllset('LMISS', true)
+rundef = DCL::glrget('RUNDEF')
+
+#----------------------------- GRAPH -----------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+DCL::sglset('LCLIP', true)
+
+DCL::grswnd(0.0, 1.0, rundef, rundef)
+DCL::usspnt(rundef, y)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uuslni(2)
+DCL::uulin(rundef, y)
+
+DCL::uusmkt(4)
+DCL::uumrk(rundef, y)
+
+DCL::grcls
+

data/demo/grph2/uupack/uupk02.rb

@@ -0,0 +1,74 @@
+#
+# uupk02.rb
+#
+# $Id: uupk02.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 20
+m = 5
+x0 = NArray.sfloat(n)
+x1 = NArray.sfloat(n)
+x2 = 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
+
+for i in 0..n-1
+  t = dt*i*2*PI
+  x0[i] = dt*i
+  y1[i] = a[0]*cos(t)
+  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
+
+  y0[i] = (y1[i]+y2[i])/2.0
+  dy = (y1[i] - y2[i]).abs
+  x1[i] = x0[i] - dy/5
+  x2[i] = x0[i] + dy/5
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+DCL::sglset('LCLIP', true)
+
+DCL::usspnt(x0, y0)
+DCL::usspnt(x1, y1)
+DCL::usspnt(x2, y2)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uverb(x0, y1, y2)
+
+DCL::uusebt(3)
+DCL::uusebi(33)
+DCL::uherb(x1, x2, y0)
+
+DCL::uusmkt(4)
+DCL::uusmki(21)
+DCL::uumrk(x0, y0)
+
+DCL::grcls
+

data/demo/grph2/uupack/uupk03.rb

@@ -0,0 +1,56 @@
+#
+# uupk03.rb
+#
+# $Id: uupk03.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
+
+y[n/2-1] = 999.0
+DCL::gllset('LMISS', true)
+rundef = DCL::glrget('RUNDEF')
+
+#----------------------------- GRAPH -----------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(0.0, 1.0, rundef, rundef)
+DCL::usspnt(rundef, y)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uvbxl(rundef, y)
+
+DCL::grfrm
+
+DCL::grswnd(rundef, rundef, 0.0, 1.0)
+DCL::usspnt(y, rundef)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uhbxl(y, rundef)
+
+DCL::grcls
+

data/demo/grph2/uupack/uupk04.rb

@@ -0,0 +1,84 @@
+#
+# uupk04.rb
+#
+# $Id: uupk04.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 20
+m = 5
+y1 = NArray.sfloat(n)
+y2 = NArray.sfloat(n)
+a = NArray.sfloat(m)
+rundef = DCL::glrget('RUNDEF')
+
+#-----------------------------------------------------------------------
+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
+
+for i in 0..n-1
+  t = dt*i*2*PI
+  y1[i] = a[0]*cos(t)
+  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
+
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+DCL::sglset('LCLIP', true)
+
+DCL::usspnt(rundef, y1)
+DCL::usspnt(rundef, y2)
+DCL::grswnd(0.0, 1.0, rundef, rundef)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uusarp(2999, 3999)
+DCL::uvbra(rundef, y1, y2)
+DCL::uvbrf(rundef, y1, y2)
+
+DCL::uuslni(41)
+DCL::uvbrl(rundef, y1)
+DCL::uuslnt(3)
+DCL::uvbrl(rundef, y2)
+
+DCL::grfrm
+
+DCL::usspnt(y1, rundef)
+DCL::usspnt(y1, rundef)
+DCL::grswnd(rundef, rundef, 0.0, 1.0)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uusarp(2999, 3999)
+DCL::uhbra(y1, y2, rundef)
+DCL::uhbrf(y1, y2, rundef)
+
+DCL::uuslni(41)
+DCL::uhbrl(y1, rundef)
+DCL::uuslnt(3)
+DCL::uhbrl(y2, rundef)
+
+DCL::grcls
+

data/demo/grph2/uupack/uupk05.rb

@@ -0,0 +1,74 @@
+#
+# uupk05.rb
+#
+# $Id: uupk05.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+n = 20
+m = 5
+y1 = NArray.sfloat(n)
+y2 = NArray.sfloat(n)
+a = NArray.sfloat(m)
+
+rundef = DCL::glrget('RUNDEF')
+
+#-----------------------------------------------------------------------
+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
+
+for i in 0..n-1
+  t = dt*i*2*PI
+  y1[i] = a[0]*cos(t)
+  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
+end
+
+#-----------------------------------------------------------------------
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sglset('LCLIP', true)
+
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::usspnt(rundef, y1)
+DCL::usspnt(rundef, y2)
+DCL::grswnd(0.0, 1.0, rundef, rundef)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uusarp(2999, 4999)
+DCL::uvbxa(rundef, rundef, y1)
+DCL::uvbxf(rundef, rundef, y1)
+
+DCL::grfrm
+
+DCL::usspnt(y1, rundef)
+DCL::usspnt(y2, rundef)
+DCL::grswnd(rundef, rundef, 0.0, 1.0)
+DCL::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::uhbxa(y1, y2, rundef)
+DCL::uhbxf(y1, y2, rundef)
+
+DCL::grcls
+