ruby-dcl 1.6.2

data/demo/grph2/udegpk/u2df04.rb

@@ -0,0 +1,52 @@
+#
+# u2df04.rb
+#
+# $Id: u2df04.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::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::usdaxs
+
+DCL::udcntr(p)
+DCL::uetone(p)
+
+DCL::grcls
+

data/demo/grph2/udegpk/u2df05.rb

@@ -0,0 +1,70 @@
+#
+# u2df05.rb
+#
+# $Id: u2df05.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 = 20
+dx2 = 60
+ymin = -90
+ymax = +90
+dy1 = 10
+dy2 = 30
+drad = PI/180
+dz = 0.05
+dp = 0.2
+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
+
+rmiss = DCL::glrget('RMISS')
+DCL::sglset('LSOFTF', false)
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uestlv(rmiss, -dp, 201)
+DCL::uestlv(dp, dp*2, 401)
+DCL::uestlv(dp*2, rmiss, 402)
+DCL::uetone(p)
+
+DCL::uxaxdv('B', dx1, dx2)
+DCL::uxaxdv('T', dx1, dx2)
+DCL::uxsttl('B', 'LONGITUDE', 0.0)
+
+DCL::uyaxdv('L', dy1, dy2)
+DCL::uyaxdv('R', dy1, dy2)
+DCL::uysttl('L', 'LATITUDE', 0.0)
+
+DCL::udgclb(p, dp)
+DCL::udcntr(p)
+
+DCL::grcls
+

data/demo/grph2/udegpk/u2df06.rb

@@ -0,0 +1,59 @@
+#
+# u2df06.rb
+#
+# $Id: u2df06.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 ---
+rmiss = DCL::glrget('RMISS')
+DCL::gllset('LMISS', true)
+
+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)
+    if ((i == 2 && j == 14) || ((8 <= i && i <= 11) && (6 <= j && j <= 10)))
+      p[i,j]=rmiss
+    else
+      p[i,j] = cos(alon) * (1-slat**2) * sin(2*PI*slat) + dz
+    end
+  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::usdaxs
+
+DCL::udcntr(p)
+DCL::uetone(p)
+
+DCL::grcls
+

data/demo/grph2/udegpk/u2df07.rb

@@ -0,0 +1,68 @@
+#
+# u2df07.rb
+#
+# $Id: u2df07.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)
+u = NArray.sfloat(nx, ny)
+v = 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
+
+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]
+      v[i,j] = p[DCL::imod(i+1,nx-1),j] - p[DCL::imod(i-1,nx-1),j]
+    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.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::usdaxs
+
+DCL::ugvect(u, v)
+
+DCL::grcls
+

data/demo/grph2/udegpk/u2df08.rb

@@ -0,0 +1,91 @@
+#
+# u2df08.rb
+#
+# $Id: u2df08.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 = 20
+dx2 = 60
+ymin = -90
+ymax = +90
+dy1 = 10
+dy2 = 30
+drad = PI/180
+dz = 0.05
+
+p = NArray.sfloat(nx, ny)
+u = NArray.sfloat(nx, ny)
+v = 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
+
+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]
+      v[i,j] = p[DCL::imod(i+1,nx-1),j] - p[DCL::imod(i-1,nx-1),j]
+    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.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::uyaxdv('L', dy1, dy2)
+DCL::uyaxdv('R', dy1, dy2)
+DCL::uysttl('L', 'LATITUDE', 0.0)
+
+DCL::udrset('RSIZEL', 0.015)
+DCL::udrset('RSIZET', 0.015)
+DCL::udcntr(p)
+
+DCL::ugrset('RSIZET', 0.015)
+DCL::uglset('LNRMAL', false)
+DCL::ugrset('XFACT1', 0.025)
+DCL::ugrset('YFACT1', 0.050)
+DCL::uglset('LUNIT', true)
+DCL::ugsut('X', 'U')
+DCL::ugsut('Y', 'V')
+DCL::ugrset('VXUNIT', 0.05)
+DCL::ugrset('VYUNIT', 0.10)
+DCL::ugvect(u, v)
+
+DCL::grcls
+

data/demo/grph2/udegpk/u2df09b.rb

@@ -0,0 +1,88 @@
+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)
+
+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::uetonb(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::udlset('LMSG', false)
+DCL::udcntr(pi)
+DCL::uetonb(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/grph2/udegpk/u2df09c.rb

@@ -0,0 +1,89 @@
+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::sgscmn(5)
+DCL::gropn iws
+
+DCL::sglset('LSOFTF', false)
+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)
+
+dp = (pmax-pmin)/kmax.to_f
+for k in 1..kmax
+  tlev1 = (k-1)*dp
+  tlev2 = tlev1 + dp
+  ipat  = 999 + (99 - k*15)*1000
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+DCL::uetonc(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::udlset('LMSG', false)
+DCL::udcntr(pi)
+DCL::uetonc(pi)
+
+DCL::slpvpr(3)
+DCL::uzlset('LABELYR', true)
+DCL::uzfact(0.8)
+DCL::uysfmt('(F4.1)')
+DCL::uyaxdv('R', dp, dp)
+
+DCL::grcls