ruby-dcl 1.6.2

data/demo/rakuraku/map3d/map3d4.rb

@@ -0,0 +1,107 @@
+#
+# map3d4.rb
+#   $Id: map3d4.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+xmin = -50
+xmax = 50
+ymin = -50
+ymax = 50
+vxmin = 0.0
+vxmax = 0.8
+vymin = 0.0
+vymax = 0.8
+zmin = -50
+zmax = 50
+vzmin = 0.0
+vzmax = 0.8
+xvp3 = 2.5
+yvp3 = -1.0
+zvp3 = 1.5
+xfc3 = (vxmax-vxmin)/2
+yfc3 = (vymax-vymin)/2
+zfc3 = (vzmax-vzmin)/2
+
+
+def aplot(ijk)
+
+  nmax = 40
+  x = NArray.sfloat(nmax+1)
+  y = NArray.sfloat(nmax+1)
+
+  DCL::slpvpr(1)
+
+  dt = 2.0*PI / nmax
+  for n in 0..nmax
+    x[n] = 40.0*sin(n*dt)
+    y[n] = 40.0*cos(n*dt)
+  end
+
+ DCL:: sgplu(x, y)
+
+  dt = 2.0*PI / 3
+  for n in 0..3
+    x[n] = 40.0*sin(n*dt)
+    y[n] = 40.0*cos(n*dt)
+  end
+
+  DCL::sgplu(x, y)
+
+  DCL::sgstxs(0.07)
+  DCL::sgtxu(0.0, 0.0, 'DENNOU')
+  cttl= format("%2.2d",ijk)
+  DCL::sgtxu(0.0, -30.0, cttl)
+
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sgiset('IFONT', 2)
+
+DCL::sgfrm
+
+#-- X-Y PLANE ----
+DCL::sgswnd(xmin, xmax, ymin, ymax)
+DCL::sgsvpt(vxmin, vxmax, vymin, vymax)
+DCL::sgstrn(1)
+DCL::sgstrf
+
+DCL::scspln(1, 2, vzmax)
+DCL::scseye(xvp3, yvp3, zvp3)
+DCL::scsobj(xfc3, yfc3, zfc3)
+DCL::scsprj
+
+aplot(1)
+
+#-- X-Z PLANE ----
+DCL::sgswnd(xmin, xmax, zmin, zmax)
+DCL::sgsvpt(vxmin, vxmax, vzmin, vzmax)
+DCL::sgstrn(1)
+DCL::sgstrf
+
+DCL::scspln(1, 3, vymin)
+DCL::scsprj
+
+aplot(2)
+
+#-- Y-Z PLANE ----
+DCL::sgswnd(ymin, ymax, zmin, zmax)
+DCL::sgsvpt(vymin, vymax, vzmin, vzmax)
+DCL::sgstrn(1)
+DCL::sgstrf
+
+DCL::scspln(2, 3, vxmax)
+DCL::scsprj
+
+aplot(3)
+
+DCL::sgcls

data/demo/rakuraku/map3d/map3d5.rb

