ruby-dcl 1.6.2

data/demo/rakuraku/kihon/kihonb.rb

@@ -0,0 +1,73 @@
+#
+# kihonb.rb
+#   $Id: kihonb.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nmax = 9
+y = NArray.sfloat(nmax)
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sgfrm
+
+x1 = 0.1
+x2 = 0.9
+xc = 0.5
+
+#-- ���� ----
+for n in 1..nmax
+  y[n-1] = 0.1*(10-n)
+  DCL::sglnv(x1, y[n-1], x2, y[n-1])
+end
+
+DCL::sglnv(xc, 0.05, xc, 0.95)
+
+#-- �ǥե���� ----
+DCL::sgtxv(xc, y[0], 'SGTXV|SUP"RST_SUB"')
+
+#-- ź�� ----
+DCL::sglset('LCNTL', true)
+DCL::sgtxv(xc, y[1], 'SGTXV|SUP"RST_SUB"')
+
+DCL::sgrset('SMALL', 0.5)
+DCL::sgrset('SHIFT', 0.5)
+DCL::sgtxv(xc, y[2], 'SGTXV|SUP"RST_SUB"')
+
+#-- �ե���� ----
+DCL::sgtxv(xc, y[3], 'ABCDEFG abcdefg')
+
+DCL::sgiset('IFONT', 2)
+DCL::sgtxv(xc, y[4], 'ABCDEFG abcdefg')
+
+DCL::sgstxi(3)
+DCL::sgtxv(xc, y[5], 'ABCDEFG abcdefg')
+
+DCL::sgiset('IFONT', 1)
+DCL::sgstxi(1)
+
+#-- ���ꥷ��ʸ�� ----
+greek1 = DCL::csgi(128) + DCL::csgi(129) + DCL::csgi(130) + DCL::csgi(131) + DCL::csgi(132) +
+         DCL::csgi(133) + DCL::csgi(134) + DCL::csgi(135) + DCL::csgi(136) + DCL::csgi(137)
+greek2 = DCL::csgi(152) + DCL::csgi(153) + DCL::csgi(154) + DCL::csgi(155) + DCL::csgi(156) +
+         DCL::csgi(157) + DCL::csgi(158) + DCL::csgi(159) + DCL::csgi(160) + DCL::csgi(161)
+
+DCL::sgtxv(xc, y[6], greek1)
+DCL::sgtxv(xc, y[7], greek2)
+
+#-- �ü쵭�� ----
+symbol = DCL::csgi(189) + DCL::csgi(190) + DCL::csgi(191) + DCL::csgi(192) + DCL::csgi(193) +
+         DCL::csgi(210) + DCL::csgi(211) + DCL::csgi(212) + DCL::csgi(217) + DCL::csgi(218)
+
+DCL::sgtxv(xc, y[8], symbol)
+
+DCL::sgcls
+

data/demo/rakuraku/kihon/kihonc.rb

