ruby-dcl 1.6.2

data/demo/rakuraku/u2d/u2d4.rb

@@ -0,0 +1,59 @@
+#
+# u2d4.rb
+#   $Id: u2d4.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
+kxmn = 9
+kxmx = 33
+kymn = 7
+kymx = 31
+
+p = 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
+  slat = sin(alat[j]*drad)
+  for i in 0..nx-1
+    p[i,j] = 3*sqrt(1-slat**2)*slat*cos(alon[i]*drad) - 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(alon[kxmn-1], alon[kxmx-1], alat[kymn-1], alat[kymx-1])
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::usdaxs
+DCL::udcntr(p[(kxmn-1)..(kxmx-1), (kymn-1)..(kymx-1)])
+
+DCL::grcls
+

data/demo/rakuraku/u2d/u2d5.rb

@@ -0,0 +1,61 @@
+#
+# u2d5.rb
+#   $Id: u2d5.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.new(Float, 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', true)
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+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::usdaxs
+DCL::udgclb(p, 0.4)
+DCL::udcntr(p)
+DCL::uetone(p)
+
+DCL::grcls
+

data/demo/rakuraku/u2d/u2d6.rb

@@ -0,0 +1,57 @@
+#
+# u2d6.rb
+#   $Id: u2d6.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nx = 11
+ny = 11
+xmin = -1
+xmax = 1
+ymin = -1
+ymax = 1
+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).to_f
+    y = ymin + (ymax-ymin)*j/(ny-1).to_f
+    u[i,j] =  x * 0.1
+    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::uglset('LNRMAL', false)
+DCL::ugrset('XFACT1', 0.5)
+DCL::ugrset('YFACT1', 0.05)
+DCL::uglset('LUNIT', true)
+DCL::ugrset('VXUNIT', 0.1)
+DCL::ugrset('VYUNIT', 0.1)
+DCL::ugrset('VXUOFF', 0.06)
+DCL::ugsut('X', 'U')
+DCL::ugsut('Y', 'V')
+DCL::ugvect(u, v)
+
+DCL::grcls
+

data/demo/rakuraku/u2d/u2d7.rb

@@ -0,0 +1,103 @@
+#
+# u2d7.rb
+#   $Id: u2d7.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
+dx1 = 1
+dx2 = 5
+dy1 = 1
+dy2 = 5
+kmax = 5
+pmin = 0
+pmax = 1
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+p = NArray.sfloat(nx, ny)
+pi = NArray.sfloat(2, kmax+1)
+
+#-- 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::gropn iws
+
+DCL::sglset('LSOFTF', true)
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+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', 'X-axis', 0.0)
+DCL::uxmttl('T', 'FIGURE TITLE', 0.0)
+
+DCL::uyaxdv('L', dy1, dy2)
+DCL::uyaxdv('R', dy1, dy2)
+DCL::uysttl('L', 'Y-axis', 0.0)
+
+DCL::uglset('LUNIT', true)
+DCL::ugrset('VXUOFF', 0.04)
+DCL::ugsut('X', 'U')
+DCL::ugsut('Y', 'V')
+DCL::ugvect(u, v)
+
+DCL::udlset('LABEL', false)
+DCL::udlset('LMSG', false)
+DCL::udgclb(p, 0.1)
+DCL::udcntr(p)
+
+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::grswnd(0.0, 1.0, pmin, pmax)
+DCL::grsvpt(0.85, 0.9, 0.45, 0.75)
+DCL::grstrn(1)
+DCL::grstrf
+
+for k in 1..kmax+1
+  pi[0,k-1] = pmin + (k-1)*dp
+  pi[1,k-1] = pmin + (k-1)*dp
+end
+
+DCL::udcntr(pi)
+DCL::uetone(pi)
+
+DCL::slpvpr(3)
+DCL::uzlset('LABELYR', true)
+DCL::uzfact(0.8)
+DCL::uysfmt('(F4.1)')
+DCL::uyaxdv('R', dp, dp)
+
+DCL::grcls

data/demo/rakuraku/uspack/uspac1.rb