@@ -0,0 +1,118 @@
+#
+# map3d5.rb
+#   $Id: map3d5.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 = -10
+xmax = 10
+ymin = -10
+ymax = 10
+vxmin = 0.2
+vxmax = 0.8
+vymin = 0.1
+vymax = 0.5
+zmin = 0.0
+zmax = 20
+vzmin = 0.0
+vzmax = 0.6
+xvp3 = -0.7
+yvp3 = -0.7
+zvp3 = 1.2
+xfc3 = (vxmax-vxmin)/2
+yfc3 = (vymax-vymin)/2
+zfc3 = (vzmax-vzmin)/2
+dx1 = 1
+dx2 = 5
+dy1 = 1
+dy2 = 4
+dz1 = 1
+dz2 = 5
+kmax = 5
+pmin = 0
+pmax = 1
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+p = NArray.sfloat(nx, ny)
+
+#-- data ---
+for j in 0..ny-1
+  for i in 0..nx-1
+    x = xmin + (xmax-xmin)*i/(nx-1).to_f
+    y = ymin + (ymax-ymin)*j/(ny-1).to_f
+    u[i,j] =  x
+    v[i,j] = -y
+    p[i,j] = exp(-x**2/64 -y**2/25)
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sgfrm
+
+#-- X-Y ʿ��: ����٥� ----
+DCL::sgswnd(xmin, xmax, ymin, ymax)
+DCL::sgsvpt(vxmin, vxmax, vymin, vymax)
+DCL::sgstrn(1)
+DCL::sgstrf
+
+DCL::scspln(1, 2, vzmin)
+DCL::scseye(xvp3, yvp3, zvp3)
+DCL::scsobj(xfc3, yfc3, zfc3)
+DCL::scsprj
+
+DCL::uxaxdv('B', dx1, dx2)
+DCL::uxaxdv('T', dx1, dx2)
+DCL::uxsttl('B', 'X-axis', 0.0)
+
+DCL::uyaxdv('L', dy1, dy2)
+DCL::uyaxdv('R', dy1, dy2)
+DCL::uysttl('L', 'Y-axis', 0.0)
+
+DCL::ugrset('RSIZET', 0.014)
+DCL::ugvect(u, v)
+
+#-- X-Y ʿ��: ���٥� ----
+vzlev = vzmin + (vzmax-vzmin)*0.6
+DCL::scspln(1, 2, vzlev)
+DCL::scsprj
+
+dp = (pmax-pmin)/kmax.to_f
+for k in 1..kmax
+  tlev1 = (k-1)*dp
+  tlev2 = tlev1 + dp
+  ipat  = 600 + k - 1
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+DCL::uetone(p)
+
+DCL::udlset('LMSG', false)
+DCL::udgclb(p, 0.1)
+DCL::udcntr(p)
+DCL::slpvpr(1)
+
+#-- X-Z ʿ�� ----
+DCL::sgswnd(xmin, xmax, zmin, zmax)
+DCL::sgsvpt(vxmin, vxmax, vzmin, vzmax)
+DCL::sgstrn(1)
+DCL::sgstrf
+
+DCL::scspln(1, 3, vymax)
+DCL::scsprj
+
+DCL::uzinit
+DCL::uyaxdv('L', dz1, dz2)
+DCL::uysttl('L', 'Z-axis', 0.0)
+
+DCL::sgcls
+

data/demo/rakuraku/map3d/map3d6.rb

