figii6.eps

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

    /stroke { oldstroke } bind def
    /show   { oldshow } bind def
    /ashow  { oldashow } bind def
    /widthshow { oldwidthshow } bind def
    /awidthshow { oldawidthshow } bind def
    /kshow  { oldkshow } bind def
    } ifelse
  end
} bind def
false PATredef
% Conditionally define setcmykcolor if not available
/setcmykcolor where { pop } {
  /setcmykcolor {
    1 sub 4 1 roll
    3 {
	3 index add neg dup 0 lt { pop 0 } if 3 1 roll
    } repeat
    setrgbcolor - pop
  } bind def
} ifelse
/PATsc {		% colorarray
  aload length		% c1 ... cn length
    dup 1 eq { pop setgray } { 3 eq { setrgbcolor } { setcmykcolor
  } ifelse } ifelse
} bind def
/PATsg {		% dict
  begin
    lw setlinewidth
    lc setlinecap
    lj setlinejoin
    ml setmiterlimit
    ds aload pop setdash
    cc aload pop setrgbcolor
    cm setmatrix
  end
} bind def

/PATDict 3 dict def
/PATsp {
  true PATredef
  PATDict begin
    /CurrentPattern exch def
    % If it's an uncolored pattern, save the color
    CurrentPattern /PaintType get 2 eq {
      /PColor exch def
    } if
    /CColor [ currentrgbcolor ] def
  end
} bind def
% PATstroke - stroke with the current pattern
/PATstroke {
  countdictstack
  save
  mark
  {
    currentpoint strokepath moveto
    PATpcalc				% proc nw nh px py
    clip newpath PATfill
    } stopped {
	(*** PATstroke Warning: Path is too complex, stroking
	  with gray) =
    cleartomark
    restore
    countdictstack exch sub dup 0 gt
	{ { end } repeat } { pop } ifelse
    gsave 0.5 setgray oldstroke grestore
  } { pop restore pop } ifelse
  newpath
} bind def
/PATtcalc {		% modmtx tilingtype PATtcalc tilematrix
  % Note: tiling types 2 and 3 are not supported
  gsave
    exch concat					% tilingtype
    matrix currentmatrix exch			% cmtx tilingtype
    % Tiling type 1 and 3: constant spacing
    2 ne {
	% Distort the pattern so that it occupies
	% an integral number of device pixels
	dup 4 get exch dup 5 get exch		% tx ty cmtx
	XStep 0 dtransform
	round exch round exch			% tx ty cmtx dx.x dx.y
	XStep div exch XStep div exch		% tx ty cmtx a b
	0 YStep dtransform
	round exch round exch			% tx ty cmtx a b dy.x dy.y
	YStep div exch YStep div exch		% tx ty cmtx a b c d
	7 -3 roll astore			% { a b c d tx ty }
    } if
  grestore
} bind def
/PATusp {
  false PATredef
  PATDict begin
    CColor PATsc
  end
} bind def

% right30
11 dict begin
/PaintType 1 def
/PatternType 1 def
/TilingType 1 def
/BBox [0 0 1 1] def
/XStep 1 def
/YStep 1 def
/PatWidth 1 def
/PatHeight 1 def
/Multi 2 def
/PaintData [
  { clippath } bind
  { 32 16 true [ 32 0 0 -16 0 16 ]
	{<00030003000c000c0030003000c000c0030003000c000c00
	30003000c000c00000030003000c000c0030003000c000c0
	030003000c000c0030003000c000c000>}
     imagemask } bind
] def
/PaintProc {
	pop
	exec fill
} def
currentdict
end
/P2 exch def
1.1111 1.1111 scale %restore scale

/cp {closepath} bind def
/ef {eofill} bind def
/gr {grestore} bind def
/gs {gsave} bind def
/sa {save} bind def
/rs {restore} bind def
/l {lineto} bind def
/m {moveto} bind def
/rm {rmoveto} bind def
/n {newpath} bind def
/s {stroke} bind def
/sh {show} bind def
/slc {setlinecap} bind def
/slj {setlinejoin} bind def
/slw {setlinewidth} bind def
/srgb {setrgbcolor} bind def
/rot {rotate} bind def
/sc {scale} bind def
/sd {setdash} bind def
/ff {findfont} bind def
/sf {setfont} bind def
/scf {scalefont} bind def
/sw {stringwidth} bind def
/tr {translate} bind def
/tnt {dup dup currentrgbcolor
  4 -2 roll dup 1 exch sub 3 -1 roll mul add
  4 -2 roll dup 1 exch sub 3 -1 roll mul add
  4 -2 roll dup 1 exch sub 3 -1 roll mul add srgb}
  bind def
