ruby-dcl 1.6.2

data/demo/grph2/grpack/grpk01.rb

@@ -0,0 +1,65 @@
+#
+# grpk01.rb
+#
+# $Id: grpk01.rb,v 1.1.1.1 2011-02-23 07:21:27 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::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(0.0, 1.0, -0.7, 0.7)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.6)
+DCL::grstrn(1)
+DCL::grstrf
+DCL::ussttl('X-AXIS', ' ', 'Y1-AXIS', ' ')
+DCL::usdaxs
+DCL::sgtxzu(0.5, 0.0, 'PAGE:1, FRAME:1', 0.02, 0, 0, 3)
+
+DCL::grfig
+
+DCL::grswnd(0.0, 1.0, -0.5, 0.5)
+DCL::grsvpt(0.2, 0.8, 0.62, 0.82)
+DCL::grstrn(1)
+DCL::grstrf
+DCL::uzlset('LABELXB', false)
+DCL::ussttl('X-AXIS', ' ', 'Y2-AXIS', ' ')
+DCL::usdaxs
+DCL::sgtxzu(0.5, 0.0, 'PAGE:1, FRAME:2', 0.02, 0, 0, 3)
+DCL::uzlset('LABELXB', true)
+
+DCL::grfrm
+
+DCL::uzfact(0.5)
+
+DCL::grswnd(0.0, 1.0, -0.7, 0.7)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.6)
+DCL::grstrn(1)
+DCL::grstrf
+DCL::ussttl('X-AXIS', ' ', 'Y1-AXIS', ' ')
+DCL::usdaxs
+DCL::sgtxzu(0.5, 0.0, 'PAGE:2, FRAME:1', 0.02, 0, 0, 3)
+
+DCL::grfig
+
+DCL::grswnd(0.0, 1.0, -0.5, 0.5)
+DCL::grsvpt(0.2, 0.8, 0.62, 0.82)
+DCL::grstrn(1)
+DCL::grstrf
+DCL::uzlset('LABELXB', false)
+DCL::ussttl('X-AXIS', ' ', 'Y2-AXIS', ' ')
+DCL::usdaxs
+DCL::sgtxzu(0.5, 0.0, 'PAGE:2, FRAME:2', 0.02, 0, 0, 3)
+DCL::uzlset('LABELXB', true)
+
+DCL::grcls
+

data/demo/grph2/ucpack/ucpk01.rb

