ruby-dcl 1.6.2

data/demo/gokuraku/step/step2.rb

@@ -0,0 +1,62 @@
+#
+# step2.rb
+#   $Id: step2.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nmax = 100
+
+x = NArray.sfloat(nmax)
+y = NArray.sfloat( nmax)
+
+#-- data ----
+iseed = 0
+x[0] = 2*(rand(iseed)-0.5)
+y[nmax-1] = x[0]
+
+for i in 1..nmax-1
+  x[i] = 2*(rand(iseed)-0.5)
+  y[i-1] = x[i]
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::usspnt(x, y)
+DCL::uspfit
+DCL::grstrf
+
+DCL::ussttl( 'X-TITLE', 'x-unit', 'Y-TITLE', 'y-unit' )
+DCL::usdaxs
+
+cnt = nmax/4
+st = 0
+ed = cnt-1
+DCL::uusmki(5)
+DCL::uusmks(0.015)
+DCL::uumrk(x[st..ed], y[st..ed])
+
+st += cnt
+ed += cnt
+DCL::uusmkt(2)
+DCL::uumrk(x[st..ed], y[st..ed])
+
+st += cnt
+ed += cnt
+DCL::uusmkt(3)
+DCL::uumrk(x[st..ed], y[st..ed])
+
+st += cnt
+ed += cnt
+DCL::uusmkt(4)
+DCL::uumrk(x[st..ed], y[st..ed])
+
+DCL::grcls
+

data/demo/gokuraku/u1d/u1d1.rb

@@ -0,0 +1,37 @@
+#
+# u1d1.rb
+#   $Id: u1d1.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nmax = 51
+xmin = 1945
+xmax = 1995
+
+y = NArray.sfloat(nmax)
+
+#-- data ----
+y0 = 0.5
+for n in 0..nmax-1
+  y[n] = 5.0*y0 + 10.0
+  y0 = 3.7*y0*(1.0-y0)
+end
+
+rundef = DCL::glrget('RUNDEF')
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, rundef, rundef)
+DCL::ussttl('TIME', 'YEAR', 'TEMPERATURE', 'DEG')
+DCL::usgrph(rundef, y)
+
+DCL::grcls
+

data/demo/gokuraku/u1d/u1d2.rb

@@ -0,0 +1,50 @@
+#
+# u1d2.rb
+#   $Id: u1d2.rb,v 1.1.1.1 2011-02-23 07:21:27 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
+
+rundef = DCL::glrget('RUNDEF')
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+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
+

data/demo/gokuraku/u2d/u2d1.rb

@@ -0,0 +1,46 @@
+#
+# u2d1.rb
+#   $Id: u2d1.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include NMath
+
+nt = 51
+nz = 21
+tmin = 0.0
+tmax = 5.0
+zmin = 20.0
+zmax = 50.0
+dt = (tmax-tmin)/(nt-1)
+dz = (zmax-zmin)/(nz-1)
+
+u = NArray.sfloat(nt, nz)
+
+#-- data ----
+zj = NArray.sfloat(nz).indgen! * dz
+uz = exp(-0.2*zj)*(zj**0.5)
+tz = -2.0*exp(-0.1*zj)
+for i in 0..nt-1
+  u[i,true] = uz*sin(3.0*(tz+i*dt))
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+DCL::grfrm
+
+DCL::grswnd(tmin, tmax, zmin, zmax)
+DCL::uspfit
+DCL::grstrf
+
+DCL::ussttl('TIME', 'YEAR', 'HEIGHT', 'km')
+DCL::usdaxs
+
+DCL::udcntr(u)
+
+DCL::grcls
+

data/demo/gokuraku/u2d/u2d2.rb

@@ -0,0 +1,49 @@
+#
+# u2d2.rb
+#   $Id: u2d2.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include NMath
+
+nt = 51
+nz = 21
+tmin = 0.0
+tmax = 5.0
+zmin = 20.0
+zmax = 50.0
+dt = (tmax-tmin)/(nt-1)
+dz = (zmax-zmin)/(nz-1)
+
+u = NArray.sfloat(nt, nz)
+
+#-- data ----
+zj = NArray.sfloat(nz).indgen! * dz
+uz = exp(-0.2*zj)*(zj**0.5)
+tz = -2.0*exp(-0.1*zj)
+for i in 0..nt-1
+  u[i,true] = uz*sin(3.0*(tz+i*dt))
+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(tmin, tmax, zmin, zmax)
+DCL::grsvpt(0.2, 0.9, 0.4, 0.8)
+DCL::uspfit
+DCL::grstrf
+
+DCL::ussttl('TIME', 'YEAR', 'HEIGHT', 'km')
+DCL::usdaxs
+
+DCL::udcntr(u)
+DCL::uetone(u)
+
+DCL::grcls
+

data/demo/gokuraku/u2d/u2d3.rb

@@ -0,0 +1,45 @@
+#
+# u2d3.rb
+#   $Id: u2d3.rb,v 1.1.1.1 2011-02-23 07:21:27 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
+dx = (xmax-xmin)/(nx-1)
+dy = (ymax-ymin)/(ny-1)
+
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+
+#-- data ----
+xi = xmin + NArray.sfloat(nx).indgen! * dx
+yj = ymin + NArray.sfloat(ny).indgen! * dy
+for j in 0..ny-1; u[true,j] = xi; end
+for i in 0..nx-1; v[i,true] = -yj; 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::uspfit
+DCL::grstrf
+
+DCL::ussttl('X', 'km', 'Y', 'km')
+DCL::usdaxs
+
+DCL::ugvect(u, v)
+
+DCL::grcls
+