@@ -0,0 +1,96 @@
+# map3d6.rb
+
+require "narray"
+require "numru/dcl"
+include NMath
+include NumRu
+
+nx=37; ny=37
+xmin=0.0;  xmax=360.0;  ymin= -90.0;  ymax= 90.0
+vxmin=-0.4; vxmax=0.4; vymin=-0.3; vymax=0.3
+zmin=-1.9;  zmax=1.9
+vzmin=-0.2; vzmax=0.2
+drad=PI/180
+
+alon = xmin + (xmax-xmin)/(nx-1) * NArray.sfloat(nx).indgen!
+alat = ymin + (ymax-ymin)/(ny-1) * NArray.sfloat(ny).indgen!
+
+slat = sin(alat*drad).reshape!(1,ny)
+p = 3*sqrt(1.0-slat**2)*slat*cos(alon*drad) - 0.5*(3*slat**2-1)
+
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+DCL::sgfrm
+
+#-- X-Y PLANE ----
+
+DCL::sgswnd(  xmin,  xmax,  ymin,  ymax )
+DCL::sgsvpt( vxmin, vxmax, vymin, vymax )
+DCL::sgstrn( 1 )
+DCL::sgstrf
+
+DCL::scspln( 1, 2, vzmin)
+DCL::scseye( -1.1, -1.1, 2.5 )
+DCL::scsobj(  0.0,  0.0, 0.0 )
+DCL::scsprj
+
+DCL::uxaxdv( 'B', 10.0, 60.0 )
+DCL::uxaxdv( 'T', 10.0, 60.0 )
+DCL::uxsttl( 'B', 'LONGITUDE', 0.0 )
+DCL::uyaxdv( 'L', 10.0, 30.0 )
+DCL::uyaxdv( 'R', 10.0, 30.0 )
+DCL::uysttl( 'L', 'LATITUDE', 0.0 )
+
+#-- X-Z PLANE ----
+
+DCL::sgswnd(  xmin,  xmax,  zmin,  zmax )
+DCL::sgsvpt( vxmin, vxmax, vzmin, vzmax )
+DCL::sgstrn( 1 )
+DCL::sgstrf
+
+DCL::scspln( 1, 3, vymax)
+DCL::scsprj
+
+DCL::uzinit
+DCL::uxaxdv( 'T', 10.0, 60.0 )
+DCL::uzlset( 'LABELXB', false )
+DCL::uxaxdv( 'B', 10.0, 60.0 )
+DCL::uyaxdv( 'L', 0.2, 1.0 )
+DCL::uyaxdv( 'R', 0.2, 1.0 )
+DCL::uysttl( 'L', 'AMPLITUDE', 0.0 )
+
+#-- Y-Z PLANE ----
+
+DCL::sgswnd(  ymin,  ymax,  zmin,  zmax )
+DCL::sgsvpt( vymin, vymax, vzmin, vzmax )
+DCL::sgstrn( 1 )
+DCL::sgstrf
+
+DCL::scspln( 2, 3, vxmax)
+DCL::scsprj
+
+DCL::uzinit
+DCL::uxaxdv( 'T', 10.0, 30.0 )
+DCL::uxaxdv( 'B', 10.0, 30.0 )
+DCL::uzlset( 'LABELYL', false )
+DCL::uyaxdv( 'L', 0.2, 1.0 )
+DCL::uyaxdv( 'R', 0.2, 1.0 )
+
+#---------------- 3-D ------------------
+
+DCL::scsvpt(vxmin, vxmax, vymin, vymax, vzmin, vzmax)
+DCL::scswnd( xmin,  xmax,  ymin,  ymax,  zmin,  zmax)
+DCL::scslog(false, false, false)
+DCL::scstrn(1)
+DCL::scstrf
+
+(ny-2).downto(0){|j|
+  (nx-2).downto(0){|i|
+    xp = [ alon[i], alon[i],   alon[i+1],  alon[i+1], alon[i] ]
+    yp = [ alat[j], alat[j+1], alat[j+1],  alat[j],   alat[j] ]
+    zp = [ p[i,j],  p[i,j+1],  p[i+1,j+1], p[i+1,j],  p[i,j] ]
+    DCL::scplu(xp, yp, zp)
+ }
+}
+
+DCL::sgcls

data/demo/rakuraku/map3d/map3d7.rb