@@ -0,0 +1,110 @@
+#
+# kihonc.rb
+#   $Id: kihonc.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sgfrm
+
+DCL::sgswnd(0.0, 10.0, 0.0, 10.0)
+DCL::sgsvpt(0.0, 1.0, 0.0, 1.0)
+DCL::sgstrn(1)
+DCL::sgstrf
+
+#-- �ǥե���� ----
+y1 = 9.0
+for i in 1..8
+  x1 = i
+  x2 = x1 + 0.1*i
+  y2 = y1 + 0.1*i
+  DCL::sglau(x1, y1, x2, y2)
+end
+
+#-- ��ʬ�Υ饤�󥿥��� ----
+y1 = 8.0
+y2 = 8.6
+for i in 1..4
+  x1 = 2*i - 1
+  x2 = x1 + 0.6
+  DCL::sgslat(i)
+  DCL::sglau(x1, y1, x2, y2)
+end
+
+DCL::sgslat(1)
+
+#-- ��ʬ�Υ饤�󥤥�ǥ��� ----
+y1 = 7.0
+y2 = 7.6
+for i in 1..4
+  x1 = 2*i
+  x2 = x1 + 0.6
+  DCL::sgslai(i)
+  DCL::sglau(x1, y1, x2, y2)
+end
+
+DCL::sgslai(1)
+
+#-- �𤸤���ʬ��Ĺ�� ----
+DCL::sglset('LPROP', false)
+DCL::sgrset('CONST', 0.03)
+
+y1 = 6.0
+for i in 1..8
+  x1 = i
+  x2 = x1 + 0.1*i
+  y2 = y1 + 0.1*i
+  DCL::sglau(x1, y1, x2, y2)
+end
+
+DCL::sglset('LPROP', true)
+
+#-- �𤸤���ʬ�γ��� ----
+y1 = 5.0
+y2 = 5.6
+for i in 1..8
+  x1 = i
+  x2 = x1 + 0.6
+  DCL::sgrset('ANGLE', 10.0*i)
+  DCL::sglau(x1, y1, x2, y2)
+end
+
+#-- �𤸤���ʬ�Τ̤�Ĥ֤� ----
+DCL::sglset('LSOFTF', true)
+DCL::sglset('LATONE', true)
+DCL::sgiset('IATONE', 655)
+
+y1 = 4.0
+y2 = 4.6
+for i in 1..8
+  x1 = i
+  x2 = x1 + 0.6
+  DCL::sgrset('ANGLE', 10.0*i)
+  DCL::sglau(x1, y1, x2, y2)
+end
+
+#-- �饤�󥵥֥ץ�ߥƥ��� ----
+#-- �ǥե���� ----
+DCL::sglnu(0.0, 3.0, 10.0, 3.0)
+
+#-- ��ʬ�Υ饤�󥤥�ǥ��� ----
+y1 = 0.5
+y2 = 2.5
+for i in 1..4
+  x1 = 2*i - 1
+  x2 = x1 + 2.0
+  DCL::sgslni(i)
+  DCL::sglnu(x1, y1, x2, y2)
+end
+
+DCL::sgcls
+

data/demo/rakuraku/layout/lay1.rb

@@ -0,0 +1,29 @@
+#
+# lay1.rb
+#   $Id: lay1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::slmgn(0.1, 0.1, 0.05, 0.05)
+DCL::sldiv('Y', 3, 2)
+DCL::slmgn(0.05, 0.05, 0.05, 0.05)
+
+for i in 1..12
+  DCL::sgfrm
+  DCL::slpvpr(1)
+  ctxt = format("FRAME%2.2d", i)
+  DCL::sgtxv(0.5, 0.5, ctxt)
+end
+
+DCL::sgcls
+

data/demo/rakuraku/layout/lay2.rb

@@ -0,0 +1,28 @@
+#
+# lay2.rb
+#   $Id: lay2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::slmgn(0.1, 0.1, 0.1, 0.1)
+DCL::slsttl( 'FIGURE TITLE',  'T', 0.0, 0.0, 0.03, 1)
+DCL::slsttl( 'PROGRAM.NAME',  'B', -1.0, 1.0, 0.02, 2)
+DCL::slsttl( '#DATE #TIME',   'B', 0.0, 0.0, 0.02, 3)
+DCL::slsttl( 'page:#PAGE',    'B', 1.0, -1.0, 0.02, 4)
+
+DCL::sgfrm
+DCL::slpvpr(1)
+DCL::sgtxv(0.5, 0.5,  'FIGURE' )
+
+DCL::sgcls
+

data/demo/rakuraku/layout/lay3.rb

@@ -0,0 +1,29 @@
+#
+# lay3.rb
+#   $Id: lay3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sglset('LFULL', true)
+DCL::slmgn(0.0, 0.0, 0.08, 0.08)
+DCL::slrat(1.0, 0.6)
+DCL::slsttl('FIGURE TITLE', 'T', 0.0, 0.0, 0.03, 1)
+DCL::slsttl('PROGRAM.NAME', 'B', -1.0, -1.0, 0.02, 2)
+DCL::slsttl('page:#PAGE', 'B', 1.0, -1.0, 0.02, 3)
+
+DCL::sgfrm
+DCL::slpvpr(1)
+DCL::sgtxv(0.5, 0.5, 'FIGURE')
+
+DCL::sgcls
+

data/demo/rakuraku/map3d/map3d1.rb

