ruby-dcl 1.6.2

data/demo/grph2/umpack/test02.rb

@@ -0,0 +1,64 @@
+#
+# test02.rb
+#
+# $Id: test02.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
+dp = 0.2
+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] = 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
+
+rmiss = DCL::glrget('RMISS')
+DCL::sglset('LSOFTF', false)
+
+DCL::grfrm
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.1, 0.9, 0.1, 0.9)
+DCL::grssim(0.4, 0.0, 0.0)
+DCL::grsmpl(165.0, 60.0, 0.0)
+DCL::grstxy(-180.0, 180.0, 0.0, 90.0)
+DCL::grstrn(30)
+DCL::grstrf
+DCL::sglset('LCLIP', true)
+
+DCL::uestlv(rmiss, -dp, 201)
+DCL::uestlv(dp, dp*2, 401)
+DCL::uestlv(dp*2, rmiss, 402)
+DCL::uetone(p)
+
+DCL::udgclb(p, dp)
+DCL::udcntr(p)
+
+DCL::umpmap('coast_world')
+DCL::umpglb
+
+DCL::grcls
+

data/demo/grph2/umpack/test03.rb

@@ -0,0 +1,84 @@
+#
+# test03.rb
+#
+# $Id: test03.rb,v 1.1.1.1 2011-02-23 07:21:27 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
+dz = 0.05
+fact = 10
+p = NArray.sfloat(nx, ny)
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+alon = NArray.sfloat(nx)
+alat = NArray.sfloat(ny)
+
+#-- data ---
+alon = alon.indgen(xmin, (xmax-xmin)/(nx-1).to_f)
+alat = alat.indgen(ymin, (ymax-ymin)/(ny-1).to_f)
+for j in 0..ny-1
+  for i in 0..nx-1
+    slat = sin(alat[j]*drad)
+    p[i,j] = cos(alon[i]*drad) * (1-slat**2) * sin(2*PI*slat) + dz
+  end
+end
+
+rmiss = DCL::glrget('RMISS')
+DCL::gllset('LMISS', true)
+
+for j in 0..ny-1
+  for i in 0..nx-1
+    if (j == 0 || j == ny-1)
+      u[i,j] = rmiss
+      v[i,j] = rmiss
+    else
+      u[i,j] = (p[i,j-1] - p[i,j+1]) * fact
+      v[i,j] = (p[DCL::imod(i+1,nx-1),j] - p[DCL::imod(i-1,nx-1),j]) * fact
+    end
+  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.1, 0.9, 0.1, 0.9)
+DCL::grssim(0.4, 0.0, 0.0)
+DCL::grsmpl(165.0, 60.0, 0.0)
+DCL::grstxy(-180.0, 180.0, 0.0, 90.0)
+DCL::grstrn(30)
+DCL::grstrf
+DCL::sglset('LCLIP', true)
+
+DCL::umlset('LGRIDMJ', false)
+DCL::umrset('DGRIDMN', 30.0)
+DCL::umpmap('coast_world')
+DCL::umpglb
+
+DCL::udcntr(p)
+
+for j in 0..ny-1
+  for i in 0..nx-1
+    if (!(u[i,j] == rmiss || v[i,j] == rmiss))
+      DCL::sglazu(alon[i],alat[j],alon[i]+u[i,j],alat[j]+v[i,j], 1, 3)
+    end
+  end
+end
+
+DCL::grcls
+

data/demo/grph2/umpack/test04.rb

@@ -0,0 +1,54 @@
+#
+# test04.rb
+#
+# $Id: test04.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+np = 14
+
+#-- data ---
+ctr = ['cyl','mer','mwd','hmr','ek6','ktd','con',
+       'coa','coc','bon','otg','pst','azm','aza']
+fct = NArray[0.12, 0.12, 0.14, 0.14, 0.14, 0.14, 0.11,
+               0.16, 0.12, 0.12, 0.40, 0.12, 0.12, 0.17]
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+DCL::sgrset('STLAT1', 45.0)
+DCL::sgrset('STLAT2', 30.0)
+
+for i in 0..np-1
+  DCL::sgfrm
+  DCL::sgssim(fct[i], 0.0, 0.0)
+  DCL::sgsmpl(0.0, 90.0, 0.0)
+  DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
+
+  if (ctr[i] == 'otg')
+    DCL::sgstxy(-180.0, 180.0, 0.0, 90.0)
+  else
+    DCL::sgstxy(-180.0, 180.0, -90.0, 90.0)
+  end
+  DCL::sgstrn(DCL::isgtrc(ctr[i]))
+  DCL::sgstrf
+
+  DCL::sglset('LCLIP', true)
+  DCL::slpwwr(1)
+  DCL::slpvpr(1)
+  cttl = DCL::sgtrnl(DCL::isgtrc(ctr[i]))
+  DCL::sgtxzr(0.5, 0.95, cttl, 0.03, 0, 0, 3)
+
+  DCL::umpmap('coast_world')
+  DCL::umpglb
+end
+
+DCL::sgcls
+