@@ -0,0 +1,98 @@
+# map3d6.rb
+
+require "narray"
+require "numru/dcl"
+include NMath
+include NumRu
+
+nx=37; ny=37
+xmin=0.0;  xmax=360.0;  ymin= -90.0;  ymax= 90.0
+vxmin=-0.4; vxmax=0.4; vymin=-0.3; vymax=0.3
+zmin=-1.9;  zmax=1.9
+vzmin=-0.2; vzmax=0.2
+drad=PI/180
+
+alon = xmin + (xmax-xmin)/(nx-1) * NArray.sfloat(nx).indgen!
+alat = ymin + (ymax-ymin)/(ny-1) * NArray.sfloat(ny).indgen!
+
+slat = sin(alat*drad).reshape!(1,ny)
+p = 3*sqrt(1.0-slat**2)*slat*cos(alon*drad) - 0.5*(3*slat**2-1)
+
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+DCL::sgfrm
+
+#-- X-Y PLANE ----
+
+DCL::sgswnd(  xmin,  xmax,  ymin,  ymax )
+DCL::sgsvpt( vxmin, vxmax, vymin, vymax )
+DCL::sgstrn( 1 )
+DCL::sgstrf
+
+DCL::scspln( 1, 2, vzmin)
+DCL::scseye( -2.1, -1.1, 0.7 )
+DCL::scsobj(  0.0,  0.0, 0.0 )
+DCL::scsprj
+
+DCL::uxaxdv( 'B', 10.0, 60.0 )
+DCL::uxaxdv( 'T', 10.0, 60.0 )
+DCL::uxsttl( 'B', 'LONGITUDE', 0.0 )
+DCL::uyaxdv( 'L', 10.0, 30.0 )
+DCL::uyaxdv( 'R', 10.0, 30.0 )
+DCL::uysttl( 'L', 'LATITUDE', 0.0 )
+
+#-- X-Z PLANE ----
+
+DCL::sgswnd(  xmin,  xmax,  zmin,  zmax )
+DCL::sgsvpt( vxmin, vxmax, vzmin, vzmax )
+DCL::sgstrn( 1 )
+DCL::sgstrf
+
+DCL::scspln( 1, 3, vymax)
+DCL::scsprj
+
+DCL::uzinit
+DCL::uxaxdv( 'T', 10.0, 60.0 )
+DCL::uzlset( 'LABELXB', false )
+DCL::uxaxdv( 'B', 10.0, 60.0 )
+DCL::uyaxdv( 'L', 0.2, 1.0 )
+DCL::uyaxdv( 'R', 0.2, 1.0 )
+DCL::uysttl( 'L', 'AMPLITUDE', 0.0 )
+
+#-- Y-Z PLANE ----
+
+DCL::sgswnd(  ymin,  ymax,  zmin,  zmax )
+DCL::sgsvpt( vymin, vymax, vzmin, vzmax )
+DCL::sgstrn( 1 )
+DCL::sgstrf
+
+DCL::scspln( 2, 3, vxmax)
+DCL::scsprj
+
+DCL::uzinit
+DCL::uxaxdv( 'T', 10.0, 30.0 )
+DCL::uxaxdv( 'B', 10.0, 30.0 )
+DCL::uzlset( 'LABELYL', false )
+DCL::uyaxdv( 'L', 0.2, 1.0 )
+DCL::uyaxdv( 'R', 0.2, 1.0 )
+
+#---------------- 3-D ------------------
+
+DCL::scsvpt(vxmin, vxmax, vymin, vymax, vzmin, vzmax)
+DCL::scswnd( xmin,  xmax,  ymin,  ymax,  zmin,  zmax)
+DCL::scslog(false, false, false)
+DCL::scstrn(1)
+DCL::scstrf
+
+(ny-2).downto(0){|j|
+  (nx-2).downto(0){|i|
+    xp = [ alon[i], alon[i],   alon[i+1],  alon[i+1], alon[i] ]
+    yp = [ alat[j], alat[j+1], alat[j+1],  alat[j],   alat[j] ]
+    zp = [ p[i,j],  p[i,j+1],  p[i+1,j+1], p[i+1,j],  p[i,j] ]
+    DCL::sctnzu(xp[0..2], yp[0..2], zp[0..2], 10999, 50999)
+    DCL::sctnzu(xp[2..4], yp[2..4], zp[2..4], 10999, 50999)
+    DCL::scplu(xp, yp, zp)
+ }
+}
+
+DCL::sgcls

data/demo/rakuraku/miss/miss1.rb

@@ -0,0 +1,76 @@
+#
+# miss1.rb
+#   $Id: miss1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include NMath
+
+
+nmax = 40
+xmin = 0.0
+xmax = 4*PI
+ymin = -1.0
+ymax = 1.0
+
+x = NArray.sfloat(nmax+1)
+y = NArray.sfloat(nmax+1)
+
+#-- data ----
+dt = xmax/nmax
+x = x.indgen * dt
+y = sin(x)
+
+n1 = nmax/4
+rmiss = DCL::glrget('RMISS')
+y[n1-1] = rmiss
+y[n1]   = rmiss
+y[n1+1] = rmiss
+
+n2 = n1*3
+y[n2-1] = rmiss
+y[n2+1] = rmiss
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+DCL::grfrm
+
+#-- default ----
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt( 0.1,  0.9,  0.7,  0.9)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::sgplu(x, y)
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt( 0.1,  0.9,  0.5,  0.7)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::sgspmt(3)
+DCL::sgpmu(x, y)
+
+#-- cover loss data ----
+DCL::gllset('LMISS', true)
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt( 0.1,  0.9,  0.3,  0.5)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::sgplu(x, y)
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt( 0.1,  0.9,  0.1,  0.3)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::sgpmu(x, y)
+
+DCL::grcls
+