ruby-dcl 1.6.2

data/demo/rakuraku/miss/miss2.rb

@@ -0,0 +1,60 @@
+#
+# miss2.rb
+#   $Id: miss2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nx = 37
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI/180
+p = NArray.sfloat(nx, ny)
+
+rmiss = DCL::glrget('RMISS')
+DCL::gllset('LMISS', true)
+
+#-- data ---
+for j in 0..ny-1
+  for i in 0..nx-1
+    alon = (xmin + (xmax-xmin)*i/(nx-1)) * drad
+    alat = (ymin + (ymax-ymin)*j/(ny-1)) * drad
+    slat = sin(alat)
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon) - 0.5*(3*slat**2-1)
+    if ( i == 5 && j == 5 )
+      p[i,j] = rmiss
+    end
+    if ( ( 7 <= i && i <= 23) && j == 29 )
+      p[i,j] = rmiss
+    end
+    if( (21 <= i && i <= 29) && (11 <= j && j <= 19))
+      p[i,j] = rmiss
+    end
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uetone(p)
+DCL::usdaxs
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/rakuraku/quick/quick1.rb

@@ -0,0 +1,35 @@
+#
+# quick1.rb
+#   $Id: quick1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nmax = 50
+
+x = NArray.sfloat(nmax+1)
+y = NArray.sfloat(nmax+1)
+
+#-- data ----
+R    = 3.7
+x[0] = 0.0
+y[0] = 0.5
+
+for n in 0..nmax-1
+  x[n+1] = x[n] + 1.0
+  y[n+1] = R*y[n]*(1.0-y[n])
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::usgrph(x, y)
+
+DCL::grcls
+

data/demo/rakuraku/quick/quick2.rb

@@ -0,0 +1,33 @@
+#
+# quick2.rb
+#   $Id: quick2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nmax = 400
+dt = PI/(nmax-1)
+
+x = NArray.sfloat(nmax)
+y = NArray.sfloat(nmax)
+
+#-- data ----
+for n in 0..nmax-1
+  t = dt*n
+  x[n] = 1.0e5*sin(6.0*t)
+  y[n] = 1.0e-4*cos(14.0*t) + 1.0
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::usgrph(x, y)
+
+DCL::grcls
+

data/demo/rakuraku/quick/quick3.rb

@@ -0,0 +1,46 @@
+#
+# quick3.rb
+#   $Id: quick3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nx = 37
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI/180
+
+p = NArray.sfloat(nx, ny)
+
+#-- data ----
+for j in 0..ny-1
+  for i in 0..nx-1
+    alon = (xmin + (xmax-xmin)*i/(nx-1))*drad
+    alat = (ymin + (ymax-ymin)*j/(ny-1))*drad
+    slat = sin(alat)
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon) - 0.5*(3*slat**2-1)
+  end
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::usdaxs
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/rakuraku/quick/quick4.rb

@@ -0,0 +1,47 @@
+#
+# quick4.rb
+#   $Id: quick4.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nx = 37
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI/180
+
+p = NArray.sfloat(nx, ny)
+
+#-- data ----
+for j in 0..ny-1
+  for i in 0..nx-1
+    alon = (xmin + (xmax-xmin)*i/(nx-1))*drad
+    alat = (ymin + (ymax-ymin)*j/(ny-1))*drad
+    slat = sin(alat)
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon) - 0.5*(3*slat**2-1)
+  end
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uetone(p)
+DCL::usdaxs
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/rakuraku/quick/quick5.rb

@@ -0,0 +1,46 @@
+#
+# quick5.rb
+#   $Id: quick5.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nx = 21
+ny = 21
+xmin = -1.0
+xmax = 1.0
+ymin = -1.0
+ymax = 1.0
+
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+
+#-- data ----
+for j in 0..ny-1
+  for i in 0..nx-1
+    x = xmin + (xmax-xmin)*i/(nx-1)
+    y = ymin + (ymax-ymin)*j/(ny-1)
+    u[i,j] = x
+    v[i,j] = -y
+  end
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::usdaxs
+DCL::ugvect(u, v)
+
+DCL::grcls
+

data/demo/rakuraku/u2d/u2d1.rb

@@ -0,0 +1,50 @@
+#
+# u2d1.rb
+#   $Id: u2d1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+nx = 37
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI/180
+p = NArray.sfloat(nx, ny)
+
+#-- data ---
+for j in 0..ny-1
+  for i in 0..nx-1
+    alon = (xmin + (xmax-xmin)*i/(nx-1)) * drad
+    alat = (ymin + (ymax-ymin)*j/(ny-1)) * drad
+    slat = sin(alat)
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon) - 0.5*(3*slat**2-1)
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::usdaxs
+DCL::udsfmt('(F4.1)')
+DCL::udgcla(0.0, 1.4, 0.2)
+DCL::udsclv(-1.0, 3, 4, 'abc', 0.028)
+DCL::uddclv(0.6)
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/rakuraku/u2d/u2d2.rb

@@ -0,0 +1,51 @@
+#
+# u2d2.rb
+#   $Id: u2d2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nx = 37
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI / 180
+p = NArray.sfloat(nx, ny)
+
+#-- data ---
+for j in 0..ny-1
+  for i in 0..nx-1
+    alon = (xmin + (xmax-xmin)*i/(nx-1)) * drad
+    alat = (ymin + (ymax-ymin)*j/(ny-1)) * drad
+    slat = sin(alat)
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon) - 0.5*(3*slat**2-1)
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::usdaxs
+DCL::udlset('LABEL', false)
+DCL::udlset('LDASH', false)
+DCL::udiset('isolid', 2)    
+DCL::udiset('icycle', 4)    
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/rakuraku/u2d/u2d3.rb

@@ -0,0 +1,70 @@
+#
+# u2d3.rb
+#   $Id: u2d3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nx = 37
+ny = 37
+my = 7
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI / 180
+p = NArray.sfloat(nx, ny)
+uy1 = NArray.sfloat(ny)
+uy2 = NArray.sfloat(my)
+
+ch = ['SP ', '60S', '30S', 'EQ ', '30N', '60N', 'NP ']
+
+#-- data ---
+for j in 0..ny-1
+  for i in 0..nx-1
+    alon = (xmin + (xmax-xmin)*i/(nx-1)) * drad
+    alat = (ymin + (ymax-ymin)*j/(ny-1)) * drad
+    slat = sin(alat)
+    uy1[j] = slat
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon) - 0.5*(3*slat**2-1)
+  end
+end
+
+for j in 0..my-1
+  uy2[j] = sin((ymin + (ymax-ymin)*j/(my-1)) * drad )
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, -1.0, 1.0)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uxaxdv('B', 10.0, 60.0)
+DCL::uxaxdv('T', 10.0, 60.0)
+DCL::uxsttl('B', 'LONGITUDE', 0.0)
+
+DCL::uyaxdv('L', 0.1, 0.5)
+DCL::uysttl('L', 'SINE LATITUDE', 0.0)
+
+DCL::uzlset('LABELYR', true)
+DCL::uyaxlb('R', uy1, uy2, ch, 3)
+DCL::uysttl('R', 'LATITUDE', 0.0)
+
+DCL::uwsgxb(xmin, xmax, nx)
+DCL::uwsgya(uy1)
+
+DCL::udcntr(p)
+
+DCL::grcls
+