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,62 @@
1
+ #
2
+ # step2.rb
3
+ # $Id: step2.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include Math
11
+
12
+ nmax = 100
13
+
14
+ x = NArray.sfloat(nmax)
15
+ y = NArray.sfloat( nmax)
16
+
17
+ #-- data ----
18
+ iseed = 0
19
+ x[0] = 2*(rand(iseed)-0.5)
20
+ y[nmax-1] = x[0]
21
+
22
+ for i in 1..nmax-1
23
+ x[i] = 2*(rand(iseed)-0.5)
24
+ y[i-1] = x[i]
25
+ end
26
+
27
+ #-- graph ----
28
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
29
+ DCL::gropn iws
30
+ DCL::grfrm
31
+
32
+ DCL::usspnt(x, y)
33
+ DCL::uspfit
34
+ DCL::grstrf
35
+
36
+ DCL::ussttl( 'X-TITLE', 'x-unit', 'Y-TITLE', 'y-unit' )
37
+ DCL::usdaxs
38
+
39
+ cnt = nmax/4
40
+ st = 0
41
+ ed = cnt-1
42
+ DCL::uusmki(5)
43
+ DCL::uusmks(0.015)
44
+ DCL::uumrk(x[st..ed], y[st..ed])
45
+
46
+ st += cnt
47
+ ed += cnt
48
+ DCL::uusmkt(2)
49
+ DCL::uumrk(x[st..ed], y[st..ed])
50
+
51
+ st += cnt
52
+ ed += cnt
53
+ DCL::uusmkt(3)
54
+ DCL::uumrk(x[st..ed], y[st..ed])
55
+
56
+ st += cnt
57
+ ed += cnt
58
+ DCL::uusmkt(4)
59
+ DCL::uumrk(x[st..ed], y[st..ed])
60
+
61
+ DCL::grcls
62
+
@@ -0,0 +1,37 @@
1
+ #
2
+ # u1d1.rb
3
+ # $Id: u1d1.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include Math
11
+
12
+ nmax = 51
13
+ xmin = 1945
14
+ xmax = 1995
15
+
16
+ y = NArray.sfloat(nmax)
17
+
18
+ #-- data ----
19
+ y0 = 0.5
20
+ for n in 0..nmax-1
21
+ y[n] = 5.0*y0 + 10.0
22
+ y0 = 3.7*y0*(1.0-y0)
23
+ end
24
+
25
+ rundef = DCL::glrget('RUNDEF')
26
+
27
+ #-- graph ----
28
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
29
+ DCL::gropn iws
30
+ DCL::grfrm
31
+
32
+ DCL::grswnd(xmin, xmax, rundef, rundef)
33
+ DCL::ussttl('TIME', 'YEAR', 'TEMPERATURE', 'DEG')
34
+ DCL::usgrph(rundef, y)
35
+
36
+ DCL::grcls
37
+
@@ -0,0 +1,50 @@
1
+ #
2
+ # u1d2.rb
3
+ # $Id: u1d2.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include Math
11
+
12
+ nmax = 50
13
+ ymin = 0.0
14
+ ymax = 50.0
15
+
16
+ x1 = NArray.sfloat(nmax+1)
17
+ x2 = NArray.sfloat(nmax+1)
18
+
19
+ #-- data ----
20
+ iseed = 0
21
+ for n in 0..nmax
22
+ y = ymin + (ymax-ymin)*n/nmax
23
+ x1[n] = 10.0*(exp(-y/20))**2 * exp((rand(iseed)-0.5)*2)**2
24
+ x2[n] = 10.0*(exp(-y/20))**2
25
+ end
26
+
27
+ rundef = DCL::glrget('RUNDEF')
28
+
29
+ #-- graph ----
30
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
31
+ DCL::gropn iws
32
+ DCL::grfrm
33
+
34
+ DCL::grswnd(rundef, rundef, ymin, ymax)
35
+ DCL::usspnt(x1, rundef)
36
+ DCL::usspnt(x2, rundef)
37
+ DCL::grstrn(3)
38
+ DCL::uspfit
39
+ DCL::grstrf
40
+
41
+ DCL::ussttl('MIXING RATIO', 'ppmv', 'HEIGHT', 'km')
42
+ DCL::usdaxs
43
+
44
+ DCL::uusmkt(2)
45
+ DCL::uumrk(x1, rundef)
46
+
47
+ DCL::uulin(x2, rundef)
48
+
49
+ DCL::grcls
50
+
@@ -0,0 +1,46 @@
1
+ #
2
+ # u2d1.rb
3
+ # $Id: u2d1.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include NMath
11
+
12
+ nt = 51
13
+ nz = 21
14
+ tmin = 0.0
15
+ tmax = 5.0
16
+ zmin = 20.0
17
+ zmax = 50.0
18
+ dt = (tmax-tmin)/(nt-1)
19
+ dz = (zmax-zmin)/(nz-1)
20
+
21
+ u = NArray.sfloat(nt, nz)
22
+
23
+ #-- data ----
24
+ zj = NArray.sfloat(nz).indgen! * dz
25
+ uz = exp(-0.2*zj)*(zj**0.5)
26
+ tz = -2.0*exp(-0.1*zj)
27
+ for i in 0..nt-1
28
+ u[i,true] = uz*sin(3.0*(tz+i*dt))
29
+ end
30
+
31
+ #-- graph ----
32
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
33
+ DCL::gropn iws
34
+ DCL::grfrm
35
+
36
+ DCL::grswnd(tmin, tmax, zmin, zmax)
37
+ DCL::uspfit
38
+ DCL::grstrf
39
+
40
+ DCL::ussttl('TIME', 'YEAR', 'HEIGHT', 'km')
41
+ DCL::usdaxs
42
+
43
+ DCL::udcntr(u)
44
+
45
+ DCL::grcls
46
+
@@ -0,0 +1,49 @@
1
+ #
2
+ # u2d2.rb
3
+ # $Id: u2d2.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include NMath
11
+
12
+ nt = 51
13
+ nz = 21
14
+ tmin = 0.0
15
+ tmax = 5.0
16
+ zmin = 20.0
17
+ zmax = 50.0
18
+ dt = (tmax-tmin)/(nt-1)
19
+ dz = (zmax-zmin)/(nz-1)
20
+
21
+ u = NArray.sfloat(nt, nz)
22
+
23
+ #-- data ----
24
+ zj = NArray.sfloat(nz).indgen! * dz
25
+ uz = exp(-0.2*zj)*(zj**0.5)
26
+ tz = -2.0*exp(-0.1*zj)
27
+ for i in 0..nt-1
28
+ u[i,true] = uz*sin(3.0*(tz+i*dt))
29
+ end
30
+
31
+ #-- graph ----
32
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
33
+ DCL::gropn iws
34
+ DCL::sglset('LSOFTF', true)
35
+ DCL::grfrm
36
+
37
+ DCL::grswnd(tmin, tmax, zmin, zmax)
38
+ DCL::grsvpt(0.2, 0.9, 0.4, 0.8)
39
+ DCL::uspfit
40
+ DCL::grstrf
41
+
42
+ DCL::ussttl('TIME', 'YEAR', 'HEIGHT', 'km')
43
+ DCL::usdaxs
44
+
45
+ DCL::udcntr(u)
46
+ DCL::uetone(u)
47
+
48
+ DCL::grcls
49
+
@@ -0,0 +1,45 @@
1
+ #
2
+ # u2d3.rb
3
+ # $Id: u2d3.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include Math
11
+
12
+ nx = 21
13
+ ny = 21
14
+ xmin = -1.0
15
+ xmax = 1.0
16
+ ymin = -1.0
17
+ ymax = 1.0
18
+ dx = (xmax-xmin)/(nx-1)
19
+ dy = (ymax-ymin)/(ny-1)
20
+
21
+ u = NArray.sfloat(nx, ny)
22
+ v = NArray.sfloat(nx, ny)
23
+
24
+ #-- data ----
25
+ xi = xmin + NArray.sfloat(nx).indgen! * dx
26
+ yj = ymin + NArray.sfloat(ny).indgen! * dy
27
+ for j in 0..ny-1; u[true,j] = xi; end
28
+ for i in 0..nx-1; v[i,true] = -yj; end
29
+
30
+ #-- graph ----
31
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
32
+ DCL::gropn iws
33
+ DCL::grfrm
34
+
35
+ DCL::grswnd(xmin, xmax, ymin, ymax)
36
+ DCL::uspfit
37
+ DCL::grstrf
38
+
39
+ DCL::ussttl('X', 'km', 'Y', 'km')
40
+ DCL::usdaxs
41
+
42
+ DCL::ugvect(u, v)
43
+
44
+ DCL::grcls
45
+
@@ -0,0 +1,61 @@
1
+ #
2
+ # u2d4.rb
3
+ # $Id: u2d4.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include NMath
11
+
12
+ nx = 21
13
+ ny = 21
14
+ xmin = -10.0
15
+ xmax = 10.0
16
+ ymin = -10.0
17
+ ymax = 10.0
18
+ kmax = 5
19
+ pmin = 0.0
20
+ pmax = 1.0
21
+ dx = (xmax-xmin)/(nx-1)
22
+ dy = (ymax-ymin)/(ny-1)
23
+ dp = (pmax-pmin)/kmax
24
+
25
+ u = NArray.sfloat(nx, ny)
26
+ v = NArray.sfloat(nx, ny)
27
+ p = NArray.sfloat(nx, ny)
28
+
29
+
30
+ #-- data ----
31
+ xi = xmin + NArray.sfloat(nx).indgen! * dx
32
+ yj = ymin + NArray.sfloat(ny).indgen! * dy
33
+ for j in 0..ny-1; u[true,j] = xi; end
34
+ for i in 0..nx-1; v[i,true] = -yj; end
35
+ p = exp(-xi.reshape(nx,1)**2/64 - yj.reshape(1,ny)**2/25)
36
+
37
+ #-- graph ----
38
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
39
+ DCL::gropn iws
40
+ DCL::sglset('LSOFTF', true)
41
+ DCL::grfrm
42
+
43
+ DCL::grswnd(xmin, xmax, ymin, ymax)
44
+ DCL::uspfit
45
+ DCL::grstrf
46
+
47
+ DCL::usdaxs
48
+ DCL::ugvect(u, v)
49
+ DCL::udcntr(p)
50
+
51
+ for k in 0..kmax-1
52
+ tlev1 = k*dp
53
+ tlev2 = tlev1 + dp
54
+ ipat = 600 + k
55
+ DCL::uestlv(tlev1, tlev2, ipat)
56
+ end
57
+
58
+ DCL::uetone(p)
59
+
60
+ DCL::grcls
61
+
@@ -0,0 +1,40 @@
1
+ #
2
+ # scpkt2.rb
3
+ #
4
+ # $Id: scpkt2.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ #-- graph ---
15
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
16
+
17
+ DCL::sglset('L2TO3', true)
18
+ DCL::sgiset('IFONT', 2)
19
+
20
+ DCL::sgopn iws
21
+
22
+ x0 = -3
23
+ dx = 0.3
24
+ for i in 1..20
25
+ DCL::sgfrm
26
+ DCL::sgswnd(0.0, 1.0, 0.0, 1.0)
27
+ DCL::sgsvpt(0.0, 1.0, 0.0, 1.0)
28
+ DCL::sgstrn(1)
29
+ DCL::sgstrf
30
+
31
+ x = x0 + dx*i
32
+ DCL::sgrset('XEYE3', x)
33
+ DCL::scsprj
34
+
35
+ DCL::slpvpr(1)
36
+ DCL::sgtxzv(0.5, 0.5, 'DENNOU', 0.1, 0, 0, 3)
37
+ end
38
+
39
+ DCL::sgcls
40
+
@@ -0,0 +1,102 @@
1
+ #
2
+ # scpkt3.rb
3
+ #
4
+ # $Id: scpkt3.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ nx = 73
15
+ ny = 37
16
+ xmin = 0
17
+ xmax = 360
18
+ ymin = -90
19
+ ymax = 90
20
+
21
+ p1 = NArray.sfloat(nx*ny)
22
+ p2 = NArray.sfloat(nx*ny)
23
+
24
+ #-- data ---
25
+ i = 0
26
+ File.foreach("t810630.dat") do |line|
27
+ line.split(" ").each do |data|
28
+ p1[i] = data.to_f
29
+ i += 1
30
+ end
31
+ end
32
+ i = 0
33
+ File.foreach("t811231.dat") do |line|
34
+ line.split(" ").each do |data|
35
+ p2[i] = data.to_f
36
+ i += 1
37
+ end
38
+ end
39
+ p1.reshape!(nx, ny)
40
+ p2.reshape!(nx, ny)
41
+
42
+ 190.step(245, 5) do |r|
43
+ amin=r
44
+ amax=r+5
45
+ # IDX=(R-180)*1.4*1000+999
46
+ idx=((r-170)*1.25).to_i*1000+999
47
+ DCL::uestlv(amin, amax, idx)
48
+ end
49
+
50
+ #-- graph ---
51
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
52
+ DCL::sgopn iws
53
+
54
+ DCL::sglset('LSOFTF', false)
55
+ DCL::umlset('LGRIDMJ', false)
56
+ DCL::umrset('DGRIDMN', 30.0)
57
+ DCL::umiset('ITYPEMN', 1)
58
+ DCL::umiset('ITYPEOUT', 1)
59
+ DCL::sgfrm
60
+
61
+ DCL::sglset('L2TO3', true)
62
+ # CALL SCSEYE(-1.0, -1.5, 2.0 )
63
+ DCL::scsobj(0.5, 0.5, 0.0)
64
+ DCL::scspln(1, 2, -0.3)
65
+ DCL::scsprj
66
+
67
+ DCL::sgssim(0.12, 0.0, 0.0)
68
+ DCL::sgsmpl(180.0, 90.0, 0.0)
69
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.9)
70
+ DCL::sgswnd(xmin, xmax, ymin, ymax)
71
+ DCL::sgstxy(-180.0, 180.0, -90.0, 90.0)
72
+ DCL::sgstrn(DCL::isgtrc('MER'))
73
+ DCL::sgstrf
74
+
75
+ DCL::sglset('LCLIP', true)
76
+ DCL::slpwwr(1)
77
+
78
+ DCL::uetone(p1)
79
+ DCL::umpmap('coast_world')
80
+ DCL::umpglb
81
+ DCL::slpvpr(1)
82
+
83
+ DCL::sglset('LCLIP', false)
84
+ DCL::sgtxzv(0.5, 0.95, 'Temp. : 10 hPa : Jun 30, 1981', 0.025, 0, 0, 3)
85
+
86
+ DCL::scspln(1, 2, 0.4)
87
+ # CALL SGRSET( 'YOFF3', -0.25 )
88
+ DCL::scsprj
89
+
90
+ DCL::sglset('LCLIP', true)
91
+ DCL::slpwwr(1)
92
+
93
+ DCL::uetone(p2)
94
+ DCL::umpmap('coast_world')
95
+ DCL::umpglb
96
+ DCL::slpvpr(1)
97
+
98
+ DCL::sglset('LCLIP', false)
99
+ DCL::sgtxzv(0.5, 0.95, 'Temp. : 10 hPa : Dec 31, 1981', 0.025, 0, 0, 3)
100
+
101
+ DCL::sgcls
102
+