| 1 |
%!PS-Adobe-3.0 |
| 2 |
%%BoundingBox: 38 24 574 768 |
| 3 |
%%Title: Enscript Output |
| 4 |
%%For: Andrea Molod |
| 5 |
%%Creator: GNU enscript 1.6.4 |
| 6 |
%%CreationDate: Mon Jun 4 17:11:38 2007 |
| 7 |
%%Orientation: Landscape |
| 8 |
%%Pages: (atend) |
| 9 |
%%DocumentMedia: Letter 612 792 0 () () |
| 10 |
%%DocumentNeededResources: (atend) |
| 11 |
%%EndComments |
| 12 |
%%BeginProlog |
| 13 |
%%BeginResource: procset Enscript-Prolog 1.6 4 |
| 14 |
% |
| 15 |
% Procedures. |
| 16 |
% |
| 17 |
|
| 18 |
/_S { % save current state |
| 19 |
/_s save def |
| 20 |
} def |
| 21 |
/_R { % restore from saved state |
| 22 |
_s restore |
| 23 |
} def |
| 24 |
|
| 25 |
/S { % showpage protecting gstate |
| 26 |
gsave |
| 27 |
showpage |
| 28 |
grestore |
| 29 |
} bind def |
| 30 |
|
| 31 |
/MF { % fontname newfontname -> - make a new encoded font |
| 32 |
/newfontname exch def |
| 33 |
/fontname exch def |
| 34 |
|
| 35 |
/fontdict fontname findfont def |
| 36 |
/newfont fontdict maxlength dict def |
| 37 |
|
| 38 |
fontdict { |
| 39 |
exch |
| 40 |
dup /FID eq { |
| 41 |
% skip FID pair |
| 42 |
pop pop |
| 43 |
} { |
| 44 |
% copy to the new font dictionary |
| 45 |
exch newfont 3 1 roll put |
| 46 |
} ifelse |
| 47 |
} forall |
| 48 |
|
| 49 |
newfont /FontName newfontname put |
| 50 |
|
| 51 |
% insert only valid encoding vectors |
| 52 |
encoding_vector length 256 eq { |
| 53 |
newfont /Encoding encoding_vector put |
| 54 |
} if |
| 55 |
|
| 56 |
newfontname newfont definefont pop |
| 57 |
} def |
| 58 |
|
| 59 |
/MF_PS { % fontname newfontname -> - make a new font preserving its enc |
| 60 |
/newfontname exch def |
| 61 |
/fontname exch def |
| 62 |
|
| 63 |
/fontdict fontname findfont def |
| 64 |
/newfont fontdict maxlength dict def |
| 65 |
|
| 66 |
fontdict { |
| 67 |
exch |
| 68 |
dup /FID eq { |
| 69 |
% skip FID pair |
| 70 |
pop pop |
| 71 |
} { |
| 72 |
% copy to the new font dictionary |
| 73 |
exch newfont 3 1 roll put |
| 74 |
} ifelse |
| 75 |
} forall |
| 76 |
|
| 77 |
newfont /FontName newfontname put |
| 78 |
|
| 79 |
newfontname newfont definefont pop |
| 80 |
} def |
| 81 |
|
| 82 |
/SF { % fontname width height -> - set a new font |
| 83 |
/height exch def |
| 84 |
/width exch def |
| 85 |
|
| 86 |
findfont |
| 87 |
[width 0 0 height 0 0] makefont setfont |
| 88 |
} def |
| 89 |
|
| 90 |
/SUF { % fontname width height -> - set a new user font |
| 91 |
/height exch def |
| 92 |
/width exch def |
| 93 |
|
| 94 |
/F-gs-user-font MF |
| 95 |
/F-gs-user-font width height SF |
| 96 |
} def |
| 97 |
|
| 98 |
/SUF_PS { % fontname width height -> - set a new user font preserving its enc |
| 99 |
/height exch def |
| 100 |
/width exch def |
| 101 |
|
| 102 |
/F-gs-user-font MF_PS |
| 103 |
/F-gs-user-font width height SF |
| 104 |
} def |
| 105 |
|
| 106 |
/M {moveto} bind def |
| 107 |
/s {show} bind def |
| 108 |
|
| 109 |
/Box { % x y w h -> - define box path |
| 110 |
/d_h exch def /d_w exch def /d_y exch def /d_x exch def |
| 111 |
d_x d_y moveto |
| 112 |
d_w 0 rlineto |
| 113 |
0 d_h rlineto |
| 114 |
d_w neg 0 rlineto |
| 115 |
closepath |
| 116 |
} def |
| 117 |
|
| 118 |
/bgs { % x y height blskip gray str -> - show string with bg color |
| 119 |
/str exch def |
| 120 |
/gray exch def |
| 121 |
/blskip exch def |
| 122 |
/height exch def |
| 123 |
/y exch def |
| 124 |
/x exch def |
| 125 |
|
| 126 |
gsave |
| 127 |
x y blskip sub str stringwidth pop height Box |
| 128 |
gray setgray |
| 129 |
fill |
| 130 |
grestore |
| 131 |
x y M str s |
| 132 |
} def |
| 133 |
|
| 134 |
/bgcs { % x y height blskip red green blue str -> - show string with bg color |
| 135 |
/str exch def |
| 136 |
/blue exch def |
| 137 |
/green exch def |
| 138 |
/red exch def |
| 139 |
/blskip exch def |
| 140 |
/height exch def |
| 141 |
/y exch def |
| 142 |
/x exch def |
| 143 |
|
| 144 |
gsave |
| 145 |
x y blskip sub str stringwidth pop height Box |
| 146 |
red green blue setrgbcolor |
| 147 |
fill |
| 148 |
grestore |
| 149 |
x y M str s |
| 150 |
} def |
| 151 |
|
| 152 |
% Highlight bars. |
| 153 |
/highlight_bars { % nlines lineheight output_y_margin gray -> - |
| 154 |
gsave |
| 155 |
setgray |
| 156 |
/ymarg exch def |
| 157 |
/lineheight exch def |
| 158 |
/nlines exch def |
| 159 |
|
| 160 |
% This 2 is just a magic number to sync highlight lines to text. |
| 161 |
0 d_header_y ymarg sub 2 sub translate |
| 162 |
|
| 163 |
/cw d_output_w cols div def |
| 164 |
/nrows d_output_h ymarg 2 mul sub lineheight div cvi def |
| 165 |
|
| 166 |
% for each column |
| 167 |
0 1 cols 1 sub { |
| 168 |
cw mul /xp exch def |
| 169 |
|
| 170 |
% for each rows |
| 171 |
0 1 nrows 1 sub { |
| 172 |
/rn exch def |
| 173 |
rn lineheight mul neg /yp exch def |
| 174 |
rn nlines idiv 2 mod 0 eq { |
| 175 |
% Draw highlight bar. 4 is just a magic indentation. |
| 176 |
xp 4 add yp cw 8 sub lineheight neg Box fill |
| 177 |
} if |
| 178 |
} for |
| 179 |
} for |
| 180 |
|
| 181 |
grestore |
| 182 |
} def |
| 183 |
|
| 184 |
% Line highlight bar. |
| 185 |
/line_highlight { % x y width height gray -> - |
| 186 |
gsave |
| 187 |
/gray exch def |
| 188 |
Box gray setgray fill |
| 189 |
grestore |
| 190 |
} def |
| 191 |
|
| 192 |
% Column separator lines. |
| 193 |
/column_lines { |
| 194 |
gsave |
| 195 |
.1 setlinewidth |
| 196 |
0 d_footer_h translate |
| 197 |
/cw d_output_w cols div def |
| 198 |
1 1 cols 1 sub { |
| 199 |
cw mul 0 moveto |
| 200 |
0 d_output_h rlineto stroke |
| 201 |
} for |
| 202 |
grestore |
| 203 |
} def |
| 204 |
|
| 205 |
% Column borders. |
| 206 |
/column_borders { |
| 207 |
gsave |
| 208 |
.1 setlinewidth |
| 209 |
0 d_footer_h moveto |
| 210 |
0 d_output_h rlineto |
| 211 |
d_output_w 0 rlineto |
| 212 |
0 d_output_h neg rlineto |
| 213 |
closepath stroke |
| 214 |
grestore |
| 215 |
} def |
| 216 |
|
| 217 |
% Do the actual underlay drawing |
| 218 |
/draw_underlay { |
| 219 |
ul_style 0 eq { |
| 220 |
ul_str true charpath stroke |
| 221 |
} { |
| 222 |
ul_str show |
| 223 |
} ifelse |
| 224 |
} def |
| 225 |
|
| 226 |
% Underlay |
| 227 |
/underlay { % - -> - |
| 228 |
gsave |
| 229 |
0 d_page_h translate |
| 230 |
d_page_h neg d_page_w atan rotate |
| 231 |
|
| 232 |
ul_gray setgray |
| 233 |
ul_font setfont |
| 234 |
/dw d_page_h dup mul d_page_w dup mul add sqrt def |
| 235 |
ul_str stringwidth pop dw exch sub 2 div ul_h_ptsize -2 div moveto |
| 236 |
draw_underlay |
| 237 |
grestore |
| 238 |
} def |
| 239 |
|
| 240 |
/user_underlay { % - -> - |
| 241 |
gsave |
| 242 |
ul_x ul_y translate |
| 243 |
ul_angle rotate |
| 244 |
ul_gray setgray |
| 245 |
ul_font setfont |
| 246 |
0 0 ul_h_ptsize 2 div sub moveto |
| 247 |
draw_underlay |
| 248 |
grestore |
| 249 |
} def |
| 250 |
|
| 251 |
% Page prefeed |
| 252 |
/page_prefeed { % bool -> - |
| 253 |
statusdict /prefeed known { |
| 254 |
statusdict exch /prefeed exch put |
| 255 |
} { |
| 256 |
pop |
| 257 |
} ifelse |
| 258 |
} def |
| 259 |
|
| 260 |
% Wrapped line markers |
| 261 |
/wrapped_line_mark { % x y charwith charheight type -> - |
| 262 |
/type exch def |
| 263 |
/h exch def |
| 264 |
/w exch def |
| 265 |
/y exch def |
| 266 |
/x exch def |
| 267 |
|
| 268 |
type 2 eq { |
| 269 |
% Black boxes (like TeX does) |
| 270 |
gsave |
| 271 |
0 setlinewidth |
| 272 |
x w 4 div add y M |
| 273 |
0 h rlineto w 2 div 0 rlineto 0 h neg rlineto |
| 274 |
closepath fill |
| 275 |
grestore |
| 276 |
} { |
| 277 |
type 3 eq { |
| 278 |
% Small arrows |
| 279 |
gsave |
| 280 |
.2 setlinewidth |
| 281 |
x w 2 div add y h 2 div add M |
| 282 |
w 4 div 0 rlineto |
| 283 |
x w 4 div add y lineto stroke |
| 284 |
|
| 285 |
x w 4 div add w 8 div add y h 4 div add M |
| 286 |
x w 4 div add y lineto |
| 287 |
w 4 div h 8 div rlineto stroke |
| 288 |
grestore |
| 289 |
} { |
| 290 |
% do nothing |
| 291 |
} ifelse |
| 292 |
} ifelse |
| 293 |
} def |
| 294 |
|
| 295 |
% EPSF import. |
| 296 |
|
| 297 |
/BeginEPSF { |
| 298 |
/b4_Inc_state save def % Save state for cleanup |
| 299 |
/dict_count countdictstack def % Count objects on dict stack |
| 300 |
/op_count count 1 sub def % Count objects on operand stack |
| 301 |
userdict begin |
| 302 |
/showpage { } def |
| 303 |
0 setgray 0 setlinecap |
| 304 |
1 setlinewidth 0 setlinejoin |
| 305 |
10 setmiterlimit [ ] 0 setdash newpath |
| 306 |
/languagelevel where { |
| 307 |
pop languagelevel |
| 308 |
1 ne { |
| 309 |
false setstrokeadjust false setoverprint |
| 310 |
} if |
| 311 |
} if |
| 312 |
} bind def |
| 313 |
|
| 314 |
/EndEPSF { |
| 315 |
count op_count sub { pos } repeat % Clean up stacks |
| 316 |
countdictstack dict_count sub { end } repeat |
| 317 |
b4_Inc_state restore |
| 318 |
} bind def |
| 319 |
|
| 320 |
% Check PostScript language level. |
| 321 |
/languagelevel where { |
| 322 |
pop /gs_languagelevel languagelevel def |
| 323 |
} { |
| 324 |
/gs_languagelevel 1 def |
| 325 |
} ifelse |
| 326 |
%%EndResource |
| 327 |
%%BeginResource: procset Enscript-Encoding-88591 1.6 4 |
| 328 |
/encoding_vector [ |
| 329 |
/.notdef /.notdef /.notdef /.notdef |
| 330 |
/.notdef /.notdef /.notdef /.notdef |
| 331 |
/.notdef /.notdef /.notdef /.notdef |
| 332 |
/.notdef /.notdef /.notdef /.notdef |
| 333 |
/.notdef /.notdef /.notdef /.notdef |
| 334 |
/.notdef /.notdef /.notdef /.notdef |
| 335 |
/.notdef /.notdef /.notdef /.notdef |
| 336 |
/.notdef /.notdef /.notdef /.notdef |
| 337 |
/space /exclam /quotedbl /numbersign |
| 338 |
/dollar /percent /ampersand /quoteright |
| 339 |
/parenleft /parenright /asterisk /plus |
| 340 |
/comma /hyphen /period /slash |
| 341 |
/zero /one /two /three |
| 342 |
/four /five /six /seven |
| 343 |
/eight /nine /colon /semicolon |
| 344 |
/less /equal /greater /question |
| 345 |
/at /A /B /C |
| 346 |
/D /E /F /G |
| 347 |
/H /I /J /K |
| 348 |
/L /M /N /O |
| 349 |
/P /Q /R /S |
| 350 |
/T /U /V /W |
| 351 |
/X /Y /Z /bracketleft |
| 352 |
/backslash /bracketright /asciicircum /underscore |
| 353 |
/quoteleft /a /b /c |
| 354 |
/d /e /f /g |
| 355 |
/h /i /j /k |
| 356 |
/l /m /n /o |
| 357 |
/p /q /r /s |
| 358 |
/t /u /v /w |
| 359 |
/x /y /z /braceleft |
| 360 |
/bar /braceright /tilde /.notdef |
| 361 |
/.notdef /.notdef /.notdef /.notdef |
| 362 |
/.notdef /.notdef /.notdef /.notdef |
| 363 |
/.notdef /.notdef /.notdef /.notdef |
| 364 |
/.notdef /.notdef /.notdef /.notdef |
| 365 |
/.notdef /.notdef /.notdef /.notdef |
| 366 |
/.notdef /.notdef /.notdef /.notdef |
| 367 |
/.notdef /.notdef /.notdef /.notdef |
| 368 |
/.notdef /.notdef /.notdef /.notdef |
| 369 |
/space /exclamdown /cent /sterling |
| 370 |
/currency /yen /brokenbar /section |
| 371 |
/dieresis /copyright /ordfeminine /guillemotleft |
| 372 |
/logicalnot /hyphen /registered /macron |
| 373 |
/degree /plusminus /twosuperior /threesuperior |
| 374 |
/acute /mu /paragraph /bullet |
| 375 |
/cedilla /onesuperior /ordmasculine /guillemotright |
| 376 |
/onequarter /onehalf /threequarters /questiondown |
| 377 |
/Agrave /Aacute /Acircumflex /Atilde |
| 378 |
/Adieresis /Aring /AE /Ccedilla |
| 379 |
/Egrave /Eacute /Ecircumflex /Edieresis |
| 380 |
/Igrave /Iacute /Icircumflex /Idieresis |
| 381 |
/Eth /Ntilde /Ograve /Oacute |
| 382 |
/Ocircumflex /Otilde /Odieresis /multiply |
| 383 |
/Oslash /Ugrave /Uacute /Ucircumflex |
| 384 |
/Udieresis /Yacute /Thorn /germandbls |
| 385 |
/agrave /aacute /acircumflex /atilde |
| 386 |
/adieresis /aring /ae /ccedilla |
| 387 |
/egrave /eacute /ecircumflex /edieresis |
| 388 |
/igrave /iacute /icircumflex /idieresis |
| 389 |
/eth /ntilde /ograve /oacute |
| 390 |
/ocircumflex /otilde /odieresis /divide |
| 391 |
/oslash /ugrave /uacute /ucircumflex |
| 392 |
/udieresis /yacute /thorn /ydieresis |
| 393 |
] def |
| 394 |
%%EndResource |
| 395 |
%%EndProlog |
| 396 |
%%BeginSetup |
| 397 |
%%IncludeResource: font Courier-Bold |
| 398 |
%%IncludeResource: font Courier |
| 399 |
/HFpt_w 10 def |
| 400 |
/HFpt_h 10 def |
| 401 |
/Courier-Bold /HF-gs-font MF |
| 402 |
/HF /HF-gs-font findfont [HFpt_w 0 0 HFpt_h 0 0] makefont def |
| 403 |
/Courier /F-gs-font MF |
| 404 |
/F-gs-font 10 10 SF |
| 405 |
/#copies 1 def |
| 406 |
% Pagedevice definitions: |
| 407 |
gs_languagelevel 1 gt { |
| 408 |
<< |
| 409 |
/PageSize [612 792] |
| 410 |
>> setpagedevice |
| 411 |
} if |
| 412 |
/d_page_w 744 def |
| 413 |
/d_page_h 536 def |
| 414 |
/d_header_x 0 def |
| 415 |
/d_header_y 536 def |
| 416 |
/d_header_w 744 def |
| 417 |
/d_header_h 0 def |
| 418 |
/d_footer_x 0 def |
| 419 |
/d_footer_y 0 def |
| 420 |
/d_footer_w 744 def |
| 421 |
/d_footer_h 0 def |
| 422 |
/d_output_w 744 def |
| 423 |
/d_output_h 536 def |
| 424 |
/cols 1 def |
| 425 |
%%EndSetup |
| 426 |
%%Page: (1) 1 |
| 427 |
%%BeginPageSetup |
| 428 |
_S |
| 429 |
90 rotate |
| 430 |
24 -574 translate |
| 431 |
/pagenum 1 def |
| 432 |
/fname (read510.m) def |
| 433 |
/fdir () def |
| 434 |
/ftail (read510.m) def |
| 435 |
/user_header_p false def |
| 436 |
/user_footer_p false def |
| 437 |
%%EndPageSetup |
| 438 |
5 523 M |
| 439 |
(% ) s |
| 440 |
5 512 M |
| 441 |
(% Read in CS 510 input fields) s |
| 442 |
5 501 M |
| 443 |
(xc1 = zeros\(511,511,6\);) s |
| 444 |
5 490 M |
| 445 |
(yc1 = zeros\(511,511,6\);) s |
| 446 |
5 479 M |
| 447 |
(dxf1 = zeros\(511,511,6\);) s |
| 448 |
5 468 M |
| 449 |
(dyf1 = zeros\(511,511,6\);) s |
| 450 |
5 457 M |
| 451 |
(ra1 = zeros\(511,511,6\);) s |
| 452 |
5 446 M |
| 453 |
(xg1 = zeros\(511,511,6\);) s |
| 454 |
5 435 M |
| 455 |
(yg1 = zeros\(511,511,6\);) s |
| 456 |
5 424 M |
| 457 |
(dxv1 = zeros\(511,511,6\);) s |
| 458 |
5 413 M |
| 459 |
(dyu1 = zeros\(511,511,6\);) s |
| 460 |
5 402 M |
| 461 |
(raz1 = zeros\(511,511,6\);) s |
| 462 |
5 391 M |
| 463 |
(dxc1 = zeros\(511,511,6\);) s |
| 464 |
5 380 M |
| 465 |
(dyc1 = zeros\(511,511,6\);) s |
| 466 |
5 369 M |
| 467 |
(raw1 = zeros\(511,511,6\);) s |
| 468 |
5 358 M |
| 469 |
(ras1 = zeros\(511,511,6\);) s |
| 470 |
5 347 M |
| 471 |
(dxg1 = zeros\(511,511,6\);) s |
| 472 |
5 336 M |
| 473 |
(dyg1 = zeros\(511,511,6\);) s |
| 474 |
5 325 M |
| 475 |
(for iface=1:6) s |
| 476 |
5 314 M |
| 477 |
(fid = fopen\(['tile00' num2str\(iface\) '.mitgrid'],'r','b'\);) s |
| 478 |
5 303 M |
| 479 |
(xc1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 480 |
5 292 M |
| 481 |
(yc1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 482 |
5 281 M |
| 483 |
(dxf1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 484 |
5 270 M |
| 485 |
(dyf1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 486 |
5 259 M |
| 487 |
(ra1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 488 |
5 248 M |
| 489 |
(xg1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 490 |
5 237 M |
| 491 |
(yg1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 492 |
5 226 M |
| 493 |
(dxv1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 494 |
5 215 M |
| 495 |
(dyu1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 496 |
5 204 M |
| 497 |
(raz1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 498 |
5 193 M |
| 499 |
(dxc1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 500 |
5 182 M |
| 501 |
(dyc1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 502 |
5 171 M |
| 503 |
(raw1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 504 |
5 160 M |
| 505 |
(ras1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 506 |
5 149 M |
| 507 |
(dxg1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 508 |
5 138 M |
| 509 |
(dyg1\(:,:,iface\) = fread\(fid,[511 511],'real*8'\);) s |
| 510 |
5 127 M |
| 511 |
(end) s |
| 512 |
5 116 M |
| 513 |
(%) s |
| 514 |
5 105 M |
| 515 |
(% Now output \(approx 1 deg\) grid) s |
| 516 |
5 94 M |
| 517 |
(xc2 = zeros\(103,103,6\);) s |
| 518 |
5 83 M |
| 519 |
(yc2 = zeros\(103,103,6\);) s |
| 520 |
5 72 M |
| 521 |
(dxf2 = zeros\(103,103,6\);) s |
| 522 |
5 61 M |
| 523 |
(dyf2 = zeros\(103,103,6\);) s |
| 524 |
5 50 M |
| 525 |
(ra2 = zeros\(103,103,6\);) s |
| 526 |
5 39 M |
| 527 |
(xg2 = zeros\(103,103,6\);) s |
| 528 |
5 28 M |
| 529 |
(yg2 = zeros\(103,103,6\);) s |
| 530 |
5 17 M |
| 531 |
(dxv2 = zeros\(103,103,6\);) s |
| 532 |
5 6 M |
| 533 |
(dyu2 = zeros\(103,103,6\);) s |
| 534 |
_R |
| 535 |
S |
| 536 |
%%Page: (2) 2 |
| 537 |
%%BeginPageSetup |
| 538 |
_S |
| 539 |
90 rotate |
| 540 |
24 -574 translate |
| 541 |
/pagenum 2 def |
| 542 |
/fname (read510.m) def |
| 543 |
/fdir () def |
| 544 |
/ftail (read510.m) def |
| 545 |
/user_header_p false def |
| 546 |
/user_footer_p false def |
| 547 |
%%EndPageSetup |
| 548 |
5 523 M |
| 549 |
(raz2 = zeros\(103,103,6\);) s |
| 550 |
5 512 M |
| 551 |
(dxc2 = zeros\(103,103,6\);) s |
| 552 |
5 501 M |
| 553 |
(dyc2 = zeros\(103,103,6\);) s |
| 554 |
5 490 M |
| 555 |
(raw2 = zeros\(103,103,6\);) s |
| 556 |
5 479 M |
| 557 |
(ras2 = zeros\(103,103,6\);) s |
| 558 |
5 468 M |
| 559 |
(dxg2 = zeros\(103,103,6\);) s |
| 560 |
5 457 M |
| 561 |
(dyg2 = zeros\(103,103,6\);) s |
| 562 |
5 446 M |
| 563 |
(% Create left, x-center, bottom and y-center indeces into 510 grid) s |
| 564 |
5 435 M |
| 565 |
(ratio = 5;) s |
| 566 |
5 424 M |
| 567 |
(ileft=[1:5:510];) s |
| 568 |
5 413 M |
| 569 |
(icent=ileft+2;) s |
| 570 |
5 402 M |
| 571 |
(jbott=[1:5:510];) s |
| 572 |
5 391 M |
| 573 |
(jcent=jbott+2;) s |
| 574 |
5 380 M |
| 575 |
(% First do the interior of each face) s |
| 576 |
5 369 M |
| 577 |
(for iface=1:6) s |
| 578 |
5 358 M |
| 579 |
(xc2\(1:102,1:102,iface\)=xc1\(icent,jcent,iface\);) s |
| 580 |
5 347 M |
| 581 |
(yc2\(1:102,1:102,iface\)=yc1\(icent,jcent,iface\);) s |
| 582 |
5 336 M |
| 583 |
(xg2\(1:102,1:102,iface\)=xg1\(ileft,jcent,iface\);) s |
| 584 |
5 325 M |
| 585 |
(yg2\(1:102,1:102,iface\)=yg1\(icent,jbott,iface\);) s |
| 586 |
5 314 M |
| 587 |
(for ipnt=0:ratio-1) s |
| 588 |
5 303 M |
| 589 |
(dxf2\(1:102,1:102,iface\)=dxf2\(1:102,1:102,iface\) + dxf1\(ileft+ipnt,jcent,iface\);) s |
| 590 |
5 292 M |
| 591 |
(dxg2\(1:102,1:102,iface\)=dxg2\(1:102,1:102,iface\) + dxg1\(ileft+ipnt,jbott,iface\);) s |
| 592 |
5 281 M |
| 593 |
(dxv2\(2:102,2:102,iface\)=dxv2\(2:102,2:102,iface\) + dxv1\(icent\(2:102\)+ipnt+1,jbott\(2:102\),iface\);) s |
| 594 |
5 270 M |
| 595 |
(dxc2\(2:102,2:102,iface\)=dxc2\(2:102,2:102,iface\) + dxc1\(icent\(2:102\)+ipnt+1,jcent\(2:102\),iface\);) s |
| 596 |
5 259 M |
| 597 |
(end) s |
| 598 |
5 248 M |
| 599 |
(for ipnt=0:ratio-1) s |
| 600 |
5 237 M |
| 601 |
(dyf2\(1:102,1:102,iface\)=dyf2\(1:102,1:102,iface\) + dyf1\(icent,jbott+ipnt,iface\);) s |
| 602 |
5 226 M |
| 603 |
(dyg2\(1:102,1:102,iface\)=dyg2\(1:102,1:102,iface\) + dyg1\(ileft,jbott+ipnt,iface\);) s |
| 604 |
5 215 M |
| 605 |
(dyu2\(2:102,2:102,iface\)=dyu2\(2:102,2:102,iface\) + dyu1\(ileft\(2:102\),jcent\(2:102\)+jpnt+1,iface\);) s |
| 606 |
5 204 M |
| 607 |
(dyc2\(2:102,2:102,iface\)=dyc2\(2:102,2:102,iface\) + dyc1\(icent\(2:102\),jcent\(2:102\)+jpnt+1,iface\);) s |
| 608 |
5 193 M |
| 609 |
(end) s |
| 610 |
5 182 M |
| 611 |
(for ipnt=0:ratio-1) s |
| 612 |
5 171 M |
| 613 |
(for jpnt=0:ratio-1) s |
| 614 |
5 160 M |
| 615 |
(ra2\(1:102,1:102,iface\)=ra2\(1:102,1:102,iface\) + ra1\(ileft+ipnt,jbott+jpnt,iface\);) s |
| 616 |
5 149 M |
| 617 |
(ras2\(2:102,2:102,iface\)=ras2\(2:102,2:102,iface\) + ras1\(ileft\(2:102\)+ipnt+1,jbott\(2:102\)+jpnt,iface\);) s |
| 618 |
5 138 M |
| 619 |
(raw2\(2:102,2:102,iface\)=raw2\(2:102,2:102,iface\) + raw1\(ileft\(2:102\)+ipnt,jbott\(2:102\)+jpnt+1,iface\);) s |
| 620 |
5 127 M |
| 621 |
(raz2\(2:102,2:102,iface\)=raz2\(2:102,2:102,iface\) + raz1\(ileft\(2:102\)+ipnt,jbott\(2:102\)+jpnt,iface\);) s |
| 622 |
5 116 M |
| 623 |
(end) s |
| 624 |
5 105 M |
| 625 |
(end) s |
| 626 |
5 94 M |
| 627 |
(% Now the edges, corners and exchanges) s |
| 628 |
5 83 M |
| 629 |
(% Exchange for lon and lat's at center and edges - no directions) s |
| 630 |
5 72 M |
| 631 |
(xc2\(103,1:102,[1 3 5]\)=xc2\(1,1:102,[2 4 6]\);) s |
| 632 |
5 61 M |
| 633 |
(xc2\(1:102,103,[1 3 5]\)=xc2\(1,1:102,[3 5 1]\);) s |
| 634 |
5 50 M |
| 635 |
(xc2\(103,103,[1 2 3 4 5 6]\)=xc2\(1,1,[3 4 5 6 1 2]\);) s |
| 636 |
5 39 M |
| 637 |
(yc2\(103,1:102,[1 3 5]\)=yc2\(1,1:102,[2 4 6]\);) s |
| 638 |
5 28 M |
| 639 |
(yc2\(103,1:102,[1 3 5]\)=yc2\(1,1:102,[2 4 6]\);) s |
| 640 |
5 17 M |
| 641 |
(yc2\(103,103,[1 2 3 4 5 6]\)=yc2\(1,1,[3 4 5 6 1 2]\);) s |
| 642 |
5 6 M |
| 643 |
(xg2\(103,1:102,[1 3 5]\)=xg2\(1,1:102,[2 4 6]\);) s |
| 644 |
_R |
| 645 |
S |
| 646 |
%%Page: (3) 3 |
| 647 |
%%BeginPageSetup |
| 648 |
_S |
| 649 |
90 rotate |
| 650 |
24 -574 translate |
| 651 |
/pagenum 3 def |
| 652 |
/fname (read510.m) def |
| 653 |
/fdir () def |
| 654 |
/ftail (read510.m) def |
| 655 |
/user_header_p false def |
| 656 |
/user_footer_p false def |
| 657 |
%%EndPageSetup |
| 658 |
5 523 M |
| 659 |
(xg2\(103,1:102,[1 3 5]\)=xg2\(1,1:102,[2 4 6]\);) s |
| 660 |
5 512 M |
| 661 |
(xg2\(103,103,[1 2 3 4 5 6]\)=xg2\(1,1,[3 4 5 6 1 2]\);) s |
| 662 |
5 501 M |
| 663 |
(yg2\(103,1:102,[1 3 5]\)=yg2\(1,1:102,[2 4 6]\);) s |
| 664 |
5 490 M |
| 665 |
(yg2\(103,1:102,[1 3 5]\)=yg2\(1,1:102,[2 4 6]\);) s |
| 666 |
5 479 M |
| 667 |
(yg2\(103,103,[1 2 3 4 5 6]\)=yg2\(1,1,[3 4 5 6 1 2]\);) s |
| 668 |
5 468 M |
| 669 |
(% Exchange for dx's and dy's at center - directions but no signs) s |
| 670 |
5 457 M |
| 671 |
(dxf2\(103,1:102,[1 3 5]\)=dxf2\(1,1:102,[2 4 6]\);) s |
| 672 |
5 446 M |
| 673 |
(dxf2\(103,1:102,[1 3 5]\)=dxf2\(1,1:102,[2 4 6]\);) s |
| 674 |
5 435 M |
| 675 |
(dxf2\(103,103,[1 2 3 4 5 6]\)=dxf2\(1,1,[3 4 5 6 1 2]\);) s |
| 676 |
_R |
| 677 |
S |
| 678 |
%%Trailer |
| 679 |
%%Pages: 3 |
| 680 |
%%DocumentNeededResources: font Courier-Bold Courier |
| 681 |
%%EOF |
Parent Directory
| 
Revision Log
| 
Revision Graph