data/demo/gokuraku/u2d/u2d4.rb

@@ -0,0 +1,61 @@
+#
+# u2d4.rb
+#   $Id: u2d4.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include NMath
+
+nx = 21
+ny = 21
+xmin = -10.0
+xmax = 10.0
+ymin = -10.0
+ymax = 10.0
+kmax = 5
+pmin = 0.0
+pmax = 1.0
+dx = (xmax-xmin)/(nx-1)
+dy = (ymax-ymin)/(ny-1)
+dp = (pmax-pmin)/kmax
+
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+p = NArray.sfloat(nx, ny)
+
+
+#-- data ----
+xi = xmin + NArray.sfloat(nx).indgen! * dx
+yj = ymin + NArray.sfloat(ny).indgen! * dy
+for j in 0..ny-1; u[true,j] = xi; end
+for i in 0..nx-1; v[i,true] = -yj; end
+p = exp(-xi.reshape(nx,1)**2/64 - yj.reshape(1,ny)**2/25)
+
+#-- 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::uspfit
+DCL::grstrf
+
+DCL::usdaxs
+DCL::ugvect(u, v)
+DCL::udcntr(p)
+
+for k in 0..kmax-1
+  tlev1 = k*dp
+  tlev2 = tlev1 + dp
+  ipat  = 600 + k
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+
+DCL::uetone(p)
+
+DCL::grcls
+

data/demo/grph1/scpack/scpkt2.rb

@@ -0,0 +1,40 @@
+#
+# scpkt2.rb
+#
+# $Id: scpkt2.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::sglset('L2TO3', true)
+DCL::sgiset('IFONT', 2)
+
+DCL::sgopn iws
+
+x0 = -3
+dx = 0.3
+for i in 1..20
+  DCL::sgfrm
+  DCL::sgswnd(0.0, 1.0, 0.0, 1.0)
+  DCL::sgsvpt(0.0, 1.0, 0.0, 1.0)
+  DCL::sgstrn(1)
+  DCL::sgstrf
+
+  x = x0 + dx*i
+  DCL::sgrset('XEYE3', x)
+  DCL::scsprj
+
+  DCL::slpvpr(1)
+  DCL::sgtxzv(0.5, 0.5, 'DENNOU', 0.1, 0, 0, 3)
+end
+
+DCL::sgcls
+

data/demo/grph1/scpack/scpkt3.rb

@@ -0,0 +1,102 @@
+#
+# scpkt3.rb
+#
+# $Id: scpkt3.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nx = 73
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+
+p1 = NArray.sfloat(nx*ny)
+p2 = NArray.sfloat(nx*ny)
+
+#-- data ---
+i = 0
+File.foreach("t810630.dat") do |line|
+  line.split(" ").each do |data|
+    p1[i] = data.to_f
+    i += 1
+  end
+end
+i = 0
+File.foreach("t811231.dat") do |line|
+  line.split(" ").each do |data|
+    p2[i] = data.to_f
+    i += 1
+  end
+end
+p1.reshape!(nx, ny)
+p2.reshape!(nx, ny)
+
+190.step(245, 5) do |r|
+  amin=r
+  amax=r+5
+#       IDX=(R-180)*1.4*1000+999
+  idx=((r-170)*1.25).to_i*1000+999
+  DCL::uestlv(amin, amax, idx)
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sglset('LSOFTF', false)
+DCL::umlset('LGRIDMJ', false)
+DCL::umrset('DGRIDMN', 30.0)
+DCL::umiset('ITYPEMN', 1)
+DCL::umiset('ITYPEOUT', 1)
+DCL::sgfrm
+
+DCL::sglset('L2TO3', true)
+#     CALL SCSEYE(-1.0, -1.5,  2.0 )
+DCL::scsobj(0.5, 0.5, 0.0)
+DCL::scspln(1, 2, -0.3)
+DCL::scsprj
+
+DCL::sgssim(0.12, 0.0, 0.0)
+DCL::sgsmpl(180.0, 90.0, 0.0)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
+DCL::sgswnd(xmin, xmax, ymin, ymax)
+DCL::sgstxy(-180.0, 180.0, -90.0, 90.0)
+DCL::sgstrn(DCL::isgtrc('MER'))
+DCL::sgstrf
+
+DCL::sglset('LCLIP', true)
+DCL::slpwwr(1)
+
+DCL::uetone(p1)
+DCL::umpmap('coast_world')
+DCL::umpglb
+DCL::slpvpr(1)
+
+DCL::sglset('LCLIP', false)
+DCL::sgtxzv(0.5, 0.95, 'Temp. : 10 hPa : Jun 30, 1981', 0.025, 0, 0, 3)
+
+DCL::scspln(1, 2, 0.4)
+#     CALL SGRSET( 'YOFF3', -0.25 )
+DCL::scsprj
+
+DCL::sglset('LCLIP', true)
+DCL::slpwwr(1)
+
+DCL::uetone(p2)
+DCL::umpmap('coast_world')
+DCL::umpglb
+DCL::slpvpr(1)
+
+DCL::sglset('LCLIP', false)
+DCL::sgtxzv(0.5, 0.95, 'Temp. : 10 hPa : Dec 31, 1981', 0.025, 0, 0, 3)
+
+DCL::sgcls
+