@@ -0,0 +1,68 @@
+#
+# map3d1.rb
+#   $Id: map3d1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+np = 14
+
+#-- data ---
+itr =         [10, 11, 12, 13, 14, 15, 20, 21, 22, 23, 30, 31, 32, 33]
+fct = NArray[12.0, 12, 14, 14, 14, 14, 11, 16, 12, 12, 40, 12, 12, 17]
+fct = fct/100.0
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn(-(iws).abs)
+
+DCL::slrat(2.0, 3.0)
+DCL::sldiv('Y', 2, 3)
+
+DCL::sglset('LCLIP', true)
+DCL::sgrset('STLAT1', 45.0)
+DCL::sgrset('STLAT2', 30.0)
+
+DCL::umlset('LGRIDMN', false)
+DCL::umiset('INDEXMJ', 1)
+
+DCL::sgstxs(0.035)
+DCL::sgstxi(3)
+
+for i in 0..np-1
+  DCL::sgfrm
+
+  DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
+  DCL::sgssim(fct[i], 0.0, 0.0)
+  DCL::sgsmpl(0.0, 90.0, 0.0)
+  if ( itr[i] == 30 )
+    DCL::sgstxy(-180.0, 180.0, 0.0, 90.0)
+  else
+    DCL::sgstxy(-180.0, 180.0, -90.0, 90.0)
+  end
+  DCL::sgstrn(itr[i])
+  DCL::sgstrf
+
+  DCL::slpwwr(1)
+  DCL::slpvpr(1)
+  cttl = DCL::sgtrnl(itr[i])
+  DCL::sgtxr(0.5, 0.95, cttl)
+  cttl = format("ITR=%2d", itr[i]);
+  DCL::sgtxr(0.5, 0.05, cttl)
+
+  DCL::umpmap('coast_world')
+  DCL::umpglb
+
+  if ( itr[i] == 23 )
+    DCL::sgfrm
+    DCL::sgfrm
+  end
+end
+
+DCL::sgcls
+

data/demo/rakuraku/map3d/map3d2.rb

@@ -0,0 +1,67 @@
+#
+# map3d2.rb
+#   $Id: map3d2.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::sglset('LSOFTF', false)
+DCL::grfrm
+
+DCL::sgswnd(xmin, xmax, ymin, ymax)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
+DCL::sgssim(0.4, 0.0, 0.0)
+DCL::sgsmpl(135.0, 35.0, 0.0)
+DCL::sgstxy(-180.0, 180.0, 0.0, 90.0)
+DCL::sgstrn(30)
+DCL::sgstrf
+
+for k in -5..3
+  tlev1 = 0.4*k
+  tlev2 = tlev1 + 0.4
+  if (k <= -1)
+    ipat = 600 + (k+1).abs
+  else
+    ipat =  30 + k
+  end
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+
+DCL::uetone(p)
+
+DCL::udgclb(p, 0.4)
+DCL::udcntr(p)
+
+DCL::umpmap('coast_world')
+DCL::umpglb
+
+DCL::grcls
+

data/demo/rakuraku/map3d/map3d3.rb

@@ -0,0 +1,66 @@
+#
+# map3d3.rb
+#   $Id: map3d3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+nx = 16
+ny = 16
+xmin = 120
+xmax = 150
+ymin = 20
+ymax = 50
+fact = 0.21
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+alon = NArray.sfloat(nx)
+alat = NArray.sfloat(ny)
+
+#-- data ---
+for i in 0..nx-1
+  alon[i] = xmin + (xmax-xmin)*i/(nx-1).to_f
+end
+for j in 0..ny-1
+  alat[j] = ymin + (ymax-ymin)*j/(ny-1).to_f
+end
+
+for j in 0..ny-1
+  for i in 0..nx-1
+    u[i,j] =   (j-(ny-1)/2.0) * fact
+    v[i,j] = - (i-(nx-1)/2.0) * fact
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+DCL::sglset('LCLIP', true)
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.1, 0.9, 0.1, 0.9)
+DCL::grstrn(10)
+DCL::umpfit
+DCL::grstrf
+
+DCL::slpvpr(3)
+DCL::umrset('DGRIDMJ', 10.0)
+DCL::umrset('DGRIDMN', 2.0)
+DCL::umpmap('coast_world')
+DCL::umpglb
+
+DCL::sgslai(3)
+for j in 0..ny-1
+  for i in 0..nx-1
+    DCL::sglau(alon[i], alat[j], alon[i]+u[i,j], alat[j]+v[i,j])
+  end
+end
+
+DCL::grcls
+