[cig-commits] r12086 - in mc/2D: . ConMan ConMan/trunk ConMan/trunk/doc ConMan/trunk/src ConMan/trunk/src/blankenbach ConMan/trunk/src/blankenbach/bm1a ConMan/trunk/src/blankenbach/bm1a/output ConMan/trunk/src/blankenbach/bm1b ConMan/trunk/src/blankenbach/bm1b/output ConMan/trunk/src/blankenbach/bm1c ConMan/trunk/src/blankenbach/bm1c/output ConMan/trunk/src/grid.src ConMan/trunk/src/mm.src ConMan/trunk/src/solver.src ConMan/trunk/src.p
wei at geodynamics.org
wei at geodynamics.org
Wed Jun 4 15:36:19 PDT 2008
Author: wei
Date: 2008-06-04 15:36:13 -0700 (Wed, 04 Jun 2008)
New Revision: 12086
Added:
mc/2D/ConMan/
mc/2D/ConMan/branches/
mc/2D/ConMan/tags/
mc/2D/ConMan/trunk/
mc/2D/ConMan/trunk/doc/
mc/2D/ConMan/trunk/doc/users_guide.pdf
mc/2D/ConMan/trunk/src.p/
mc/2D/ConMan/trunk/src.p/src.p.tar.gz
mc/2D/ConMan/trunk/src/
mc/2D/ConMan/trunk/src/ConMan.F
mc/2D/ConMan/trunk/src/Makefile
mc/2D/ConMan/trunk/src/adjbc.f
mc/2D/ConMan/trunk/src/batchelor.f
mc/2D/ConMan/trunk/src/bc.h
mc/2D/ConMan/trunk/src/blankenbach/
mc/2D/ConMan/trunk/src/blankenbach/Makefile
mc/2D/ConMan/trunk/src/blankenbach/bench_table.bash
mc/2D/ConMan/trunk/src/blankenbach/blankenbach.aux
mc/2D/ConMan/trunk/src/blankenbach/blankenbach.dvi
mc/2D/ConMan/trunk/src/blankenbach/blankenbach.log
mc/2D/ConMan/trunk/src/blankenbach/blankenbach.ps
mc/2D/ConMan/trunk/src/blankenbach/blankenbach.tex
mc/2D/ConMan/trunk/src/blankenbach/bm1a/
mc/2D/ConMan/trunk/src/blankenbach/bm1a/geom.128
mc/2D/ConMan/trunk/src/blankenbach/bm1a/geom.32
mc/2D/ConMan/trunk/src/blankenbach/bm1a/geom.64
mc/2D/ConMan/trunk/src/blankenbach/bm1a/geom.96
mc/2D/ConMan/trunk/src/blankenbach/bm1a/in.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/in.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/in.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/in.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/geoid.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/geoid.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/geoid.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/geoid.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/mean.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/mean.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/mean.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/mean.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/out.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/out.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/out.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/out.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/temp.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/temp.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/temp.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/temp.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/ts.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/ts.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/ts.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/output/ts.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/r.64.4i
mc/2D/ConMan/trunk/src/blankenbach/bm1a/r.64.4o
mc/2D/ConMan/trunk/src/blankenbach/bm1a/run.128.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/run.32.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/run.64.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/run.96.4
mc/2D/ConMan/trunk/src/blankenbach/bm1a/stdout
mc/2D/ConMan/trunk/src/blankenbach/bm1b/
mc/2D/ConMan/trunk/src/blankenbach/bm1b/geoid.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/geom.128
mc/2D/ConMan/trunk/src/blankenbach/bm1b/geom.32
mc/2D/ConMan/trunk/src/blankenbach/bm1b/geom.64
mc/2D/ConMan/trunk/src/blankenbach/bm1b/geom.96
mc/2D/ConMan/trunk/src/blankenbach/bm1b/in.128.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/in.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/in.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/in.96.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/mean.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/out.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/geoid.128.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/geoid.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/geoid.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/geoid.96.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/mean.128.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/mean.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/mean.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/mean.96.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/out.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/out.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/temp.128.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/temp.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/temp.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/temp.96.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/ts.128.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/ts.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/ts.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/output/ts.96.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/r.64.5i
mc/2D/ConMan/trunk/src/blankenbach/bm1b/r.64.5o
mc/2D/ConMan/trunk/src/blankenbach/bm1b/run.128.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/run.32.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/run.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/run.96.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/stdout
mc/2D/ConMan/trunk/src/blankenbach/bm1b/temp.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1b/ts.64.5
mc/2D/ConMan/trunk/src/blankenbach/bm1c/
mc/2D/ConMan/trunk/src/blankenbach/bm1c/geoid.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/geom.128
mc/2D/ConMan/trunk/src/blankenbach/bm1c/geom.32
mc/2D/ConMan/trunk/src/blankenbach/bm1c/geom.64
mc/2D/ConMan/trunk/src/blankenbach/bm1c/geom.96
mc/2D/ConMan/trunk/src/blankenbach/bm1c/in.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/in.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/in.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/in.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/mean.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/out.128.63
mc/2D/ConMan/trunk/src/blankenbach/bm1c/out.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/geoid.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/geoid.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/geoid.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/geoid.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/mean.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/mean.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/mean.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/mean.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/out.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/out.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/out.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/out.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/temp.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/temp.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/temp.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/temp.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/ts.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/ts.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/ts.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/output/ts.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/r.64.6i
mc/2D/ConMan/trunk/src/blankenbach/bm1c/r.64.6o
mc/2D/ConMan/trunk/src/blankenbach/bm1c/r.96.6o2
mc/2D/ConMan/trunk/src/blankenbach/bm1c/run.128.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/run.32.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/run.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/run.96.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/stdout
mc/2D/ConMan/trunk/src/blankenbach/bm1c/temp.64.6
mc/2D/ConMan/trunk/src/blankenbach/bm1c/ts.64.6
mc/2D/ConMan/trunk/src/blankenbach/do_plot
mc/2D/ConMan/trunk/src/blankenbach/dotimeseries
mc/2D/ConMan/trunk/src/blankenbach/geom.128
mc/2D/ConMan/trunk/src/blankenbach/geom.32
mc/2D/ConMan/trunk/src/blankenbach/geom.64
mc/2D/ConMan/trunk/src/blankenbach/geom.96
mc/2D/ConMan/trunk/src/blankenbach/mypost
mc/2D/ConMan/trunk/src/blankenbach/tv.ps.cmd
mc/2D/ConMan/trunk/src/common.h
mc/2D/ConMan/trunk/src/differ.f
mc/2D/ConMan/trunk/src/eg2.F
mc/2D/ConMan/trunk/src/eglib.f
mc/2D/ConMan/trunk/src/elminp.F
mc/2D/ConMan/trunk/src/error.f
mc/2D/ConMan/trunk/src/f_tlhs.F
mc/2D/ConMan/trunk/src/f_tmres.F
mc/2D/ConMan/trunk/src/f_trhsimp.F
mc/2D/ConMan/trunk/src/f_vrres.F
mc/2D/ConMan/trunk/src/f_vstf.F
mc/2D/ConMan/trunk/src/fe.h
mc/2D/ConMan/trunk/src/fluxke.F
mc/2D/ConMan/trunk/src/g2assm.f
mc/2D/ConMan/trunk/src/genien.F
mc/2D/ConMan/trunk/src/genmat.F
mc/2D/ConMan/trunk/src/genshg.F
mc/2D/ConMan/trunk/src/genshp.F
mc/2D/ConMan/trunk/src/geom.12
mc/2D/ConMan/trunk/src/geom.bm1
mc/2D/ConMan/trunk/src/geom.slava
mc/2D/ConMan/trunk/src/grid.src/
mc/2D/ConMan/trunk/src/grid.src/Makefile
mc/2D/ConMan/trunk/src/grid.src/clear.f
mc/2D/ConMan/trunk/src/grid.src/common.h
mc/2D/ConMan/trunk/src/grid.src/genBC.f
mc/2D/ConMan/trunk/src/grid.src/genNBC.f
mc/2D/ConMan/trunk/src/grid.src/genfl.f
mc/2D/ConMan/trunk/src/grid.src/genfl12.f
mc/2D/ConMan/trunk/src/grid.src/geniBC.f
mc/2D/ConMan/trunk/src/grid.src/iclear_mov.f
mc/2D/ConMan/trunk/src/grid.src/move.f
mc/2D/ConMan/trunk/src/grid.src/smove.f
mc/2D/ConMan/trunk/src/in.12
mc/2D/ConMan/trunk/src/in.all
mc/2D/ConMan/trunk/src/in.bm1
mc/2D/ConMan/trunk/src/in.bm5
mc/2D/ConMan/trunk/src/in.bm6
mc/2D/ConMan/trunk/src/in.compress1
mc/2D/ConMan/trunk/src/in.slava5
mc/2D/ConMan/trunk/src/in.test
mc/2D/ConMan/trunk/src/input.F
mc/2D/ConMan/trunk/src/intemp.F
mc/2D/ConMan/trunk/src/libs/
mc/2D/ConMan/trunk/src/materials.h
mc/2D/ConMan/trunk/src/mean.bm1
mc/2D/ConMan/trunk/src/mm.src/
mc/2D/ConMan/trunk/src/mm.src/Makefile
mc/2D/ConMan/trunk/src/mm.src/Makefile.fc
mc/2D/ConMan/trunk/src/mm.src/Makefile.fce
mc/2D/ConMan/trunk/src/mm.src/driver
mc/2D/ConMan/trunk/src/mm.src/driver.f
mc/2D/ConMan/trunk/src/mm.src/heapexit.f
mc/2D/ConMan/trunk/src/mm.src/icharln.f
mc/2D/ConMan/trunk/src/mm.src/mm2000.f
mc/2D/ConMan/trunk/src/mm.src/mm2000.h
mc/2D/ConMan/trunk/src/mm.src/mm2000_32.h
mc/2D/ConMan/trunk/src/mm.src/mm2000_64.h
mc/2D/ConMan/trunk/src/mm.src/mmcomprs.f
mc/2D/ConMan/trunk/src/mm.src/mmgetblk.F
mc/2D/ConMan/trunk/src/mm.src/mmgetindex.f
mc/2D/ConMan/trunk/src/mm.src/mmprint.f
mc/2D/ConMan/trunk/src/mm.src/mmrelblk.f
mc/2D/ConMan/trunk/src/mm.src/mmsc.c
mc/2D/ConMan/trunk/src/mm.src/mmsc_orig.c
mc/2D/ConMan/trunk/src/mm.src/mmverify.f
mc/2D/ConMan/trunk/src/mydate.f
mc/2D/ConMan/trunk/src/mytime.f
mc/2D/ConMan/trunk/src/out.bm1
mc/2D/ConMan/trunk/src/pickdt.F
mc/2D/ConMan/trunk/src/print.f
mc/2D/ConMan/trunk/src/print_reg.F
mc/2D/ConMan/trunk/src/printd.f
mc/2D/ConMan/trunk/src/prj4q.f
mc/2D/ConMan/trunk/src/process.F
mc/2D/ConMan/trunk/src/projct.F
mc/2D/ConMan/trunk/src/prtstr_reg.F
mc/2D/ConMan/trunk/src/r.12in
mc/2D/ConMan/trunk/src/r.12out
mc/2D/ConMan/trunk/src/r.bm1in
mc/2D/ConMan/trunk/src/r.bm1out
mc/2D/ConMan/trunk/src/rdflux.F
mc/2D/ConMan/trunk/src/rheol.F
mc/2D/ConMan/trunk/src/rheol.newt.F
mc/2D/ConMan/trunk/src/rheol.nn.sdk.F
mc/2D/ConMan/trunk/src/run.12
mc/2D/ConMan/trunk/src/run.test
mc/2D/ConMan/trunk/src/runbm1
mc/2D/ConMan/trunk/src/solver.h
mc/2D/ConMan/trunk/src/solver.src/
mc/2D/ConMan/trunk/src/solver.src/Makefile
mc/2D/ConMan/trunk/src/solver.src/back.f
mc/2D/ConMan/trunk/src/solver.src/coldot.f
mc/2D/ConMan/trunk/src/solver.src/factor.f
mc/2D/ConMan/trunk/src/solver.src/unback.f
mc/2D/ConMan/trunk/src/solver.src/unfact.f
mc/2D/ConMan/trunk/src/solver.src/untest.f
mc/2D/ConMan/trunk/src/solver.src/zfactor.f
mc/2D/ConMan/trunk/src/stress.F
mc/2D/ConMan/trunk/src/stress.h
mc/2D/ConMan/trunk/src/temp.12
mc/2D/ConMan/trunk/src/temp.bm1
mc/2D/ConMan/trunk/src/timdrv.F
mc/2D/ConMan/trunk/src/timdrv.nn.F
mc/2D/ConMan/trunk/src/timer.f
mc/2D/ConMan/trunk/src/tmass.F
mc/2D/ConMan/trunk/src/vadd.F
mc/2D/ConMan/trunk/src/variables.h
Log:
Initial import of ConMan code package
Added: mc/2D/ConMan/trunk/doc/users_guide.pdf
===================================================================
(Binary files differ)
Property changes on: mc/2D/ConMan/trunk/doc/users_guide.pdf
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: mc/2D/ConMan/trunk/src/ConMan.F
===================================================================
--- mc/2D/ConMan/trunk/src/ConMan.F (rev 0)
+++ mc/2D/ConMan/trunk/src/ConMan.F 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,134 @@
+ program conman
+
+c----------------------------------------------------------------------
+c *********************************************************************
+c *********************************************************************
+c *********************************************************************
+c *********************************************************************
+c ** **
+c ** ***** ==> C O N M A N <== ***** **
+c ** **
+c ** Convective Mantle Flow code **
+c ** **
+c ** A Petrov Galerkin Finite Element Program to Analyze **
+c ** Incompressible Thermal Convective Flows. **
+c ** **
+c ** **
+c ** written by: Arthur Raefsky and Scott King . **
+c ** **
+c ** **
+c *********************************************************************
+c *********************************************************************
+c----------------------------------------------------------------------
+c
+ implicit double precision (a-h,o-z)
+c
+c.... remove above card for single-precision operation
+c
+ character*80 name1,name2,name3,name4,name6,name7,name8
+ common /temp1 / x(50000)
+ include 'common.h'
+c
+c.... set up the input/output files
+c
+ read(5,1001) name1
+ read(5,1001) name2
+ read(5,1001) name3
+ read(5,1001) name4
+ read(5,1001) name5
+ read(5,1001) name6
+ read(5,1001) name7
+ read(5,1001) name8
+ 1001 format( a80 )
+ open (unit=iin, file=name1 , status='old')
+ open (unit=igeom, file=name2 , status='old')
+ open (unit=iout, file=name3 , status='unknown')
+ open (unit=irsin, file=name4 , status='unknown')
+c open (unit=irsout, file=name5 , status='unknown')
+ open (unit=itsout, file=name6 , status='unknown' )
+ open (unit=itout, file=name7 , status='unknown')
+ open (unit=imout, file=name8 , status='unknown' )
+c
+c.... Input phase
+c
+ call timer ('begin ')
+ call input
+ close (iin)
+ close (igeom)
+c
+c.... Execution phase
+c
+ call proces
+ call timer ('end ')
+
+ write(iout,2000)
+ write(iout,2001) ntotal
+ rndpbyte=ntotal*4.e-6
+ write(iout,2002) rndpbyte
+
+ call mmprint(iout)
+
+ 2000 format(
+ & ' M e m o r y A l l o c a t i o n S t a t i s t i c s',//)
+ 2001 format('ntotal = ',i12 ,' Total Words of Memory Program ')
+ 2002 format('ndpbyte = ',f12.2,' Mb Memory Allocated')
+
+c
+c.... Close files
+c
+ close (iout)
+ close (irsin)
+ close (itsout)
+ close (itout)
+ close (imout)
+c
+c.... end of the program
+c
+ call exit (istat)
+ end
+ block data
+c----------------------------------------------------------------------
+c
+ implicit double precision (a-h,o-z)
+c
+c.... remove above card for single-precision operation
+c
+c.... dynamic storage pointer management data
+c
+c common /point / mbeg, mend, mprec, mdct, ndict
+c
+c.... io channels
+ common /io / iin , igeom, iout , itsout , itout ,
+ & imout , irsin , irsout
+c
+c.... useful constants
+ common /const / zero , pt25, pt33, pt5, pt66, one,
+ & onept5, two , three, four, sixten, eps7
+c
+c----------------------------------------------------------------------
+c----------------------------------------------------------------------
+c
+c
+c data mbeg, mend, mprec, mdct, ndict
+c & / 1,130000000, 2, 1, 0 /
+c
+ data iin , igeom , iout , itsout , itout , imout ,
+ & irsin , irsout
+ & / 10, 11 , 12, 13, 14, 15,
+ & 16, 17 /
+c
+ data zero, pt25,
+ & pt33, pt5,
+ & pt66, one,
+ & onept5, two,
+ & three, four,
+ & sixten, eps7
+ & / 0.00000000000000d+0, 2.50000000000000d-1,
+ & 3.33333333333333d-1, 5.00000000000000d-1,
+ & 6.66666666666666d-1, 1.00000000000000d+0,
+ & 1.50000000000000d+0, 2.00000000000000d+0,
+ & 3.00000000000000d+0, 4.00000000000000d+0,
+ & 1.60000000000000d+1, 1.00000000000000d-7 /
+c
+c----------------------------------------------------------------------
+ end
Added: mc/2D/ConMan/trunk/src/Makefile
===================================================================
--- mc/2D/ConMan/trunk/src/Makefile (rev 0)
+++ mc/2D/ConMan/trunk/src/Makefile 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,62 @@
+#
+# makefile for the program conman on the Sun
+#
+COMMON=ConMan.o adjbc.o eg2.o eglib.o elminp.o error.o \
+ f_vstf.o f_vrres.o fluxke.o g2assm.o \
+ genien.o genmat.o genshg.o genshp.o \
+ input.o intemp.o pickdt.o print.o printd.o process.o \
+ rdflux.o stress.o timer.o tmass.o \
+ vadd.o mydate.o mytime.o projct.o prj4q.o print_reg.o prtstr_reg.o \
+ batchelor.o
+
+SRC_HOME=/home/scott/ConMan02
+
+IMPLICIT=f_tlhs.o f_trhsimp.o
+PICARD=f_tlhs.o f_trhsimp.o
+EXPLICIT=f_tmres.o
+
+LIBS= libs/libsolver.a libs/libmm.a libs/libgrid.a
+
+PROGRAM=conman
+#
+# switch for flags around for picard and implicit versions
+# use -m64 -r8 -i8 for 64-bit architecture (also see changed in mm.src)
+#
+FC=ifort
+FFLAGS = -O2 -DIMPLICIT -I.
+FFLAGS = -O2 -DIMPLICIT -DPICARD -I. -m64 -r8 -i8
+FFLAGS = -O2 -I. -m64 -r8 -i8
+
+#
+explicit: $(COMMON) $(EXPLICIT) rheol.newt.o timdrv.o
+ (cd solver.src ; make)
+ (cd mm.src ; make)
+ (cd grid.src ; make)
+ $(FC) -I$(SRC_HOME)/libs -o conman.exp rheol.newt.o timdrv.o $(COMMON) $(EXPLICIT) $(LIBS)
+
+picard: $(COMMON) $(PICARD) rheol.newt.o timdrv.o
+ (cd solver.src ; make)
+ (cd mm.src ; make)
+ (cd grid.src ; make)
+ $(FC) $(FFLAGS) -I$(SRC_HOME)/libs -o conman.pic rheol.newt.o timdrv.o $(COMMON) $(PICARD) $(LIBS)
+
+picardnn: $(COMMON) $(PICARD) rheol.nn.sdk.o timdrv.nn.o
+ (cd solver.src ; make)
+ (cd mm.src ; make)
+ (cd grid.src ; make)
+ $(FC) -I$(SRC_HOME)/libs -o conman.nn.pic rheol.nn.sdk.o timdrv.nn.o $(COMMON) $(PICARD) $(LIBS)
+
+implicit: $(COMMON) $(IMPLICIT)
+ (cd solver.src ; make)
+ (cd mm.src ; make)
+ (cd grid.src ; make)
+ $(FC) -o conman.imp $(COMMON) $(IMPLICIT) $(LIBS)
+
+clean:
+ rm -rf core $(COMMON) $(PICARD) $(EXPLICIT) $(IMPLICIT) rheol.newt.o timdrv.o libs/lib*
+ (cd solver.src ; make clean)
+ (cd mm.src ; make clean)
+ (cd grid.src ; make clean)
+
+veryclean:
+ rm -rf core $(COMMON) $(PICARD) $(EXPLICIT) $(IMPLICIT) conman.???
Added: mc/2D/ConMan/trunk/src/adjbc.f
===================================================================
--- mc/2D/ConMan/trunk/src/adjbc.f (rev 0)
+++ mc/2D/ConMan/trunk/src/adjbc.f 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,41 @@
+ subroutine adjbc(lmv, stiff, ivel, vbcl, vbcr)
+c
+c
+ implicit double precision (a-h,o-z)
+c
+c.... remove above card for single-precision operation
+c
+ include 'common.h'
+c
+ dimension lmv(numel,*),stiff(*),vbcl(*),vbcr(*)
+ integer idiag(8)
+c
+ idiag(1) = 1
+ idiag(2) = 3
+ idiag(3) = 6
+ idiag(4) = 10
+ idiag(5) = 15
+ idiag(6) = 21
+ idiag(7) = 28
+ idiag(8) = 36
+c
+ do j = 1, 8
+ do i = 1, 8
+ if ( lmv(ivel,i) .ne. 0 ) then
+c
+ if ( j .ge. i) then
+ igadd = idiag( j ) - j + i
+ else
+ igadd = idiag( i ) - i + j
+ endif
+ vbcr( lmv(ivel,i) ) = vbcr( lmv(ivel,i) )
+ & - stiff(igadd) * vbcl(j)
+c
+ end if
+ enddo
+ enddo
+c
+c.... return
+c
+ return
+ end
Added: mc/2D/ConMan/trunk/src/batchelor.f
===================================================================
--- mc/2D/ConMan/trunk/src/batchelor.f (rev 0)
+++ mc/2D/ConMan/trunk/src/batchelor.f 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,67 @@
+ subroutine batchelor(xorig, yorig, u, v)
+c
+c A program to implement the Batchelor cornerflow solution
+c See Batchelor (1967) pp. 224-225
+c
+c Written: S. D. King November 1, 2002
+c
+c constants
+c
+ implicit double precision (a-h,o-z)
+ common /io / iin, igeom, iout , itsout , itout , imout,
+ & irsin , irsout
+
+ two = 2.0e0
+ four = 4.0e0
+ pi = four*atan(1.0)
+c
+c velocity of the downgoing slab
+c
+c Uplate = 1.538330
+ Uplate = 1.54
+ Uplate = 1.53
+ visc = 1.0e3
+c
+c the general expression for the streamfunction is given by:
+c psi = Ax + By + (Cx + Dy)*arctan(y/x)
+c
+c and
+c u = - d psi/dy = -B -x*(Cx + Dy)/(x^2+y^2) - D arctan(y/x)
+c
+c v = d psi/dx = A - y*(Cx + Dy)/(x^2+y^2) + C arctan(y/x)
+c
+c for the arc corner with a slab at theta = pi/4 where the downgoing
+c slab velocity is U*sqrt(two)/two, we have
+c
+ A = 0.0e0
+ B = pi*Uplate/(two-(pi**2)/four)
+ C = -pi*Uplate/(two-(pi**2)/four)
+ D = -Uplate*two*(two-pi/two)/(two-(pi**2)/four)
+c
+c works for my grids
+ x = 1.0e3*(610.0e0-xorig)
+ y = 1.0e3*(550.0e0-yorig)
+c works for changyeol's grids
+c x = 1.0e3*(xorig-50.0e0)
+c y = 1.0e3*(550.0e0-yorig)
+ press = -2*visc*(C*x + D*y)/((x**2)+(y**2))
+ psi = A*x + B*y + (C*x + D*y)*atan2(y,x)
+ u = -B - x*(C*x + D*y)/((x**2)+(y**2))
+ & - D*atan2(y,x)
+ v = A - y*(C*x + D*y)/((x**2)+(y**2))
+ & + C*atan2(y,x)
+c if (u .le. Uplate) u = Uplate
+c if (v .le. Uplate) v = Uplate
+c needed for sdk grids
+ u = -u
+ v = -v
+ if (y .le. 0) then
+ u = 0.0d0
+ v = 0.0d0
+ endif
+c write(iout,*) xorig,yorig,x,y,u,v
+c
+c end
+c
+ return
+ end
Added: mc/2D/ConMan/trunk/src/bc.h
===================================================================
--- mc/2D/ConMan/trunk/src/bc.h (rev 0)
+++ mc/2D/ConMan/trunk/src/bc.h 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,19 @@
+ common /bccom/ ipidv , ipidt , ipvbc , iptbc , ipisuf , ipsuf
+
+ pointer (ipidv,idv)
+ dimension idv(ndof,numnp)
+
+ pointer (ipvbc,vbc)
+ dimension vbc(ndof,numnp)
+
+ pointer (ipidt,idt)
+ dimension idt(numnp)
+
+ pointer (iptbc,tbc)
+ dimension tbc(numnp)
+
+ pointer (ipisuf,isuf)
+ dimension isuf(2,numsuf)
+
+ pointer (ipsuf,suf)
+ dimension suf(3,numsuf)
Added: mc/2D/ConMan/trunk/src/blankenbach/Makefile
===================================================================
--- mc/2D/ConMan/trunk/src/blankenbach/Makefile (rev 0)
+++ mc/2D/ConMan/trunk/src/blankenbach/Makefile 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,177 @@
+#
+# Makefile to test ConMan using 2D Blankenbach Benchmarks
+# Written 3/24/99 Scott King
+#
+
+CONMANDIR = ~/Desktop/Temp_Dep_Rheology_Project/src.production
+LOCDIR = ${CONMANDIR}/blankenbach
+CONMAN = ${CONMANDIR}/conman.exp
+
+
+blankenbach_1: blankenbach_1a blankenbach_1b blankenbach_1c
+blankenbach_1a: blank_1a_32 blank_1a_64 blank_1a_96
+blankenbach_1b: blank_1b_32 blank_1b_64 blank_1b_96
+blankenbach_1c: blank_1c_32 blank_1c_64 blank_1c_96
+blankenbach_64: blank_1a_64 blank_1b_64 blank_1c_64
+blankenbach_96: blank_1a_96 blank_1b_96 blank_1c_96
+blankenbach_128: blank_1b_128 blank_1b_128 blank_1c_128
+
+
+blank_1a_32:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1A 32 by 32 GRID"
+ - at cd ${LOCDIR}/bm1a; \
+ ${CONMAN} < run.32.4 > stdout; \
+ if (diff output/ts.32.4 ts.32.4) then true; \
+ else echo "Possible problem with TIMESERIES bm1a 32 by 32, diffs above"; fi ; \
+ if (diff output/geoid.32.4 geoid.32.4) then true; \
+ else echo "Possible problem with GEOID bm1a 32 by 32, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1A 32 by 32 GRID"
+ - at echo "========================================="
+
+blank_1a_64:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1A 64 by 64 GRID"
+ - at cd ${LOCDIR}/bm1a; \
+ ${CONMAN} < run.64.4 > stdout; \
+ if (diff output/ts.64.4 ts.64.4) then true; \
+ else echo "Possible problem with TIMESERIES bm1a 64 by 64, diffs above"; fi ; \
+ if (diff output/geoid.64.4 geoid.64.4) then true; \
+ else echo "Possible problem with GEOID bm1a 64 by 64, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1A 64 by 64 GRID"
+ - at echo "========================================="
+
+blank_1a_96:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1A 96 by 96 GRID"
+ - at cd ${LOCDIR}/bm1a; \
+ ${CONMAN} < run.96.4 > stdout; \
+ if (diff output/ts.96.4 ts.96.4) then true; \
+ else echo "Possible problem with TIMESERIES bm1a 96 by 96, diffs above"; fi ; \
+ if (diff output/geoid.96.4 geoid.96.4) then true; \
+ else echo "Possible problem with GEOID bm1a 96 by 96, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1A 96 by 96 GRID"
+ - at echo "========================================="
+
+blank_1a_128:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1A 128 by 128 GRID"
+ - at cd ${LOCDIR}/bm1a; \
+ ${CONMAN} < run.128.4 > stdout; \
+ if (diff output/ts.128.4 ts.128.4) then true; \
+ else echo "Possible problem with TIMESERIES bm1a 128 by 128, diffs above"; fi ; \
+ if (diff output/geoid.128.4 geoid.128.4) then true; \
+ else echo "Possible problem with GEOID bm1a 128 by 128, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1A 128 by 128 GRID"
+ - at echo "========================================="
+
+blank_1b_32:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1B 32 by 32 GRID"
+ - at cd ${LOCDIR}/bm1b; \
+ ${CONMAN} < run.32.5 > stdout; \
+ if (diff output/ts.32.5 ts.32.5) then true; \
+ else echo "Possible problem with TIMESERIES bm1b 32 by 32, diffs above"; fi ; \
+ if (diff output/geoid.32.5 geoid.32.5) then true; \
+ else echo "Possible problem with GEOID bm1b 32 by 32, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1B 32 by 32 GRID"
+ - at echo "========================================="
+
+blank_1b_64:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1B 64 by 64 GRID"
+ - at cd ${LOCDIR}/bm1b; \
+ ${CONMAN} < run.64.5 > stdout; \
+ if (diff output/ts.64.5 ts.64.5) then true; \
+ else echo "Possible problem with TIMESERIES bm1b 64 by 64, diffs above"; fi ; \
+ if (diff output/geoid.64.5 geoid.64.5) then true; \
+ else echo "Possible problem with GEOID bm1b 64 by 64, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1B 64 by 64 GRID"
+ - at echo "========================================="
+
+blank_1b_96:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1B 96 by 96 GRID"
+ - at cd ${LOCDIR}/bm1b; \
+ ${CONMAN} < run.96.5 > stdout; \
+ if (diff output/ts.96.5 ts.96.5) then true; \
+ else echo "Possible problem with TIMESERIES bm1b 96 by 96, diffs above"; fi ; \
+ if (diff output/geoid.96.5 geoid.96.5) then true; \
+ else echo "Possible problem with GEOID bm1b 96 by 96, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1B 96 by 96 GRID"
+ - at echo "========================================="
+
+blank_1b_128:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1B 128 by 128 GRID"
+ - at cd ${LOCDIR}/bm1b; \
+ ${CONMAN} < run.128.5 > stdout; \
+ if (diff output/ts.128.5 ts.128.5) then true; \
+ else echo "Possible problem with TIMESERIES bm1b 128 by 128, diffs above"; fi ; \
+ if (diff output/geoid.128.5 geoid.128.5) then true; \
+ else echo "Possible problem with GEOID bm1b 128 by 128, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1B 128 by 128 GRID"
+ - at echo "========================================="
+
+blank_1c_32:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1C 32 by 32 GRID"
+ - at cd ${LOCDIR}/bm1c; \
+ ${CONMAN} < run.32.6 > stdout; \
+ if (diff output/ts.32.6 ts.32.6) then true; \
+ else echo "Possible problem with TIMESERIES bm1c 32 by 32, diffs above"; fi ; \
+ if (diff output/geoid.32.6 geoid.32.6) then true; \
+ else echo "Possible problem with GEOID bm1c 32 by 32, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1C 32 by 32 GRID"
+ - at echo "========================================="
+
+blank_1c_64:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1C 64 by 64 GRID"
+ - at cd ${LOCDIR}/bm1c; \
+ ${CONMAN} < run.64.6 > stdout; \
+ if (diff output/ts.64.6 ts.64.6) then true; \
+ else echo "Possible problem with TIMESERIES bm1c 64 by 64, diffs above"; fi ; \
+ if (diff output/geoid.64.6 geoid.64.6) then true; \
+ else echo "Possible problem with GEOID bm1c 64 by 64, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1C 64 by 64 GRID"
+ - at echo "========================================="
+
+blank_1c_96:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1C 96 by 96 GRID"
+ - at cd ${LOCDIR}/bm1c; \
+ ${CONMAN} < run.96.6 > stdout; \
+ if (diff output/ts.96.6 ts.96.6) then true; \
+ else echo "Possible problem with TIMESERIES bm1c 96 by 96, diffs above"; fi ; \
+ if (diff output/geoid.96.6 geoid.96.6) then true; \
+ else echo "Possible problem with GEOID bm1c 96 by 96, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1C 96 by 96 GRID"
+ - at echo "========================================="
+
+blank_1c_128:
+ - at echo "BEGINNING BLANKENBACH BENCHMARK 1C 128 by 128 GRID"
+ - at cd ${LOCDIR}/bm1c; \
+ ${CONMAN} < run.128.6 > stdout; \
+ if (diff output/ts.128.6 ts.128.6) then true; \
+ else echo "Possible problem with TIMESERIES bm1c 128 by 128, diffs above"; fi ; \
+ if (diff output/geoid.128.6 geoid.128.6) then true; \
+ else echo "Possible problem with GEOID bm1c 128 by 128, diffs above"; fi
+ - at echo "COMPLETED BLANKENBACH BENCHMARK 1C 128 by 128 GRID"
+ - at echo "========================================="
+
+
+clean:
+ - at cd ${LOCDIR}/bm1a; \
+ ${RM} -f mean.32.4 mean.64.4 mean.96.4 mean.128.4 ; \
+ ${RM} -f temp.32.4 temp.64.4 temp.96.4 temp.128.4 ; \
+ ${RM} -f ts.32.4 ts.64.4 ts.96.4 ts.128.4 ; \
+ ${RM} -f out.32.4 out.64.4 out.96.4 out.128.4 ; \
+ ${RM} -f geoid.32.4 geoid.64.4 geoid.96.4 geoid.128.4 ; \
+ ${RM} -f r.32.4i r.64i.4i r.96.4i r.128.4i ; \
+ ${RM} -f r.32.4o r.64o.4o r.96.4o r.128.4o
+ - at cd ${LOCDIR}/bm1b; \
+ ${RM} -f mean.32.5 mean.65.5 mean.96.5 mean.128.5 ; \
+ ${RM} -f temp.32.5 temp.65.5 temp.96.5 temp.128.5 ; \
+ ${RM} -f ts.32.5 ts.65.5 ts.96.5 ts.128.5 ; \
+ ${RM} -f out.32.5 out.65.5 out.96.5 out.128.5 ; \
+ ${RM} -f geoid.32.5 geoid.65.5 geoid.96.5 geoid.128.5 ; \
+ ${RM} -f r.32.5i r.65i.5i r.96.5i r.128.5i ; \
+ ${RM} -f r.32.5o r.65o.5o r.96.5o r.128.5o
+ - at cd ${LOCDIR}/bm1c; \
+ ${RM} -f mean.32.6 mean.66.6 mean.96.6 mean.128.6 ; \
+ ${RM} -f temp.32.6 temp.66.6 temp.96.6 temp.128.6 ; \
+ ${RM} -f ts.32.6 ts.66.6 ts.96.6 ts.128.6 ; \
+ ${RM} -f out.32.6 out.66.6 out.96.6 out.128.6 ; \
+ ${RM} -f geoid.32.6 geoid.66.6 geoid.96.6 geoid.128.6 ; \
+ ${RM} -f r.32.6i r.66i.6i r.96.6i r.128.6i ; \
+ ${RM} -f r.32.6o r.66o.6o r.96.6o r.128.6o
Added: mc/2D/ConMan/trunk/src/blankenbach/bench_table.bash
===================================================================
--- mc/2D/ConMan/trunk/src/blankenbach/bench_table.bash (rev 0)
+++ mc/2D/ConMan/trunk/src/blankenbach/bench_table.bash 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,211 @@
+# Make a LaTeX table for Benchmarks
+
+#########################################
+# User Adjustable Parameters #
+#########################################
+FILE=blankenbach
+LATEX='latex'
+DVIPS='/sw/bin/dvips'
+GHOST='/sw/bin/gv'
+export FILE LATEX DVIPS GHOST
+
+#########################################
+# End User Adjustable Parameters #
+#########################################
+
+TEXFILE=$FILE.tex
+DVIFILE=$FILE.dvi
+PSFILE=$FILE.ps
+export TEXFILE DVIFILE PSFILE
+
+/bin/rm -f $TEXFILE
+echo '\documentclass[11pt]{article}' >> $TEXFILE
+echo '\begin{document}' >> $TEXFILE
+#########################################
+# Rayleigh Number 10,000 Blankenbach #
+#########################################
+echo '\begin{table}' >> $TEXFILE
+echo '\centering' >> $TEXFILE
+echo '\begin{tabular} {| r c c c c c c r |} \hline ' >> $TEXFILE
+echo '\multicolumn{1}{| c }{Grid} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{V$_{rms}$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Nusselt No.} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Topo$_L$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Topo$_R$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Geoid$_L$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Geoid$_R$} & ' >> $TEXFILE
+echo '\multicolumn{1}{ c |}{Run Time (sec)} \\ ' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+ vrms=`tail -1 bm1a/ts.32.4 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1a/ts.32.4 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==34 {printf "%7.3f", $2}' bm1a/geoid.32.4`
+ geoidl=`awk 'NR==34 {printf "%5.3f", $4}' bm1a/geoid.32.4`
+ topor=`awk 'NR==66 {printf "%7.3f", $2}' bm1a/geoid.32.4`
+ geoidr=`awk 'NR==66 {printf "%5.3f", $4}' bm1a/geoid.32.4`
+ time=`awk '$1=="total" {print $2}' bm1a/out.32.4`
+ export vrms nu topol geoidl topor geoidr time
+echo '32 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1a/ts.64.4 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1a/ts.64.4 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==66 {printf "%7.3f", $2}' bm1a/geoid.64.4`
+ geoidl=`awk 'NR==66 {printf "%5.3f", $4}' bm1a/geoid.64.4`
+ topor=`awk 'NR==130 {printf "%7.3f", $2}' bm1a/geoid.64.4`
+ geoidr=`awk 'NR==130 {printf "%5.3f", $4}' bm1a/geoid.64.4`
+ time=`awk '$1=="total" {print $2}' bm1a/out.64.4`
+ export vrms nu topol geoidl topor geoidr time
+echo '64 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1a/ts.96.4 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1a/ts.96.4 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==98 {printf "%7.3f", $2}' bm1a/geoid.96.4`
+ geoidl=`awk 'NR==98 {printf "%5.3f", $4}' bm1a/geoid.96.4`
+ topor=`awk 'NR==194 {printf "%7.3f", $2}' bm1a/geoid.96.4`
+ geoidr=`awk 'NR==194 {printf "%5.3f", $4}' bm1a/geoid.96.4`
+ time=`awk '$1=="total" {print $2}' bm1a/out.96.4`
+ export vrms nu topol geoidl topor geoidr time
+echo '96 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1a/ts.128.4 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1a/ts.128.4 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==130 {printf "%7.3f", $2}' bm1a/geoid.128.4`
+ geoidl=`awk 'NR==130 {printf "%5.3f", $4}' bm1a/geoid.128.4`
+ topor=`awk 'NR==258 {printf "%7.3f", $2}' bm1a/geoid.128.4`
+ geoidr=`awk 'NR==258 {printf "%5.3f", $4}' bm1a/geoid.128.4`
+ time=`awk '$1=="total" {print $2}' bm1a/out.128.4`
+ export vrms nu topol geoidl topor geoidr time
+echo '128 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\dag C$_{ext}$ & 42.865 & 4.884 & 2254.021 & -2903.221 & 54.822 & -62.622 & \null \\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\multicolumn{8}{|l|}{\dag Christensen$^{\prime}$s extrapolated values.} \\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\end{tabular}' >> $TEXFILE
+echo '\caption{Blankenbach (1989) Benchmark 1a: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^4$ using ConMan} ' >> $TEXFILE
+echo '\end{table}' >> $TEXFILE
+
+#########################################
+# Rayleigh Number 100,000 Blankenbach #
+#########################################
+echo '\begin{table}' >> $TEXFILE
+echo '\centering' >> $TEXFILE
+echo '\begin{tabular} {| r c c c c c c r |} \hline ' >> $TEXFILE
+echo '\multicolumn{1}{| c }{Grid} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{V$_{rms}$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Nusselt No.} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Topo$_L$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Topo$_R$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Geoid$_L$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Geoid$_R$} & ' >> $TEXFILE
+echo '\multicolumn{1}{ c |}{Run Time (sec)} \\ ' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+ vrms=`tail -1 bm1b/ts.32.5 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1b/ts.32.5 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==1 {printf "%7.3f", $2}' bm1b/geoid.32.5`
+ geoidl=`awk 'NR==1 {printf "%5.3f", $4}' bm1b/geoid.32.5`
+ topor=`awk 'NR==33 {printf "%7.3f", $2}' bm1b/geoid.32.5`
+ geoidr=`awk 'NR==33 {printf "%5.3f", $4}' bm1b/geoid.32.5`
+ time=`awk '$1=="total" {print $2}' bm1b/out.32.5`
+ export vrms nu topol geoidl topor geoidr time
+echo '32 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1b/ts.64.5 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1b/ts.64.5 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==1 {printf "%7.3f", $2}' bm1b/geoid.64.5`
+ geoidl=`awk 'NR==1 {printf "%5.3f", $4}' bm1b/geoid.64.5`
+ topor=`awk 'NR==65 {printf "%7.3f", $2}' bm1b/geoid.64.5`
+ geoidr=`awk 'NR==65 {printf "%5.3f", $4}' bm1b/geoid.64.5`
+ time=`awk '$1=="total" {print $2}' bm1b/out.64.5`
+ export vrms nu topol geoidl topor geoidr time
+echo '64 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1b/ts.96.5 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1b/ts.96.5 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==1 {printf "%7.3f", $2}' bm1b/geoid.96.5`
+ geoidl=`awk 'NR==1 {printf "%5.3f", $4}' bm1b/geoid.96.5`
+ topor=`awk 'NR==97 {printf "%7.3f", $2}' bm1b/geoid.96.5`
+ geoidr=`awk 'NR==97 {printf "%5.3f", $4}' bm1b/geoid.96.5`
+ time=`awk '$1=="total" {print $2}' bm1b/out.96.5`
+ export vrms nu topol geoidl topor geoidr time
+echo '96 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1b/ts.128.5 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1b/ts.128.5 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==1 {printf "%7.3f", $2}' bm1b/geoid.128.5`
+ geoidl=`awk 'NR==1 {printf "%5.3f", $4}' bm1b/geoid.128.5`
+ topor=`awk 'NR==129 {printf "%7.3f", $2}' bm1b/geoid.128.5`
+ geoidr=`awk 'NR==129 {printf "%5.3f", $4}' bm1b/geoid.128.5`
+ time=`awk '$1=="total" {print $2}' bm1b/out.128.5`
+ export vrms nu topol geoidl topor geoidr time
+echo '128 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\dag C$_{ext}$ & 193.214 & 10.534 & 1460.986 & -2004.205 & 27.703 & -32.016 & \null \\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\multicolumn{8}{|l|}{\dag Christensen$^{\prime}$s extrapolated values.} \\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\end{tabular}' >> $TEXFILE
+echo '\caption{Blankenbach (1989) Benchmark 1b: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^5$ using ConMan} ' >> $TEXFILE
+echo '\end{table}' >> $TEXFILE
+
+#########################################
+# Rayleigh Number 1,000,000 Blankenbach #
+#########################################
+echo '\begin{table}' >> $TEXFILE
+echo '\centering' >> $TEXFILE
+echo '\begin{tabular} {| r c c c c c c r |} \hline ' >> $TEXFILE
+echo '\multicolumn{1}{| c }{Grid} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{V$_{rms}$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Nusselt No.} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Topo$_L$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Topo$_R$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Geoid$_L$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c }{Geoid$_R$} &' >> $TEXFILE
+echo '\multicolumn{1}{ c |}{Run Time (sec)} \\ ' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+ vrms=`tail -1 bm1c/ts.32.6 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1c/ts.32.6 | awk '{printf "%5.3f", $4}'`
+ topor=`awk 'NR==1 {printf "%7.3f", $2}' bm1c/geoid.32.6`
+ geoidr=`awk 'NR==1 {printf "%5.3f", $4}' bm1c/geoid.32.6`
+ topol=`awk 'NR==33 {printf "%7.3f", $2}' bm1c/geoid.32.6`
+ geoidl=`awk 'NR==33 {printf "%5.3f", $4}' bm1c/geoid.32.6`
+ time=`awk '$1=="total" {print $2}' bm1c/out.32.6`
+ export vrms nu topol geoidl topor geoidr time
+echo '32 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1c/ts.64.6 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1c/ts.64.6 | awk '{printf "%5.3f", $4}'`
+ topor=`awk 'NR==1 {printf "%7.3f", $2}' bm1c/geoid.64.6`
+ geoidr=`awk 'NR==1 {printf "%5.3f", $4}' bm1c/geoid.64.6`
+ topol=`awk 'NR==65 {printf "%7.3f", $2}' bm1c/geoid.64.6`
+ geoidl=`awk 'NR==65 {printf "%5.3f", $4}' bm1c/geoid.64.6`
+ time=`awk '$1=="total" {print $2}' bm1c/out.64.6`
+ export vrms nu topol geoidl topor geoidr time
+echo '64 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1c/ts.96.6 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1c/ts.96.6 | awk '{printf "%5.3f", $4}'`
+ topor=`awk 'NR==1 {printf "%7.3f", $2}' bm1c/geoid.96.6`
+ geoidr=`awk 'NR==1 {printf "%5.3f", $4}' bm1c/geoid.96.6`
+ topol=`awk 'NR==97 {printf "%7.3f", $2}' bm1c/geoid.96.6`
+ geoidl=`awk 'NR==97 {printf "%5.3f", $4}' bm1c/geoid.96.6`
+ time=`awk '$1=="total" {print $2}' bm1c/out.96.6`
+ export vrms nu topol geoidl topor geoidr time
+echo '96 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+ vrms=`tail -1 bm1c/ts.128.6 | awk '{printf "%5.3f", sqrt(2*$5)}'`
+ nu=`tail -1 bm1c/ts.128.6 | awk '{printf "%5.3f", $4}'`
+ topol=`awk 'NR==1 {printf "%7.3f", $2}' bm1c/geoid.128.6`
+ geoidl=`awk 'NR==1 {printf "%5.3f", $4}' bm1c/geoid.128.6`
+ topor=`awk 'NR==129 {printf "%7.3f", $2}' bm1c/geoid.128.6`
+ geoidr=`awk 'NR==129 {printf "%5.3f", $4}' bm1c/geoid.128.6`
+ time=`awk '$1=="total" {print $2}' bm1c/out.128.6`
+ export vrms nu topol geoidl topor geoidr time
+echo '128 &' $vrms '&' $nu '&' $topol '&' $topor '&' $geoidl '&' $geoidr '&' $time '\\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\dag C$_{ext}$ & 833.989 & 21.997 & 931.962 & -1283.813 & 13.4520& -15.034& \null \\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\multicolumn{8}{|l|}{\dag Christensen$^{\prime}$s extrapolated values.} \\' >> $TEXFILE
+echo '\hline ' >> $TEXFILE
+echo '\end{tabular}' >> $TEXFILE
+echo '\caption{Blankenbach (1989) Benchmark 1c: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^6$ using ConMan.} ' >> $TEXFILE
+echo '\end{table}' >> $TEXFILE
+
+echo '\end{document}' >> $TEXFILE
+
+#########################################
+# Run LaTeX and make Postscript file #
+#########################################
+$LATEX $TEXFILE
+$DVIPS $DVIFILE -o $PSFILE
+$GHOST $PSFILE
Property changes on: mc/2D/ConMan/trunk/src/blankenbach/bench_table.bash
___________________________________________________________________
Name: svn:executable
+
Added: mc/2D/ConMan/trunk/src/blankenbach/blankenbach.aux
===================================================================
--- mc/2D/ConMan/trunk/src/blankenbach/blankenbach.aux (rev 0)
+++ mc/2D/ConMan/trunk/src/blankenbach/blankenbach.aux 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,4 @@
+\relax
+\@writefile{lot}{\contentsline {table}{\numberline {1}{\ignorespaces Blankenbach (1989) Benchmark 1a: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^4$ using ConMan}}{1}}
+\@writefile{lot}{\contentsline {table}{\numberline {2}{\ignorespaces Blankenbach (1989) Benchmark 1b: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^5$ using ConMan}}{1}}
+\@writefile{lot}{\contentsline {table}{\numberline {3}{\ignorespaces Blankenbach (1989) Benchmark 1c: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^6$ using ConMan.}}{1}}
Added: mc/2D/ConMan/trunk/src/blankenbach/blankenbach.dvi
===================================================================
(Binary files differ)
Property changes on: mc/2D/ConMan/trunk/src/blankenbach/blankenbach.dvi
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: mc/2D/ConMan/trunk/src/blankenbach/blankenbach.log
===================================================================
--- mc/2D/ConMan/trunk/src/blankenbach/blankenbach.log (rev 0)
+++ mc/2D/ConMan/trunk/src/blankenbach/blankenbach.log 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,84 @@
+This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6) (format=latex 2007.5.18) 18 APR 2008 12:14
+entering extended mode
+ %&-line parsing enabled.
+**blankenbach.tex
+(./blankenbach.tex
+LaTeX2e <2005/12/01>
+Babel <v3.8h> and hyphenation patterns for english, usenglishmax, dumylang, noh
+yphenation, arabic, basque, bulgarian, coptic, welsh, czech, slovak, german, ng
+erman, danish, esperanto, spanish, catalan, galician, estonian, farsi, finnish,
+ french, greek, monogreek, ancientgreek, croatian, hungarian, interlingua, ibyc
+us, indonesian, icelandic, italian, latin, mongolian, dutch, norsk, polish, por
+tuguese, pinyin, romanian, russian, slovenian, uppersorbian, serbian, swedish,
+turkish, ukenglish, ukrainian, loaded.
+(/usr/local/texlive/2007/texmf-dist/tex/latex/base/article.cls
+Document Class: article 2005/09/16 v1.4f Standard LaTeX document class
+(/usr/local/texlive/2007/texmf-dist/tex/latex/base/size11.clo
+File: size11.clo 2005/09/16 v1.4f Standard LaTeX file (size option)
+)
+\c at part=\count79
+\c at section=\count80
+\c at subsection=\count81
+\c at subsubsection=\count82
+\c at paragraph=\count83
+\c at subparagraph=\count84
+\c at figure=\count85
+\c at table=\count86
+\abovecaptionskip=\skip41
+\belowcaptionskip=\skip42
+\bibindent=\dimen102
+)
+(./blankenbach.aux)
+\openout1 = `blankenbach.aux'.
+
+LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 2.
+LaTeX Font Info: ... okay on input line 2.
+LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 2.
+LaTeX Font Info: ... okay on input line 2.
+LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 2.
+LaTeX Font Info: ... okay on input line 2.
+LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 2.
+LaTeX Font Info: ... okay on input line 2.
+LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 2.
+LaTeX Font Info: ... okay on input line 2.
+LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 2.
+LaTeX Font Info: ... okay on input line 2.
+LaTeX Font Info: External font `cmex10' loaded for size
+(Font) <10.95> on input line 5.
+LaTeX Font Info: External font `cmex10' loaded for size
+(Font) <8> on input line 5.
+LaTeX Font Info: External font `cmex10' loaded for size
+(Font) <6> on input line 5.
+LaTeX Font Info: Try loading font information for OMS+cmr on input line 20.
+
+(/usr/local/texlive/2007/texmf-dist/tex/latex/base/omscmr.fd
+File: omscmr.fd 1999/05/25 v2.5h Standard LaTeX font definitions
+)
+LaTeX Font Info: Font shape `OMS/cmr/m/n' in size <10.95> not available
+(Font) Font shape `OMS/cmsy/m/n' tried instead on input line 20.
+
+Overfull \hbox (77.74152pt too wide) in paragraph at lines 5--25
+ [][]
+ []
+
+
+Overfull \hbox (83.21652pt too wide) in paragraph at lines 29--49
+ [][]
+ []
+
+
+Overfull \hbox (79.29317pt too wide) in paragraph at lines 53--73
+ [][]
+ []
+
+[1] (./blankenbach.aux) )
+Here is how much of TeX's memory you used:
+ 236 strings out of 94074
+ 2709 string characters out of 1163007
+ 54804 words of memory out of 1500000
+ 3599 multiletter control sequences out of 10000+50000
+ 6379 words of font info for 23 fonts, out of 1200000 for 2000
+ 645 hyphenation exceptions out of 8191
+ 21i,9n,19p,253b,251s stack positions out of 5000i,500n,6000p,200000b,5000s
+
+Output written on blankenbach.dvi (1 page, 4364 bytes).
Added: mc/2D/ConMan/trunk/src/blankenbach/blankenbach.ps
===================================================================
--- mc/2D/ConMan/trunk/src/blankenbach/blankenbach.ps (rev 0)
+++ mc/2D/ConMan/trunk/src/blankenbach/blankenbach.ps 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,972 @@
+%!PS-Adobe-2.0
+%%Creator: dvips(k) 5.92b Copyright 2002 Radical Eye Software
+%%Title: blankenbach.dvi
+%%Pages: 1
+%%PageOrder: Ascend
+%%BoundingBox: 0 0 596 842
+%%DocumentFonts: CMR10 CMMI8 CMSY10 CMSY8 CMR8
+%%EndComments
+%DVIPSWebPage: (www.radicaleye.com)
+%DVIPSCommandLine: /sw/bin/dvips blankenbach.dvi -o blankenbach.ps
+%DVIPSParameters: dpi=600, compressed
+%DVIPSSource: TeX output 2004.08.18:2052
+%%BeginProcSet: texc.pro
+%!
+/TeXDict 300 dict def TeXDict begin/N{def}def/B{bind def}N/S{exch}N/X{S
+N}B/A{dup}B/TR{translate}N/isls false N/vsize 11 72 mul N/hsize 8.5 72
+mul N/landplus90{false}def/@rigin{isls{[0 landplus90{1 -1}{-1 1}ifelse 0
+0 0]concat}if 72 Resolution div 72 VResolution div neg scale isls{
+landplus90{VResolution 72 div vsize mul 0 exch}{Resolution -72 div hsize
+mul 0}ifelse TR}if Resolution VResolution vsize -72 div 1 add mul TR[
+matrix currentmatrix{A A round sub abs 0.00001 lt{round}if}forall round
+exch round exch]setmatrix}N/@landscape{/isls true N}B/@manualfeed{
+statusdict/manualfeed true put}B/@copies{/#copies X}B/FMat[1 0 0 -1 0 0]
+N/FBB[0 0 0 0]N/nn 0 N/IEn 0 N/ctr 0 N/df-tail{/nn 8 dict N nn begin
+/FontType 3 N/FontMatrix fntrx N/FontBBox FBB N string/base X array
+/BitMaps X/BuildChar{CharBuilder}N/Encoding IEn N end A{/foo setfont}2
+array copy cvx N load 0 nn put/ctr 0 N[}B/sf 0 N/df{/sf 1 N/fntrx FMat N
+df-tail}B/dfs{div/sf X/fntrx[sf 0 0 sf neg 0 0]N df-tail}B/E{pop nn A
+definefont setfont}B/Cw{Cd A length 5 sub get}B/Ch{Cd A length 4 sub get
+}B/Cx{128 Cd A length 3 sub get sub}B/Cy{Cd A length 2 sub get 127 sub}
+B/Cdx{Cd A length 1 sub get}B/Ci{Cd A type/stringtype ne{ctr get/ctr ctr
+1 add N}if}B/id 0 N/rw 0 N/rc 0 N/gp 0 N/cp 0 N/G 0 N/CharBuilder{save 3
+1 roll S A/base get 2 index get S/BitMaps get S get/Cd X pop/ctr 0 N Cdx
+0 Cx Cy Ch sub Cx Cw add Cy setcachedevice Cw Ch true[1 0 0 -1 -.1 Cx
+sub Cy .1 sub]/id Ci N/rw Cw 7 add 8 idiv string N/rc 0 N/gp 0 N/cp 0 N{
+rc 0 ne{rc 1 sub/rc X rw}{G}ifelse}imagemask restore}B/G{{id gp get/gp
+gp 1 add N A 18 mod S 18 idiv pl S get exec}loop}B/adv{cp add/cp X}B
+/chg{rw cp id gp 4 index getinterval putinterval A gp add/gp X adv}B/nd{
+/cp 0 N rw exit}B/lsh{rw cp 2 copy get A 0 eq{pop 1}{A 255 eq{pop 254}{
+A A add 255 and S 1 and or}ifelse}ifelse put 1 adv}B/rsh{rw cp 2 copy
+get A 0 eq{pop 128}{A 255 eq{pop 127}{A 2 idiv S 128 and or}ifelse}
+ifelse put 1 adv}B/clr{rw cp 2 index string putinterval adv}B/set{rw cp
+fillstr 0 4 index getinterval putinterval adv}B/fillstr 18 string 0 1 17
+{2 copy 255 put pop}for N/pl[{adv 1 chg}{adv 1 chg nd}{1 add chg}{1 add
+chg nd}{adv lsh}{adv lsh nd}{adv rsh}{adv rsh nd}{1 add adv}{/rc X nd}{
+1 add set}{1 add clr}{adv 2 chg}{adv 2 chg nd}{pop nd}]A{bind pop}
+forall N/D{/cc X A type/stringtype ne{]}if nn/base get cc ctr put nn
+/BitMaps get S ctr S sf 1 ne{A A length 1 sub A 2 index S get sf div put
+}if put/ctr ctr 1 add N}B/I{cc 1 add D}B/bop{userdict/bop-hook known{
+bop-hook}if/SI save N @rigin 0 0 moveto/V matrix currentmatrix A 1 get A
+mul exch 0 get A mul add .99 lt{/QV}{/RV}ifelse load def pop pop}N/eop{
+SI restore userdict/eop-hook known{eop-hook}if showpage}N/@start{
+userdict/start-hook known{start-hook}if pop/VResolution X/Resolution X
+1000 div/DVImag X/IEn 256 array N 2 string 0 1 255{IEn S A 360 add 36 4
+index cvrs cvn put}for pop 65781.76 div/vsize X 65781.76 div/hsize X}N
+/p{show}N/RMat[1 0 0 -1 0 0]N/BDot 260 string N/Rx 0 N/Ry 0 N/V{}B/RV/v{
+/Ry X/Rx X V}B statusdict begin/product where{pop false[(Display)(NeXT)
+(LaserWriter 16/600)]{A length product length le{A length product exch 0
+exch getinterval eq{pop true exit}if}{pop}ifelse}forall}{false}ifelse
+end{{gsave TR -.1 .1 TR 1 1 scale Rx Ry false RMat{BDot}imagemask
+grestore}}{{gsave TR -.1 .1 TR Rx Ry scale 1 1 false RMat{BDot}
+imagemask grestore}}ifelse B/QV{gsave newpath transform round exch round
+exch itransform moveto Rx 0 rlineto 0 Ry neg rlineto Rx neg 0 rlineto
+fill grestore}B/a{moveto}B/delta 0 N/tail{A/delta X 0 rmoveto}B/M{S p
+delta add tail}B/b{S p tail}B/c{-4 M}B/d{-3 M}B/e{-2 M}B/f{-1 M}B/g{0 M}
+B/h{1 M}B/i{2 M}B/j{3 M}B/k{4 M}B/w{0 rmoveto}B/l{p -4 w}B/m{p -3 w}B/n{
+p -2 w}B/o{p -1 w}B/q{p 1 w}B/r{p 2 w}B/s{p 3 w}B/t{p 4 w}B/x{0 S
+rmoveto}B/y{3 2 roll p a}B/bos{/SS save N}B/eos{SS restore}B end
+
+%%EndProcSet
+%%BeginProcSet: f7b6d320.enc
+% Thomas Esser, Dec 2002. public domain
+%
+% Encoding for:
+% cmb10 cmbx10 cmbx12 cmbx5 cmbx6 cmbx7 cmbx8 cmbx9 cmbxsl10
+% cmdunh10 cmr10 cmr12 cmr17cmr6 cmr7 cmr8 cmr9 cmsl10 cmsl12 cmsl8
+% cmsl9 cmss10cmss12 cmss17 cmss8 cmss9 cmssbx10 cmssdc10 cmssi10
+% cmssi12 cmssi17 cmssi8cmssi9 cmssq8 cmssqi8 cmvtt10
+%
+/TeXf7b6d320Encoding [
+/Gamma /Delta /Theta /Lambda /Xi /Pi /Sigma /Upsilon /Phi /Psi /Omega
+/ff /fi /fl /ffi /ffl /dotlessi /dotlessj /grave /acute /caron /breve
+/macron /ring /cedilla /germandbls /ae /oe /oslash /AE /OE /Oslash
+/suppress /exclam /quotedblright /numbersign /dollar /percent /ampersand
+/quoteright /parenleft /parenright /asterisk /plus /comma /hyphen
+/period /slash /zero /one /two /three /four /five /six /seven /eight
+/nine /colon /semicolon /exclamdown /equal /questiondown /question /at
+/A /B /C /D /E /F /G /H /I /J /K /L /M /N /O /P /Q /R /S /T /U /V /W /X
+/Y /Z /bracketleft /quotedblleft /bracketright /circumflex /dotaccent
+/quoteleft /a /b /c /d /e /f /g /h /i /j /k /l /m /n /o /p /q /r /s /t /u
+/v /w /x /y /z /endash /emdash /hungarumlaut /tilde /dieresis /suppress
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /space
+/Gamma /Delta /Theta /Lambda /Xi /Pi /Sigma /Upsilon /Phi /Psi /.notdef
+/.notdef /Omega /ff /fi /fl /ffi /ffl /dotlessi /dotlessj /grave /acute
+/caron /breve /macron /ring /cedilla /germandbls /ae /oe /oslash /AE
+/OE /Oslash /suppress /dieresis /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+] def
+
+%%EndProcSet
+%%BeginProcSet: aae443f0.enc
+% Thomas Esser, Dec 2002. public domain
+%
+% Encoding for:
+% cmmi10 cmmi12 cmmi5 cmmi6 cmmi7 cmmi8 cmmi9 cmmib10
+%
+/TeXaae443f0Encoding [
+/Gamma /Delta /Theta /Lambda /Xi /Pi /Sigma /Upsilon /Phi /Psi /Omega
+/alpha /beta /gamma /delta /epsilon1 /zeta /eta /theta /iota /kappa
+/lambda /mu /nu /xi /pi /rho /sigma /tau /upsilon /phi /chi /psi
+/omega /epsilon /theta1 /pi1 /rho1 /sigma1 /phi1 /arrowlefttophalf
+/arrowleftbothalf /arrowrighttophalf /arrowrightbothalf /arrowhookleft
+/arrowhookright /triangleright /triangleleft /zerooldstyle /oneoldstyle
+/twooldstyle /threeoldstyle /fouroldstyle /fiveoldstyle /sixoldstyle
+/sevenoldstyle /eightoldstyle /nineoldstyle /period /comma /less /slash
+/greater /star /partialdiff /A /B /C /D /E /F /G /H /I /J /K /L /M /N
+/O /P /Q /R /S /T /U /V /W /X /Y /Z /flat /natural /sharp /slurbelow
+/slurabove /lscript /a /b /c /d /e /f /g /h /i /j /k /l /m /n /o /p
+/q /r /s /t /u /v /w /x /y /z /dotlessi /dotlessj /weierstrass /vector
+/tie /psi /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/space /Gamma /Delta /Theta /Lambda /Xi /Pi /Sigma /Upsilon /Phi /Psi
+/.notdef /.notdef /Omega /alpha /beta /gamma /delta /epsilon1 /zeta /eta
+/theta /iota /kappa /lambda /mu /nu /xi /pi /rho /sigma /tau /upsilon
+/phi /chi /psi /tie /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef
+] def
+
+%%EndProcSet
+%%BeginProcSet: bbad153f.enc
+% Thomas Esser, Dec 2002. public domain
+%
+% Encoding for:
+% cmsy10 cmsy5 cmsy6 cmsy7 cmsy8 cmsy9
+%
+/TeXbbad153fEncoding [
+/minus /periodcentered /multiply /asteriskmath /divide /diamondmath
+/plusminus /minusplus /circleplus /circleminus /circlemultiply
+/circledivide /circledot /circlecopyrt /openbullet /bullet
+/equivasymptotic /equivalence /reflexsubset /reflexsuperset /lessequal
+/greaterequal /precedesequal /followsequal /similar /approxequal
+/propersubset /propersuperset /lessmuch /greatermuch /precedes /follows
+/arrowleft /arrowright /arrowup /arrowdown /arrowboth /arrownortheast
+/arrowsoutheast /similarequal /arrowdblleft /arrowdblright /arrowdblup
+/arrowdbldown /arrowdblboth /arrownorthwest /arrowsouthwest /proportional
+/prime /infinity /element /owner /triangle /triangleinv /negationslash
+/mapsto /universal /existential /logicalnot /emptyset /Rfractur /Ifractur
+/latticetop /perpendicular /aleph /A /B /C /D /E /F /G /H /I /J /K
+/L /M /N /O /P /Q /R /S /T /U /V /W /X /Y /Z /union /intersection
+/unionmulti /logicaland /logicalor /turnstileleft /turnstileright
+/floorleft /floorright /ceilingleft /ceilingright /braceleft /braceright
+/angbracketleft /angbracketright /bar /bardbl /arrowbothv /arrowdblbothv
+/backslash /wreathproduct /radical /coproduct /nabla /integral
+/unionsq /intersectionsq /subsetsqequal /supersetsqequal /section
+/dagger /daggerdbl /paragraph /club /diamond /heart /spade /arrowleft
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/minus /periodcentered /multiply /asteriskmath /divide /diamondmath
+/plusminus /minusplus /circleplus /circleminus /.notdef /.notdef
+/circlemultiply /circledivide /circledot /circlecopyrt /openbullet
+/bullet /equivasymptotic /equivalence /reflexsubset /reflexsuperset
+/lessequal /greaterequal /precedesequal /followsequal /similar
+/approxequal /propersubset /propersuperset /lessmuch /greatermuch
+/precedes /follows /arrowleft /spade /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+/.notdef /.notdef /.notdef /.notdef /.notdef /.notdef /.notdef
+] def
+
+%%EndProcSet
+%%BeginProcSet: texps.pro
+%!
+TeXDict begin/rf{findfont dup length 1 add dict begin{1 index/FID ne 2
+index/UniqueID ne and{def}{pop pop}ifelse}forall[1 index 0 6 -1 roll
+exec 0 exch 5 -1 roll VResolution Resolution div mul neg 0 0]FontType 0
+ne{/Metrics exch def dict begin Encoding{exch dup type/integertype ne{
+pop pop 1 sub dup 0 le{pop}{[}ifelse}{FontMatrix 0 get div Metrics 0 get
+div def}ifelse}forall Metrics/Metrics currentdict end def}{{1 index type
+/nametype eq{exit}if exch pop}loop}ifelse[2 index currentdict end
+definefont 3 -1 roll makefont/setfont cvx]cvx def}def/ObliqueSlant{dup
+sin S cos div neg}B/SlantFont{4 index mul add}def/ExtendFont{3 -1 roll
+mul exch}def/ReEncodeFont{CharStrings rcheck{/Encoding false def dup[
+exch{dup CharStrings exch known not{pop/.notdef/Encoding true def}if}
+forall Encoding{]exch pop}{cleartomark}ifelse}if/Encoding exch def}def
+end
+
+%%EndProcSet
+%%BeginFont: CMR8
+%!PS-AdobeFont-1.1: CMR8 1.0
+%%CreationDate: 1991 Aug 20 16:39:40
+% Copyright (C) 1997 American Mathematical Society. All Rights Reserved.
+11 dict begin
+/FontInfo 7 dict dup begin
+/version (1.0) readonly def
+/Notice (Copyright (C) 1997 American Mathematical Society. All Rights Reserved) readonly def
+/FullName (CMR8) readonly def
+/FamilyName (Computer Modern) readonly def
+/Weight (Medium) readonly def
+/ItalicAngle 0 def
+/isFixedPitch false def
+end readonly def
+/FontName /CMR8 def
+/PaintType 0 def
+/FontType 1 def
+/FontMatrix [0.001 0 0 0.001 0 0] readonly def
+/Encoding 256 array
+0 1 255 {1 index exch /.notdef put} for
+dup 0 /.notdef put
+readonly def
+/FontBBox{-36 -250 1070 750}readonly def
+/UniqueID 5000791 def
+currentdict end
+currentfile eexec
+D9D66F633B846A97B686A97E45A3D0AA052A014267B7904EB3C0D3BD0B83D891
+016CA6CA4B712ADEB258FAAB9A130EE605E61F77FC1B738ABC7C51CD46EF8171
+9098D5FEE67660E69A7AB91B58F29A4D79E57022F783EB0FBBB6D4F4EC35014F
+D2DECBA99459A4C59DF0C6EBA150284454E707DC2100C15B76B4C19B84363758
+469A6C558785B226332152109871A9883487DD7710949204DDCF837E6A8708B8
+2BDBF16FBC7512FAA308A093FE5CF4E9D2405B169CD5365D6ECED5D768D66D6C
+68618B8C482B341F8CA38E9BB9BAFCFAAD9C2F3FD033B62690986ED43D9C9361
+3645B82392D5CAE11A7CB49D7E2E82DCD485CBA1772CE422BB1D7283AD675B65
+48A7EA0069A883EC1DAA3E1F9ECE7586D6CF0A128CD557C7E5D7AA3EA97EBAD3
+9619D1BFCF4A6D64768741EDEA0A5B0EFBBF347CDCBE2E03D756967A16B613DB
+0FC45FA2A3312E0C46A5FD0466AB097C58FFEEC40601B8395E52775D0AFCD7DB
+8AB317333110531E5C44A4CB4B5ACD571A1A60960B15E450948A5EEA14DD330F
+EA209265DB8E1A1FC80DCD3860323FD26C113B041A88C88A21655878680A4466
+FA10403D24BB97152A49B842C180E4D258C9D48F21D057782D90623116830BA3
+9902B3C5F2F2DD01433B0D7099C07DBDE268D0FFED5169BCD03D48B2F058AD62
+D8678C626DC7A3F352152C99BA963EF95F8AD11DB8B0D351210A17E4C2C55AD8
+9EB64172935D3C20A398F3EEEEC31551966A7438EF3FEE422C6D4E05337620D5
+ACC7B52BED984BFAAD36EF9D20748B05D07BE4414A63975125D272FAD83F76E6
+10FFF8363014BE526D580873C5A42B70FA911EC7B86905F13AFE55EB0273F582
+83158793B8CC296B8DE1DCCF1250FD57CB0E035C7EDA3B0092ED940D37A05493
+2EC54E09B984FCA4AB7D2EA182BCF1263AA244B07EC0EA901C077A059F709F30
+4384CB5FA748F2054FAD9A7A43D4EA427918BD414F766531136B60C3477C6632
+BEFE3897B58C19276A301926C2AEF2756B367319772C9B201C49B4D935A8267B
+041D6F1783B6AEA4DAC4F5B3507D7032AA640AAB12E343A4E9BDCF419C04A721
+3888B25AF4E293AACED9A6BDC78E61DA1C424C6503CC1885F762BFD4563ABA4F
+D926227FAD2B2D4975835541E32E74F60B94D706659522BF4CD9E74EB5C62BAA
+D70DFE70BB4A82A00FF1248158C21082280E3370CB98291962F8082C963BFBCE
+9030B5F7FCB39E22EE6D254CC6204D30DCDA58496CD7587E6B571CFD6D346EB5
+76F3F4C8D42E88652652FD7A342F9424E7A3571337B2FCDA1D187290C5A18680
+FD13A692DB9FC79D736B7FF4FA62676763FC1CF35315AD3CA290FB06B3364955
+CAEBD2D7309223E16EC7F6E04FBEC7EB812154CBE37C1C22FF6AE5824ED4E33D
+2058385AB4EDE485ACFA4A78E8AC3700B38F52327846D0D3FE5A6E1BFBF9F631
+833AECE2DBD389C0A75E31BD74304C34B2E0A5D5681755ADACCF769723C77662
+844C222533FF950890A097CAC1108B8AFD1023FC9C4317BD8C5385F8D3B36633
+E894B01DF1A962F43C899FA685BF244A553ACA90FF7F72A78EC4A44976FCA087
+E048945DAED1F11D946951180CD7E8FF593E5DA359A69B08EB85D14A35628E48
+701FEB250728B290A58CF989CE74E4A3DCBC2A7CA588545AF7E826865567AA66
+591182C6447C9BB4E279064714FDBEA796124CF703D0AF2E51D51BEACD5F7EE1
+8796B6ECF00E5C792D74F203E1029523BA8D267CAE420CB1D8071F086DF5E1D8
+6C77FC2751DD98F5451C0E48707FD1E32AA34B71FD920882FD3415AEA1A867B0
+7410AF7994652691F7AAECBE6BBEC44EF07B3E980814E1E7A7BE79D74B9B75EE
+155E90F95ACDD4009635B2FFFC613505E0850E2EB61DF43A5701763E7656F59A
+5AF46CD0FCE5297617AC0BBD0791E1B04C0BF7B8EA3B7DE73F714A5453DDF843
+1E4809C1945FD7AA617B22DC9C98261FFFDD1447F80AE125664AFD655FCC94CA
+D596ECE719774096D2CC1752B0333AC5A8C2B197C22CF8BB852CF1746F4903F0
+BFAE9FE8DADCC5A6B946ABCF4365912D32BF0CCDCBD7826BBEECAEBEC6FCA436
+B0338EB9206C52255AB88DCDD65B9E2E67
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+cleartomark
+%%EndFont
+%%BeginFont: CMSY8
+%!PS-AdobeFont-1.1: CMSY8 1.0
+%%CreationDate: 1991 Aug 15 07:22:10
+% Copyright (C) 1997 American Mathematical Society. All Rights Reserved.
+11 dict begin
+/FontInfo 7 dict dup begin
+/version (1.0) readonly def
+/Notice (Copyright (C) 1997 American Mathematical Society. All Rights Reserved) readonly def
+/FullName (CMSY8) readonly def
+/FamilyName (Computer Modern) readonly def
+/Weight (Medium) readonly def
+/ItalicAngle -14.035 def
+/isFixedPitch false def
+end readonly def
+/FontName /CMSY8 def
+/PaintType 0 def
+/FontType 1 def
+/FontMatrix [0.001 0 0 0.001 0 0] readonly def
+/Encoding 256 array
+0 1 255 {1 index exch /.notdef put} for
+dup 0 /.notdef put
+readonly def
+/FontBBox{-30 -955 1185 779}readonly def
+/UniqueID 5000818 def
+currentdict end
+currentfile eexec
+D9D66F633B846A97B686A97E45A3D0AA052F09F9C8ADE9D907C058B87E9B6964
+7D53359E51216774A4EAA1E2B58EC3176BD1184A633B951372B4198D4E8C5EF4
+A213ACB58AA0A658908035BF2ED8531779838A960DFE2B27EA49C37156989C85
+E21B3ABF72E39A89232CD9F4237FC80C9E64E8425AA3BEF7DED60B122A52922A
+221A37D9A807DD01161779DDE7D5FC1B2109839E5B52DFBB2A7C1B5D8E7E8AA0
+5B10EA43D6A8ED61AF5B23D49920D8F79DAB6A59062134D84AC0100187A6CD1F
+80F5DDD9D222ACB1C23326A7656A635C4A241CCD32CBFDF8363206B8AA36E107
+1477F5496111E055C7491002AFF272E46ECC46422F0380D093284870022523FB
+DA1716CC4F2E2CCAD5F173FCBE6EDDB874AD255CD5E5C0F86214393FCB5F5C20
+9C3C2BB5886E36FC3CCC21483C3AC193485A46E9D22BD7201894E4D45ADD9BF1
+CC5CF6A5010B5654AC0BE0DA903DB563B13840BA3015F72E51E3BC80156388BA
+F83C7D393392BCBC227771CDCB976E93302531886DDA73EBC9178917EFD0C20B
+133F1E59BBEF76037FEFAD9A067DABC0F3B4064CC55152A17D667AC39B66A4E6
+1C2F8E36A56B569BB14584735C1DFDA557D5DD1D22D96D617D824C864B00E027
+E17B8C2A768E15D69B2EEBDFEB468E287818A7C802A13E0329FDD3D256EA8EA7
+30CAA928BC525604E92B2B78132360EE7025D54F427D59D81B8BAF192A6F79C7
+18B98500457845855BD50E25B93FF6ABD58898CF45746265DDD4D038BFF2B804
+A46D8AD01D35DCB23B7E35F7139260EF998E51E1E93EC3AA295F7BCEB8196803
+6B334A5E472169EA7D3E0CE6B2DBEE318B12F20BBAE8A7F314F2A2
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+cleartomark
+%%EndFont
+%%BeginFont: CMSY10
+%!PS-AdobeFont-1.1: CMSY10 1.0
+%%CreationDate: 1991 Aug 15 07:20:57
+% Copyright (C) 1997 American Mathematical Society. All Rights Reserved.
+11 dict begin
+/FontInfo 7 dict dup begin
+/version (1.0) readonly def
+/Notice (Copyright (C) 1997 American Mathematical Society. All Rights Reserved) readonly def
+/FullName (CMSY10) readonly def
+/FamilyName (Computer Modern) readonly def
+/Weight (Medium) readonly def
+/ItalicAngle -14.035 def
+/isFixedPitch false def
+end readonly def
+/FontName /CMSY10 def
+/PaintType 0 def
+/FontType 1 def
+/FontMatrix [0.001 0 0 0.001 0 0] readonly def
+/Encoding 256 array
+0 1 255 {1 index exch /.notdef put} for
+dup 0 /.notdef put
+readonly def
+/FontBBox{-29 -960 1116 775}readonly def
+/UniqueID 5000820 def
+currentdict end
+currentfile eexec
+D9D66F633B846A97B686A97E45A3D0AA052F09F9C8ADE9D907C058B87E9B6964
+7D53359E51216774A4EAA1E2B58EC3176BD1184A633B951372B4198D4E8C5EF4
+A213ACB58AA0A658908035BF2ED8531779838A960DFE2B27EA49C37156989C85
+E21B3ABF72E39A89232CD9F4237FC80C9E64E8425AA3BEF7DED60B122A52922A
+221A37D9A807DD01161779DDE7D31FF2B87F97C73D63EECDDA4C49501773468A
+27D1663E0B62F461F6E40A5D6676D1D12B51E641C1D4E8E2771864FC104F8CBF
+5B78EC1D88228725F1C453A678F58A7E1B7BD7CA700717D288EB8DA1F57C4F09
+0ABF1D42C5DDD0C384C7E22F8F8047BE1D4C1CC8E33368FB1AC82B4E96146730
+DE3302B2E6B819CB6AE455B1AF3187FFE8071AA57EF8A6616B9CB7941D44EC7A
+71A7BB3DF755178D7D2E4BB69859EFA4BBC30BD6BB1531133FD4D9438FF99F09
+4ECC068A324D75B5F696B8688EEB2F17E5ED34CCD6D047A4E3806D000C199D7C
+515DB70A8D4F6146FE068DC1E5DE8BC5703711DA090312BA3FC00A08C453C609
+C627A8B8802AA63AE69E2D766D27C2FB99E94FC4217B8395DE944F6BCE9E7765
+0DCCDB37CE4DCF1D7A9B6BAC985A1B72DF48ADDF3B23EBFC45889DB2A2B565A9
+6BB248468F054DA16C281CD76150241A0023E20C240C81C71C92674F2798015C
+49C57C7DD4CBEDB4FFCB7D86EBF82683D0EAB5F47991CF336250AA4786BF720E
+BD814E1A71628524C0DCDCB0F2E11BBB289A6889A4F1261AB79D920851B9FBBE
+325D615F381960510B1E8C77D4FED56BDF56D60BA5982EC39419F4097A7E2E43
+33EC6A080FAD9326E38AC3B2017CD6BD013118F4E74C78D0CF9BEA973493A004
+1F0F0539742DBDB319B8587855AC8A934DC3A5C86EC454E5E1B902183B750A56
+923649C81294F8DE99FCEAA794A4A5186C4BADCE68A4278328D76CAE58C4037B
+3392013DB55E57AC99E17FA0F68343BD518CFE48174CB1A9A31F70ADAC42E967
+4C29C1A4B1C565B094F20DEF0B883B377B5C727542CA79C688F0FD0F80DAB18C
+6F2073A1A4E0A37B33
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+cleartomark
+%%EndFont
+%%BeginFont: CMMI8
+%!PS-AdobeFont-1.1: CMMI8 1.100
+%%CreationDate: 1996 Jul 23 07:53:54
+% Copyright (C) 1997 American Mathematical Society. All Rights Reserved.
+11 dict begin
+/FontInfo 7 dict dup begin
+/version (1.100) readonly def
+/Notice (Copyright (C) 1997 American Mathematical Society. All Rights Reserved) readonly def
+/FullName (CMMI8) readonly def
+/FamilyName (Computer Modern) readonly def
+/Weight (Medium) readonly def
+/ItalicAngle -14.04 def
+/isFixedPitch false def
+end readonly def
+/FontName /CMMI8 def
+/PaintType 0 def
+/FontType 1 def
+/FontMatrix [0.001 0 0 0.001 0 0] readonly def
+/Encoding 256 array
+0 1 255 {1 index exch /.notdef put} for
+dup 0 /.notdef put
+readonly def
+/FontBBox{-24 -250 1110 750}readonly def
+/UniqueID 5087383 def
+currentdict end
+currentfile eexec
+D9D66F633B846A97B686A97E45A3D0AA0529731C99A784CCBE85B4993B2EEBDE
+3B12D472B7CF54651EF21185116A69AB1096ED4BAD2F646635E019B6417CC77B
+532F85D811C70D1429A19A5307EF63EB5C5E02C89FC6C20F6D9D89E7D91FE470
+B72BEFDA23F5DF76BE05AF4CE93137A219ED8A04A9D7D6FDF37E6B7FCDE0D90B
+986423E5960A5D9FBB4C956556E8DF90CBFAEC476FA36FD9A5C8175C9AF513FE
+D919C2DDD26BDC0D99398B9F4D03D6A8F05B47AF95EF28A9C561DBDC98C47CF5
+5250011D19E9366EB6FD153D3A100CAA6212E3D5D93990737F8D326D347B7EDC
+4391C9DF440285B8FC159D0E98D4258FC57892DDF753642CD526A96ACEDA4120
+788F22B1D09F149794E66DD1AC2C2B3BC6FEC59D626F427CD5AE9C54C7F78F62
+C36F49B3C2E5E62AFB56DCEE87445A12A942C14AE618D1FE1B11A9CF9FAA1F32
+617B598CE5058715EF3051E228F72F651040AD99A741F247C68007E68C84E9D1
+D0BF99AA5D777D88A7D3CED2EA67F4AE61E8BC0495E7DA382E82DDB2B009DD63
+532C74E3BE5EC555A014BCBB6AB31B8286D7712E0E926F8696830672B8214E9B
+5D0740C16ADF0AFD47C4938F373575C6CA91E46D88DE24E682DEC44B57EA8AF8
+4E57D45646073250D82C4B50CBBB0B369932618301F3D4186277103B53B3C9E6
+DB42D6B30115F67B9D078220D5752644930643BDF9FACF684EBE13E39B65055E
+B1BD054C324962025EC79E1D155936FE32D9F2224353F2A46C3558EF216F6BB2
+A304BAF752BEEC36C4440B556AEFECF454BA7CBBA7537BCB10EBC21047333A89
+8936419D857CD9F59EBA20B0A3D9BA4A0D3395336B4CDA4BA6451B6E4D1370FA
+D9BDABB7F271BC1C6C48D9DF1E5A6FAE788F5609DE3C48D47A67097C547D9817
+AD389FD24494D365D62017182ABEBCC79E7127AC417197B6109301EFEFD60FA3
+94298CD8A0CFA8861BE3814E05368B4CDD4ED901F4CCAA145D0C7A5A6080618A
+96D5B23FFA9F70DCA3BC0C0C9B9469A96F50117D57FAFACDC9F059CFD81B50B9
+A1654593E5FFA2E906F3C0B6861BC85CEBAB83464F26C049EC0D70F25DC9E399
+FCBA99C4E026953949028A6B7E92E0A9F922D3DBDBBD07DA1B1C76BD1D05ADA2
+A309199A8C278B6CE76AA34ECA5FAC0370469626C0114ADE1B875C14839D7B00
+BBB91354CA5E4575EC2E27F559F1BDD6305D043AE923E00E99172BDDA92882CB
+A8D0DC8C3FF65AD700C8FEAA5C33D10656001BA8F59EAC3C63F7E7C747266AD2
+D2DB772220F044DA33C6E301019580EDA7668E75737975DD62412C16BD556FC0
+CAFE5B885ED44D38ED3BEF582F545E8DA514B6F0DBBE65F9F7DEB94807293733
+9C51915FB1C1D57A42C76237CB6435913CA8BF70928EACF1470D14C4434B6DA0
+DFEFB41A77E69DD471EBD8618CC2F64A25F18287E155D99C397141CC77C09D32
+01F5D18558934FB948ABB429489210677DED577F2D5CAE994CEAB43B66BBD5F6
+D5EAB350F0B99686ECA3A0C47B90647E998CB8EBC49362B399005BEC78A94C3F
+8CDBAE64E3C837FF7442AE2D3D9611456875E78341C2ABD7587BFD2022D63C3B
+EF519ACE2D321793F9F19FAFD9D3A5A0C7434FB8AA6283D0EB17C92119D8BC5A
+3B77B6F837DC96D745360219E1DDCE6D5290B19243C13ABA68BE2BBF7D97296C
+2D90263DA77894200097B92F3E773FD650471F4375075CD1A5E70F7419DA5B12
+77AD99261D9C9625317FEB04A47392F73A2FC3CC6F6ACA7E98CD202750AA4BEE
+6F5CF24BB2F33E14AC2E6B9226E485CFB4DC4080CF81D8F7BA6E2EEF9057A57B
+9B3B6AA7FB65D2E02E0B48D893979D5BB5948BB6444D93A69AF4DABB918F7936
+259D18E1B2E096CD5D193F834C613D4F079573ACF2976533D0359F755042D453
+7863161521A303DF3E06A7A591BB71940FE92A1E9440E9797126F9F8C5A3ACC4
+C14CF13AE4BAB28DEC6DCA0C49FB95CD0A9719AB8AC5CD125C7B781A8A4F9220
+A5BD7DB60CD7F3E4C1FC1F175134CF64E543D71777C2C02A08C0819F185CBB7F
+E96044035DD20237C8222E2709934A4E7F598095337A1F29E560A803768B6691
+0A18271D3AFCCA2C175AA4F8FD0341CE507D9BAEF5E4E6A0024BC67A78A24315
+8AFB9F754EDD0949FC67724697634242E72B3B112A12F05F472347F32FDA23A6
+97081409E76E6F0F21EA333B2DE4EDEE693E9D32DEAAA24DC3F8BC48BA906B45
+5A4F3B3E29950AFEDBA755F991E76EC75C0A51138C3F9AF5529327D1C7F01889
+88C9DD21CC94F5B7CBE10E8E0BBD79851F0B3F0CAAC11C03485273B95032B6F4
+ABB7ED24E278371AADCF66838BA0185EC84370725D34190791C1204CFCD5480E
+9DF4BF9C034C59204191D1214CB6986F976B7ECFC90A04135A8027617B590176
+7A9400A651FD507A059FE4BD6BCD2BCB2E7A841F7C50732C65054F48FE6343F5
+DE44D7FBAEAF694222FB9552C5D88BC123E2F8DA39D92850D494FC67A7144DEF
+37BA0DFA1EEEC1146EE4583AE6838C3E33884CC59410D8E4083AAE6679C0DD1A
+3AD3A1B703594393A97C63E4675F22C9D8BBEA96E98152D7F28D162C8F9F9B55
+D3888ED78EDF9ABE16174B5F254DD42326D3A00DEB1D60BD07F2388091440FED
+3A7CBD2CE63496C090CBD67A3DA312F75CCAEE26DC6B2C544557D1DF8A24F861
+9ED8EA50327E2072820B2A91D019AD1942D1323A6F7434B82D271D62F5D77CDB
+4CF8DD29DC6535E86C7605B95C73CED36A185217B7EE1AF01145C9DD73C76299
+EF74CA08F9132564C113EBD8E68FFC47DA0A56BAB15521648455963E14EC0595
+B31F24CA320093F551F5101235C694D907EA3A7D1B8B13F97E5E502ABFA33047
+858935BA1BB0A5CAE4E05BE9EAFC0D592601A9DDB795610F8823E9E33E86DD26
+51F76304653850732DCDBB0D680FD42B601511BF76530F4C1D64B8305A4D31AB
+03529DFCC2720AD5664A89E77F56816DD9A93D44466E309C6745A8A3516A5595
+9DE56D15878125AD7E9FD18DFEAE2BFEAA59C073038E46FB163D45A79FA4511B
+1147B60130FE628F33010467EF17AE8DDEF413C00FC5A718CE1F22710BCDA86D
+556617F34D03D2E085D1DDDC467AE797BA5C2E8ACE5403FF69AD3AE0AE9B3075
+FCD88B5CBECF697887F90F14AFA78A79EE40D1D5E75417C0E51E3B738881A801
+EC96CA7B08CF4C401FE7B10AF9D078793A8B8A12DB3D3DB445D90B3643C34510
+40D1713741925F49D6E083FDAB4F9CF0CDD1E33A6049E803392FDFC637E80C55
+01E0C40EBA0A2426F07CE76EA9713478D1B135236D7C28163AD2C4E9B513D171
+D01EA9719FF103318B0DA7DF377E06EB8F09F73C17B59011B1680D426BA7FF7F
+82B707224D68D7D1FC002820B229458982F2EBD4B4338BF52D557C660E8B1203
+5AEAA3BDF74DFEDF90288A8FDD4A6C65E11A9F2B98F1746404385710CAFDBD1E
+886CA42B3817D664EDB38D1538F8F11CD447FB038712625B5DDB416AA1C60708
+104FFEBD26AA1B3F851E06C6E32D60F1BA475EC0655BD42EFFB691A8D78E105F
+5C7A868E48239BE2A08A35998421F78DC062FCD09867C7AFF1B38DF5A27622A8
+9CC6BB52DFE1B6C2D99617726D5B972F2FD4A115338F52404D333958083B90DD
+3DDB37289143C58E42A071FDF8BB018C09A27134E73B99E9A418389BAEE12D22
+4D9D06070B006D1C3CA08E9E2EB85FD25CD21E74D520894F5BC78420CBAFA156
+FD2059117128AE517ACBC705FD868F17FE7E7C9D5F8B3BDD77F3089BEEAC4ADC
+2251879926D5F61B6DBB69A8ACCF402C1816D40040F80CCB0C72E3E69392E58B
+4D03086CCF84BCD45A03A5D9A9C879F734874946AF81
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+cleartomark
+%%EndFont
+%%BeginFont: CMR10
+%!PS-AdobeFont-1.1: CMR10 1.00B
+%%CreationDate: 1992 Feb 19 19:54:52
+% Copyright (C) 1997 American Mathematical Society. All Rights Reserved.
+11 dict begin
+/FontInfo 7 dict dup begin
+/version (1.00B) readonly def
+/Notice (Copyright (C) 1997 American Mathematical Society. All Rights Reserved) readonly def
+/FullName (CMR10) readonly def
+/FamilyName (Computer Modern) readonly def
+/Weight (Medium) readonly def
+/ItalicAngle 0 def
+/isFixedPitch false def
+end readonly def
+/FontName /CMR10 def
+/PaintType 0 def
+/FontType 1 def
+/FontMatrix [0.001 0 0 0.001 0 0] readonly def
+/Encoding 256 array
+0 1 255 {1 index exch /.notdef put} for
+dup 0 /.notdef put
+readonly def
+/FontBBox{-251 -250 1009 969}readonly def
+/UniqueID 5000793 def
+currentdict end
+currentfile eexec
+D9D66F633B846A97B686A97E45A3D0AA052A014267B7904EB3C0D3BD0B83D891
+016CA6CA4B712ADEB258FAAB9A130EE605E61F77FC1B738ABC7C51CD46EF8171
+9098D5FEE67660E69A7AB91B58F29A4D79E57022F783EB0FBBB6D4F4EC35014F
+D2DECBA99459A4C59DF0C6EBA150284454E707DC2100C15B76B4C19B84363758
+469A6C558785B226332152109871A9883487DD7710949204DDCF837E6A8708B8
+2BDBF16FBC7512FAA308A093FE5CF7158F1163BC1F3352E22A1452E73FECA8A4
+87100FB1FFC4C8AF409B2067537220E605DA0852CA49839E1386AF9D7A1A455F
+D1F017CE45884D76EF2CB9BC5821FD25365DDEA6E45F332B5F68A44AD8A530F0
+92A36FAC8D27F9087AFEEA2096F839A2BC4B937F24E080EF7C0F9374A18D565C
+295A05210DB96A23175AC59A9BD0147A310EF49C551A417E0A22703F94FF7B75
+409A5D417DA6730A69E310FA6A4229FC7E4F620B0FC4C63C50E99E179EB51E4C
+4BC45217722F1E8E40F1E1428E792EAFE05C5A50D38C52114DFCD24D54027CBF
+2512DD116F0463DE4052A7AD53B641A27E81E481947884CE35661B49153FA19E
+0A2A860C7B61558671303DE6AE06A80E4E450E17067676E6BBB42A9A24ACBC3E
+B0CA7B7A3BFEA84FED39CCFB6D545BB2BCC49E5E16976407AB9D94556CD4F008
+24EF579B6800B6DC3AAF840B3FC6822872368E3B4274DD06CA36AF8F6346C11B
+43C772CC242F3B212C4BD7018D71A1A74C9A94ED0093A5FB6557F4E0751047AF
+D72098ECA301B8AE68110F983796E581F106144951DF5B750432A230FDA3B575
+5A38B5E7972AABC12306A01A99FCF8189D71B8DBF49550BAEA9CF1B97CBFC7CC
+96498ECC938B1A1710B670657DE923A659DB8757147B140A48067328E7E3F9C3
+7D1888B284904301450CE0BC15EEEA00E48CCD6388F3FC3BEFD8D9C400015B65
+0F2F536D035626B1FF0A69D732C7A1836D635C30C06BED4327737029E5BA5830
+B9E88A4024C3326AD2F34F47B54739B48825AD6699F7D117EA4C4AEC4440BF6D
+AA0099DEFD326235965C63647921828BF269ECC87A2B1C8CAD6C78B6E561B007
+97BE2BC7CA32B4534075F6491BE959D1F635463E71679E527F4F456F774B2AF8
+FEF3D8C63B283796A9AD847424B4E6508546C36223A3B17EB82A56592F27FC27
+F1D49D5FF4BBC0E16231807AF7E195AA7D0D01C7566243448B222D57B811EAE6
+DE9370F84E207DC9BEC731AD6040FD9B804FA14CA264B73136F9AF34390319F6
+A543D5D4D7FDDDF2F76651E557683614521110DEC1CCAC426117DDA7D6CF1B8B
+7879B21FDC78BAB3C944BCDCD6A65B67F3692F0A8D5E36FB783A63D4FBC9842E
+2CBC2720A7206F42A99AEC79FBBA92A27965AF40A71E05E4BA8D7FC58C828226
+7A8B47241F73590CE626EDA930A1D18D758636250D9D6A09E133378CD725F586
+F40320EF186DDF87AF8C2F5470EADFF39C60384EB4DD859794CDD6E4AA8E8A7A
+18E246579B5FD876A7C0AFB4D935B09F36A74412957D9BE1DA54DAC6239E5574
+A9AB77DC4A6CFBAF23EC30579BA578299EA68B0BFC5B9B4EAA4D3CADB581AC57
+6236C14F8F8D56406DC427BBE7913BA2990D22F36D646D3D31A6DB69512AACF4
+0A46538C40FE5F6D369871E017069547B6CD9793652C2EE8EF391B88DC95DC58
+DD90F27F47D996223903C9A02FB30DB43963F1DD6EAAE43F5E43F19B54FCA3EB
+114CFCC68A6B5106C697C5876011047F6679201540AD102C571BF5DD26AA443B
+6B73E80D855F5EBFF6A5BCEC484EFF7634D31618ED926015CF883D323F319D95
+D79BEA4091291DAE4ED9091359807685D980DDE8915E8718CD2CEA9D364B2E4C
+3DB1B8C4238F51C1E63C26E8A0D90EF1D4D2C9B4D52AFE11892F2CCF26DEBEF2
+20ECC90AD1EC92C7A2F4D1A838F1F3750B8A5E51E58568B2DE12D36CAB055DFC
+7D065C746CE7DAE1F2834FAF9CBE539EC06CABC6498495C913CE0EFA7FA87088
+F7376129F3F8A533334B624CF62D9013101548D93B43BA92D6533CA56769149F
+2A1B8CA40154EE515385DEC0F987AFA8E058647D2E3FF551AD08932C59B0BB2E
+AE0EB03749B40B6A4B3B94F24FB13CB7D83968F3F151438223CB1115379D36AC
+F35CC99FD1C9E059678560648FC797A6D22500F8CA1DD90BB7149E3C32ADAB69
+D9AB35CE944B6ADE05782CF731A759F68E9791D19CECFB26421B070877E6BC65
+EC82F51B81F2F939694DA32E41557DE6B7D055BA949C01CD34FC9B62B5E71B22
+81BBED290647C1B584FD784D9CDFFF89F65EF025BDC7838AFABD55154B54F8AF
+803F70491A0641ADE28AE83760C79C08C1265830384BDF4FF32CC5BA82DFAEC4
+C878C506E54F4678652D5639AFD082E4051871DACBE3E5A2569F2AB22CF952C5
+D4AFAE065DC9E9EB650262CCC5BAD17F6A43B60AF8645D78AE47A18CB6B14E6F
+7D9D5891B53F6532D0FB83540033FFA39DB14BFCE9EC0CF8E1487F47415380B4
+36D141090BCD9BCCB084B53031803F82772A752012896F122727523CB8088760
+48521F937A458A63396A578F0213154792B5712012EDEC8689DB24B872A9241E
+57321D17F92AB30E532ECBE5BED8990B04B6F7E0A7244A1C02889BBCE8F8A40E
+DFF04FD5CC893C2048EA9E49E7DA9443173E0163A0B6D795D64EFAE6C9BBAEA9
+6C74613CFDFE1D91C694230E836DF2E0B819B37FBFA55654144437DE084A959A
+B21031C540D15A1C1B29BE2725F851EC9C25C25B4BB10A8FEBCBC8630A6C1EA8
+30E3D1A60C87065936608BBC90F2819CA08A087D084DB86E6D6A17E38C3F6211
+323C27A6D8454A9D7C305668391BE3D1AB1D77AC83CFC32D636ECF5A30BE9654
+86F4996A56D946FE49AC0143B8286D771430F14E17E897166513559D7AD5FD38
+2AE132EDB60884004902D2533CB1DCBA729DCA229757DB91C5A66C3376DF806A
+B04BBF94DDE9B389B57B5A4408B16FF20E3A3FB3C9B26A785B2C951C021A5448
+23CB1BFBC0574EE8F79226965BAEA48452056A06FE7A3CB6B1814E7E5C35A3F2
+9C53C075238C067EC20AD746802ECA45F0A56988BE871C215FBEA7D023D1CD16
+13812A4DDB84B475D5C8A53CBFFD7CA5FFAE9E305849E32688C051DD752AB420
+5249D7457A1C0246665E4723F0314D25B4CF113D413997A8F3342F679A98D700
+CAA06A94034801EB398C22B881A32C0937775E7272F3CB9545B8B853BF33312B
+DB31DB13BB2612C588CA6226E5E83A60336E890B41CDD8AC4904B331A71CA54F
+388A28DC39EAF3DE36B93E94F8379C657EEAD2B35CD15EE02DCB85EF39D0659D
+A518D6CC1B1B0D5A8E8D38253F2B9F70C0974E97A6B7EDA9BE9DB17E883C4144
+CD79F0365860AFAED12A79FC2C2813504CCEBAE832674780D63035505049D180
+E94A58BAB54F30A3E657723DD05976639A8C93552EF8694DB21193EC54DCAECF
+101A97FD8E032F0B0A5EC7CC6436E53057BAD3D8B1E6ED7DB5B36511262EDA93
+DB7D954FCA90D8F8ECBB7360F2CFC91DCAF1ED221B5C3C3C59DD3726B02E9EFD
+8DFA05D203517366A1490A8DDB93FDA0FF4117B6770A26B9A128664F3D9F3AB4
+7828A74D81841C650934B313CA4CF61F37312725966E5F2B00D8E8600898D4ED
+F191FAF4D7332B664EFF2D47F923C160B613E907889ED30DA31B34E6C3264024
+CC5BD30F93D3760EABE6327EC7B915AD718CA093A1C64EB99B096CF8C60ED3CB
+4FACE8CB62E5248E5AD5EB26A9CC36621A2D5033DF7D10D57D13429E754CF893
+F3FA71593680E72CD076222EE06CA691D85A3F2AFD3D7488E4E80F4EBADF104D
+C2EF604E0A0BAB4051D78BDB663145DC30A97BE8AEF15825FF7A7BDAAC8D819D
+D65DCBB143EDF84C3A6752DED245E482A1BD2BB21C7839F80663B8D678C4CBA6
+D8E04C33759E9969D59AA054F848052AA2BDBFE4A5E4FCF21C8627A8F40C8570
+787551EEB264FB59E43EF463367320EA5FDF14AE28645F806C4FA26C82A1ECAE
+A8430CE67571A49302F752663938F40654A2A7BEF101F9D848C974473D27E170
+C39247A1268DC9CE96995F4228A594B7202144B2637647FAF4EDCB9332B4009F
+F29FD1B050974D0D2CA2DB143831A3A94B8B979AEB8FA6BB4C72C424850A1ECA
+26BD6FD46E5F868C7DE6D3C2A08ADDB899D3F8C3A49F2BB0A7A527A0FFC3853A
+E9CCE928E7DE3E4AA20FC99BC1BB4616C7E9172F04DD5F9DBB4840F522049091
+EBD97B58F0D307E0E800ABBACF357DFDF4D70B80AE43C8AB548BAA11641F5C2D
+AC6591432063B78AA494C8CFE11E23E78D7D6DD3A66D17560291812858AF10BB
+761440D9276FAEC0EB2CE7050BC17FCC10063DE84FD4E63F09F8BF4EE610EC6C
+3932E3EB5DF5189720A713F9787C58AC13610E497B875F3028654A88D310FD78
+F25E3DD5986A7CA7A028BE7429B1DF7112E20D34377487DB3BD19068B6AF09C5
+274FBBE02CB49AD514F8FA56127096D9B8C12FED24D66ABEE35498BCDD20B8F7
+13A8B7B37D202FCC044B7456D7C57B4C932B6B544EBBADD969D6E6F5B72FA922
+3ECEC5CE61603A0CBF60644C8E785FF198BD841B80602802BBB6535A1F2E7A5B
+7E49603065742DBBE98C625782DB5CB5D46A0667096940955B92DBCB8B1ABC4C
+93FDAD380D3ED3AAFAACCAB25D36695A84088111560E00442AABDAD53E380378
+45E906D8B1D77D3D906DFDB779E7EB82562C58F8FD66DE979157A1105953392E
+D536BAAE92E1FE4BEAA08F3888D285EFC950C008B768AC6CBF042E005A02F3CB
+FB291509248DADB0E8F014438FC70323F6719F163979D46E4F9C1AC6825C7DE5
+F28446027A7507B35A2C688407E7B0A867A6E129A1A76C94C26A2EEC166BD295
+599C46169D93A44FA5FD6A8C59E2B0ABB544CBA2C027E42603203F0792DFECFB
+A5B10489743573277616CEDB6FA54755E97D7ACEBBAB0E999046F750CDDBCF1A
+27D6942E8667E87800A559DD282787B6E50159E2175507919C70C7E63F9A18DB
+F48D85507D3E54D40FDED041729B5AA9307B3292B8D43D22BCDA1348EFEC35AF
+B1FBC6E2704446DC891B0968D965AFD0989F740D6675922092D4689F7536AA2C
+722AF4415C52036F661AB2E46D8230192C3A9768ACA7F2C4B4D70DA4C1899655
+A929754D47093E7B6486350468CC58B15261A9746F355D05A0B70C767D48558D
+FA1A315D66C900F554B9F6F65E471E8BEE6CFD754911DAE67C31D5985020155A
+2974344A459903712FE7569CCB5E98B990E2FB9A54D9A4250FD647E876A22D64
+A4B9ACBA7F0DFB147CA6511CAE5F140A239E2C853BDD89EE35B4A70C17F70DD4
+4BA59681BABF0A0F5109CDAB4AF885A75F3AA82C271370974758E87D183BE87E
+40EC7511C33AC6175B1CE915692017227EBC355165AB9A565F27E9FED1942332
+654748D05EED73A71E78D0DF232A4932E2684EDC0C57D45A6C52E2F06F02879B
+7FA12BF33A0F12238E0B5642B544C64A1FBC9A45DC6406D66FA68AEBE2A52F8F
+CD5D5D4EDB96F6F2FD09A41930F3D55FF1142A12D7EC21905DDE29D14D21F34B
+2705D67E7570E9A45B532B6BB01FDB44761CA1E425FDE93C745509FA3660376A
+653DF00F66B7C0F5277E1C84A4E7485D57E6BB19AFFD114A23AB610A9D8F839B
+92C15F70FE3EDF0AD8D8BF4EF80F10CF799EB263898CF78F45857ABD2EEC5E80
+C38A941E3031B3EF6A13A52847460E00572F80D9A3CA9879B27ED7374D8116B0
+F1D613DB684184606C943B2979188E73295236F9576586F87AD88B49029497B0
+F4DB22869C8E74E1F10EE39497032E2E5DCF21BC3AE097FF46F62D078F9AC5F1
+0B00F189E3AEE9FA8679E3384F2FBB3C94853E2EFFBF6F7998A1F170464595B4
+26EF9336F1A5FD9585CF0CA864B1066CA3328B512FC46A558E2DBD7CDB712AEF
+B4094D75489D04BB86042C4362B0AC33ABDE75726F2E6869CF4166D5289CEAB0
+1F3186EB47784D60019DA83A130574B5DD0FBB13EB94E4FC8F1FB16A797D971A
+50CDA76D62DAF52F45A76FA61623E44811A3F60011989431A1EC79D8ACE329E5
+4FA9A674D02467963C067A01AD1F4574FB3F17DC73A1DDBBA3BC40AA85C94F1E
+EEC27E46C5C622A54ABD5A4B28ACE0B0A370E164EB9E47FCD995561271C33631
+D7492FC20B7DF3EB97243BF78772E30BE766FF0B0BF89D88F1AA47CBF791B1B6
+7FD49C9D7E0D54C57F0B5506F028E6B565DB3D49ACED992047131D19B8BE0BB0
+D1163DAD0F75300C7EBEB20A11594904A96E1A32B93ED079BFF5DE5176ED291C
+660A1FEB86D1BD71AD8007351A246AA66679DA7721195389454EF48F376A577F
+2AA0EF25719652CA40F422856051E0DC758F98430C95E44DEE661505928024A7
+554D77EAE5B816A15110D7683369DBC01AC5521EBE2F68E248FCBDC0B1D57FAC
+2751FE18FD6D3B4C10D46FD0E937F46D96FB427744F75C1245A2B3E841361C78
+A092901C1FD123D389F964462DA4D57CFE9BA9C0E3D0131490ADCA047471B53E
+7C7F995CF04C2EE072569A327CC57FF9DEA7EB5D3CF61491ABB30645B57A932A
+290408D0B8420A9ED5FDCC2B737C6A3954262DE9273EBBDE0ECF242B213A26C9
+80C4A75F69458E6CA85D8270DEF8C2B31BC4A719E396F4686FBDC0CF1F4559CB
+B665E0E70635B3BE649E97A13B7F6EA1C24A0AF12A22AD23B94BEE23FB7171B2
+815438004C2A46F9EBE537B9BDDD28DC296671D7670F90DF443C6C2FAD8ED3A9
+9BEABF360C5DB0E7AB2446D000D80B8CC23501075785DAEEDBE4110A3CDCBF92
+893C55A0F1B055E50A7AE41711B96BBEC1AB4284EC5A5EF8624B751D73BABA75
+8037B962732C51A3C9EDB98E76DC53A86FE013AAE661BAF841F89E3D746CF0E9
+6CC00F1CB7F788EBA06F5531CD640E020244BFAAEAB9D6DD8C61B58A3B22DFDC
+C69E8DE14EB456ECD543C3F83F99099BB3C98660694F9180013268FE523038E5
+B2E42B916FEE5C92F2A2A0C60BB9698F08DA41265B1E5ABDA07AFDA1300CB09F
+6CA0A48DE9713A4B33B458454E7CC2412A360E436EA7F64FB2D7E4AAB8792DD8
+418AF92560CDB5C3BFE602EAB66C419D752C3F075F079C47BC940E7D54311671
+6E61C78A51D4C4B4FF6A610FADC205CE21427BA853CD3C86877B11442A9A7043
+5325EF9AD60B3D8FC73258280E8CB1B7B8DB5A9236BAF20DF2792FD440192730
+303B771B88639F74E552F3EC79BDF5F41A631BD6AA2C471099E18D76316C230B
+8F0BDE94A8224EE6342AA2F9E1DD5F62A26FA371D0C37726990F6DB6FBA02684
+F9E5E47C9AB584837358EBA22B0175E261857C86956B640AB9E0708B17ACBB73
+46F39CABB5E760211090C075E30422545CFE7AA4E8F5F51042377B0F62998F2D
+C1A6D273684D8449B6EF23F3238B46F38DC35E4D327225038F24A74ECB98E303
+AC61A5FD6A14B07B65009FD1FE0A3AD09B4F7843924C9E6B4EAB1638840F2988
+030FACB508EEC2060AD5A21589D970ED04AC60E3C5A9087572AD4DCCAA2DF6C2
+26C8383B593F691F2FC31D9F772EC08F6651C103B8DD65519727C575C54B5CC7
+7106AF7881B2AF1AFED1E190120175EE4D03F8EDC694CF2C1B287F6DC922C2D7
+C2BF29A7363132E94D3610FCEF6E6052CF02D9B4820E65EE7139FED9E1B12838
+60271147FB4E111C4ADE2686926A064C4ACE571677E8DDBAA07C49EB19DFA136
+3253AAD0E07AD2CEE04AAE2F9C2B7D457A3299B7C08A993952F12780C5E2B205
+8A5FFDC75C350797AA706EBFB4FB03FC29AED5B0295749647FFF45E972F57D34
+79DB6B3F9CBA848048E612F9F74DB6556D83B00A9C8ED38F40B0551B4149A313
+A6885250DABFBACE66A7AC1202274B9BD1C1EFB759DB251C5DEC15E39427BB74
+ADE7425816FC33EFEDBE1204AA4AADF92AB2B960692A45434EE6A2581E80B50A
+2704BDB69360CDAB86EF25FBC01DDD3750BBBC2005A1499783947A1900A46F61
+100E95E98CB8AAB39CA83FB7CEC5C732976375F90188A0E076D0B9EDC0145C55
+D59A83B8AD0546D72A638008B3738D5F0030F7D71A5B003E8E8B50BF34EF23E0
+DA2EA2077EBC110E8D4F210129716F2D01B189826C71BA9ADA11D8537C72C41B
+72F4B210AC1C0E860334374C11D16580D0B2096C9D8C072EB83178D4C9C21EF4
+68A32F8DCAA06E09C1DFEC37B17FBB34F45F79E8AB62EA123E0E46D276E0A5C9
+9EFDA220DC6EB5DADCB28F1AEBE8FECFE14E21DF91469E501775D2D219BAE384
+551BDD29026EB92165790C92417292BCA4D730EA5622C4607E87C37661389D95
+9AD454B4F90940B53B953A62BB6DA6F464F86EAD550065F6B70113499B0C6742
+521F7D2C3E238B8CD51DB109280CCB5F8E41CC4B8FDBE3D86A448659B5EB4B52
+FCE83FDD5080F29087EBF3F89108FFF1CA198591D41E4BFD77579D1BA2820909
+CFEEAC2ABEDEC2BB4F7ABF3563B54D4DCA3BDA55B242F8BEACBCEF43BE3871F8
+B6BB23139940FF3C6B0A9467DFA051C635CBC525B5F9CA95B3DFA360EBF330D3
+6E6B67F1FE652B08C793EFCB34E7CC85CB335402499A96EAE18BF3BC8FA15BB5
+CEBC40AD7D40E982F6315B5ACBCF8742A5228C7DF893D1A3DFCD728EDF84450B
+3A542B4560EB393D5A494F1EC01B6E5625973CB579E753EE7BD211C7E4AC0233
+74C91B3893C16391FCB3746BA4DED80FE9101D38A2D16A5449194E17C65859D3
+5829055C9597815DB9F18661EAD6179CCBA48426063EF7D74E2624CEA9C15347
+36125D962DCBCAD8342B751D4BFBCE023743B2464CA7452D8C412642B3B77153
+15B0651E51E1C0871BAEE1596E79746C31E6A5214126AD1EF36967B403D55AC3
+9F1A832590ED77867691C40EDF5798CA30AC596F9BB9F7CEB00ACBBAD4721C93
+9081EE528DC797806B8DFB063CE70B4D213BA105B76016D3CF1AA90572881823
+E9ADD4DE33CD09DAF42DF4114F919229AA7D9AC84246F585FD485254C01A57A0
+EAE48BC04501B50A57E2C384078FF934F38164176072F7296F04AD51B2D2ECEF
+2E6AA4D2933215A7E98B21107D368E3F8C1272FF4014BA7AB385C69A816FFD25
+326FC1E1507C0B6EF5ADCF378186E36AE252C069CAA9439093F215F6CEFAFE59
+7FA835F2BD581713B86AE39216B32E93A69C1D73561B357FFAAD0FD1D2A89EF5
+03CF5B46068143B69D5A2A013E1A9C6C6F803F532D41DCE42BFD2B920D9FAF4B
+370CE0EE5E6F0BA284EA154215E4CC7F449EA606332AD5CAE782622142BE69F5
+6461065B6B1AFC2E3C1BDDFF1F23E29C695886E3C37A2C271DC495DF5517BA83
+AC3D99CA3DF068409C340A80B029FE2581ADCC6241FC9980B448799D5FCFA469
+061BFE6B74275C7D3E329A4BB9F2F5CE15D99C479A855DA8B4297B75E75C6018
+BC5F19D64511D727BDAFBEE079F1B2E54BD0978E81ACE2DA389A6BF4C308D984
+83489D05A06296E6756B4748BB8652ABA4D953EF892FED112FBC521AFD44B4B5
+F0E9B225BEC3EC80B520A5424278E1F52FE3F297974CE4A112FDE1D62DEA6597
+D3AFF8E059C237B0AF2500740FBB0ED7C03B11F481731A7CB1960E83E4B50EF6
+0844CD3A747C9A0C48321C83F67B208D784DA0AB0342644ACCC97F438A84F279
+0D0A56A371A2A9D65E25527A7C11CDB4C0BA6A7840CB01520170E9C110EC8863
+86F46332F057466DEF7417B2AE1CD0F7F0CA2ABEA50C26E2F5323ED7C4CAD040
+274C87EB3BFDEAD0EB4097A098F5ECD85A52C54410D71680BDA91BF39E1E514A
+7C7FA83969A39E2F9AD38C105C59FCAC4663B964F6674077E0B2971C2CC68726
+D2982547B187880C64C77C9ED778AFFC9FF20BC2357EDCE553E72F455CCD76A9
+7525FE42177BB64A93CDC312EE96886222CA5DA55CD1FDE165B7BBDC01207B42
+4F5A804B106C3F342FE9D15F45972C49513720520852D071B6D48ED8118C3602
+843CF86C64284DFC2966898C2A36D2D0938D02FF899B6380CAC5670A9548C2B1
+40B193458E778E0DBD52F58A03BA0618C8299AC36C9636ABADD02C943F6A0264
+E7E9B383A50E125157638F78245B1B0FE00F541BB9AAE6731AA57A3FBABB81F4
+6BA0A25ABCE226200CCEE74ABDBD90B3F649DBE616D9A13B63BAA453DF3B5FE6
+028F999A2B1B9115CACC40991A3A0E22663E3A4CDD45C009CBD50588623CB855
+0556A619E6BB4F58B9B5BA347E4EF6FB6B63259BB4F13219E0BA243F980BA167
+C6875F60AE769471D866B5F7E4EF628761278399213F99CFB3D507E9CDC764E5
+9739A74E83FDE8C9198F51626676F06EDB233E535ACF2FC4E893F748E61C8246
+4D59959D9CAF25AFAA2DA85073092FE312BEBFE2D915D1E3BE8F5AFFD2EDDF29
+9621B62DB047F75C76625136B1FB986F68FD835EE6C1D1DFBE2798BA3A766336
+45424556C0CF0B0EFC6F7CC2386FF195B73E0A8AFF7060083DB424373F645D5E
+901A9D4BDFA7EF4AA96B19002E8E88CF34E6010FC6726A482827D20BADBC156D
+1450E46B7858BB3AE1878DDF06A0022D2D5FC59474637201411855596DE5A37E
+38C2BCCD94019AF30D7BCD936DD128207AC47316CAFB76A9E48377CF61CAAEA3
+CBE8BAFAD4A4CB9B19316CE5991E9887F838299AEEBE85E25945C94AFAA5548D
+FA3D0131256E59D78CC0A630308637522C4AB984E278B37B5F9A973C1924D64A
+BD94EC12F23EC16D9A6F17D781BDE886D2D5A927B7EF1CD25251692A576D4E83
+25930790E46CC506DB373C0A4A54D577F35868A3A4F74BCA45C9DC952436F7B7
+E4906EE55C86EDD230D798A5E0313B8CF43C7E637F5A95BEF3DD385FEDCF8343
+B08753F02AB9DF4A7050263B3931FAD938836E8F04A941E296DD6407B5059757
+45E6DBA31AC57F790F35CB8E57DB3381EC4909515BAB8FD6F01C5E50E4AE641F
+79826BACD9ABCCDF01D5D737828638B606D4E1EDF9289DC3A75E20398F325901
+EFBE09187507DF6C754F93DB49867757135735A88AC460004B9CAC0B8A308BDD
+48DAD3B81618D357E5BEDFA16A7E5812FA3A8448B87AFBB9D4AD1A793DE94DDA
+0541D75D494DCA4C6933F2D884B96D9EF031AB6F2B4220DD0ABA23C3B2D9CFA9
+63B68E3D6CFE830B852F0CAEA9EC709B67A7D28A601A10ADC0855046F3928182
+3E7B59410EF943A8F8DFE78CF1985D38E773F0E95682140BACBC44B9AA541378
+DB8B6AE642FD940890C079A8544166B0267F81A6B50FFE65400EAAE45D07B962
+0BB63680EEFCF2E028D48DAE5C251F9FFC9011659680F3AEF4419944A6FAB49C
+06BC461F92584173433FA0263B6B4A5C1FAF29635CC9CED3C85D2C687009E776
+3A8B2231C3089E59AE6A716A3ED883FDB149B7465B0240E44BF87A8A3700DA04
+E404C2DD91FE885E6242CA076F9A6A1C2221A703122D7297D0BBBB51C4FCE753
+00369D0904530727FCE8B184AE2DA1A39547A712A4F80B1CA5D68253AE731BD2
+BB2146441AC8D7FCFCE5626F7179ABC6829AC6BCD48B252BD57CD26A66367E55
+10B9BB604DB86506F9F92E7F70AC5CED92D0321EC1DA1BA102EE7999830CD109
+BEF289B44B940A6010B4121A161B5493E37C94A6359C6AE9823A06D3BF1A3349
+9AB1A5EA2CE1023DDB17838264047E518C6FBDDB43D72856FB03925250C79173
+EF9CD0D6F9D0ADAA52A425E0205EC26AF1E423F3402DE6ABFEBEF332838CFA0E
+1BBA5DB818F8602BC2DC13166F2FFEFD5416FBA6E1E5C1CB78BB7CD6D813EB91
+3D4F60B72B0BC7B0D511633EB7D0C9961B2729E08844ACD961F6D705567FEBD2
+46D604CBD98C898559B10C1F2E3EB80571EA48A6FC05EFE3030391BC9913361A
+13F4ED8862F26B2BB762AE060932E17531A1BB1558FCC40CFCE2D6AAAB3F636C
+151D5A29280ECCEE0A6C5B5FBE8FAAEDFB7B4C6BCD492E755F5A59575D3F5A36
+A843CAE18D4307ECBE2D85763396782FE4B393F44D6E7623463EA56D8AB7D4D2
+0CA51D1ACD9E58D6ACC60E85A0D3264A1C48BA4E54207CD3485F853232D91396
+5502E44E4A443A4410789E25DFB5862E0DFB2CD47A9C137034992AA7ADA4BB60
+C31CCAB6C226E97779CF2F3DE2A4D523A6228998CB4A930594D171E96C6D4D3D
+B56197B941DF5B677E9A1BBDAFD52BADF156315B39F8911A0B1A57885ABC720B
+179CE1FADC1024085D96745D64678EA9B2FEC47EA2A6E1A5F714581D2888DC75
+06AE25730E75E8837C2AD1A4E64A9B7B61FC22F5663429AC87895880F6D83C95
+1C63FB76623C2340367B38F7B3DD52CD698239B73AD5F0585F18002A39AD3D1E
+161BAFBF6895D34720649D237BDAC66D78CEFAE151926F520D78FF8FDAB248E4
+2CC4F426E03D6C35F2BE2E46E05AC5320301D7B1DDBFEC26DFAE95108258E514
+15E481D96EE15BB4FE9FC1E0D3828DC9355F383A69D4A4E5DF6469416C896494
+3E683B442CA498BD87B5C24B120EE9C612072EE84EE3957074C349E1E1042104
+F2A59307CB8A25EBA004FAE2D1CD6837782B089D69149CAC2A73C2CEA173AA5C
+14F2BDB160D3E320069C55F7613AB1AA162E60C4F429990FFD2765A864DF60B0
+B12B1A06EBF3D7AE83B6C4640975AF45C31705B73641F4938C902904560C86D3
+9CDC690E69DE8896150B4946A0940B87F722838C993087F9D72D09349C506E97
+62615FBE54A4DA129FCD621FC384A16040A4F438648C6891A926C552515EF14E
+FF01BE86DFA631B93489ED4967A6EC9E5E51DBC27481DE45F6AC92C93AEAA566
+1801D40B16B08E7924F539DB1F271FBA5A553B69D1F39564CD4D5341FDE31386
+83E8BABF0E91B9113F140B30D0ABB9D6CFCBAD7D36E302B07728D38783EF13B1
+A5F2514CD568D6DABA691332EF7B2DC1FEE6D2B280A77F721B245288E542FBAF
+0FA83D8DDD47DA2B7ECC2B820AF1C82A8977442112B79458D0DB791AD787F1D9
+6D144BF530AB44A10D8C2481B778D72C70B32A5F4688916D46770F606FDCBA5E
+083E3C699AC9682FCA3C451020846B804EE36C3029E975647814938537722BD9
+9B0136E68804E3BC9BCA06FCC59B4161887864D06C078C42F5872D9452A8D3A0
+0AC19FBA80B1BAD7422FD74973A34DFC78BEF63F08480E515A4FD51CE7E24AB2
+494FE559C0C3AE5553ECDC24245782C66A3DD3629EBF1CD3950BF0C1245BB6DA
+A52BA2518DE38E406850CB642ED26A0321255EF6FC30A98136B8F0760235BE86
+981FD80BDA2CEC800B0C08E0D546CD38398D9F936368E4FFD0F7E2E8B9D68A80
+34ADB6F8AB35DF361170B74E5CF963C96BCC54BE8A7D35C5AA45A9720A66C7F0
+3963B6802141BB300AE9A14A64DB4FDFBC26BF56223F140A2DE200FB8F1DA206
+148C247F44404F50DCBCE9B0444A91FCE717B128ACE82F2AE881C15AD14B15C9
+1A30AE2924E36ECAAC8690DD2C438EE3FC24204C51C61D0E17A3090EBA6E4157
+A6CA6683DC25D965E86E013F1E1855C7D7992B92882C7B781410AE4ECA7CDDE0
+B2DF0912E55AAFA01D0615371FEDD1B60AB5736E1330E5F660BC9A31855C38B4
+DEB3ECD7CEBBD78D4DBC7D6FC6E2FEF2900759F4427CC247C6134FE4504E9CD3
+F4D060F266203911036342E308E6235109995F2551F71A42F5D89C82D0908BFB
+E2B984BA35AC3E8412CB90DB173BCF035A50C24D1CE0264DD061BD6E4B1130B4
+A685F53E839714BC6051110C82FB71960EAE88E3AD6EC655E84A8B5DE1B271F2
+C9F5D76F46AB74BFB17885CA16E7A0326B42A605D206797B45FE513314C1ED3A
+BA91974D03A445BFEC3820144C1A2EA1344A8A28F3F02D652B48079FBBAC34F8
+439B0170901BAAA186042B7D079292C98EE7F7AA53B5F60E4EE0C57C1A93E559
+3DC6BE4D186C6ADAC4048BA3834E685308E229CA60DAD94B8E15D7E1C7C2A840
+CE7D85AA7E38D3FBB6D703702AD5F52098C6AD16E97E3775160E4A6DD7C987
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+0000000000000000000000000000000000000000000000000000000000000000
+cleartomark
+%%EndFont
+TeXDict begin 39158280 55380996 1000 600 600 (blankenbach.dvi)
+ at start /Fa 201[35 35 35 52[{ TeXf7b6d320Encoding ReEncodeFont }3
+66.4176 /CMR8 rf /Fb 207[19 48[{ TeXbbad153fEncoding ReEncodeFont }1
+66.4176 /CMSY8 rf /Fc 134[40 121[{ TeXbbad153fEncoding ReEncodeFont }1
+90.9091 /CMSY10 rf /Fd 135[40 3[25 33 32 4[62 7[33 18[53
+5[48 76[{ TeXaae443f0Encoding ReEncodeFont }8 66.4176
+/CMMI8 rf /Fe 134[48 48 66 48 51 35 36 36 1[51 45 51
+76 25 48 1[25 51 45 1[40 51 40 51 45 10[68 1[66 51 67
+3[68 83 5[71 2[69 66 64 7[25 45 45 45 45 45 45 45 45
+45 45 1[25 30 25 2[35 35 40[{ TeXf7b6d320Encoding ReEncodeFont }48
+90.9091 /CMR10 rf end
+%%EndProlog
+%%BeginSetup
+%%Feature: *Resolution 600dpi
+TeXDict begin
+%%PaperSize: A4
+ end
+%%EndSetup
+%%Page: 1 1
+TeXDict begin 1 0 bop 448 647 3635 4 v 446 760 4 113
+v 511 726 a Fe(Grid)138 b(V)900 740 y Fd(r)r(ms)1158
+726 y Fe(Nusselt)31 b(No.)145 b(T)-8 b(op)s(o)1965 740
+y Fd(L)2218 726 y Fe(T)g(op)s(o)2420 740 y Fd(R)2635
+726 y Fe(Geoid)2867 740 y Fd(L)3019 726 y Fe(Geoid)3251
+740 y Fd(R)3410 726 y Fe(Run)29 b(Time)h(\(sec\))p 4081
+760 V 448 764 3635 4 v 446 877 4 113 v 615 843 a(32)101
+b(42.964)229 b(4.892)f(2238.711)103 b(-2879.020)120 b(54.074)i(-61.715)
+477 b(218.35)p 4081 877 V 446 989 V 615 956 a(64)101
+b(42.890)229 b(4.886)f(2250.099)103 b(-2897.102)120 b(54.630)i(-62.392)
+432 b(1628.02)p 4081 989 V 446 1102 V 615 1069 a(96)101
+b(42.876)229 b(4.885)f(2252.273)103 b(-2900.498)120 b(54.736)i(-62.520)
+432 b(5526.62)p 4081 1102 V 446 1215 V 570 1181 a(128)101
+b(42.871)229 b(4.885)f(2253.039)103 b(-2901.690)120 b(54.774)i(-62.565)
+386 b(15922.29)p 4081 1215 V 448 1219 3635 4 v 446 1332
+4 113 v 498 1298 a Fc(y)p Fe(C)604 1312 y Fd(ext)806
+1298 y Fe(42.865)229 b(4.884)f(2254.021)103 b(-2903.221)120
+b(54.822)i(-62.622)p 4081 1332 V 448 1335 3635 4 v 446
+1448 4 113 v 498 1414 a Fc(y)p Fe(Christensen)1005 1381
+y Fb(0)1028 1414 y Fe(s)30 b(extrap)s(olated)h(v)-5 b(alues.)p
+4081 1448 V 448 1451 3635 4 v 448 1611 a(T)d(able)41
+b(1:)60 b(Blank)m(en)m(bac)m(h)41 b(\(1989\))h(Benc)m(hmark)e(1a:)61
+b(Steady)40 b(State,)j(2D,)e(constan)m(t)448 1724 y(viscosit)m(y)24
+b(con)m(v)m(ection)g(in)e(a)g(1)g(b)m(y)g(1)g(b)s(o)m(x)g(with)f(Ra)m
+(yleigh)i(n)m(um)m(b)s(er)e(10)2805 1691 y Fa(4)2867
+1724 y Fe(using)g(ConMan)p 448 2136 3680 4 v 446 2249
+4 113 v 511 2215 a(Grid)160 b(V)922 2229 y Fd(r)r(ms)1204
+2215 y Fe(Nusselt)30 b(No.)145 b(T)-8 b(op)s(o)2010 2229
+y Fd(L)2264 2215 y Fe(T)g(op)s(o)2466 2229 y Fd(R)2680
+2215 y Fe(Geoid)2912 2229 y Fd(L)3065 2215 y Fe(Geoid)3297
+2229 y Fd(R)3455 2215 y Fe(Run)29 b(Time)i(\(sec\))p
+4126 2249 V 448 2252 3680 4 v 446 2365 4 113 v 615 2331
+a(32)101 b(194.237)207 b(10.556)f(1430.652)103 b(-1939.314)119
+b(26.308)k(-30.355)477 b(124.46)p 4126 2365 V 446 2478
+V 615 2444 a(64)101 b(193.422)207 b(10.540)f(1453.292)103
+b(-1986.864)119 b(27.298)k(-31.543)431 b(1200.03)p 4126
+2478 V 446 2591 V 615 2557 a(96)101 b(193.303)207 b(10.539)f(1457.178)
+103 b(-1997.129)119 b(27.515)k(-31.805)431 b(5607.44)p
+4126 2591 V 446 2704 V 570 2670 a(128)101 b(193.277)207
+b(10.539)f(1458.659)103 b(-2000.774)119 b(27.596)k(-31.902)386
+b(17010.54)p 4126 2704 V 448 2707 3680 4 v 446 2820 4
+113 v 498 2786 a Fc(y)p Fe(C)604 2800 y Fd(ext)806 2786
+y Fe(193.214)207 b(10.534)f(1460.986)103 b(-2004.205)119
+b(27.703)k(-32.016)p 4126 2820 V 448 2824 3680 4 v 446
+2936 4 113 v 498 2903 a Fc(y)p Fe(Christensen)1005 2870
+y Fb(0)1028 2903 y Fe(s)30 b(extrap)s(olated)h(v)-5 b(alues.)p
+4126 2936 V 448 2940 3680 4 v 448 3099 a(T)d(able)40
+b(2:)59 b(Blank)m(en)m(bac)m(h)41 b(\(1989\))h(Benc)m(hmark)e(1b:)58
+b(Steady)40 b(State,)i(2D,)f(constan)m(t)448 3212 y(viscosit)m(y)24
+b(con)m(v)m(ection)g(in)e(a)g(1)g(b)m(y)g(1)g(b)s(o)m(x)g(with)f(Ra)m
+(yleigh)i(n)m(um)m(b)s(er)e(10)2805 3179 y Fa(5)2867
+3212 y Fe(using)g(ConMan)p 448 3625 3648 4 v 446 3738
+4 113 v 511 3704 a(Grid)160 b(V)922 3718 y Fd(r)r(ms)1204
+3704 y Fe(Nusselt)30 b(No.)123 b(T)-8 b(op)s(o)1988 3718
+y Fd(L)2218 3704 y Fe(T)g(op)s(o)2420 3718 y Fd(R)2641
+3704 y Fe(Geoid)2873 3718 y Fd(L)3032 3704 y Fe(Geoid)3264
+3718 y Fd(R)3423 3704 y Fe(Run)29 b(Time)h(\(sec\))p
+4094 3738 V 448 3741 3648 4 v 446 3854 4 113 v 615 3820
+a(32)101 b(850.318)207 b(20.428)f(818.862)102 b(-1189.116)126
+b(11.778)j(-13.033)477 b(325.17)p 4094 3854 V 446 3967
+V 615 3933 a(64)101 b(838.021)207 b(22.016)f(905.007)102
+b(-1253.356)126 b(12.664)j(-14.275)432 b(2866.17)p 4094
+3967 V 446 4080 V 615 4046 a(96)101 b(835.751)207 b(22.026)f(921.826)
+102 b(-1269.228)126 b(13.049)j(-14.659)386 b(17217.40)p
+4094 4080 V 446 4193 V 570 4159 a(128)101 b(834.961)207
+b(22.002)f(926.278)102 b(-1275.250)126 b(13.214)j(-14.814)386
+b(52139.65)p 4094 4193 V 448 4196 3648 4 v 446 4309 4
+113 v 498 4275 a Fc(y)p Fe(C)604 4289 y Fd(ext)806 4275
+y Fe(833.989)207 b(21.997)f(931.962)102 b(-1283.813)i(13.4520)i
+(-15.034)p 4094 4309 V 448 4312 3648 4 v 446 4425 4 113
+v 498 4391 a Fc(y)p Fe(Christensen)1005 4358 y Fb(0)1028
+4391 y Fe(s)30 b(extrap)s(olated)h(v)-5 b(alues.)p 4094
+4425 V 448 4428 3648 4 v 448 4588 a(T)d(able)41 b(3:)61
+b(Blank)m(en)m(bac)m(h)42 b(\(1989\))g(Benc)m(hmark)f(1c:)61
+b(Steady)41 b(State,)i(2D,)f(constan)m(t)448 4701 y(viscosit)m(y)35
+b(con)m(v)m(ection)h(in)d(a)h(1)f(b)m(y)g(1)h(b)s(o)m(x)f(with)g(Ra)m
+(yleigh)i(n)m(um)m(b)s(er)d(10)2931 4668 y Fa(6)3004
+4701 y Fe(using)h(Con-)448 4814 y(Man.)1920 5225 y(1)p
+eop end
+%%Trailer
+
+userdict /end-hook known{end-hook}if
+%%EOF
Added: mc/2D/ConMan/trunk/src/blankenbach/blankenbach.tex
===================================================================
--- mc/2D/ConMan/trunk/src/blankenbach/blankenbach.tex (rev 0)
+++ mc/2D/ConMan/trunk/src/blankenbach/blankenbach.tex 2008-06-04 22:36:13 UTC (rev 12086)
@@ -0,0 +1,75 @@
+\documentclass[11pt]{article}
+\begin{document}
+\begin{table}
+\centering
+\begin{tabular} {| r c c c c c c r |} \hline
+\multicolumn{1}{| c }{Grid} &
+\multicolumn{1}{ c }{V$_{rms}$} &
+\multicolumn{1}{ c }{Nusselt No.} &
+\multicolumn{1}{ c }{Topo$_L$} &
+\multicolumn{1}{ c }{Topo$_R$} &
+\multicolumn{1}{ c }{Geoid$_L$} &
+\multicolumn{1}{ c }{Geoid$_R$} &
+\multicolumn{1}{ c |}{Run Time (sec)} \\
+\hline
+32 & nan & 0.685 & 2238.711 & -2879.020 & 54.074 & -61.715 & 0.00 \\
+64 & nan & 0.671 & 2250.099 & -2897.102 & 54.630 & -62.392 & 0.00 \\
+96 & nan & 0.650 & 2252.273 & -2900.498 & 54.736 & -62.520 & 0.00 \\
+128 & 42.871 & 4.885 & 2253.039 & -2901.690 & 54.774 & -62.565 & 15922.29 \\
+\hline
+\dag C$_{ext}$ & 42.865 & 4.884 & 2254.021 & -2903.221 & 54.822 & -62.622 & \null \\
+\hline
+\multicolumn{8}{|l|}{\dag Christensen$^{\prime}$s extrapolated values.} \\
+\hline
+\end{tabular}
+\caption{Blankenbach (1989) Benchmark 1a: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^4$ using ConMan}
+\end{table}
+\begin{table}
+\centering
+\begin{tabular} {| r c c c c c c r |} \hline
+\multicolumn{1}{| c }{Grid} &
+\multicolumn{1}{ c }{V$_{rms}$} &
+\multicolumn{1}{ c }{Nusselt No.} &
+\multicolumn{1}{ c }{Topo$_L$} &
+\multicolumn{1}{ c }{Topo$_R$} &
+\multicolumn{1}{ c }{Geoid$_L$} &
+\multicolumn{1}{ c }{Geoid$_R$} &
+\multicolumn{1}{ c |}{Run Time (sec)} \\
+\hline
+32 & nan & 0.137 & 1430.652 & -1939.314 & 26.308 & -30.355 & 0.00 \\
+64 & nan & 0.135 & 1453.292 & -1986.864 & 27.298 & -31.543 & 0.00 \\
+96 & nan & 0.178 & 1457.178 & -1997.129 & 27.515 & -31.805 & 0.00 \\
+128 & 193.277 & 10.539 & 1458.659 & -2000.774 & 27.596 & -31.902 & 17010.54 \\
+\hline
+\dag C$_{ext}$ & 193.214 & 10.534 & 1460.986 & -2004.205 & 27.703 & -32.016 & \null \\
+\hline
+\multicolumn{8}{|l|}{\dag Christensen$^{\prime}$s extrapolated values.} \\
+\hline
+\end{tabular}
+\caption{Blankenbach (1989) Benchmark 1b: Steady State, 2D, constant viscosity convection in a 1 by 1 box with Rayleigh number $10^5$ using ConMan}
+\end{table}
+\begin{table}
+\centering
+\begin{tabular} {| r c c c c c c r |} \hline
+\multicolumn{1}{| c }{Grid} &
+\multicolumn{1}{ c }{V$_{rms}$} &
+\multicolumn{1}{ c }{Nusselt No.} &
+\multicolumn{1}{ c }{Topo$_L$} &
+\multicolumn{1}{ c }{Topo$_R$} &
+\multicolumn{1}{ c }{Geoid$_L$} &
+\multicolumn{1}{ c }{Geoid$_R$} &
+\multicolumn{1}{ c |}{Run Time (sec)} \\
+\hline
+32 & 3.110 & 0.093 & 818.862 & -1189.116 & 11.778 & -13.033 & 0.00 \\
+64 & 1.577 & 0.088 & 905.007 & -1253.356 & 12.664 & -14.275 & 0.00 \\
+96 & 1.055 & 0.145 & 921.826 & -1269.228 & 13.049 & -14.659 & 0.00 \\
+128 & 834.961 & 22.002 & 926.278 & -1275.250 & 13.214 & -14.814 & 52139.65 \\
+\hline
+\dag C$_{ext}$ & 833.989 & 21.997 & 931.962 & -1283.813 & 13.4520& -15.034&