ruby-dcl 1.6.2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (319) hide show
  1. data/ChangeLog +321 -0
  2. data/GenWrapper/Makefile +14 -0
  3. data/GenWrapper/cproto2init.rb +41 -0
  4. data/GenWrapper/dcl_narrayed_funcs.rb +181 -0
  5. data/GenWrapper/dcl_rb_footing +29 -0
  6. data/GenWrapper/dcl_rb_heading +129 -0
  7. data/GenWrapper/def.rb +20 -0
  8. data/GenWrapper/elim_ary_size.rb +398 -0
  9. data/GenWrapper/etc/Mk_proto +375 -0
  10. data/GenWrapper/etc/dcl_gen +13 -0
  11. data/GenWrapper/etc/p_header +63 -0
  12. data/GenWrapper/etc/p_init +410 -0
  13. data/GenWrapper/pparse.rb +137 -0
  14. data/GenWrapper/proto/Makefile +35 -0
  15. data/GenWrapper/proto/grph1_00 +17 -0
  16. data/GenWrapper/proto/grph1_csgi.fp +27 -0
  17. data/GenWrapper/proto/grph1_scpack.fp +442 -0
  18. data/GenWrapper/proto/grph1_sgpack.fp +1406 -0
  19. data/GenWrapper/proto/grph1_slpack.fp +316 -0
  20. data/GenWrapper/proto/grph1_stpack.fp +411 -0
  21. data/GenWrapper/proto/grph1_swpack.fp +746 -0
  22. data/GenWrapper/proto/grph2_00 +16 -0
  23. data/GenWrapper/proto/grph2_grpack.fp +139 -0
  24. data/GenWrapper/proto/grph2_ucpack.fp +522 -0
  25. data/GenWrapper/proto/grph2_udpack.fp +692 -0
  26. data/GenWrapper/proto/grph2_uepack.fp +601 -0
  27. data/GenWrapper/proto/grph2_ugpack.fp +436 -0
  28. data/GenWrapper/proto/grph2_uhpack.fp +246 -0
  29. data/GenWrapper/proto/grph2_ulpack.fp +516 -0
  30. data/GenWrapper/proto/grph2_umpack.fp +676 -0
  31. data/GenWrapper/proto/grph2_uspack.fp +989 -0
  32. data/GenWrapper/proto/grph2_uupack.fp +663 -0
  33. data/GenWrapper/proto/grph2_uvpack.fp +246 -0
  34. data/GenWrapper/proto/grph2_uwpack.fp +203 -0
  35. data/GenWrapper/proto/grph2_uxpack.fp +295 -0
  36. data/GenWrapper/proto/grph2_uypack.fp +295 -0
  37. data/GenWrapper/proto/grph2_uzpack.fp +574 -0
  38. data/GenWrapper/proto/math1_00 +12 -0
  39. data/GenWrapper/proto/math1_blklib.fp +58 -0
  40. data/GenWrapper/proto/math1_chrlib.fp +83 -0
  41. data/GenWrapper/proto/math1_fnclib.fp +79 -0
  42. data/GenWrapper/proto/math1_gnmlib.fp +118 -0
  43. data/GenWrapper/proto/math1_gt2dlib.fp +144 -0
  44. data/GenWrapper/proto/math1_ifalib.fp +123 -0
  45. data/GenWrapper/proto/math1_indxlib.fp +222 -0
  46. data/GenWrapper/proto/math1_intlib.fp +46 -0
  47. data/GenWrapper/proto/math1_lrllib.fp +276 -0
  48. data/GenWrapper/proto/math1_maplib.fp +24 -0
  49. data/GenWrapper/proto/math1_oslib.fp +28 -0
  50. data/GenWrapper/proto/math1_rfalib.fp +420 -0
  51. data/GenWrapper/proto/math1_rfblib.fp +51 -0
  52. data/GenWrapper/proto/math1_sublib.fp +60 -0
  53. data/GenWrapper/proto/math1_syslib.fp +767 -0
  54. data/GenWrapper/proto/math1_vialib.fp +339 -0
  55. data/GenWrapper/proto/math1_viblib.fp +264 -0
  56. data/GenWrapper/proto/math1_vralib.fp +339 -0
  57. data/GenWrapper/proto/math1_vrblib.fp +264 -0
  58. data/GenWrapper/proto/math1_xfclib.fp +68 -0
  59. data/GenWrapper/proto/math2_00 +12 -0
  60. data/GenWrapper/proto/math2_fftlib.fp +803 -0
  61. data/GenWrapper/proto/math2_intrlib.fp +32 -0
  62. data/GenWrapper/proto/math2_odelib.fp_notused +568 -0
  63. data/GenWrapper/proto/math2_rnmlib.fp +54 -0
  64. data/GenWrapper/proto/math2_shtlib.fp +1292 -0
  65. data/GenWrapper/proto/math2_shtlib.fp_old +1294 -0
  66. data/GenWrapper/proto/math2_vstlib.fp +84 -0
  67. data/GenWrapper/proto/misc1_00 +17 -0
  68. data/GenWrapper/proto/misc1_chnlib.fp +47 -0
  69. data/GenWrapper/proto/misc1_datelib.fp +352 -0
  70. data/GenWrapper/proto/misc1_fmtlib.fp +19 -0
  71. data/GenWrapper/proto/misc1_misclib.fp +26 -0
  72. data/GenWrapper/proto/misc1_randlib.fp +36 -0
  73. data/GenWrapper/proto/misc1_timelib.fp +151 -0
  74. data/GenWrapper/proto2c.rb +13 -0
  75. data/GenWrapper/prototype.rb +209 -0
  76. data/GenWrapper/util.rb +9 -0
  77. data/GenWrapper/variable.rb +671 -0
  78. data/README +29 -0
  79. data/Rakefile +44 -0
  80. data/ToDo +2 -0
  81. data/dcl_cary2obj.c +354 -0
  82. data/dcl_narrayed_funcs.c +518 -0
  83. data/dcl_obj2cary.c +512 -0
  84. data/dcl_rubydcloriginal.c +129 -0
  85. data/demo/gokuraku/hop/hop.rb +26 -0
  86. data/demo/gokuraku/jump/jump1.rb +64 -0
  87. data/demo/gokuraku/jump/jump2.rb +50 -0
  88. data/demo/gokuraku/layout/lay1.rb +46 -0
  89. data/demo/gokuraku/layout/lay2.rb +44 -0
  90. data/demo/gokuraku/step/step0.rb +34 -0
  91. data/demo/gokuraku/step/step1.rb +84 -0
  92. data/demo/gokuraku/step/step2.rb +62 -0
  93. data/demo/gokuraku/u1d/u1d1.rb +37 -0
  94. data/demo/gokuraku/u1d/u1d2.rb +50 -0
  95. data/demo/gokuraku/u2d/u2d1.rb +46 -0
  96. data/demo/gokuraku/u2d/u2d2.rb +49 -0
  97. data/demo/gokuraku/u2d/u2d3.rb +45 -0
  98. data/demo/gokuraku/u2d/u2d4.rb +61 -0
  99. data/demo/grph1/scpack/scpkt2.rb +40 -0
  100. data/demo/grph1/scpack/scpkt3.rb +102 -0
  101. data/demo/grph1/scpack/scpkt6.rb +84 -0
  102. data/demo/grph1/scpack/scpkt7.rb +59 -0
  103. data/demo/grph1/scpack/scpkt8.rb +60 -0
  104. data/demo/grph1/scpack/t810630.dat +296 -0
  105. data/demo/grph1/scpack/t811231.dat +296 -0
  106. data/demo/grph1/sgpack/sgfont.rb +52 -0
  107. data/demo/grph1/sgpack/sgfonz.rb +53 -0
  108. data/demo/grph1/sgpack/sgksx1.rb +44 -0
  109. data/demo/grph1/sgpack/sgksx2.rb +90 -0
  110. data/demo/grph1/sgpack/sgksx3.rb +75 -0
  111. data/demo/grph1/sgpack/sglidx.rb +40 -0
  112. data/demo/grph1/sgpack/sgltyp.rb +48 -0
  113. data/demo/grph1/sgpack/sgpk01.rb +39 -0
  114. data/demo/grph1/sgpack/sgpk02.rb +79 -0
  115. data/demo/grph1/sgpack/sgpk03.rb +107 -0
  116. data/demo/grph1/sgpack/sgpk04.rb +86 -0
  117. data/demo/grph1/sgpack/sgpk05.rb +66 -0
  118. data/demo/grph1/sgpack/sgpk06.rb +91 -0
  119. data/demo/grph1/sgpack/sgpk07.rb +66 -0
  120. data/demo/grph1/sgpack/sgpk08.rb +67 -0
  121. data/demo/grph1/sgpack/sgpk09.rb +61 -0
  122. data/demo/grph1/sgpack/sgpk10.rb +58 -0
  123. data/demo/grph1/sgpack/sgtclr.rb +59 -0
  124. data/demo/grph1/sgpack/sgtone.rb +70 -0
  125. data/demo/grph1/sgpack/sgtonz.rb +66 -0
  126. data/demo/grph1/slpack/slpk01.rb +37 -0
  127. data/demo/grph1/slpack/slpk02.rb +29 -0
  128. data/demo/grph1/slpack/slpk03.rb +36 -0
  129. data/demo/grph1/slpack/slpk04.rb +36 -0
  130. data/demo/grph2/g2pack/g2pk01.rb +78 -0
  131. data/demo/grph2/g2pack/g2pk02.rb +75 -0
  132. data/demo/grph2/grpack/grpk01.rb +65 -0
  133. data/demo/grph2/ucpack/ucpk01.rb +95 -0
  134. data/demo/grph2/ucpack/ucpk02.rb +95 -0
  135. data/demo/grph2/udegpk/u2df01.rb +50 -0
  136. data/demo/grph2/udegpk/u2df02.rb +76 -0
  137. data/demo/grph2/udegpk/u2df03.rb +80 -0
  138. data/demo/grph2/udegpk/u2df04.rb +52 -0
  139. data/demo/grph2/udegpk/u2df05.rb +70 -0
  140. data/demo/grph2/udegpk/u2df06.rb +59 -0
  141. data/demo/grph2/udegpk/u2df07.rb +68 -0
  142. data/demo/grph2/udegpk/u2df08.rb +91 -0
  143. data/demo/grph2/udegpk/u2df09b.rb +88 -0
  144. data/demo/grph2/udegpk/u2df09c.rb +89 -0
  145. data/demo/grph2/udegpk/u2df09e.rb +88 -0
  146. data/demo/grph2/udegpk/u2df09f.rb +89 -0
  147. data/demo/grph2/ulpack/ulpk01.rb +111 -0
  148. data/demo/grph2/ulpack/ulpk01n.rb +111 -0
  149. data/demo/grph2/ulpack/ulpk02.rb +111 -0
  150. data/demo/grph2/ulpack/ulpk02n.rb +111 -0
  151. data/demo/grph2/umpack/t811231.dat +296 -0
  152. data/demo/grph2/umpack/test01.rb +69 -0
  153. data/demo/grph2/umpack/test02.rb +64 -0
  154. data/demo/grph2/umpack/test03.rb +84 -0
  155. data/demo/grph2/umpack/test04.rb +54 -0
  156. data/demo/grph2/umpack/test05.rb +85 -0
  157. data/demo/grph2/umpack/test06.rb +44 -0
  158. data/demo/grph2/umpack/test07.rb +43 -0
  159. data/demo/grph2/umpack/test08.rb +37 -0
  160. data/demo/grph2/umpack/test09.rb +86 -0
  161. data/demo/grph2/umpack/test10.rb +52 -0
  162. data/demo/grph2/umpack/umpk01.rb +65 -0
  163. data/demo/grph2/umpack/umpk02.rb +54 -0
  164. data/demo/grph2/umpack/umpk03.rb +46 -0
  165. data/demo/grph2/umpack/umpk04.rb +65 -0
  166. data/demo/grph2/umpack/umpk05.rb +84 -0
  167. data/demo/grph2/uspack/uspk01.rb +39 -0
  168. data/demo/grph2/uspack/uspk02.rb +47 -0
  169. data/demo/grph2/uspack/uspk03.rb +38 -0
  170. data/demo/grph2/uspack/uspk04.rb +64 -0
  171. data/demo/grph2/uspack/uspk05.rb +68 -0
  172. data/demo/grph2/uspack/uspk06.rb +43 -0
  173. data/demo/grph2/uspack/uspk07.rb +62 -0
  174. data/demo/grph2/uspack/uspk08.rb +46 -0
  175. data/demo/grph2/uspack/uspk09.rb +89 -0
  176. data/demo/grph2/uspack/uspk10.rb +48 -0
  177. data/demo/grph2/uspack/uspk11.rb +71 -0
  178. data/demo/grph2/uspack/uspk12.rb +51 -0
  179. data/demo/grph2/uupack/uupk01.rb +50 -0
  180. data/demo/grph2/uupack/uupk02.rb +74 -0
  181. data/demo/grph2/uupack/uupk03.rb +56 -0
  182. data/demo/grph2/uupack/uupk04.rb +84 -0
  183. data/demo/grph2/uupack/uupk05.rb +74 -0
  184. data/demo/grph2/uupack/uupk06.rb +77 -0
  185. data/demo/grph2/uupack/uupk07.rb +88 -0
  186. data/demo/grph2/uxyzpk/uxyz01.rb +38 -0
  187. data/demo/grph2/uxyzpk/uxyz02.rb +36 -0
  188. data/demo/grph2/uxyzpk/uxyz03.rb +37 -0
  189. data/demo/grph2/uxyzpk/uxyz04.rb +54 -0
  190. data/demo/grph2/uxyzpk/uxyz05.rb +41 -0
  191. data/demo/grph2/uxyzpk/uxyz06.rb +42 -0
  192. data/demo/grph2/uxyzpk/uxyz07.rb +38 -0
  193. data/demo/grph2/uxyzpk/uxyz08.rb +60 -0
  194. data/demo/grph2/uxyzpk/uxyz09.rb +52 -0
  195. data/demo/grph2/uxyzpk/uxyz10.rb +57 -0
  196. data/demo/grph2/ximage/MEMO +23 -0
  197. data/demo/grph2/ximage/tomsclm.dat +1344 -0
  198. data/demo/grph2/ximage/ximg01.rb +145 -0
  199. data/demo/math1/gt2dlib/gt2d01.rb +33 -0
  200. data/demo/math2/fftlib/fftl01.rb +28 -0
  201. data/demo/math2/fftlib/fftl02.rb +224 -0
  202. data/demo/math2/intrlib/intr01.rb +22 -0
  203. data/demo/math2/rnmlib/rnml01.rb +21 -0
  204. data/demo/math2/vstlib/vstl01.rb +29 -0
  205. data/demo/math2/vstlib/vstl02.rb +38 -0
  206. data/demo/rakuraku/color/color1.rb +61 -0
  207. data/demo/rakuraku/color/color2.rb +57 -0
  208. data/demo/rakuraku/color/color3.rb +81 -0
  209. data/demo/rakuraku/kihon/kihon1.rb +39 -0
  210. data/demo/rakuraku/kihon/kihon2.rb +56 -0
  211. data/demo/rakuraku/kihon/kihon3.rb +57 -0
  212. data/demo/rakuraku/kihon/kihon4.rb +70 -0
  213. data/demo/rakuraku/kihon/kihon5.rb +56 -0
  214. data/demo/rakuraku/kihon/kihon6.rb +40 -0
  215. data/demo/rakuraku/kihon/kihon7.rb +98 -0
  216. data/demo/rakuraku/kihon/kihon8.rb +93 -0
  217. data/demo/rakuraku/kihon/kihon9.rb +77 -0
  218. data/demo/rakuraku/kihon/kihona.rb +123 -0
  219. data/demo/rakuraku/kihon/kihonb.rb +73 -0
  220. data/demo/rakuraku/kihon/kihonc.rb +110 -0
  221. data/demo/rakuraku/layout/lay1.rb +29 -0
  222. data/demo/rakuraku/layout/lay2.rb +28 -0
  223. data/demo/rakuraku/layout/lay3.rb +29 -0
  224. data/demo/rakuraku/map3d/map3d1.rb +68 -0
  225. data/demo/rakuraku/map3d/map3d2.rb +67 -0
  226. data/demo/rakuraku/map3d/map3d3.rb +66 -0
  227. data/demo/rakuraku/map3d/map3d4.rb +107 -0
  228. data/demo/rakuraku/map3d/map3d5.rb +118 -0
  229. data/demo/rakuraku/map3d/map3d6.rb +96 -0
  230. data/demo/rakuraku/map3d/map3d7.rb +98 -0
  231. data/demo/rakuraku/miss/miss1.rb +76 -0
  232. data/demo/rakuraku/miss/miss2.rb +60 -0
  233. data/demo/rakuraku/quick/quick1.rb +35 -0
  234. data/demo/rakuraku/quick/quick2.rb +33 -0
  235. data/demo/rakuraku/quick/quick3.rb +46 -0
  236. data/demo/rakuraku/quick/quick4.rb +47 -0
  237. data/demo/rakuraku/quick/quick5.rb +46 -0
  238. data/demo/rakuraku/u2d/u2d1.rb +50 -0
  239. data/demo/rakuraku/u2d/u2d2.rb +51 -0
  240. data/demo/rakuraku/u2d/u2d3.rb +70 -0
  241. data/demo/rakuraku/u2d/u2d4.rb +59 -0
  242. data/demo/rakuraku/u2d/u2d5.rb +61 -0
  243. data/demo/rakuraku/u2d/u2d6.rb +57 -0
  244. data/demo/rakuraku/u2d/u2d7.rb +103 -0
  245. data/demo/rakuraku/uspack/uspac1.rb +35 -0
  246. data/demo/rakuraku/uspack/uspac2.rb +66 -0
  247. data/demo/rakuraku/uspack/uspac3.rb +37 -0
  248. data/demo/rakuraku/uspack/uspac4.rb +49 -0
  249. data/demo/rakuraku/uxyz/uxyz1.rb +37 -0
  250. data/demo/rakuraku/uxyz/uxyz2.rb +35 -0
  251. data/demo/rakuraku/uxyz/uxyz3.rb +38 -0
  252. data/demo/rakuraku/uxyz/uxyz4.rb +36 -0
  253. data/demo/rakuraku/uxyz/uxyz5.rb +57 -0
  254. data/demo/rakuraku/uxyz/uxyz6.rb +43 -0
  255. data/demo/rakuraku/uxyz/uxyz7.rb +39 -0
  256. data/demo/rakuraku/uxyz/uxyz8.rb +56 -0
  257. data/demo/rubydcloriginal/uemrkz1.rb +37 -0
  258. data/dummy.c +9 -0
  259. data/extconf.rb +133 -0
  260. data/grph1_csgi.c +119 -0
  261. data/grph1_scpack.c +1388 -0
  262. data/grph1_sgpack.c +3876 -0
  263. data/grph1_slpack.c +432 -0
  264. data/grph1_stpack.c +1044 -0
  265. data/grph1_swpack.c +1922 -0
  266. data/grph1_zgpack.c.org +141 -0
  267. data/grph2_grpack.c +368 -0
  268. data/grph2_ucpack.c +1232 -0
  269. data/grph2_udpack.c +1404 -0
  270. data/grph2_uepack.c +1443 -0
  271. data/grph2_ugpack.c +1073 -0
  272. data/grph2_uhpack.c +943 -0
  273. data/grph2_ulpack.c +1304 -0
  274. data/grph2_umpack.c +1540 -0
  275. data/grph2_uspack.c +2481 -0
  276. data/grph2_uupack.c +1689 -0
  277. data/grph2_uvpack.c +943 -0
  278. data/grph2_uwpack.c +552 -0
  279. data/grph2_uxpack.c +555 -0
  280. data/grph2_uypack.c +555 -0
  281. data/grph2_uzpack.c +1455 -0
  282. data/init.c.default +187 -0
  283. data/init.c.gtk +189 -0
  284. data/lib/dcl.rb +5084 -0
  285. data/math1_blklib.c +227 -0
  286. data/math1_chrlib.c +239 -0
  287. data/math1_fnclib.c +254 -0
  288. data/math1_gnmlib.c +356 -0
  289. data/math1_gt2dlib.c +514 -0
  290. data/math1_ifalib.c +427 -0
  291. data/math1_indxlib.c +709 -0
  292. data/math1_intlib.c +167 -0
  293. data/math1_lrllib.c +817 -0
  294. data/math1_maplib.c +172 -0
  295. data/math1_oslib.c +111 -0
  296. data/math1_rfalib.c +1417 -0
  297. data/math1_rfblib.c +232 -0
  298. data/math1_sublib.c +221 -0
  299. data/math1_syslib.c +2025 -0
  300. data/math1_vialib.c +1156 -0
  301. data/math1_viblib.c +1027 -0
  302. data/math1_vralib.c +1156 -0
  303. data/math1_vrblib.c +1027 -0
  304. data/math1_xfclib.c +142 -0
  305. data/math2_fftlib.c +833 -0
  306. data/math2_intrlib.c +138 -0
  307. data/math2_rnmlib.c +226 -0
  308. data/math2_shtlib.c +1568 -0
  309. data/math2_vstlib.c +366 -0
  310. data/misc1_chnlib.c +179 -0
  311. data/misc1_datelib.c +953 -0
  312. data/misc1_fmtlib.c +99 -0
  313. data/misc1_misclib.c +118 -0
  314. data/misc1_randlib.c +142 -0
  315. data/misc1_timelib.c +380 -0
  316. data/obsolete/README +2 -0
  317. data/obsolete/dcl_ext.html +112 -0
  318. data/obsolete/dcl_ext.rb +313 -0
  319. metadata +433 -0
