home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 October / usenetsourcesnewsgroupsinfomagicoctober1994disk2.iso / misc / volume18 / planet / part04 < prev    next >
Internet Message Format  |  1991-04-10  |  54KB
  1. From: allen@viewlogic.com (Dave Allen)
  2. Newsgroups: comp.sources.misc
  3. Subject: REPOST: v18i004:  planet - planet generation simulator, Part04/04
  4. Message-ID: <1991Apr9.041825.8147@sparky.IMD.Sterling.COM>
  5. Date: 9 Apr 91 04:18:25 GMT
  6. Approved: kent@sparky.imd.sterling.com
  7. X-Checksum-Snefru: b2104cf2 01fc0da8 88f3c561 60378f00
  8.  
  9. Submitted-by: Dave Allen <allen@viewlogic.com>
  10. Posting-number: Volume 18, Issue 4
  11. Archive-name: planet/part04
  12. Supersedes: tec: Volume 10, Issue 77-78
  13.  
  14. #! /bin/sh
  15. # into a shell via "sh file" or similar.  To overwrite existing files,
  16. # type "sh file -c".
  17. # The tool that generated this appeared in the comp.sources.unix newsgroup;
  18. # send mail to comp-sources-unix@uunet.uu.net if you want that tool.
  19. # Contents:  ./README ./doc/params.tec ./doc/readme.clim
  20. #   ./doc/readme.globe ./doc/readme.tec ./example/clim.in
  21. #   ./example/clim.out ./example/tec.out.1 ./example/tec.out.2
  22. #   ./example/tec.out.3 ./src/clim.h ./src/climate.c ./src/const.h
  23. #   ./src/globe.c ./src/tec.h ./src/unix.c ./src/wind.c
  24. # Wrapped by kent@sparky on Mon Apr  8 22:39:15 1991
  25. PATH=/bin:/usr/bin:/usr/ucb ; export PATH
  26. echo If this archive is complete, you will see the following message:
  27. echo '          "shar: End of archive 4 (of 4)."'
  28. if test -f './README' -a "${1}" != "-c" ; then 
  29.   echo shar: Will not clobber existing file \"'./README'\"
  30. else
  31.   echo shar: Extracting \"'./README'\" \(849 characters\)
  32.   sed "s/^X//" >'./README' <<'END_OF_FILE'
  33. XThis is version 1 of a scientific planet generator (as opposed to a
  34. Xdie-rolling planet generator).  It is written in C and runs in graphics
  35. Xmode under Xwindows, or in text-only mode under UN*X.  There are three
  36. Xprograms: 
  37. X
  38. X   - tec: a continental drift simulator
  39. X   - clim: a climate generator
  40. X   - globe: computes and displays a rotating planet
  41. X
  42. XCopyright (c) 1991 by David Allen.  You may distribute this freely as long
  43. Xas you leave my name and copyright notice intact.
  44. X
  45. XI'd really like your comments and feedback; send e-mail to
  46. Xallen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore Ave,
  47. XMaynard, MA 01754.
  48. X
  49. XThere are three subdirectories:
  50. X
  51. X   - src: contains all source and a UN*X Makefile
  52. X   - example: contains some sample input and output files
  53. X   - doc: contains details about how to use the programs and how they work.
  54. X
  55. XHave fun!
  56. END_OF_FILE
  57.   if test 849 -ne `wc -c <'./README'`; then
  58.     echo shar: \"'./README'\" unpacked with wrong size!
  59.   fi
  60.   # end of './README'
  61. fi
  62. if test -f './doc/params.tec' -a "${1}" != "-c" ; then 
  63.   echo shar: Will not clobber existing file \"'./doc/params.tec'\"
  64. else
  65.   echo shar: Extracting \"'./doc/params.tec'\" \(5857 characters\)
  66.   sed "s/^X//" >'./doc/params.tec' <<'END_OF_FILE'
  67. XIf you give tec a filename as a command line argument, it will read 
  68. Xparameters from that file.  If you hand it '-' as an argument, it will
  69. Xread them from stdin.  For example, you could type "tec - < foo", or
  70. X"tec foo".  The '-' option is handy for pipes.
  71. X
  72. XThe parameter file is optional; all of the parameters you can change
  73. Xhave defaults.  Thus, the parameter file need contain only the parameters
  74. Xyou want to change.  A parameter file looks like LISP; for example,
  75. Xto change the XSIZE parameter, the file would have the one line
  76. X
  77. X(XSIZE 40)
  78. X
  79. XParameters can also be vectors; for example
  80. X
  81. X(MOVERATE (1 1 1))
  82. X
  83. XThe parameters are either easy or hard.  This rating concerns how much
  84. Xyou have to know about the simulation to change the parameters.  Easy
  85. Xones just affect run time, resolution, and so on.  Hard ones should
  86. Xonly be changed if you understand the code pretty well.
  87. X
  88. XEasy parameters:
  89. XPRINTMODE - 0 produces no output, 1 produces a long text file (default),
  90. X            3 produces a summary text file and 4 produces PostScript
  91. X            grayscale drawings.
  92. XDRAWEVERY - How many steps should be iterated for each output.  Default 1,
  93. X            but this produces an awful lot of text.  Reasonable values
  94. X            for non-interactive runs are 5 or 10.  Supercontinents appear
  95. X            every 25 steps or so.  A value of zero implies that only the
  96. X            final map, at the last timestep, should be output.
  97. XMAXSTEP   - How many steps to run before terminating.  Default 100.
  98. XXSIZE     - Horizontal size of arrays.  Default 90.  To produce summaries
  99. X            that fit on a normal page, I used 72.
  100. XYSIZE     - Vertical size of arrays.  Default 90.  To produce the summaries
  101. X            included in tec.out.[1-3], I used 48.
  102. XHYDROPCT -  The percentage of the grid that should be covered by water.
  103. X            Default is 70.  A threshhold is found for the initial blob to meet
  104. X            this goal within about 2-3%.  If DRAWEVERY is also set to 0 (see
  105. X            above), then at the last step the ocean level is forced up or
  106. X            down to satisfy the requirement.  This may cause a PANIC error if
  107. X            there just isn't enough land to do it.
  108. XUNDERSCAN - If a potential rift gets within this many squares of the edge, it
  109. X            is thrown away.  This keeps new fragments from appearing near the
  110. X            edge, which in turn means that far less land falls off the edge
  111. X            of the map.  This is a hack; the default is 0, which is no change;
  112. X            for 100 x 100 maps, 10 is the best choice.
  113. X
  114. XSo to produce a text file with 10 pictures, you could use 
  115. X
  116. X(PRINTMODE 2) (DRAWEVERY 10) (XSIZE 72) (YSIZE 48)
  117. X
  118. XTo produce a 20-page PostScript output with a little more detail, use
  119. X
  120. X(PRINTMODE 4) (DRAWEVERY 5)
  121. X
  122. XHard Parameters:
  123. XZINIT     - Initial altitude of supercontinent.  Default 22.
  124. XZSUBSUME  - Altitude added to leading edge of drifting plate.  Default 16.
  125. XZCOAST    - Sea level.  Everything below this is under water.  Default 16.
  126. XZSHELF    - Minimum altitude for continental shelf.  Default 8.
  127. XZMOUNTAIN - Used in PRINTMODE 2; altitude above which a square is counted
  128. X            as mountain instead of land.  Default 48.
  129. XZERODE    - Minimum altitude used in erosion.  Default 16.  Decrease this
  130. X            number and continents grow faster at the margins; increase it
  131. X            and they grow more slowly.  Increase it too much and weird
  132. X            things happen to your mountains.
  133. XMAXBUMP   - If two plates overlap on more than this number of squares,
  134. X            the plates must be merged together.  Default 50. 
  135. XBUMPTOL   - If the relative velocity of two touching plates is less than this
  136. X            value, the plates will be merged together.  Default 50.
  137. XMOVERATE  - A vector of real numbers describing the rate at which plates
  138. X            move apart, slow down, and then drift back together.  Each
  139. X            plate has an age, measured in steps from its creation.  The age
  140. X            is used as an offset into the vector.  The result is a real
  141. X            between -1.0 and +1.0 which is multiplied by the plate's original
  142. X            movement vector.  The default vector is
  143. X 
  144. X            (1.0 1.0 1.0 0.7 0.4 0.1 -0.2 -0.5 -0.8 -1.0)
  145. X
  146. X            It can be interpreted this way.  The first three steps a plate is
  147. X            in motion, it moves at its normal velocity.  During the 4-6
  148. X            steps, it slows down.  During the 7-10 steps, it moves the
  149. X            opposite way from the way it started, at increasing speeds.
  150. X            Thereafter, it keeps moving at its final velocity.
  151. XRIFTPCT   - Percent chance that a rift will occur in a given step.  Default 40.
  152. XMAXCTRTRY - The number of tries the rift routine will make to find an
  153. X            acceptable place to start a rift.  Default 50.
  154. XRIFTDIST  - The minimum allowable distance between a proposed rift center
  155. X            and the nearest coast.  Default 5.
  156. XDOERODE   - Whether or not to compute erosion.  Default 1; if 0, the program
  157. X            will run about twice as fast but will generate way too many
  158. X            mountains, since the only thing that gets rid of mountains
  159. X            is erosion.
  160. XERODERND  - The rounding factor used in erosion.  Default 4.  As this 
  161. X            parameter goes from 0 to 7, erosion occurs faster.
  162. XBENDEVERY - Each time a rift grows by this many squares, it is allowed to
  163. X            bend a little bit.  Default 6.
  164. XBENDBY    - When a rift bends, it bends by a random angle in the range
  165. X            -BENDBY ... +BENDBY multiplied by pi/2000.  Default value is 100.
  166. X            If BENDBY is too high or BENDEVERY is too low, very curvy rifts
  167. X            will be generated, with unrealistic results.
  168. XSPPEDBASE - The minimum speed at which new plates will travel.  Default 200.
  169. XSPEEDRNG  - The actual speed at which a plate travels is SPEEDBASE +
  170. X            rnd (SPEEDRNG).  Default value is 300.
  171. END_OF_FILE
  172.   if test 5857 -ne `wc -c <'./doc/params.tec'`; then
  173.     echo shar: \"'./doc/params.tec'\" unpacked with wrong size!
  174.   fi
  175.   # end of './doc/params.tec'
  176. fi
  177. if test -f './doc/readme.clim' -a "${1}" != "-c" ; then 
  178.   echo shar: Will not clobber existing file \"'./doc/readme.clim'\"
  179. else
  180.   echo shar: Extracting \"'./doc/readme.clim'\" \(3183 characters\)
  181.   sed "s/^X//" >'./doc/readme.clim' <<'END_OF_FILE'
  182. XThis is the README file for Version 1 of clim, a program which computes
  183. Xthe climate for a planet: temperature, pressure zones, wind patterns, rain
  184. Xpatterns, and finally climate based on temperature and rain.  It is
  185. Xwritten in C and runs on any Xwindows system (X11R4).  Copyright (c) 1991
  186. Xby David Allen.  You may distribute this freely as long as you leave my
  187. Xname and copyright notice intact.
  188. X
  189. XI'd really like your comments and feedback; send e-mail to
  190. Xallen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore Ave,
  191. XMaynard, MA 01754.
  192. X
  193. XWHAT DOES CLIM DO?
  194. X
  195. XGiven a topographical map, for example, one produced by TEC, this program
  196. Xcomputes temperature, pressure, wind, rain and climate.  The information
  197. Xused in these computations comes from any college-level meterology text.
  198. XOutput is produced as an animation on the graphics screen, on in text
  199. Xfiles including PostScript output.
  200. X
  201. XI have included a sample text map showing climate for one world; it is
  202. Xin example/clim.out.
  203. X
  204. XHOW DO I RUN IT?
  205. X
  206. XTo compile clim, go to the src directory and type "make clim".  To run it,
  207. Xgo to the example directory and type "../src/clim clim.in".  This directs
  208. Xthe program to read the map stored in the file clim.in.
  209. X
  210. XWhen clim runs, it opens a graphics screen.  It prints about 40 lines of
  211. Xstatus information to stderr, then begins to animate the results.
  212. X
  213. XAt any time, you can type "q" into the graphics window to end the program.
  214. XYou can flip between the types of display with these keys:
  215. X   h: display temperature; default.  Temperature is displayed using the
  216. X      full spectrum; purple is coldest, then blue, and so on; red is
  217. X      hottest.  There are about 30 visibly different shades.  Black
  218. X      lines indicate coastline; white lines surround mountains.
  219. X   p: display pressure.  Zones of low pressure are surrounded by white
  220. X      lines; zones of high pressure are surrounded by black lines.  Ocean
  221. X      is blue, land is red, and mountains are grey.
  222. X   w: display wind patterns.  Winds run in circles, especially over the
  223. X      ocean.  Black horizontal lines are westbound; white horizontals are
  224. X      eastbound; black verticals are southbound; white verticals are 
  225. X      northbound.  Ocean is blue, land is red, mountains are grey.
  226. X   r: display rainfall.  Amount of rainfall is displayed using the full
  227. X      spectrum; purple is driest, then blue; red is wettest.  Black lines
  228. X      indicate coastline; white lines surround mountains.
  229. X   c: display climate.  The following colors are used:
  230. X      dark blue:  Ocean           orange:        Savannah
  231. X      light grey: Icebergs        dark green:    Decidous forest 
  232. X      dark grey:  Tundra          bright green:  Jungle
  233. X      brown:      Steppe          bright blue:   Swamp
  234. X      yellow:     Desert
  235. X
  236. XTo produce your own maps for use with clim, you can use tec.  For example,
  237. Xcreate a tec input file with this single line:
  238. X
  239. X(XSIZE 60) (YSIZE 30) (PRINTMODE 3) (MAXSTEP 20)
  240. X
  241. XThen type "tec filename > tec.out" (where filename is the name of the
  242. Xinput file you just created).  Tec should run for a couple of seconds
  243. Xand then produce a file called tec.out.  Then type "clim tec.out" and
  244. Xyou will see the climate for your own planet.
  245. END_OF_FILE
  246.   if test 3183 -ne `wc -c <'./doc/readme.clim'`; then
  247.     echo shar: \"'./doc/readme.clim'\" unpacked with wrong size!
  248.   fi
  249.   # end of './doc/readme.clim'
  250. fi
  251. if test -f './doc/readme.globe' -a "${1}" != "-c" ; then 
  252.   echo shar: Will not clobber existing file \"'./doc/readme.globe'\"
  253. else
  254.   echo shar: Extracting \"'./doc/readme.globe'\" \(1735 characters\)
  255.   sed "s/^X//" >'./doc/readme.globe' <<'END_OF_FILE'
  256. XThis is the README file for Version 1 of globe, a program which takes
  257. Xa square map, wraps it around a globe, and displays an animation of the
  258. Xglobe rotating.  It is written in C and runs on any Xwindows system (X11R4).
  259. XCopyright (c) 1991 by David Allen.  You may distribute this freely as long
  260. Xas you leave my name and copyright notice intact.
  261. X
  262. XI'd really like your comments and feedback; send e-mail to
  263. Xallen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore Ave,
  264. XMaynard, MA 01754.
  265. X
  266. XWHAT DOES THE GLOBE PROGRAM DO?
  267. X
  268. XGiven a square map, for example, one produced by TEC, this program shows
  269. Xwhat the map would look like from space.  Since a square is flat, and
  270. Xa sphere is round, some distortion is involved.  
  271. X
  272. XHOW DO I RUN IT?
  273. X
  274. XTo compile globe, go to the src directory and type "make globe".  To run it,
  275. Xgo to the example directory and type "../src/globe globe.in".  This directs
  276. Xthe program to read the map stored in the file globe.in.
  277. X
  278. XWhen globe runs, it opens a graphics screen.  After a couple of seconds,
  279. Xa purple globe pops up.  There are 18 frames in the animation; the colors
  280. Xindicate altitudes.  Unlike tec, the ocean is displayed as bright purple
  281. Xrather than black.
  282. X
  283. XAt any time, you can type "q" into the graphics window to end the program.
  284. XIt will stop automatically after 10000 steps (a fairly long time).
  285. X
  286. XYou can display your own maps created using tec.  For example, create
  287. Xa tec input file with this single line:
  288. X
  289. X(MAXSTEP 20) (UNDERSCAN 10) (PRINTMODE 1)
  290. X
  291. XThen type "tec filename > tec.out" (where filename is the name of the
  292. Xinput file you just created).  Tec should run for a couple of seconds
  293. Xand then produce a file called tec.out.  Then type "globe tec.out" and
  294. Xyou will see your own globe rotating.
  295. END_OF_FILE
  296.   if test 1735 -ne `wc -c <'./doc/readme.globe'`; then
  297.     echo shar: \"'./doc/readme.globe'\" unpacked with wrong size!
  298.   fi
  299.   # end of './doc/readme.globe'
  300. fi
  301. if test -f './doc/readme.tec' -a "${1}" != "-c" ; then 
  302.   echo shar: Will not clobber existing file \"'./doc/readme.tec'\"
  303. else
  304.   echo shar: Extracting \"'./doc/readme.tec'\" \(3896 characters\)
  305.   sed "s/^X//" >'./doc/readme.tec' <<'END_OF_FILE'
  306. XThis is the README file for Version 5 of tec, a plate tectonic simulator.
  307. XIt is written in C and runs on any Xwindows system (X11R4) or any UN*X
  308. Xmachine using text or PostScript output.  Copyright (c) 1991 by David Allen.
  309. XYou may distribute this freely as long as you leave my name and copyright
  310. Xnotice intact.
  311. X
  312. XI'd really like your comments and feedback; send e-mail to
  313. Xallen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore Ave,
  314. XMaynard, MA 01754.
  315. X
  316. XVersion history:
  317. X1.  Posted to rec.games.frp and rec.games.programmer on 11 / 15 / 89
  318. X    - Runs only on Amiga, produces long text summary.
  319. X    - Includes PostScript translator for any PostScript printer.
  320. X2.  (Unreleased)
  321. X    - Added text-only mode so it runs on UN*X machines.
  322. X    - Added summary text output mode that fits on 80-character lines.
  323. X    - PostScript support merged into main program.
  324. X3.  Posted to comp.sources.misc, v10i077-78, on 2 / 15 / 90.
  325. X    - Parameters can come from an input file (see params.tec).
  326. X    - Ported to Sun workstation by the author.
  327. X    - Ported to IBM PC by Peter Lind.
  328. X4.  (Unreleased)
  329. X    - Ported to Xwindows by George Ferguson.
  330. X    - Autoscaling to produce specified hydrographics percentage.
  331. X    - Underscan mode to reduce drift "off the edge of the world".
  332. X5.  Posted to comp.sources.misc on 3 / 23 / 1991.
  333. X    - Merged with globe display and climate generator
  334. X    - Scope reduced to Xwindow graphics
  335. X
  336. XCredits:
  337. XBrian Love - conversion to Amiga Lattice C 5.02, Amiga beta tester.
  338. XPeter Lind - port to IBM PC Turbo C with CGA/MCGA/EGA/VGA.
  339. XAndrew Hay - improved random number seeding (now BSD independent).
  340. XFred Smith - added support for IBM PC Hercules mono graphics card.
  341. XHoward Johnson, Dave Pare - bug fixes.
  342. XGeorge Ferguson - port to XWindows.
  343. X
  344. X
  345. XWHAT DOES THE TECTONIC SIMULATOR DO?
  346. X
  347. XTEC draws a map of some imaginary world and puts one big continent on it.
  348. XThe program then graphically simulates the breakup of this supercontinent by
  349. Xrifts and the drifting of the continents.  Mountains are built by
  350. Xsubsuming and by continental collision, and are reduced by erosion.  As
  351. Xmore time passes, the continents tend to drift back into a supercontinent
  352. Xand split up again, forming what some authors call a "supercontinent cycle."
  353. X
  354. XThe program produces a map of an imaginary world which contains realistic
  355. Xlooking mountain ranges and continents which sometimes look like they fit
  356. Xtogether (like South America and Africa on Earth).
  357. X
  358. XI've included some sample output in the files example/tec.out.[1-3]; they are
  359. Xshort text files.  0's indicate ocean, 1's indicate land, and 2's indicate
  360. Xmountains.  Take a look.
  361. X
  362. X
  363. XHOW DO I RUN IT?
  364. X
  365. XTo compile tec, go to the src directory and type "make tec".  To run tec,
  366. Xtype "tec".  You can provide many optional parameter settings; place them
  367. Xin a file and run tec with "tec filename".  Output from tec goes to stdout;
  368. Xredirect it to a file to save it, as in "tec filename > outfile".
  369. X
  370. XWhen tec runs, it opens a graphics screen.  A big purple blob pops up.
  371. XThis is the initial supercontinent.  Every second or so (depending on
  372. Xmachine speed) another "step" will be computed and drawn.
  373. X
  374. XAltitudes are represented by colors; black is the lowest (ocean), and goes
  375. Xthrough the spectrum from lowest to highest: purple, blue, green, yellow,
  376. Xorange, red and finally white.  There are about 30 visibly different altitude
  377. Xranges, since each color has several shades.
  378. X
  379. XAs time passes, you will see the supercontinent split up; green and yellow
  380. Xmountains will appear on the leading edges of drifting continents; when
  381. Xcontinents slow down, reverse course, and collide, higher orange and red
  382. Xmountains will appear.  If you watch carefully, you can see erosion on
  383. Xmountains that are not growing; they turn darker and travel down through
  384. Xthe color ranges.
  385. X
  386. XAt any time, you can type "q" into the graphics window to end the program.
  387. XIt will stop automatically after 100 steps.
  388. END_OF_FILE
  389.   if test 3896 -ne `wc -c <'./doc/readme.tec'`; then
  390.     echo shar: \"'./doc/readme.tec'\" unpacked with wrong size!
  391.   fi
  392.   # end of './doc/readme.tec'
  393. fi
  394. if test -f './example/clim.in' -a "${1}" != "-c" ; then 
  395.   echo shar: Will not clobber existing file \"'./example/clim.in'\"
  396. else
  397.   echo shar: Extracting \"'./example/clim.in'\" \(1840 characters\)
  398.   sed "s/^X//" >'./example/clim.in' <<'END_OF_FILE'
  399. X(LAND (
  400. X000000000000000000000011111000000000011111111111100000000000
  401. X000000000000000000000111111100000000011111111111110000000000
  402. X000000000000000000001111111110000000011111111111110000000000
  403. X000000000000000000011111111110000000011111111111110000000000
  404. X000000000000000000111111111110000000011111111111110000000000
  405. X000000000000001111111122111110000000011111111111100000000000
  406. X000000000000011111111222111100000000011111111111100000000000
  407. X000000000000011111112222111100000000001111111111000000000000
  408. X000000000000011111222221111000000000001111111111000000000000
  409. X000000000000111112222211111000000000000111111111000000000000
  410. X000000000000111122222211110000000000000111111111000000000000
  411. X000000000000111222221111100000000000000111111111000000000000
  412. X000000000000111222211111000000000000000111111111000000000000
  413. X000000000000111222111110000000000000000111111110000000000000
  414. X000000000000111221111100000000000000111111111110000000000000
  415. X000000000000111211111100000000000001111111111110000000000000
  416. X000000000000111111111000000000000001111111111110000000000000
  417. X000000000000111111110000000000000011111111111100000000000000
  418. X000000000000111111110000000000000011111111111000000000000000
  419. X000000000000111111100000000000000111111111110000000000000000
  420. X000000000000111111000000000000001111111111110000000000000000
  421. X000000000000111111000000000000001111111111000000000000000000
  422. X000000000000000000000000000000011111111111000000000000000000
  423. X000000000000000000000000000000111111111111000000000000000000
  424. X000000000000000000000000000000111111111111000000000000000000
  425. X000000000000000000000000000001111111111110000000000000000000
  426. X000000000000000000000000000001111111111110000000000000000000
  427. X000000000000000000000000000001111111111100000000000000000000
  428. X000000000000000000000000000000111111111100000000000000000000
  429. X000000000000000000000000000000011111110000000000000000000000))
  430. END_OF_FILE
  431.   if test 1840 -ne `wc -c <'./example/clim.in'`; then
  432.     echo shar: \"'./example/clim.in'\" unpacked with wrong size!
  433.   fi
  434.   # end of './example/clim.in'
  435. fi
  436. if test -f './example/clim.out' -a "${1}" != "-c" ; then 
  437.   echo shar: Will not clobber existing file \"'./example/clim.out'\"
  438. else
  439.   echo shar: Extracting \"'./example/clim.out'\" \(1830 characters\)
  440.   sed "s/^X//" >'./example/clim.out' <<'END_OF_FILE'
  441. X333333333333333333333333333333333333333333333333333333333333
  442. X333333333333333333333333333333333333333333333333333333333333
  443. X333333333333333333333333333334444444444444433333333333333333
  444. X333333333333333333333333333300088885555555550000003333333333
  445. X000000000000000000000000000000008888555555555888000000000000
  446. X000000000000000080000000000000000088555555588888800000000000
  447. X000000000000885550000000000000000000005555555558850000000000
  448. X000000000000887667000000000000000000006655555555557000000000
  449. X000000000000877666780000000000000000006666666667667700000000
  450. X000000000000877666788000000000000000007666666667777880000000
  451. X000000000000876666778800000000000000007766667778887880000000
  452. X000000000000866666777888000000000000006666677778888800000000
  453. X000000000000666666777779900000000000007797900000009000000000
  454. X000000000000666666777777999000000000000000000000000000000000
  455. X000000000000777766777777799990000000000000000000000000000000
  456. X000000000000999777777777779999900000000000000000000000000000
  457. X000000000000A999999999999999999A0000000000000000000000000000
  458. X000000000000A99997777777777777900000000000000000000000000000
  459. X000000000000999979777677777770000000000000000000000000000000
  460. X000000000000988888766677667700000000000000000000000000000000
  461. X000000000000008888876666677800000000000000000000000000000000
  462. X000000000000000088887767677000000000000000000000000000000000
  463. X000000000000000000080007700000000000000000000000000000000000
  464. X000000000000000000000000000000000000000000000000000000000000
  465. X000000000000000000000000000000000000000000000000000000000000
  466. X000000000000000000000000000000000000000000000000000000000000
  467. X000000000000000000000000000000000000000000000000000000000000
  468. X333333333333333333333333333333333333333333333333333333333333
  469. X333333333333333333333333333333333333333333333333333333333333
  470. X333333333333333333333333333333333333333333333333333333333333
  471. END_OF_FILE
  472.   if test 1830 -ne `wc -c <'./example/clim.out'`; then
  473.     echo shar: \"'./example/clim.out'\" unpacked with wrong size!
  474.   fi
  475.   # end of './example/clim.out'
  476. fi
  477. if test -f './example/tec.out.1' -a "${1}" != "-c" ; then 
  478.   echo shar: Will not clobber existing file \"'./example/tec.out.1'\"
  479. else
  480.   echo shar: Extracting \"'./example/tec.out.1'\" \(3504 characters\)
  481.   sed "s/^X//" >'./example/tec.out.1' <<'END_OF_FILE'
  482. X000000000000000000000000000000000000222222221110000000000000000000000000
  483. X000000000000001111110000000000000001222111111100000000000000000000000000
  484. X000000000000011111111000000000000011211111111000000000000000000000000000
  485. X000000000001111111111111000000000111111111110000001110000000000000000000
  486. X000000000011111111111111100000000111111111100000011100000000000000000000
  487. X000000000011111111111111110000001111111111000000011110000000000000000000
  488. X000000000111111111111111111000011111111100000000111110000000000000000000
  489. X000000000111111111111111111100111111111000000001111110000000000000000000
  490. X000000000111111111111111111001111111110000000011111111000000000000000000
  491. X000000000111111111111111111111111111100000000111111111000000000000000000
  492. X000000000111111111111111111111111110000000011111111111000000000000000000
  493. X000000000111111111111111111111111100000000111111111111000000000000000000
  494. X000000000111111111112211111222110000000001111111111110000000000000000000
  495. X000000000111111111122222221111100000000011111221111111000000000000000000
  496. X000000000111111111122222111000000000001111122212221111100000000000000000
  497. X000000000011111111112222100000000000011112222212111111100000000000000000
  498. X000000000011111111111121000000000011111222222211111111110000000000000000
  499. X000000000011111111111010000000000111112222222211111111111000000000000000
  500. X000000000011111111100000000001111111222221112111111111111100000000000000
  501. X000000000011111111000000000112211222222211122111111111111111000000000000
  502. X000000000011111111000000001211111222221111111111111111111110000000000000
  503. X000000000000001110000000011111122222111111111111111111111110000000000000
  504. X000000000000000000000000111111222221111111111111111111100000000000000000
  505. X000000000000000000000001111212222111111111111101111111000000000000000000
  506. X000000000000000000000001222122221111111111111001111111000000000000000000
  507. X000000000000000000000001222222111111111111111101111110000000000000000000
  508. X000000000000000000000001222111111111111111111100110000000000000000000000
  509. X000000000000000000000001211111111111111111111000000000000000000000000000
  510. X000000000000000000000001111111111111111111111000000000000000000000000000
  511. X000000000000000000000000111111111111111111110000000000000000000000000000
  512. X000000000000000000000000111111111111111100000000000000000000000000000000
  513. X000000000000000000000000111111111111111000000000000000000000000000000000
  514. X000000000000000000000000011011111111100000000000000000000000000000000000
  515. X000000000000000000000000000011111111000000000000000000000000000000000000
  516. X000000000000000000000000000000110000000000000000000000000000000000000000
  517. X000000000000000000000000000000000000000000000000000000000000000000000000
  518. X000000000000000000000000000000000000000000000000000000000000000000000000
  519. X000000000000000000000000000000000000000000000000000000000000000000000000
  520. X000000000000000000000000000000000000000000000000000000000000000000000000
  521. X000000000000000000000000000000000000000000000000000000000000000000000000
  522. X000000000000000000000000000000000000000000000000000000000000000000000000
  523. X000000000000000000000000000000000000000000000000000000000000000000000000
  524. X000000000000000000000000000000000000000000000000000000000000000000000000
  525. X000000000000000000000000000000000000000000000000000000000000000000000000
  526. X000000000000000000000000000000000000000000000000000000000000000000000000
  527. X000000000000000000000000000000000000000000000000000000000000000000000000
  528. X000000000000000000000000000000000000000000000000000000000000000000000000
  529. X000000000000000000000000000000000000000000000000000000000000000000000000
  530. END_OF_FILE
  531.   if test 3504 -ne `wc -c <'./example/tec.out.1'`; then
  532.     echo shar: \"'./example/tec.out.1'\" unpacked with wrong size!
  533.   fi
  534.   # end of './example/tec.out.1'
  535. fi
  536. if test -f './example/tec.out.2' -a "${1}" != "-c" ; then 
  537.   echo shar: Will not clobber existing file \"'./example/tec.out.2'\"
  538. else
  539.   echo shar: Extracting \"'./example/tec.out.2'\" \(3504 characters\)
  540.   sed "s/^X//" >'./example/tec.out.2' <<'END_OF_FILE'
  541. X000000000000000000000000000000000000000000000000011111111111221000000000
  542. X000000000000000000000000000000000000000000000000001111111111222000000000
  543. X000000000000000000000000000011000000000000000000001111111111122100000000
  544. X000000000000000000000000011111110000000000000000000111111111111100000000
  545. X000000000000000000000000111111111100000000000000000111111111111100000000
  546. X000000000000000000000000011111111100000000000000000011111111111001100000
  547. X000000000000000000110000000011111100000000000000000011111111111011110000
  548. X000000000000000001110000000001111100000000000000000001111111211121111000
  549. X000000000000000001110110000000011110000000000000000001111111221121111000
  550. X000000000000000001111100000000000111000000000000000000111111221121111000
  551. X000000000000000001111110000000000111000001000000000000111111221121111000
  552. X000000000000000001111111110000000111000001000000000000011111221111111000
  553. X000000000000000001111111111000001111000001000000000000001111211112111000
  554. X000000000000000001211111111100011111000001000000000000001111221112111000
  555. X000000000000000001111111111110111111000001100000000000000111122112111000
  556. X000000000000000001111111111111111111000001110000000000000111122112111000
  557. X000000000000000001111111111111111111000001111000000000000011122112111000
  558. X000000000000000001111111111111111111000001111000000000000011112111210000
  559. X000000000000000000111111111111111111000001111100000000000001112211110000
  560. X000000000000000000111111111111111111000001111110000000000000112211110000
  561. X000000000000000000111111111111111111000001111110000000000000011221110000
  562. X000000000000000000111111111111111111000001111111000000000000001111100000
  563. X000000000000000000121111111111111111000001111111000000000000000111000000
  564. X000000000000000000011111111111111111000001111111000000000000000010000000
  565. X000000000000000000001111111111111111000001111111100000000000000000000000
  566. X000000000000000000000000121111111111000001111111110000000000000000000000
  567. X000000000000000000000000111111111110000001111111111000000000000000000000
  568. X000000000000000000000000121111111100000001111111111100000000000000000000
  569. X000000000000000000000000011111111000000011111111111110000000000000000000
  570. X000000000000000000000000001111111000000011111111111110000000000000000000
  571. X000000000000000011111100000111111000000011111111111111000000000000000000
  572. X000000000000000112222110000000011000000011111111111111000000000000000000
  573. X000000000001011222222210000000001000000011111111111111100000000000000000
  574. X000000000011112221112211000000000000000001111111111111110000000000000000
  575. X000000000012222111111121000000000000000001111111111111110000000000000000
  576. X000000000012211111111121110000000000000001111111111111100000000000000000
  577. X000000000011111111111112211111000000000001111111111111000000000000000000
  578. X000000000001111111111112222221110000000001111111111110000000000000000000
  579. X000000000011111111111112222222221100000000000000000000000000000000000000
  580. X000000000011111111112222222222222110000000000000000000000000000000000000
  581. X000000000011111111112222222222222210000000000000000000000000000000000000
  582. X000000000001111001111111111111111110000000000000000000000000000000000000
  583. X000000000000000000000000000000000000000000000000000000000000000000000000
  584. X000000000000000000000000000000000000000000000000000000000000000000000000
  585. X000000000000000000000000000000000000000000000000000000000000000000000000
  586. X000000000000000000000000000000000000000000000000000000000000000000000000
  587. X000000000000000000000000000000000000000000000000000000000000000000000000
  588. X000000000000000000000000000000000000000000000000000000000000000000000000
  589. END_OF_FILE
  590.   if test 3504 -ne `wc -c <'./example/tec.out.2'`; then
  591.     echo shar: \"'./example/tec.out.2'\" unpacked with wrong size!
  592.   fi
  593.   # end of './example/tec.out.2'
  594. fi
  595. if test -f './example/tec.out.3' -a "${1}" != "-c" ; then 
  596.   echo shar: Will not clobber existing file \"'./example/tec.out.3'\"
  597. else
  598.   echo shar: Extracting \"'./example/tec.out.3'\" \(3504 characters\)
  599.   sed "s/^X//" >'./example/tec.out.3' <<'END_OF_FILE'
  600. X000000000000000000000000000000000000000000000000000000000000000000000000
  601. X000000000000000000000000000000000000000000000000000000000000000000000000
  602. X000000000000000000000000000000000000000000000000000000000000000000000000
  603. X000000000000000000000000000000000000000000000000000000000000000000000000
  604. X000000000000000000000000000001110000000000000000000011000000000000000000
  605. X000000000000000000000000000001101111111111111111111111000000000000000000
  606. X000000000000000000000000000000001111111111111111111111000000000000000000
  607. X000000000000000000000011000000000111111111111111111111000000000000000000
  608. X000000000000000000000011100000000001111111111111111111110000000000000000
  609. X000000000000000000000011111100000000011111111111111111110000000000000000
  610. X000000000000000000000011111100000000000111111111111111100000000000000000
  611. X000000000000000000000011111110000000000011111111111111100000000000000000
  612. X000000000000000000000011111111110000000001111111111110111111110000000000
  613. X000000000000000000000011111111111000000001111111111101111111111000000000
  614. X000000000000000000000011111111111110000000111111000001111111111110000000
  615. X000000000000000000000011111111111111000000011110000011111111111110000000
  616. X000000000000000000000011111111111110000000111110000111111111111000000000
  617. X000000000000000000000011111111111111100001111110001111111111110000000000
  618. X000000000000000000000011111111111111110011111111111111111111110000000000
  619. X000000000000000000000001111111111111111111111111111111111111111100000000
  620. X000000000000000000000001111111111111111111111111111111111111111110000000
  621. X000000000000000000000001111111111111111111111111111111111111111110000000
  622. X000000000000000000000001111111112221111111111111111111111111111110000000
  623. X000000000000000000000001111111112221111111111111111111111000010000000000
  624. X000000000000000000000000111111112211111111111111111111110000000000000000
  625. X000000000000000000000000000011122111111111111111111111111000000000000000
  626. X000000000000000000000000000111121111111111111111111111111100000000000000
  627. X000000000000000000000000001111111111111011111111111111111100000000000000
  628. X000000000000000000000000011112111111110011111111111111111100000000000000
  629. X000000000000000111111111111122111111110011111111111111111000000000000000
  630. X000000000000000111111111111111111111110011111111111111110000000000000000
  631. X000000000000000111111112222111111111110011111111111111000000000000000000
  632. X000000000000000111111222211111111111110011111111111111100000000000000000
  633. X000000000000000111111221111111111111111001111111111111110000000000000000
  634. X000000000000000111111111111111111111111101111111111111100000000000000000
  635. X000000000000000111111111111111111111111101111111111111000000000000000000
  636. X000000000000000111111111111111111111111001111111111110000000000000000000
  637. X000000000000000111111111111000111111100001111101101100000000000000000000
  638. X000000000000000111111111000000000000000000000000000000000000000000000000
  639. X000000000000000001111110000000000000000000000000000000000000000000000000
  640. X000000000000000000111000000000000000000000000000000000000000000000000000
  641. X000000000000000000010000000000000000000000000000000000000000000000000000
  642. X000000000000000000000000000000000000000000000000000000000000000000000000
  643. X000000000000000000000000000000000000000000000000000000000000000000000000
  644. X000000000000000000000000000000000000000000000000000000000000000000000000
  645. X000000000000000000000000000000000000000000000000000000000000000000000000
  646. X000000000000000000000000000000000000000000000000000000000000000000000000
  647. X000000000000000000000000000000000000000000000000000000000000000000000000
  648. END_OF_FILE
  649.   if test 3504 -ne `wc -c <'./example/tec.out.3'`; then
  650.     echo shar: \"'./example/tec.out.3'\" unpacked with wrong size!
  651.   fi
  652.   # end of './example/tec.out.3'
  653. fi
  654. if test -f './src/clim.h' -a "${1}" != "-c" ; then 
  655.   echo shar: Will not clobber existing file \"'./src/clim.h'\"
  656. else
  657.   echo shar: Extracting \"'./src/clim.h'\" \(675 characters\)
  658.   sed "s/^X//" >'./src/clim.h' <<'END_OF_FILE'
  659. X#define MAXB       4
  660. X#define LINE_DIAG  1
  661. X#define LINE_CORN  2
  662. X#define LINE_0V    1
  663. X#define LINE_1V    2
  664. X#define LINE_0H    4
  665. X#define LINE_1H    8
  666. X#define N          LINE_0V
  667. X#define S          LINE_1V
  668. X#define E          LINE_0H
  669. X#define W          LINE_1H
  670. X#define PR_LOW     1
  671. X#define PR_HIGH    2
  672. X#define PR_HEQ     3
  673. X#define M_MAIN     0
  674. X#define M_HEAT     1
  675. X#define M_PRESS    2
  676. X#define M_WIND     3
  677. X#define M_RAIN     4
  678. X#define M_CLIM     5
  679. X#define C_OCEAN    0
  680. X#define C_BARELAND 1
  681. X#define C_MOUNTAIN 2
  682. X#define C_OCEANICE 3
  683. X#define C_TUNDRA   4
  684. X#define C_STEPPE   5
  685. X#define C_DESERT   6
  686. X#define C_SAVANA   7
  687. X#define C_DECID    8
  688. X#define C_JUNGLE   9
  689. X#define C_SWAMP   10
  690. END_OF_FILE
  691.   if test 675 -ne `wc -c <'./src/clim.h'`; then
  692.     echo shar: \"'./src/clim.h'\" unpacked with wrong size!
  693.   fi
  694.   # end of './src/clim.h'
  695. fi
  696. if test -f './src/climate.c' -a "${1}" != "-c" ; then 
  697.   echo shar: Will not clobber existing file \"'./src/climate.c'\"
  698. else
  699.   echo shar: Extracting \"'./src/climate.c'\" \(4025 characters\)
  700.   sed "s/^X//" >'./src/climate.c' <<'END_OF_FILE'
  701. X/* This program is Copyright (c) 1991 David Allen.  It may be freely
  702. X   distributed as long as you leave my name and copyright notice on it.
  703. X   I'd really like your comments and feedback; send e-mail to
  704. X   allen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore
  705. X   Avenue, Maynard, MA 01754. */
  706. X
  707. X/* This file contains the routines that compute local climate */
  708. X
  709. X#include "const.h"
  710. X#include "clim.h"
  711. X
  712. X
  713. X/* Input arrays are l from main.c, rn from rain.c, and t from heat.c.  Note
  714. X   that the unscaled temperatures are used.  The output array is cl.  Array
  715. X   lm (from main.c) is used by the drawing routine.*/
  716. X
  717. Xextern unsigned char l[MAXX][MAXY], lm[MAXX][MAXY], rn[MAXB][MAXX][MAXY];
  718. Xextern int tt[MAXB][MAXX][MAXY];
  719. Xunsigned char cl[MAXX][MAXY];
  720. X
  721. X
  722. X/* These parameters are defined in main.c or heat.c */
  723. Xextern int BSIZE, XSIZE, YSIZE, TEMPSCALE, PRINTMODE;
  724. X
  725. X/* These parameters are used here, and are described in params.doc */
  726. Xint tempcut[] = { 0, 40, 90, 120 }, raincut[] = { 40, 60, 110, 160, 180 };
  727. Xint PRINTCLIM = 0, MTDELTA = 20, TCSIZE = 4, RCSIZE = 5;
  728. Xdouble ICEBERGK = 263.0;
  729. X
  730. X
  731. X/* This array is the heart of the climate routine; temperature increases
  732. X   going down the array, and rainfall increases going from left to right. */
  733. X
  734. Xstatic unsigned char climkey[4][5] = {
  735. X   { C_TUNDRA, C_TUNDRA, C_TUNDRA, C_TUNDRA, C_TUNDRA },
  736. X   { C_STEPPE, C_STEPPE, C_DECID,  C_DECID,  C_DECID  },
  737. X   { C_DESERT, C_SAVANA, C_DECID,  C_JUNGLE, C_SWAMP  },
  738. X   { C_DESERT, C_SAVANA, C_JUNGLE, C_SWAMP,  C_SWAMP  } };
  739. X
  740. X
  741. Xclimpar (s) char *s; {
  742. X   /* This function is called by mainpar() in main.c; it simply tests input
  743. X   parameter names to see if they are defined in this file.  Each of the
  744. X   above ints are defined in this file.  If the input string matches here,
  745. X   the function returns true. */
  746. X   if      (CMP ("ICEBERGK"))    getdbl  (&ICEBERGK,        M_CLIM);
  747. X   else if (CMP ("PRINTCLIM"))   getlng  (&PRINTCLIM,       M_CLIM);
  748. X   else if (CMP ("TEMPCUT"))     getlvec (&TCSIZE, tempcut, M_CLIM);
  749. X   else if (CMP ("RAINCUT"))     getlvec (&RCSIZE, raincut, M_CLIM);
  750. X   else if (CMP ("MTDELTA"))     getlng  (&MTDELTA,         M_CLIM);
  751. X   else if (CMP ("PRINTCLIM"))   getlng  (&PRINTCLIM,       M_CLIM);
  752. X   else return (0);
  753. X   return (1); }
  754. X
  755. X
  756. Xclimdraw () { draw (DRAW_CLIM, LINE_CORN, cl, lm); }
  757. X   /* This routine calls draw with the proper arguments to display climate */
  758. X
  759. X
  760. Xclimcomp () {
  761. X   /* The outer loop looks at each square.  If it is ocean, the climate will
  762. X   be ocean unless the temperature is below ICEBREGK degrees all year round.
  763. X   If it is land, then the average rainfall and temperature (in Farenheit) are
  764. X   computed for the square.  If the square is mountain, it is colder; the
  765. X   temperature is decreased.  These two figures are then turned into array
  766. X   indices by using the tempcut and raincut parameter vectors.  The climate
  767. X   for the square is then simply a table lookup.  Finally, the array is printed
  768. X   if desired. */
  769. X
  770. X   register int i, j, buf;
  771. X   int noice, avetemp, averain, ttt, r;
  772. X
  773. X   usermessage ("Computing climate"); checkmouse ();
  774. X   for (j=0; j<YSIZE; j++) for (i=0; i<XSIZE; i++) {
  775. X      if (!l[i][j]) { /* ocean */
  776. X         for (noice=0, buf=0; buf<BSIZE; buf++)
  777. X            noice |= (tt[buf][i][j] > TEMPSCALE * ICEBERGK);
  778. X         cl[i][j] = noice ? C_OCEAN : C_OCEANICE; }
  779. X      else { /* land or mountain */
  780. X         for (averain=0, avetemp=0, buf=0; buf<BSIZE; buf++) {
  781. X            averain += rn[buf][i][j]; avetemp += tt[buf][i][j]; }
  782. X         averain /= BSIZE; avetemp /= BSIZE;
  783. X         avetemp = ((double) (avetemp / TEMPSCALE) - 273.0) * 1.8 + 32.0;
  784. X         if (l[i][j] == 2) avetemp -= MTDELTA;
  785. X         ttt = 0; while ((avetemp > tempcut[ttt]) && (ttt < TCSIZE-1)) ttt++;
  786. X         r  = 0; while ((averain > raincut[r])  && (r  < RCSIZE-1)) r++;
  787. X         cl[i][j] = climkey[ttt][r]; } }
  788. X   if (PRINTCLIM) {
  789. X      if (PRINTMODE == PRINTMODE_GREY)
  790. X         putmat ("CLIMATE", -1, PRINTMODE_CLIM, cl, lm);
  791. X      else putmat ("CLIMATE", -1, PRINTMODE_SHORT, cl, lm); } }
  792. END_OF_FILE
  793.   if test 4025 -ne `wc -c <'./src/climate.c'`; then
  794.     echo shar: \"'./src/climate.c'\" unpacked with wrong size!
  795.   fi
  796.   # end of './src/climate.c'
  797. fi
  798. if test -f './src/const.h' -a "${1}" != "-c" ; then 
  799.   echo shar: Will not clobber existing file \"'./src/const.h'\"
  800. else
  801.   echo shar: Extracting \"'./src/const.h'\" \(467 characters\)
  802.   sed "s/^X//" >'./src/const.h' <<'END_OF_FILE'
  803. X#define MAXX 100
  804. X#define MAXY 100
  805. X#define PRINTMODE_NONE  0
  806. X#define PRINTMODE_LONG  1
  807. X#define PRINTMODE_SCALE 2
  808. X#define PRINTMODE_SHORT 3
  809. X#define PRINTMODE_GREY  4
  810. X#define PRINTMODE_CLIM  5
  811. X#define DRAW_GREY  1
  812. X#define DRAW_LAND  2
  813. X#define DRAW_CLIM  3
  814. X#define DRAW_TEC   4
  815. X#define LINE_DIAG  1
  816. X#define LINE_CORN  2
  817. X#define LINE_NONE  3
  818. X
  819. X#define ABS(xx) ((xx < 0) ? -xx: xx)
  820. X#define CMP(x) !strcmp (s, x)
  821. X
  822. Xextern double cos(), sin(), sqrt(), atof(), asin(), atan2();
  823. END_OF_FILE
  824.   if test 467 -ne `wc -c <'./src/const.h'`; then
  825.     echo shar: \"'./src/const.h'\" unpacked with wrong size!
  826.   fi
  827.   # end of './src/const.h'
  828. fi
  829. if test -f './src/globe.c' -a "${1}" != "-c" ; then 
  830.   echo shar: Will not clobber existing file \"'./src/globe.c'\"
  831. else
  832.   echo shar: Extracting \"'./src/globe.c'\" \(3176 characters\)
  833.   sed "s/^X//" >'./src/globe.c' <<'END_OF_FILE'
  834. X/* This program is Copyright (c) 1991 David Allen.  It may be freely
  835. X   distributed as long as you leave my name and copyright notice on it.
  836. X   I'd really like your comments and feedback; send e-mail to
  837. X   allen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore
  838. X   Avenue, Maynard, MA 01754. */
  839. X
  840. X/* This file contains a spinning globe animation */
  841. X
  842. X#include "const.h"
  843. X#define N ((MAXX / 2) - 1)
  844. X#define AREA (MAXX * MAXX)
  845. X#define RADIUS 45
  846. X#define PI 3.14159
  847. X#define TWOPI (2.0 * PI)
  848. X#define HALFPI (PI / 2.0)
  849. X#define SET(n,x,y,i,j) m[n][(x)+N][(y)+N] = m[0][(i)+N][(j)+N]
  850. X
  851. Xunsigned char m[21][MAXX][MAXY];
  852. Xint rx[] = {0, 0, 0, 0, 0}, ry[] = {0, N, 0, N, 0};
  853. Xint xofs, prime, xflip, sflip, sxflip, change[1], step = 0;
  854. Xint MAXSTEP = 10000, PRINTMODE = PRINTMODE_SCALE, XSIZE = MAXX, YSIZE = MAXX;
  855. Xint ZCOAST = 0;
  856. X
  857. X
  858. Xinit (s) char *s; { int i;
  859. X   fileinit ();
  860. X   if (!s || !(*s)) panic ("no valid filename"); else getparams (s);
  861. X   for (i=1, xofs=45; i<6; xofs-=10, i++) {
  862. X      prime = 16-i; xflip = i+5;  sflip = 15+i; sxflip = 6-i; compute (); } }
  863. X
  864. X
  865. Xonestep () { static int i = 0;
  866. X   draw (DRAW_TEC, LINE_NONE, m[i], 0); i++; step++;
  867. X   if (i == 5) i = 6; if (i == 16) i = 17; if (i == 21) i = 1; }
  868. X
  869. X
  870. Xmainpar (s) char *s; { int vec [AREA], index = AREA; register int i, j, k;
  871. X
  872. X   if (CMP ("LAND")) {
  873. X      getlvec (&index, vec, 0);
  874. X      for (k=0, j=0; j<MAXX; j++) for (i=0; i<MAXX; i++, k++)
  875. X         m[0][i][j] = (vec[k] > 16) ? vec[k] : 16;
  876. X      return (1); }
  877. X   return (0); }
  878. X
  879. X
  880. Xcompute () {
  881. X   register int i, j, top; int x, y;
  882. X
  883. X   rx[0] = xofs + N; rx[4] = rx[0]; rx[2] = xofs - N;
  884. X   for (j=0; j<RADIUS; j++) {
  885. X      checkmouse ();
  886. X      top = sqrt ((double) (RADIUS * RADIUS - j * j + 0.5));
  887. X      for (i=1-top; i<top; i++) {
  888. X         point (i, j, &x, &y);
  889. X         SET (prime,   i,  j,  x,  y); SET (prime,   i, -j,  x, -y);
  890. X         SET (xflip,  -i,  j, -x,  y); SET (xflip,  -i, -j, -x, -y);
  891. X         permute (x, y, &x, &y);
  892. X         SET (sxflip,  i, -j, -x,  y); SET (sxflip,  i,  j, -x, -y);
  893. X         SET (sflip,  -i, -j,  x,  y); SET (sflip,  -i,  j,  x, -y); } } }
  894. X
  895. X
  896. Xpoint (x, y, i, j) int x, y, *i, *j; {
  897. X   double r, theta, topx, topy, botx, boty; int bot, top;
  898. X
  899. X   r = asin (sqrt ((double)(x * x + y * y)) / RADIUS) * 2.0 / PI;
  900. X   if (x || y) theta = atan2 ((double)y, (double)x);
  901. X   else theta = 0;
  902. X   if (theta <= 0) theta += TWOPI;
  903. X   bot = 0; while (theta > HALFPI) { bot++; theta -= HALFPI; }
  904. X   top = bot + 1; theta = theta / HALFPI;
  905. X   mapr (bot, r, &botx, &boty); mapr (top, r, &topx, &topy);
  906. X   *i = botx + theta * (topx - botx) + 0.5;
  907. X   *j = boty + theta * (topy - boty) + 0.5;
  908. X   if (*i >= N) *i = N-1; if (*j >= N) *j = N-1; }
  909. X
  910. X
  911. Xmapr (n, r, x, y) int n; double r, *x, *y; {
  912. X   *x = r * (rx[n] - xofs) + xofs;
  913. X   *y = r * ry[n];
  914. X   if (*x >= N) { *y = *x - N; *x = N; if (n == 4) *y = -(*y); } }
  915. X
  916. X
  917. Xpermute (x, y, i, j) int x, y, *i, *j; {
  918. X   if ((x >= 0) && (y >= 0)) { *i =  N - y - 1; *j =  x - N + 0; }
  919. X   if ((x <  0) && (y >= 0)) { *i =  y - N + 0; *j = -N - x - 1; }
  920. X   if ((x >= 0) && (y <  0)) { *i =  y + N + 0; *j =  N - x - 1; }
  921. X   if ((x <  0) && (y <  0)) { *i = -N - y - 1; *j =  x + N + 0; } }
  922. X
  923. X
  924. Xgreyscale (x) int x; { }
  925. END_OF_FILE
  926.   if test 3176 -ne `wc -c <'./src/globe.c'`; then
  927.     echo shar: \"'./src/globe.c'\" unpacked with wrong size!
  928.   fi
  929.   # end of './src/globe.c'
  930. fi
  931. if test -f './src/tec.h' -a "${1}" != "-c" ; then 
  932.   echo shar: Will not clobber existing file \"'./src/tec.h'\"
  933. else
  934.   echo shar: Extracting \"'./src/tec.h'\" \(531 characters\)
  935.   sed "s/^X//" >'./src/tec.h' <<'END_OF_FILE'
  936. X#define MAXSPLINTER 10
  937. X#define MAXFRAG 100
  938. X#define MAXPLATE 40
  939. X#define MAXCHANGE 10
  940. X#define REALSCALE 100
  941. X
  942. Xstruct plate { int dx, dy, odx, ody, rx, ry, age, area, id, next; };
  943. X
  944. Xextern int XSIZE,    YSIZE,        MAXSTEP;
  945. Xextern int MAXLIFE,    MAXBUMP,    BUMPTOL;
  946. Xextern int DRAWEVERY,    HYDROPCT,    PRINTMODE;
  947. Xextern int ZINIT,    ZSUBSUME,    ZCOAST;
  948. Xextern int ZSHELF,    ZMOUNTAIN,    ZERODE;
  949. Xextern int RIFTPCT,    DOERODE,    ERODERND;
  950. Xextern int MAXCTRTRY,    RIFTDIST,    BENDEVERY;
  951. Xextern int BENDBY,    SPEEDBASE,    SPEEDRNG;
  952. Xextern int UNDERSCAN;
  953. Xextern double MR[];
  954. END_OF_FILE
  955.   if test 531 -ne `wc -c <'./src/tec.h'`; then
  956.     echo shar: \"'./src/tec.h'\" unpacked with wrong size!
  957.   fi
  958.   # end of './src/tec.h'
  959. fi
  960. if test -f './src/unix.c' -a "${1}" != "-c" ; then 
  961.   echo shar: Will not clobber existing file \"'./src/unix.c'\"
  962. else
  963.   echo shar: Extracting \"'./src/unix.c'\" \(862 characters\)
  964.   sed "s/^X//" >'./src/unix.c' <<'END_OF_FILE'
  965. X/* This program is Copyright (c) 1991 David Allen.  It may be freely
  966. X   distributed as long as you leave my name and copyright notice on it.
  967. X   I'd really like your comments and feedback; send e-mail to
  968. X   allen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore Ave,
  969. X   Maynard, MA 01754. */
  970. X
  971. X#include "const.h"
  972. Xextern int XSIZE, YSIZE, MAXSTEP, step;
  973. Xint picktype;
  974. X
  975. Xmain (argc,argv) int argc; char **argv; {
  976. X   unsigned short seed16v[3];
  977. X
  978. X   /* Initialize random number generator */
  979. X   srand (time ((long *) 0));    /* initialize rand() */
  980. X   seed16v[0] = rand(); seed16v[1] = rand(); seed16v[2] = rand();
  981. X   seed48 (seed16v);    /* initialize lrand48() */
  982. X
  983. X   init (*++argv);
  984. X   for (step=0; step<MAXSTEP; step++) onestep();
  985. X   return (0); }
  986. X
  987. Xcheckmouse () { }
  988. X 
  989. Xpanic (s) char *s; { printf ("PANIC: %s\n", s); exit (1); }
  990. X
  991. Xdraw (ctype, ltype, cra, lra) { }
  992. END_OF_FILE
  993.   if test 862 -ne `wc -c <'./src/unix.c'`; then
  994.     echo shar: \"'./src/unix.c'\" unpacked with wrong size!
  995.   fi
  996.   # end of './src/unix.c'
  997. fi
  998. if test -f './src/wind.c' -a "${1}" != "-c" ; then 
  999.   echo shar: Will not clobber existing file \"'./src/wind.c'\"
  1000. else
  1001.   echo shar: Extracting \"'./src/wind.c'\" \(5268 characters\)
  1002.   sed "s/^X//" >'./src/wind.c' <<'END_OF_FILE'
  1003. X/* This program is Copyright (c) 1991 David Allen.  It may be freely
  1004. X   distributed as long as you leave my name and copyright notice on it.
  1005. X   I'd really like your comments and feedback; send e-mail to
  1006. X   allen@viewlogic.com, or send us-mail to David Allen, 10 O'Moore
  1007. X   Avenue, Maynard, MA 01754. */
  1008. X
  1009. X/* This file contains the wind computation routines. */
  1010. X
  1011. X#include "const.h"
  1012. X#include "clim.h"
  1013. X#define MAXPRESS 255
  1014. X
  1015. X/* Input arrays are pr, from pressure.c, and l, from main.c.  The array
  1016. X   hl is used as temporary storage for highs and lows, while p is used
  1017. X   to store the smoothed pressure map.  The output array, wd, contains
  1018. X   an edge map with wind directions. */
  1019. X
  1020. Xextern unsigned char pr[MAXB][MAXX][MAXY], l[MAXX][MAXY];
  1021. Xunsigned char wd[MAXB][MAXX][MAXY];
  1022. Xstatic char hl[2][MAXX][MAXY];
  1023. Xstatic unsigned char p[MAXX][MAXY];
  1024. X
  1025. X/* The externs below are parameters defined in main.c.  The other two ints
  1026. X   are parameters described in params.doc. */
  1027. X
  1028. Xextern int XSIZE, YSIZE, BSIZE, PRINTMODE;
  1029. Xint BARSEP = 16, PRINTWIND = 0;
  1030. X
  1031. X
  1032. Xwindpar (s) char *s; {
  1033. X   /* This function is called by mainpar() in main.c; it simply checks to
  1034. X   see if the parameter is defined in this file.  If so, it returns true. */
  1035. X   if (CMP ("BARSEP")) getlng (&BARSEP, M_WIND);
  1036. X   else if (CMP ("PRINTWIND")) getlng (&PRINTWIND, M_WIND);
  1037. X   else return (0);
  1038. X   return (1); }
  1039. X
  1040. X
  1041. Xwindcomp () {
  1042. X   /* This is the main function in this file; it calls getpress() to create
  1043. X   a smoothed pressure map, then getwind() to put isobars (wind lines) on
  1044. X   the output map.  The last step makes sure that contradictory winds are
  1045. X   removed, such as N and S winds in the same square. */
  1046. X
  1047. X   register int i, j, x, buf;
  1048. X
  1049. X   for (buf=0; buf<BSIZE; buf++) {
  1050. X      status (M_WIND, buf); getpress (buf); getwind (buf);
  1051. X      for (j=0; j<YSIZE; j++) for (i=0; i<XSIZE; i++) {
  1052. X         x = wd[buf][i][j];
  1053. X         if (x & N) x &= (~S); if (x & E) x &= (~W);
  1054. X         wd[buf][i][j] = x; } }
  1055. X
  1056. X   if (PRINTWIND) for (buf=0; buf<BSIZE; buf++) {
  1057. X      if (PRINTMODE == PRINTMODE_GREY) 
  1058. X         putmat ("WIND", buf, PRINTMODE_SHORT, l, wd[buf]);
  1059. X      else putmat ("WIND", buf, PRINTMODE_SHORT, wd[buf], 0); } }
  1060. X
  1061. X
  1062. Xwindraw (n) int n; { draw (DRAW_LAND, LINE_DIAG, l, wd[n]); }
  1063. X   /* This function calls draw with the right arguments to display wind */
  1064. X
  1065. X
  1066. Xgetpress (buf) int buf; {
  1067. X   /* This function takes the high and low markings from pressure.c and creates
  1068. X   a smoothed function.  Highs turn into MAXPRESS and lows turn into 0. */
  1069. X
  1070. X   register int i, j;
  1071. X
  1072. X   for (j=0; j<YSIZE; j++) for (i=0; i<XSIZE; i++) {
  1073. X      /* Zero out the arrays to be used */
  1074. X      wd[buf][i][j] = 0; hl[0][i][j] = 0; hl[1][i][j] = 0;
  1075. X
  1076. X      /* Fill hl[0] with the low pressure zones, and hl[1] with highs */
  1077. X      if      (pr[buf][i][j] == PR_LOW)  hl[0][i][j] = -1;
  1078. X      else if (pr[buf][i][j] == PR_HIGH) hl[1][i][j] = -1;
  1079. X      else if (pr[buf][i][j] == PR_HEQ)  hl[0][i][j] = -1; }
  1080. X
  1081. X   /* Set each square in hl[0] to the distance from that square to the */
  1082. X   /* nearest low, and each square in hl[1] to the distance to a high. */
  1083. X   range (hl[0]); range (hl[1]);
  1084. X
  1085. X   /* The final pressure, in array p, is zero if a low is there and */
  1086. X   /* MAXPRESS if a high is there.  Otherwise, the pressure in a square is */
  1087. X   /* proportional to the ratio of (distance from the square to the nearest */
  1088. X   /* low) to (total of distance from nearest high and nearest low).  This */
  1089. X   /* gives a smooth curve between the extremes. */
  1090. X   for (j=0; j<YSIZE; j++) for (i=0; i<XSIZE; i++) {
  1091. X      if (hl[1][i][j] == -1) p[i][j] = MAXPRESS;
  1092. X      else if (hl[0][i][j] == -1) p[i][j] = 0;
  1093. X      else p[i][j] = (MAXPRESS*hl[0][i][j]) / (hl[0][i][j] + hl[1][i][j]); } }
  1094. X
  1095. X
  1096. Xgetwind (buf) int buf; {
  1097. X   /* This function draws isobars around the pressure map created above.  These
  1098. X   isobars are the directions of wind flow.  The isobars are given a direction
  1099. X   depending on whether the square is above or below the heat equator; north of
  1100. X   the heat equator, the winds blow counterclockwise out from a low, while
  1101. X   south of it, the opposite is true. */
  1102. X
  1103. X   register int i, j; int a, b, e, bar;
  1104. X
  1105. X   /* Step from 0 to MAXPRESS by BARSEP; bar is the pressure for which this */
  1106. X   /* isobar will be drawn. */
  1107. X   for (bar=BARSEP; bar<=MAXPRESS; bar+=BARSEP) {
  1108. X      checkmouse ();
  1109. X      for (i=0; i<XSIZE; i++) for (e=0, j=0; j<YSIZE; j++) {
  1110. X
  1111. X         /* Set e if this square is south of the heat equator */
  1112. X         a = p[i][j]; if (pr[buf][i][j] == 3) e = 1;
  1113. X
  1114. X         /* Provided the square is not at the top of the array, compare the */
  1115. X         /* pressure here to the pressure one square up.  This gives the */
  1116. X         /* direction of the wind in terms of east / west flow. */
  1117. X         if (j) {
  1118. X            b = p[i][j-1];
  1119. X            if ((a < bar) && (b >= bar)) wd[buf][i][j] |= (e ? E : W);
  1120. X            if ((a >= bar) && (b < bar)) wd[buf][i][j] |= (e ? W : E); }
  1121. X
  1122. X         /* Compare the pressure here to the pressure one square to the left */
  1123. X         /* (including wraparound); this gives the wind direction in terms */
  1124. X         /* of north / south flow. */
  1125. X         b = i ? p[i-1][j] : p[XSIZE-1][j];
  1126. X         if ((a < bar) && (b >= bar)) wd[buf][i][j] |= (e ? N : S);
  1127. X         if ((a >= bar) && (b < bar)) wd[buf][i][j] |= (e ? S : N); } } }
  1128. END_OF_FILE
  1129.   if test 5268 -ne `wc -c <'./src/wind.c'`; then
  1130.     echo shar: \"'./src/wind.c'\" unpacked with wrong size!
  1131.   fi
  1132.   # end of './src/wind.c'
  1133. fi
  1134. echo shar: End of archive 4 \(of 4\).
  1135. cp /dev/null ark4isdone
  1136. MISSING=""
  1137. for I in 1 2 3 4 ; do
  1138.     if test ! -f ark${I}isdone ; then
  1139.     MISSING="${MISSING} ${I}"
  1140.     fi
  1141. done
  1142. if test "${MISSING}" = "" ; then
  1143.     echo You have unpacked all 4 archives.
  1144.     rm -f ark[1-9]isdone
  1145. else
  1146.     echo You still must unpack the following archives:
  1147.     echo "        " ${MISSING}
  1148. fi
  1149. exit 0
  1150. exit 0 # Just in case...
  1151. -- 
  1152. Kent Landfield                   INTERNET: kent@sparky.IMD.Sterling.COM
  1153. Sterling Software, IMD           UUCP:     uunet!sparky!kent
  1154. Phone:    (402) 291-8300         FAX:      (402) 291-4362
  1155. Please send comp.sources.misc-related mail to kent@uunet.uu.net.
  1156.