@@ -0,0 +1,95 @@
+#
+# ucpk01.rb
+#
+# $Id: ucpk01.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+jd0 = 19920401
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn(-(iws.abs))
+
+DCL::sglset('LFULL', true)
+DCL::uzfact(0.7)
+DCL::slrat(0.75, 1.0)
+DCL::sldiv('T', 1, 7)
+
+DCL::sgfrm
+DCL::uzinit
+nd=30
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxacl('B', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=90
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxacl('B', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=180
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxacl('B', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=400
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxacl('B', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=60
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxady('B', jd0, nd)
+ctl = 'UCXADY (ND=60)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+nd=120
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxamn('B', jd0, nd)
+ctl = 'UCXAMN (ND=120)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+nd=2000
+DCL::sgswnd(0.0, nd.to_f, 0.0, 1.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucxayr('B', jd0, nd)
+ctl = 'UCXAYR (ND=2000)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgcls
+

data/demo/grph2/ucpack/ucpk02.rb

@@ -0,0 +1,95 @@
+#
+# ucpk02.rb
+#
+# $Id: ucpk02.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+jd0 = 19920401
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn(iws.abs)
+
+DCL::sglset('LFULL', true)
+DCL::uzfact(0.7)
+DCL::slrat(1.0, 0.75)
+DCL::sldiv('Y', 7, 1)
+
+DCL::sgfrm
+DCL::uzinit
+nd=30
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyacl('L', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=90
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyacl('L', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=180
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyacl('L', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=400
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyacl('L', jd0, nd)
+
+DCL::sgfrm
+DCL::uzinit
+nd=60
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyady('L', jd0, nd)
+ctl = 'UCYADY (ND=60)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+nd=120
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyamn('L', jd0, nd)
+ctl = 'UCYAMN (ND=120)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+nd=2000
+DCL::sgswnd(0.0, 1.0, 0.0, nd.to_f)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(1)
+DCL::sgstrf
+DCL::ucyayr('L', jd0, nd)
+ctl = 'UCYAYR (ND=2000)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgcls
+

data/demo/grph2/udegpk/u2df01.rb

@@ -0,0 +1,50 @@
+#
+# u2df01.rb
+#
+# $Id: u2df01.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nx = 19
+ny = 19
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+drad = PI/180
+dz = 0.05
+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] = cos(alon) * (1-slat**2) * sin(2*PI*slat) + dz
+  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/grph2/udegpk/u2df02.rb

@@ -0,0 +1,76 @@
+#
+# u2df02.rb
+#
+# $Id: u2df02.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nx = 19
+ny = 19
+xmin = 0
+xmax = 360
+dx1 = 10
+dx2 = 60
+ymin = -90
+ymax = 90
+my = 7
+nc = 3
+drad = PI/180
+dz = 0.05
+p = NArray.sfloat(nx, ny)
+uy1 = NArray.sfloat(ny)
+uy2 = NArray.sfloat(my)
+
+#-- data ---
+ch = ['SP ', '60S', '30S', 'EQ ', '30N', '60N', 'NP ']
+
+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] = cos(alon) * (1-slat**2) * sin(2*PI*slat) + dz
+  end
+end
+
+uy1.indgen(ymin, (ymax-ymin)/(ny-1))
+uy1 = NMath.sin(uy1*drad)
+uy2.indgen(ymin, (ymax-ymin)/(my-1))
+uy2 = NMath.sin(uy2*drad)
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, uy1[0], uy1[ny-1])
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uxaxdv('B', dx1, dx2)
+DCL::uxaxdv('T', dx1, dx2)
+DCL::uxsttl('B', 'LONGITUDE', 0.0)
+
+DCL::uyaxlb('L', uy1, uy2, ch, nc)
+DCL::uyaxlb('R', uy1, uy2, ch, nc)
+DCL::uysttl('L', 'LATITUDE', 0.0)
+
+DCL::uwsgxb(xmin, xmax, nx)
+DCL::uwsgya(uy1)
+
+DCL::udsfmt('(F6.1)')
+DCL::udgclb(p, 0.2)
+DCL::udsclv(0.1, 1, 4, ' ', 0.01)
+
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/grph2/udegpk/u2df03.rb

@@ -0,0 +1,80 @@
+#
+# u2df03.rb
+#
+# $Id: u2df03.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nx = 19
+ny = 19
+xmin = 0
+xmax = 360
+dx1 = 10
+dx2 = 60
+ymin = -90
+ymax = 90
+my = 7
+nc = 3
+drad = PI/180
+dz = 0.05
+
+XMNZ = 60
+XMXZ = 300
+NXZ = 13
+NYZ = 13
+IXZ1 = 4
+IYZ1 = 4
+
+p = NArray.sfloat(nx, ny)
+uy1 = NArray.sfloat(ny)
+uy2 = NArray.sfloat(my)
+
+#-- data ---
+ch = ['SP ', '60S', '30S', 'EQ ', '30N', '60N', 'NP ']
+
+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] = cos(alon) * (1-slat**2) * sin(2*PI*slat) + dz
+  end
+end
+
+uy1.indgen(ymin, (ymax-ymin)/(ny-1).to_f)
+uy1 = NMath.sin(uy1*drad)
+uy2.indgen(ymin, (ymax-ymin)/(my-1).to_f)
+uy2 = NMath.sin(uy2*drad)
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, uy1[0], uy1[ny-1])
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uxaxdv('B', dx1, dx2)
+DCL::uxaxdv('T', dx1, dx2)
+DCL::uxsttl('B', 'LONGITUDE', 0.0)
+
+DCL::uyaxlb('L', uy1, uy2, ch, nc)
+DCL::uyaxlb('R', uy1, uy2, ch, nc)
+DCL::uysttl('L', 'LATITUDE', 0.0)
+
+DCL::uwsgxb(XMNZ, XMXZ, NXZ)
+DCL::uwsgya(uy1[(IYZ1-1)..-1])
+
+DCL::udcntr(p[(IXZ1-1)..-1,(IYZ1-1)..-1])
+
+DCL::grcls
+