@@ -0,0 +1,107 @@
1
+ #
2
+ # map3d4.rb
3
+ # $Id: map3d4.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+ include NumRu
10
+ include Math
11
+
12
+
13
+ xmin = -50
14
+ xmax = 50
15
+ ymin = -50
16
+ ymax = 50
17
+ vxmin = 0.0
18
+ vxmax = 0.8
19
+ vymin = 0.0
20
+ vymax = 0.8
21
+ zmin = -50
22
+ zmax = 50
23
+ vzmin = 0.0
24
+ vzmax = 0.8
25
+ xvp3 = 2.5
26
+ yvp3 = -1.0
27
+ zvp3 = 1.5
28
+ xfc3 = (vxmax-vxmin)/2
29
+ yfc3 = (vymax-vymin)/2
30
+ zfc3 = (vzmax-vzmin)/2
31
+
32
+
33
+ def aplot(ijk)
34
+
35
+ nmax = 40
36
+ x = NArray.sfloat(nmax+1)
37
+ y = NArray.sfloat(nmax+1)
38
+
39
+ DCL::slpvpr(1)
40
+
41
+ dt = 2.0*PI / nmax
42
+ for n in 0..nmax
43
+ x[n] = 40.0*sin(n*dt)
44
+ y[n] = 40.0*cos(n*dt)
45
+ end
46
+
47
+ DCL:: sgplu(x, y)
48
+
49
+ dt = 2.0*PI / 3
50
+ for n in 0..3
51
+ x[n] = 40.0*sin(n*dt)
52
+ y[n] = 40.0*cos(n*dt)
53
+ end
54
+
55
+ DCL::sgplu(x, y)
56
+
57
+ DCL::sgstxs(0.07)
58
+ DCL::sgtxu(0.0, 0.0, 'DENNOU')
59
+ cttl= format("%2.2d",ijk)
60
+ DCL::sgtxu(0.0, -30.0, cttl)
61
+
62
+ end
63
+
64
+ #-- graph ---
65
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
66
+ DCL::sgopn iws
67
+
68
+ DCL::sgiset('IFONT', 2)
69
+
70
+ DCL::sgfrm
71
+
72
+ #-- X-Y PLANE ----
73
+ DCL::sgswnd(xmin, xmax, ymin, ymax)
74
+ DCL::sgsvpt(vxmin, vxmax, vymin, vymax)
75
+ DCL::sgstrn(1)
76
+ DCL::sgstrf
77
+
78
+ DCL::scspln(1, 2, vzmax)
79
+ DCL::scseye(xvp3, yvp3, zvp3)
80
+ DCL::scsobj(xfc3, yfc3, zfc3)
81
+ DCL::scsprj
82
+
83
+ aplot(1)
84
+
85
+ #-- X-Z PLANE ----
86
+ DCL::sgswnd(xmin, xmax, zmin, zmax)
87
+ DCL::sgsvpt(vxmin, vxmax, vzmin, vzmax)
88
+ DCL::sgstrn(1)
89
+ DCL::sgstrf
90
+
91
+ DCL::scspln(1, 3, vymin)
92
+ DCL::scsprj
93
+
94
+ aplot(2)
95
+
96
+ #-- Y-Z PLANE ----
97
+ DCL::sgswnd(ymin, ymax, zmin, zmax)
98
+ DCL::sgsvpt(vymin, vymax, vzmin, vzmax)
99
+ DCL::sgstrn(1)
100
+ DCL::sgstrf
101
+
102
+ DCL::scspln(2, 3, vxmax)
103
+ DCL::scsprj
104
+
105
+ aplot(3)
106
+
107
+ DCL::sgcls
@@ -0,0 +1,118 @@
1
+ #
2
+ # map3d5.rb
3
+ # $Id: map3d5.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+ include NumRu
10
+ include Math
11
+
12
+
13
+ nx = 21
14
+ ny = 21
15
+ xmin = -10
16
+ xmax = 10
17
+ ymin = -10
18
+ ymax = 10
19
+ vxmin = 0.2
20
+ vxmax = 0.8
21
+ vymin = 0.1
22
+ vymax = 0.5
23
+ zmin = 0.0
24
+ zmax = 20
25
+ vzmin = 0.0
26
+ vzmax = 0.6
27
+ xvp3 = -0.7
28
+ yvp3 = -0.7
29
+ zvp3 = 1.2
30
+ xfc3 = (vxmax-vxmin)/2
31
+ yfc3 = (vymax-vymin)/2
32
+ zfc3 = (vzmax-vzmin)/2
33
+ dx1 = 1
34
+ dx2 = 5
35
+ dy1 = 1
36
+ dy2 = 4
37
+ dz1 = 1
38
+ dz2 = 5
39
+ kmax = 5
40
+ pmin = 0
41
+ pmax = 1
42
+ u = NArray.sfloat(nx, ny)
43
+ v = NArray.sfloat(nx, ny)
44
+ p = NArray.sfloat(nx, ny)
45
+
46
+ #-- data ---
47
+ for j in 0..ny-1
48
+ for i in 0..nx-1
49
+ x = xmin + (xmax-xmin)*i/(nx-1).to_f
50
+ y = ymin + (ymax-ymin)*j/(ny-1).to_f
51
+ u[i,j] = x
52
+ v[i,j] = -y
53
+ p[i,j] = exp(-x**2/64 -y**2/25)
54
+ end
55
+ end
56
+
57
+ #-- graph ---
58
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
59
+ DCL::sgopn iws
60
+
61
+ DCL::sgfrm
62
+
63
+ #-- X-Y ʿ��: ����٥� ----
64
+ DCL::sgswnd(xmin, xmax, ymin, ymax)
65
+ DCL::sgsvpt(vxmin, vxmax, vymin, vymax)
66
+ DCL::sgstrn(1)
67
+ DCL::sgstrf
68
+
69
+ DCL::scspln(1, 2, vzmin)
70
+ DCL::scseye(xvp3, yvp3, zvp3)
71
+ DCL::scsobj(xfc3, yfc3, zfc3)
72
+ DCL::scsprj
73
+
74
+ DCL::uxaxdv('B', dx1, dx2)
75
+ DCL::uxaxdv('T', dx1, dx2)
76
+ DCL::uxsttl('B', 'X-axis', 0.0)
77
+
78
+ DCL::uyaxdv('L', dy1, dy2)
79
+ DCL::uyaxdv('R', dy1, dy2)
80
+ DCL::uysttl('L', 'Y-axis', 0.0)
81
+
82
+ DCL::ugrset('RSIZET', 0.014)
83
+ DCL::ugvect(u, v)
84
+
85
+ #-- X-Y ʿ��: ���٥� ----
86
+ vzlev = vzmin + (vzmax-vzmin)*0.6
87
+ DCL::scspln(1, 2, vzlev)
88
+ DCL::scsprj
89
+
90
+ dp = (pmax-pmin)/kmax.to_f
91
+ for k in 1..kmax
92
+ tlev1 = (k-1)*dp
93
+ tlev2 = tlev1 + dp
94
+ ipat = 600 + k - 1
95
+ DCL::uestlv(tlev1, tlev2, ipat)
96
+ end
97
+ DCL::uetone(p)
98
+
99
+ DCL::udlset('LMSG', false)
100
+ DCL::udgclb(p, 0.1)
101
+ DCL::udcntr(p)
102
+ DCL::slpvpr(1)
103
+
104
+ #-- X-Z ʿ�� ----
105
+ DCL::sgswnd(xmin, xmax, zmin, zmax)
106
+ DCL::sgsvpt(vxmin, vxmax, vzmin, vzmax)
107
+ DCL::sgstrn(1)
108
+ DCL::sgstrf
109
+
110
+ DCL::scspln(1, 3, vymax)
111
+ DCL::scsprj
112
+
113
+ DCL::uzinit
114
+ DCL::uyaxdv('L', dz1, dz2)
115
+ DCL::uysttl('L', 'Z-axis', 0.0)
116
+
117
+ DCL::sgcls
118
+
@@ -0,0 +1,96 @@
1
+ # map3d6.rb
2
+
3
+ require "narray"
4
+ require "numru/dcl"
5
+ include NMath
6
+ include NumRu
7
+
8
+ nx=37; ny=37
9
+ xmin=0.0; xmax=360.0; ymin= -90.0; ymax= 90.0
10
+ vxmin=-0.4; vxmax=0.4; vymin=-0.3; vymax=0.3
11
+ zmin=-1.9; zmax=1.9
12
+ vzmin=-0.2; vzmax=0.2
13
+ drad=PI/180
14
+
15
+ alon = xmin + (xmax-xmin)/(nx-1) * NArray.sfloat(nx).indgen!
16
+ alat = ymin + (ymax-ymin)/(ny-1) * NArray.sfloat(ny).indgen!
17
+
18
+ slat = sin(alat*drad).reshape!(1,ny)
19
+ p = 3*sqrt(1.0-slat**2)*slat*cos(alon*drad) - 0.5*(3*slat**2-1)
20
+
21
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
22
+ DCL::sgopn iws
23
+ DCL::sgfrm
24
+
25
+ #-- X-Y PLANE ----
26
+
27
+ DCL::sgswnd( xmin, xmax, ymin, ymax )
28
+ DCL::sgsvpt( vxmin, vxmax, vymin, vymax )
29
+ DCL::sgstrn( 1 )
30
+ DCL::sgstrf
31
+
32
+ DCL::scspln( 1, 2, vzmin)
33
+ DCL::scseye( -1.1, -1.1, 2.5 )
34
+ DCL::scsobj( 0.0, 0.0, 0.0 )
35
+ DCL::scsprj
36
+
37
+ DCL::uxaxdv( 'B', 10.0, 60.0 )
38
+ DCL::uxaxdv( 'T', 10.0, 60.0 )
39
+ DCL::uxsttl( 'B', 'LONGITUDE', 0.0 )
40
+ DCL::uyaxdv( 'L', 10.0, 30.0 )
41
+ DCL::uyaxdv( 'R', 10.0, 30.0 )
42
+ DCL::uysttl( 'L', 'LATITUDE', 0.0 )
43
+
44
+ #-- X-Z PLANE ----
45
+
46
+ DCL::sgswnd( xmin, xmax, zmin, zmax )
47
+ DCL::sgsvpt( vxmin, vxmax, vzmin, vzmax )
48
+ DCL::sgstrn( 1 )
49
+ DCL::sgstrf
50
+
51
+ DCL::scspln( 1, 3, vymax)
52
+ DCL::scsprj
53
+
54
+ DCL::uzinit
55
+ DCL::uxaxdv( 'T', 10.0, 60.0 )
56
+ DCL::uzlset( 'LABELXB', false )
57
+ DCL::uxaxdv( 'B', 10.0, 60.0 )
58
+ DCL::uyaxdv( 'L', 0.2, 1.0 )
59
+ DCL::uyaxdv( 'R', 0.2, 1.0 )
60
+ DCL::uysttl( 'L', 'AMPLITUDE', 0.0 )
61
+
62
+ #-- Y-Z PLANE ----
63
+
64
+ DCL::sgswnd( ymin, ymax, zmin, zmax )
65
+ DCL::sgsvpt( vymin, vymax, vzmin, vzmax )
66
+ DCL::sgstrn( 1 )
67
+ DCL::sgstrf
68
+
69
+ DCL::scspln( 2, 3, vxmax)
70
+ DCL::scsprj
71
+
72
+ DCL::uzinit
73
+ DCL::uxaxdv( 'T', 10.0, 30.0 )
74
+ DCL::uxaxdv( 'B', 10.0, 30.0 )
75
+ DCL::uzlset( 'LABELYL', false )
76
+ DCL::uyaxdv( 'L', 0.2, 1.0 )
77
+ DCL::uyaxdv( 'R', 0.2, 1.0 )
78
+
79
+ #---------------- 3-D ------------------
80
+
81
+ DCL::scsvpt(vxmin, vxmax, vymin, vymax, vzmin, vzmax)
82
+ DCL::scswnd( xmin, xmax, ymin, ymax, zmin, zmax)
83
+ DCL::scslog(false, false, false)
84
+ DCL::scstrn(1)
85
+ DCL::scstrf
86
+
87
+ (ny-2).downto(0){|j|
88
+ (nx-2).downto(0){|i|
89
+ xp = [ alon[i], alon[i], alon[i+1], alon[i+1], alon[i] ]
90
+ yp = [ alat[j], alat[j+1], alat[j+1], alat[j], alat[j] ]
91
+ zp = [ p[i,j], p[i,j+1], p[i+1,j+1], p[i+1,j], p[i,j] ]
92
+ DCL::scplu(xp, yp, zp)
93
+ }
94
+ }
95
+
96
+ DCL::sgcls
@@ -0,0 +1,98 @@
1
+ # map3d6.rb
2
+
3
+ require "narray"
4
+ require "numru/dcl"
5
+ include NMath
6
+ include NumRu
7
+
8
+ nx=37; ny=37
9
+ xmin=0.0; xmax=360.0; ymin= -90.0; ymax= 90.0
10
+ vxmin=-0.4; vxmax=0.4; vymin=-0.3; vymax=0.3
11
+ zmin=-1.9; zmax=1.9
12
+ vzmin=-0.2; vzmax=0.2
13
+ drad=PI/180
14
+
15
+ alon = xmin + (xmax-xmin)/(nx-1) * NArray.sfloat(nx).indgen!
16
+ alat = ymin + (ymax-ymin)/(ny-1) * NArray.sfloat(ny).indgen!
17
+
18
+ slat = sin(alat*drad).reshape!(1,ny)
19
+ p = 3*sqrt(1.0-slat**2)*slat*cos(alon*drad) - 0.5*(3*slat**2-1)
20
+
21
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
22
+ DCL::sgopn iws
23
+ DCL::sgfrm
24
+
25
+ #-- X-Y PLANE ----
26
+
27
+ DCL::sgswnd( xmin, xmax, ymin, ymax )
28
+ DCL::sgsvpt( vxmin, vxmax, vymin, vymax )
29
+ DCL::sgstrn( 1 )
30
+ DCL::sgstrf
31
+
32
+ DCL::scspln( 1, 2, vzmin)
33
+ DCL::scseye( -2.1, -1.1, 0.7 )
34
+ DCL::scsobj( 0.0, 0.0, 0.0 )
35
+ DCL::scsprj
36
+
37
+ DCL::uxaxdv( 'B', 10.0, 60.0 )
38
+ DCL::uxaxdv( 'T', 10.0, 60.0 )
39
+ DCL::uxsttl( 'B', 'LONGITUDE', 0.0 )
40
+ DCL::uyaxdv( 'L', 10.0, 30.0 )
41
+ DCL::uyaxdv( 'R', 10.0, 30.0 )
42
+ DCL::uysttl( 'L', 'LATITUDE', 0.0 )
43
+
44
+ #-- X-Z PLANE ----
45
+
46
+ DCL::sgswnd( xmin, xmax, zmin, zmax )
47
+ DCL::sgsvpt( vxmin, vxmax, vzmin, vzmax )
48
+ DCL::sgstrn( 1 )
49
+ DCL::sgstrf
50
+
51
+ DCL::scspln( 1, 3, vymax)
52
+ DCL::scsprj
53
+
54
+ DCL::uzinit
55
+ DCL::uxaxdv( 'T', 10.0, 60.0 )
56
+ DCL::uzlset( 'LABELXB', false )
57
+ DCL::uxaxdv( 'B', 10.0, 60.0 )
58
+ DCL::uyaxdv( 'L', 0.2, 1.0 )
59
+ DCL::uyaxdv( 'R', 0.2, 1.0 )
60
+ DCL::uysttl( 'L', 'AMPLITUDE', 0.0 )
61
+
62
+ #-- Y-Z PLANE ----
63
+
64
+ DCL::sgswnd( ymin, ymax, zmin, zmax )
65
+ DCL::sgsvpt( vymin, vymax, vzmin, vzmax )
66
+ DCL::sgstrn( 1 )
67
+ DCL::sgstrf
68
+
69
+ DCL::scspln( 2, 3, vxmax)
70
+ DCL::scsprj
71
+
72
+ DCL::uzinit
73
+ DCL::uxaxdv( 'T', 10.0, 30.0 )
74
+ DCL::uxaxdv( 'B', 10.0, 30.0 )
75
+ DCL::uzlset( 'LABELYL', false )
76
+ DCL::uyaxdv( 'L', 0.2, 1.0 )
77
+ DCL::uyaxdv( 'R', 0.2, 1.0 )
78
+
79
+ #---------------- 3-D ------------------
80
+
81
+ DCL::scsvpt(vxmin, vxmax, vymin, vymax, vzmin, vzmax)
82
+ DCL::scswnd( xmin, xmax, ymin, ymax, zmin, zmax)
83
+ DCL::scslog(false, false, false)
84
+ DCL::scstrn(1)
85
+ DCL::scstrf
86
+
87
+ (ny-2).downto(0){|j|
88
+ (nx-2).downto(0){|i|
89
+ xp = [ alon[i], alon[i], alon[i+1], alon[i+1], alon[i] ]
90
+ yp = [ alat[j], alat[j+1], alat[j+1], alat[j], alat[j] ]
91
+ zp = [ p[i,j], p[i,j+1], p[i+1,j+1], p[i+1,j], p[i,j] ]
92
+ DCL::sctnzu(xp[0..2], yp[0..2], zp[0..2], 10999, 50999)
93
+ DCL::sctnzu(xp[2..4], yp[2..4], zp[2..4], 10999, 50999)
94
+ DCL::scplu(xp, yp, zp)
95
+ }
96
+ }
97
+
98
+ DCL::sgcls
@@ -0,0 +1,76 @@
1
+ #
2
+ # miss1.rb
3
+ # $Id: miss1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include NMath
11
+
12
+
13
+ nmax = 40
14
+ xmin = 0.0
15
+ xmax = 4*PI
16
+ ymin = -1.0
17
+ ymax = 1.0
18
+
19
+ x = NArray.sfloat(nmax+1)
20
+ y = NArray.sfloat(nmax+1)
21
+
22
+ #-- data ----
23
+ dt = xmax/nmax
24
+ x = x.indgen * dt
25
+ y = sin(x)
26
+
27
+ n1 = nmax/4
28
+ rmiss = DCL::glrget('RMISS')
29
+ y[n1-1] = rmiss
30
+ y[n1] = rmiss
31
+ y[n1+1] = rmiss
32
+
33
+ n2 = n1*3
34
+ y[n2-1] = rmiss
35
+ y[n2+1] = rmiss
36
+
37
+ #-- graph ----
38
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
39
+ DCL::sgopn iws
40
+ DCL::grfrm
41
+
42
+ #-- default ----
43
+ DCL::grswnd(xmin, xmax, ymin, ymax)
44
+ DCL::grsvpt( 0.1, 0.9, 0.7, 0.9)
45
+ DCL::grstrn(1)
46
+ DCL::grstrf
47
+
48
+ DCL::sgplu(x, y)
49
+
50
+ DCL::grswnd(xmin, xmax, ymin, ymax)
51
+ DCL::grsvpt( 0.1, 0.9, 0.5, 0.7)
52
+ DCL::grstrn(1)
53
+ DCL::grstrf
54
+
55
+ DCL::sgspmt(3)
56
+ DCL::sgpmu(x, y)
57
+
58
+ #-- cover loss data ----
59
+ DCL::gllset('LMISS', true)
60
+
61
+ DCL::grswnd(xmin, xmax, ymin, ymax)
62
+ DCL::grsvpt( 0.1, 0.9, 0.3, 0.5)
63
+ DCL::grstrn(1)
64
+ DCL::grstrf
65
+
66
+ DCL::sgplu(x, y)
67
+
68
+ DCL::grswnd(xmin, xmax, ymin, ymax)
69
+ DCL::grsvpt( 0.1, 0.9, 0.1, 0.3)
70
+ DCL::grstrn(1)
71
+ DCL::grstrf
72
+
73
+ DCL::sgpmu(x, y)
74
+
75
+ DCL::grcls
76
+