/shd {dup dup currentrgbcolor 4 -2 roll mul 4 -2 roll mul
  4 -2 roll mul srgb} bind def
/$F2psBegin {$F2psDict begin /$F2psEnteredState save def} def
/$F2psEnd {$F2psEnteredState restore end} def
%%EndProlog

$F2psBegin
10 setmiterlimit
n 0 842 m 0 0 l 595 0 l 595 842 l cp clip
 0.06000 0.06000 sc
7.500 slw
% Polyline
gs  clippath
7165 4750 m 7365 4800 l 7165 4850 l 7515 4850 l 7515 4750 l  cp clip
n 1200 4800 m 7500 4800 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 80.00 320.00] PATmp PATsp ef gr PATusp gs col-1 s gr gr

% arrowhead
60.000 slw
n 7165 4750 m 7365 4800 l 7165 4850 l  col-1 s
7.500 slw
% Polyline
gs  clippath
3973 551 m 4125 412 l 4056 606 l 4250 315 l 4167 260 l  cp clip
n 1200 4800 m 4200 300 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 80.00 20.00] PATmp PATsp ef gr PATusp gs col-1 s gr gr

% arrowhead
60.000 slw
n 3973 551 m 4125 412 l 4056 606 l  col-1 s
45.000 slw
% Polyline
n 3810 3615 m 4710 3615 l 6015 3000 l 6405 2400 l 5505 2400 l 4200 3000 l
 3810 3600 l cp gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 254.00 160.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4200 3000 m 6000 3000 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 280.00 200.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 5535 2400 m 4710 3600 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 314.00 160.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4815 630 m 5490 2400 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 321.00 42.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4830 660 m 4230 2985 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 282.00 44.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4845 600 m 3825 3615 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 255.00 40.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4860 660 m 4725 3615 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 315.00 44.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4860 675 m 6000 3000 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 324.00 45.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4845 660 m 6390 2385 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 323.00 44.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
7.500 slw
% Polyline
n 6600 4800 m 9015 1215 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 440.00 81.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 5685 4800 m 8115 1200 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 379.00 80.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 3885 4785 m 6315 1200 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 259.00 80.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 3000 4800 m 5400 1200 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 200.00 80.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 2100 4800 m 4500 1200 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 140.00 80.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 4800 4800 m 7200 1200 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 320.00 80.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 2400 3000 m 7800 3000 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 160.00 200.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 1545 4200 m 6945 4200 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 103.00 280.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 1965 3600 m 7365 3600 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 131.00 240.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 2775 2400 m 8175 2400 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 185.00 160.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
n 3180 1800 m 8580 1800 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 212.00 120.00] PATmp PATsp ef gr PATusp gs col-1 s gr 
% Polyline
gs  clippath
5740 924 m 5650 1010 l 5689 892 l 5547 1122 l 5598 1154 l  cp clip
n 6090 300 m 5580 1125 l gs /PC [[1.00 1.00 1.00] [0.00 0.00 0.00]] def
15.00 15.00 sc P2 [16 0 0 -8 372.00 20.00] PATmp PATsp ef gr PATusp gs col-1 s gr gr

% arrowhead
60.000 slw
n 5740 924 m 5650 1010 l 5689 892 l  col-1 s
/Times-Italic ff 270.00 scf sf
7455 5145 m
gs 1 -1 sc (x) col-1 sh gr
/Times-Italic ff 270.00 scf sf
3795 450 m
gs 1 -1 sc (y) col-1 sh gr
/Times-Italic ff 270.00 scf sf
6165 255 m
gs 1 -1 sc (w) col-1 sh gr
/Times-Italic ff 210.00 scf sf
6420 165 m
gs 1 -1 sc (i) col-1 sh gr
$F2psEnd
rs
end