data/demo/grph2/umpack/test05.rb

@@ -0,0 +1,85 @@
+#
+# test05.rb
+#
+# $Id: test05.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+np = 14
+nx = 73
+ny = 37
+xmin = 0
+xmax = 360
+ymin = -90
+ymax = 90
+p = NArray.new(Float, nx*ny)
+
+#-- data ---
+ctr = ['cyl','mer','mwd','hmr','ek6','ktd','con',
+       'coa','coc','bon','otg','pst','azm','aza']
+fct = NArray[0.12, 0.12, 0.14, 0.14, 0.14, 0.14, 0.11,
+               0.16, 0.12, 0.12, 0.40, 0.12, 0.12, 0.17]
+
+i = 0
+File.foreach('t811231.dat') do |line|
+  line.split(" ").each do |data|
+    p[i] = data.to_f
+    i += 1
+  end
+end
+p.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::sgrset('STLAT1', 75.0)
+DCL::sgrset('STLAT2', 60.0)
+DCL::umlset('LGRIDMJ', false)
+DCL::umrset('DGRIDMN', 30.0)
+
+for i in 0..np-1
+  DCL::sgfrm
+  DCL::sgssim(fct[i], 0.0, 0.0)
+  DCL::sgsmpl(165.0, 60.0, 0.0)
+  DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
+  DCL::sgswnd(xmin, xmax, ymin, ymax)
+
+  if (ctr[i] == 'otg')
+    DCL::sgstxy(-180.0, 180.0, 0.0, 90.0)
+  else
+    DCL::sgstxy(-180.0, 180.0, -90.0, 90.0)
+  end
+  DCL::sgstrn(DCL::isgtrc(ctr[i]))
+  DCL::sgstrf
+
+  DCL::sglset('LCLIP', true)
+  DCL::slpwwr(1)
+  DCL::slpvpr(1)
+  cttl = DCL::sgtrnl(DCL::isgtrc(ctr[i]))
+  DCL::sgtxzr(0.5, 0.95, cttl, 0.03, 0, 0, 3)
+
+  DCL::uetone(p)
+  DCL::udcntr(p)
+  DCL::umpmap('coast_world')
+  DCL::umpglb
+  DCL::slpvpr(1)
+end
+
+DCL::sgcls
+

data/demo/grph2/umpack/test06.rb

@@ -0,0 +1,44 @@
+#
+# test06.rb
+#
+# $Id: test06.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+np = 14
+
+#-- data ---
+ctr = ['cyl','mer','mwd','hmr','ek6','ktd','con',
+       'coa','coc','bon','otg','pst','azm','aza']
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::umlset('LGLOBE', true)
+
+for i in 0..np-1
+  DCL::grfrm
+  DCL::grsvpt(0.1, 0.9, 0.1, 0.9)
+  DCL::grstrn(DCL::isgtrc(ctr[i]))
+  DCL::umpfit
+  DCL::grstrf
+
+  DCL::sglset('LCLIP', true)
+  DCL::slpwwr(1)
+  DCL::slpvpr(1)
+  cttl =DCL:: sgtrnl(DCL::isgtrc(ctr[i]))
+  DCL::sgtxzr(0.5, 0.95, cttl, 0.03, 0, 0, 3)
+
+  DCL::umpmap('coast_world')
+  DCL::umpglb
+end
+
+DCL::grcls
+

data/demo/grph2/umpack/test07.rb