@@ -0,0 +1,35 @@
+#
+# uspac1.rb
+#   $Id: uspac1.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 ---
+x.indgen(1950.0, 1.0) 
+r    = 3.7
+y[0] = 0.5
+for n in 0..nmax-1
+  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::ussttl('TIME', 'YEAR', 'HEAT FLUX', 'W/m|2"')
+DCL::usgrph(x, y)
+
+DCL::grcls
+

data/demo/rakuraku/uspack/uspac2.rb

@@ -0,0 +1,66 @@
+#
+# uspac2.rb
+#   $Id: uspac2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nmax = 201
+imax = 5
+x = NArray.sfloat(nmax)
+y0 = NArray.sfloat(nmax)
+y1 = NArray.sfloat(nmax)
+y2 = NArray.sfloat(nmax)
+a = NArray.sfloat(imax)
+
+#-- data ---
+for i in 0..imax-1
+  ii = 2*i + 1
+  a[i] = ((-1.0)**(i+1))*2/(ii*PI)
+end
+
+x = x.indgen * 1.0/(nmax-1)
+for n in 0..nmax-1
+  t     = 2.0*PI*x[n]
+  if (t < PI/2.0 || t >= PI*3.0/2.0)
+    y0[n] = 0.0
+  else
+    y0[n] = 1.0
+  end
+  y1[n] = 0.5 + a[0]*cos(t) 
+  y2[n] = 0.5
+  for i in 0..imax-1
+    ii = 2*i + 1
+    y2[n] = y2[n] + a[i]*cos(ii*t)
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::usspnt(x, y0)
+DCL::usspnt(x, y1)
+DCL::usspnt(x, y2)
+DCL::uspfit
+DCL::grstrf
+
+DCL::ussttl('FREQUENCY', '/DAY', 'RESPONSE', ' ')
+DCL::usdaxs
+
+DCL::uulin(x, y0)
+DCL::uuslnt(2)
+DCL::uuslni(3)
+DCL::uulin(x, y1)
+DCL::uuslnt(3)
+DCL::uulin(x, y2)
+
+DCL::grcls
+

data/demo/rakuraku/uspack/uspac3.rb

@@ -0,0 +1,37 @@
+#
+# uspac3.rb
+#   $Id: uspac3.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
+xmin = 1950
+xmax = 2000
+y = NArray.sfloat(nmax+1)
+
+#-- data ---
+r    = 3.7
+y[0] = 0.5
+for n in 0..nmax-1
+  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
+
+rundef = DCL::glrget('RUNDEF')
+DCL::grswnd(xmin, xmax, rundef, rundef)
+DCL::ussttl('TIME', 'YEAR', 'HEAT FLUX', 'W/m|2"')
+DCL::usgrph(rundef, y)
+
+DCL::grcls
+

data/demo/rakuraku/uspack/uspac4.rb

@@ -0,0 +1,49 @@
+#
+# uspac4.rb
+#   $Id: uspac4.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
+ymin = 0.0
+ymax = 50.0
+x1 = NArray.sfloat(nmax+1)
+x2 = NArray.sfloat(nmax+1)
+
+#-- data ---
+iseed = 0
+for n in 0..nmax
+  y = ymin + (ymax-ymin)*n/nmax
+  x1[n] = 10.0*(exp(-y/20))**2 * exp((rand(iseed)-0.5)*2)**2
+  x2[n] = 10.0*(exp(-y/20))**2
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+rundef = DCL::glrget('RUNDEF')
+DCL::grswnd(rundef, rundef, ymin, ymax)
+DCL::usspnt(x1, rundef)
+DCL::usspnt(x2, rundef)
+DCL::grstrn(3)
+DCL::uspfit
+DCL::grstrf
+
+DCL::ussttl('MIXING RATIO', 'ppmv', 'HEIGHT', 'km')
+DCL::usdaxs
+
+DCL::uusmkt(2)
+DCL::uumrk(x1, rundef)
+DCL::uulin(x2, rundef)
+
+DCL::grcls
+