@@ -0,0 +1,43 @@
+#
+# test07.rb
+#
+# $Id: test07.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+np = 14
+
+#-- data ---
+ctr = ['cyl','mer','mwd','hmr','ek6','ktd','con',
+       'coa','coc','bon','otg','pst','azm','aza']
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+for i in 0..np-1
+  DCL::grfrm
+  DCL::grswnd(123.0, 147.0, 30.0, 46.0)
+  DCL::grsvpt(0.1, 0.9, 0.1, 0.9)
+  DCL::grstrn(DCL::isgtrc(ctr[i]))
+  DCL::umpfit
+  DCL::grstrf
+
+  DCL::sglset('LCLIP', true)
+  DCL::slpwwr(1)
+  DCL::slpvpr(1)
+  cttl = DCL::sgtrnl(DCL::isgtrc(ctr[i]))
+  DCL::sgtxzr(0.5, 0.95, cttl, 0.03, 0, 0, 3)
+
+  DCL::umpmap('coast_japan')
+  DCL::umpglb
+end
+
+DCL::grcls
+

data/demo/grph2/umpack/test08.rb

@@ -0,0 +1,37 @@
+#
+# test08.rb
+#
+# $Id: test08.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('LFULL', true)
+DCL::gropn iws
+
+DCL::slrat(1.0, 0.6)
+DCL::grfrm
+
+DCL::grsmpl(140.0, 90.0, 0.0)
+DCL::grstrn(DCL::isgtrc('HMR'))
+DCL::umpfit
+DCL::grstrf
+DCL::umpglb
+ 
+DCL::umiset('INDEXOUT', 3)
+DCL::umpmap('coast_world')
+DCL::umiset('INDEXOUT', 30)
+DCL::umpmap('border_world')
+DCL::umiset('INDEXOUT', 20)
+DCL::umpmap('plate_world')
+
+DCL::grcls
+

data/demo/grph2/umpack/test09.rb

@@ -0,0 +1,86 @@
+#
+# test09.rb
+#
+# T Horinouchi  2004/11/11
+
+require "numru/dcl"
+include NumRu
+
+raise "This program is for DCL 5.3 or later" unless DCL::DCLVERSION >= '5.3'
+
+nx=73
+ny=37
+xmin=0.0
+xmax=360.0
+ymin=-90.0
+ymax=90.0
+
+ctr = [ 'CYL','MER','MWD','HMR','EK6','KTD','CON','COA','COC','BON',
+        'OTG','PST','AZM','AZA']
+fct = [ 0.12, 0.12, 0.14, 0.14, 0.14, 0.14,0.11, 0.16, 0.12, 0.12,
+        0.40, 0.12, 0.12, 0.17 ]
+np = ctr.length  # == fct.length
+
+p = NArray.sfloat(nx*ny)
+i = 0
+File.foreach('t811231.dat') do |line|
+  line.split(" ").each do |data|
+    p[i] = data.to_f
+    i += 1
+  end
+end
+p.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.shift || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sglset('LSOFTF', false)
+DCL::sgrset('STLAT1', 75.0)
+DCL::sgrset('STLAT2', 60.0)
+DCL::umlset('LGRIDMJ', false)
+DCL::umrset('DGRIDMN', 30.0)
+
+nn = DCL::sgqcmn
+
+for i in 0..np-1
+  icn= (i % nn) + 1
+  print "colormap ",icn,": ",DCL::swqcmd(icn),"\n"
+  DCL::sgscmn(icn)
+  DCL::sgopn(iws)
+  DCL::sgfrm
+  DCL::sgssim(fct[i], 0.0, 0.0)
+  DCL::sgsmpl(165.0, 60.0, 0.0)
+  DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
+  DCL::sgswnd(xmin, xmax, ymin, ymax)
+  if (ctr[i] == 'OTG')
+    DCL::sgstxy(-180.0, 180.0, 0.0, 90.0)
+  else
+    DCL::sgstxy(-180.0, 180.0, -90.0, 90.0)
+  end
+  DCL::sgstrn(DCL::isgtrc(ctr[i]))
+  DCL::sgstrf
+
+  DCL::sglset('LCLIP', true)
+  DCL::slpwwr(1)
+  DCL::slpvpr(1)
+  cttl = DCL::sgtrnl(DCL::isgtrc(ctr[i]))
+  DCL::sgtxzr(0.5, 0.95, cttl, 0.03, 0, 0, 3)
+
+  DCL::uetone(p)
+  DCL::udcntr(p)
+  DCL::umfmap('coast_world')
+  DCL::umpmap('coast_world')
+  DCL::umpglb
+  DCL::slpvpr(1)
+  DCL::swpcls
+end
+
+DCL::sgcls