diff --git a/ammoni/figure_ammoni/E-distance-nh3-w-bsse.eps b/ammoni/figure_ammoni/E-distance-nh3-w-bsse.eps
new file mode 100644
index 0000000..73407b6
--- /dev/null
+++ b/ammoni/figure_ammoni/E-distance-nh3-w-bsse.eps
@@ -0,0 +1,1460 @@
+%!PS-Adobe-2.0 EPSF-2.0
+%%Title: E-distance-NH3-w-bsse.eps
+%%Creator: gnuplot 5.2 patchlevel 2
+%%CreationDate: Fri Oct 18 16:12:03 2019
+%%DocumentFonts: (atend)
+%%BoundingBox: 50 50 900 900
+%%EndComments
+%%BeginProlog
+/gnudict 256 dict def
+gnudict begin
+%
+% The following true/false flags may be edited by hand if desired.
+% The unit line width and grayscale image gamma correction may also be changed.
+%
+/Color true def
+/Blacktext false def
+/Solid false def
+/Dashlength 1 def
+/Landscape false def
+/Level1 false def
+/Level3 false def
+/Rounded false def
+/ClipToBoundingBox false def
+/SuppressPDFMark false def
+/TransparentPatterns false def
+/gnulinewidth 15.000 def
+/userlinewidth gnulinewidth def
+/Gamma 1.0 def
+/BackgroundColor {-1.000 -1.000 -1.000} def
+%
+/vshift -200 def
+/dl1 {
+  10.0 Dashlength userlinewidth gnulinewidth div mul mul mul
+  Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
+} def
+/dl2 {
+  10.0 Dashlength userlinewidth gnulinewidth div mul mul mul
+  Rounded { currentlinewidth 0.75 mul add } if
+} def
+/hpt_ 31.5 def
+/vpt_ 31.5 def
+/hpt hpt_ def
+/vpt vpt_ def
+/doclip {
+  ClipToBoundingBox {
+    newpath 50 50 moveto 900 50 lineto 900 900 lineto 50 900 lineto closepath
+    clip
+  } if
+} def
+%
+% Gnuplot Prolog Version 5.1 (Oct 2015)
+%
+%/SuppressPDFMark true def
+%
+/M {moveto} bind def
+/L {lineto} bind def
+/R {rmoveto} bind def
+/V {rlineto} bind def
+/N {newpath moveto} bind def
+/Z {closepath} bind def
+/C {setrgbcolor} bind def
+/f {rlineto fill} bind def
+/g {setgray} bind def
+/Gshow {show} def   % May be redefined later in the file to support UTF-8
+/vpt2 vpt 2 mul def
+/hpt2 hpt 2 mul def
+/Lshow {currentpoint stroke M 0 vshift R 
+	Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def
+/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
+	Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def
+/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R 
+	Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def
+/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
+  /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
+/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
+ {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
+/BL {stroke userlinewidth 2 mul setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+/AL {stroke userlinewidth 2 div setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+/UL {dup gnulinewidth mul /userlinewidth exch def
+	dup 1 lt {pop 1} if 10 mul /udl exch def} def
+/PL {stroke userlinewidth setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+3.8 setmiterlimit
+% Classic Line colors (version 5.0)
+/LCw {1 1 1} def
+/LCb {0 0 0} def
+/LCa {0 0 0} def
+/LC0 {1 0 0} def
+/LC1 {0 1 0} def
+/LC2 {0 0 1} def
+/LC3 {1 0 1} def
+/LC4 {0 1 1} def
+/LC5 {1 1 0} def
+/LC6 {0 0 0} def
+/LC7 {1 0.3 0} def
+/LC8 {0.5 0.5 0.5} def
+% Default dash patterns (version 5.0)
+/LTB {BL [] LCb DL} def
+/LTw {PL [] 1 setgray} def
+/LTb {PL [] LCb DL} def
+/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
+/LT0 {PL [] LC0 DL} def
+/LT1 {PL [2 dl1 3 dl2] LC1 DL} def
+/LT2 {PL [1 dl1 1.5 dl2] LC2 DL} def
+/LT3 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC3 DL} def
+/LT4 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
+/LT5 {PL [4 dl1 2 dl2] LC5 DL} def
+/LT6 {PL [1.5 dl1 1.5 dl2 1.5 dl1 1.5 dl2 1.5 dl1 6 dl2] LC6 DL} def
+/LT7 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC7 DL} def
+/LT8 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC8 DL} def
+/SL {[] 0 setdash} def
+/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
+/Dia {stroke [] 0 setdash 2 copy vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath stroke
+  Pnt} def
+/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
+  currentpoint stroke M
+  hpt neg vpt neg R hpt2 0 V stroke
+ } def
+/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath stroke
+  Pnt} def
+/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
+  hpt2 vpt2 neg V currentpoint stroke M
+  hpt2 neg 0 R hpt2 vpt2 V stroke} def
+/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath stroke
+  Pnt} def
+/Star {2 copy Pls Crs} def
+/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath fill} def
+/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath fill} def
+/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath stroke
+  Pnt} def
+/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath fill} def
+/DiaF {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath fill} def
+/Pent {stroke [] 0 setdash 2 copy gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath stroke grestore Pnt} def
+/PentF {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath fill grestore} def
+/Circle {stroke [] 0 setdash 2 copy
+  hpt 0 360 arc stroke Pnt} def
+/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
+/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
+/C1 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C2 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C3 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C4 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 180 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C5 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc
+	2 copy moveto
+	2 copy vpt 180 270 arc closepath fill
+	vpt 0 360 arc} bind def
+/C6 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C7 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C8 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 270 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C9 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 270 450 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
+	2 copy moveto
+	2 copy vpt 90 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C11 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 180 arc closepath fill
+	2 copy moveto
+	2 copy vpt 270 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C12 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 180 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C13 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc closepath fill
+	2 copy moveto
+	2 copy vpt 180 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C14 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 360 arc closepath fill
+	vpt 0 360 arc} bind def
+/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
+	neg 0 rlineto closepath} bind def
+/Square {dup Rec} bind def
+/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
+/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
+/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
+/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
+/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
+/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
+/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
+	exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
+/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
+/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
+	2 copy vpt Square fill Bsquare} bind def
+/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
+/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
+/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
+	Bsquare} bind def
+/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
+	Bsquare} bind def
+/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
+/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
+	2 copy vpt Square fill Bsquare} bind def
+/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
+	2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
+/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
+/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
+/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
+/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
+/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
+/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
+/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
+/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
+/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
+/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
+/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
+/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
+/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
+/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
+/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
+/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
+/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
+/DiaE {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath stroke} def
+/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath stroke} def
+/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath stroke} def
+/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath stroke} def
+/PentE {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath stroke grestore} def
+/CircE {stroke [] 0 setdash 
+  hpt 0 360 arc stroke} def
+/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
+/DiaW {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V Opaque stroke} def
+/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V Opaque stroke} def
+/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V Opaque stroke} def
+/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V Opaque stroke} def
+/PentW {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  Opaque stroke grestore} def
+/CircW {stroke [] 0 setdash 
+  hpt 0 360 arc Opaque stroke} def
+/BoxFill {gsave Rec 1 setgray fill grestore} def
+/Density {
+  /Fillden exch def
+  currentrgbcolor
+  /ColB exch def /ColG exch def /ColR exch def
+  /ColR ColR Fillden mul Fillden sub 1 add def
+  /ColG ColG Fillden mul Fillden sub 1 add def
+  /ColB ColB Fillden mul Fillden sub 1 add def
+  ColR ColG ColB setrgbcolor} def
+/BoxColFill {gsave Rec PolyFill} def
+/PolyFill {gsave Density fill grestore grestore} def
+/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
+%
+% PostScript Level 1 Pattern Fill routine for rectangles
+% Usage: x y w h s a XX PatternFill
+%	x,y = lower left corner of box to be filled
+%	w,h = width and height of box
+%	  a = angle in degrees between lines and x-axis
+%	 XX = 0/1 for no/yes cross-hatch
+%
+/PatternFill {gsave /PFa [ 9 2 roll ] def
+  PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
+  PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
+  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
+  clip
+  currentlinewidth 0.5 mul setlinewidth
+  /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
+  0 0 M PFa 5 get rotate PFs -2 div dup translate
+  0 1 PFs PFa 4 get div 1 add floor cvi
+	{PFa 4 get mul 0 M 0 PFs V} for
+  0 PFa 6 get ne {
+	0 1 PFs PFa 4 get div 1 add floor cvi
+	{PFa 4 get mul 0 2 1 roll M PFs 0 V} for
+ } if
+  stroke grestore} def
+%
+/languagelevel where
+ {pop languagelevel} {1} ifelse
+dup 2 lt
+	{/InterpretLevel1 true def
+	 /InterpretLevel3 false def}
+	{/InterpretLevel1 Level1 def
+	 2 gt
+	    {/InterpretLevel3 Level3 def}
+	    {/InterpretLevel3 false def}
+	 ifelse }
+ ifelse
+%
+% PostScript level 2 pattern fill definitions
+%
+/Level2PatternFill {
+/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
+	bind def
+/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} 
+>> matrix makepattern
+/Pat1 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
+	0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
+>> matrix makepattern
+/Pat2 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
+	8 8 L 8 0 L 0 0 L fill}
+>> matrix makepattern
+/Pat3 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
+	0 12 M 12 0 L stroke}
+>> matrix makepattern
+/Pat4 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
+	0 -4 M 12 8 L stroke}
+>> matrix makepattern
+/Pat5 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
+	0 12 M 8 -4 L 4 12 M 10 0 L stroke}
+>> matrix makepattern
+/Pat6 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
+	0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
+>> matrix makepattern
+/Pat7 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
+	12 0 M -4 8 L 12 4 M 0 10 L stroke}
+>> matrix makepattern
+/Pat8 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
+	-4 0 M 12 8 L -4 4 M 8 10 L stroke}
+>> matrix makepattern
+/Pat9 exch def
+/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
+/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
+/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
+/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
+/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
+/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
+/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
+} def
+%
+%
+%End of PostScript Level 2 code
+%
+/PatternBgnd {
+  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
+} def
+%
+% Substitute for Level 2 pattern fill codes with
+% grayscale if Level 2 support is not selected.
+%
+/Level1PatternFill {
+/Pattern1 {0.250 Density} bind def
+/Pattern2 {0.500 Density} bind def
+/Pattern3 {0.750 Density} bind def
+/Pattern4 {0.125 Density} bind def
+/Pattern5 {0.375 Density} bind def
+/Pattern6 {0.625 Density} bind def
+/Pattern7 {0.875 Density} bind def
+} def
+%
+% Now test for support of Level 2 code
+%
+Level1 {Level1PatternFill} {Level2PatternFill} ifelse
+%
+/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
+dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
+currentdict end definefont pop
+%
+%
+% Encoding for ISO-8859-1 (also called Latin1)
+%
+/reencodeISO {
+dup dup findfont dup length dict begin
+{ 1 index /FID ne { def }{ pop pop } ifelse } forall
+currentdict /CharStrings known {
+	CharStrings /Idieresis known {
+		/Encoding ISOLatin1Encoding def } if
+} if
+currentdict end definefont
+} def
+/ISOLatin1Encoding [
+/.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
+/space/exclam/quotedbl/numbersign/dollar/percent/ampersand/quoteright
+/parenleft/parenright/asterisk/plus/comma/minus/period/slash
+/zero/one/two/three/four/five/six/seven/eight/nine/colon/semicolon
+/less/equal/greater/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/backslash/bracketright
+/asciicircum/underscore/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/braceleft/bar/braceright/asciitilde
+/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
+/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
+/.notdef/dotlessi/grave/acute/circumflex/tilde/macron/breve
+/dotaccent/dieresis/.notdef/ring/cedilla/.notdef/hungarumlaut
+/ogonek/caron/space/exclamdown/cent/sterling/currency/yen/brokenbar
+/section/dieresis/copyright/ordfeminine/guillemotleft/logicalnot
+/hyphen/registered/macron/degree/plusminus/twosuperior/threesuperior
+/acute/mu/paragraph/periodcentered/cedilla/onesuperior/ordmasculine
+/guillemotright/onequarter/onehalf/threequarters/questiondown
+/Agrave/Aacute/Acircumflex/Atilde/Adieresis/Aring/AE/Ccedilla
+/Egrave/Eacute/Ecircumflex/Edieresis/Igrave/Iacute/Icircumflex
+/Idieresis/Eth/Ntilde/Ograve/Oacute/Ocircumflex/Otilde/Odieresis
+/multiply/Oslash/Ugrave/Uacute/Ucircumflex/Udieresis/Yacute
+/Thorn/germandbls/agrave/aacute/acircumflex/atilde/adieresis
+/aring/ae/ccedilla/egrave/eacute/ecircumflex/edieresis/igrave
+/iacute/icircumflex/idieresis/eth/ntilde/ograve/oacute/ocircumflex
+/otilde/odieresis/divide/oslash/ugrave/uacute/ucircumflex/udieresis
+/yacute/thorn/ydieresis
+] def
+/Metrics {ExtendTextBox Gswidth} def
+/Lwidth {currentpoint stroke M 0 vshift R Metrics} def
+/Rwidth {currentpoint stroke M dup stringwidth pop neg vshift R Metrics} def
+/Cwidth {currentpoint stroke M dup stringwidth pop -2 div vshift R Metrics} def
+/GLwidth {currentpoint stroke M 0 vshift R {ExtendTextBox} forall} def
+/GRwidth {currentpoint stroke M dup Gwidth vshift R {ExtendTextBox} forall} def
+/GCwidth {currentpoint stroke M dup Gwidth 2 div vshift R {ExtendTextBox} forall} def
+/GLwidth2 {0 Gwidth AddGlyphWidth} def
+/GRwidth2 {Gwidth -1 mul 0 AddGlyphWidth} def
+/GCwidth2 {Gwidth 2 div dup -1 mul AddGlyphWidth} def
+/AddGlyphWidth { dup TBx2 gt {userdict /TBx2 3 -1 roll put} {pop} ifelse
+                 dup TBx1 lt {userdict /TBx1 3 -1 roll put} {pop} ifelse } def
+/MFshow {
+   { dup 5 get 3 ge
+     { 5 get 3 eq {gsave} {grestore} ifelse }
+     {dup dup 0 get findfont exch 1 get scalefont setfont
+     [ currentpoint ] exch dup 2 get 0 exch R dup 5 get 2 ne {dup dup 6
+     get exch 4 get {textshow} {Metrics pop 0 R} ifelse }if dup 5 get 0 eq
+     {dup 3 get {2 get neg 0 exch R pop} {pop aload pop M} ifelse} {dup 5
+     get 1 eq {dup 2 get exch dup 3 get exch 6 get Gswidth pop -2 div
+     dup 0 R} {dup 6 get Gswidth pop -2 div 0 R 6 get
+     textshow 2 index {aload pop M neg 3 -1 roll neg R pop pop} {pop pop pop
+     pop aload pop M} ifelse }ifelse }ifelse }
+     ifelse }
+   forall} def
+/Gswidth {dup type /stringtype eq {stringwidth} {pop (n) stringwidth} ifelse} def
+/MFwidth {0 exch { dup 5 get 3 ge { 5 get 3 eq { 0 } { pop } ifelse }
+ {dup 3 get{dup dup 0 get findfont exch 1 get scalefont setfont
+     6 get Gswidth pop add} {pop} ifelse} ifelse} forall} def
+/MLshow { currentpoint stroke M
+  0 exch R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/MRshow { currentpoint stroke M
+  exch dup MFwidth neg 3 -1 roll R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/MCshow { currentpoint stroke M
+  exch dup MFwidth -2 div 3 -1 roll R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/XYsave    { [( ) 1 2 true false 3 ()] } bind def
+/XYrestore { [( ) 1 2 true false 4 ()] } bind def
+/Times-New-Roman reencodeISO def
+Level1 SuppressPDFMark or 
+{} {
+/SDict 10 dict def
+systemdict /pdfmark known not {
+  userdict /pdfmark systemdict /cleartomark get put
+} if
+SDict begin [
+  /Title (E-distance-NH3-w-bsse.eps)
+  /Subject (gnuplot plot)
+  /Creator (gnuplot 5.2 patchlevel 2)
+%  /Producer (gnuplot)
+%  /Keywords ()
+  /CreationDate (Fri Oct 18 16:12:03 2019)
+  /DOCINFO pdfmark
+end
+} ifelse
+%
+% Support for boxed text - Ethan A Merritt Sep 2016
+%
+/InitTextBox { userdict /TBy2 3 -1 roll put userdict /TBx2 3 -1 roll put
+           userdict /TBy1 3 -1 roll put userdict /TBx1 3 -1 roll put
+	   /Boxing true def } def
+/ExtendTextBox { dup type /stringtype eq
+    { Boxing { gsave dup false charpath pathbbox
+      dup TBy2 gt {userdict /TBy2 3 -1 roll put} {pop} ifelse
+      dup TBx2 gt {userdict /TBx2 3 -1 roll put} {pop} ifelse
+      dup TBy1 lt {userdict /TBy1 3 -1 roll put} {pop} ifelse
+      dup TBx1 lt {userdict /TBx1 3 -1 roll put} {pop} ifelse
+      grestore } if }
+    {} ifelse} def
+/PopTextBox { newpath TBx1 TBxmargin sub TBy1 TBymargin sub M
+               TBx1 TBxmargin sub TBy2 TBymargin add L
+	       TBx2 TBxmargin add TBy2 TBymargin add L
+	       TBx2 TBxmargin add TBy1 TBymargin sub L closepath } def
+/DrawTextBox { PopTextBox stroke /Boxing false def} def
+/FillTextBox { gsave PopTextBox fill grestore /Boxing false def} def
+0 0 0 0 InitTextBox
+/TBxmargin 20 def
+/TBymargin 20 def
+/Boxing false def
+/textshow { ExtendTextBox Gshow } def
+%
+end
+%%EndProlog
+%%Page: 1 1
+gnudict begin
+gsave
+doclip
+50 50 translate
+0.050 0.050 scale
+0 setgray
+newpath
+(Times-New-Roman) findfont 600 scalefont setfont
+BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 17006.00 17006.00 BoxColFill} if
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 2520 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 2520 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-6)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 3498 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 3498 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-5)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 4476 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 4476 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-4)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 5453 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 5453 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-3)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 6431 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 6431 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-2)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 7409 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 7409 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-1)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 8387 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 8387 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 0)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 9365 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 9365 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 1)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 10343 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 10343 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 2)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 11320 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 11320 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 3)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 12298 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 12298 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 4)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 13276 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 13276 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 5)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 14254 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 14254 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 6)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 15232 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 15232 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 7)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 16209 M
+63 0 V
+13642 0 R
+-63 0 V
+stroke
+1860 16209 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 8)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2220 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+2220 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 2.5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+3466 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+3466 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 3)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+4712 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+4712 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 3.5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+5958 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+5958 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 4)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+7204 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+7204 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 4.5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+8450 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+8450 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+9695 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+9695 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 5.5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+10941 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+10941 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 6)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+12187 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+12187 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 6.5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+13433 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+13433 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 7)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+14679 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+14679 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 7.5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+15925 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+15925 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 8)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+1.000 UL
+LTB
+LCb setrgbcolor
+2220 16405 N
+0 -13885 V
+13705 0 V
+0 13885 V
+-13705 0 V
+Z stroke
+1.000 UP
+1.000 UL
+LTb
+LCb setrgbcolor
+LCb setrgbcolor
+540 9462 M
+currentpoint gsave translate -270 rotate 0 0 moveto
+[ [(Times-New-Roman) 650.0 0.0 true true 0 (Binding Energy \(kcal.mol)]
+[(Times-New-Roman) 520.0 325.0 true true 0 (-1)]
+[(Times-New-Roman) 650.0 0.0 true true 0 (\))]
+] -281.7 MCshow
+grestore
+LTb
+LCb setrgbcolor
+9072 1020 M
+[ [(Times-New-Roman) 650.0 0.0 true true 0 (Distance of N and O atoms in cluster NH)]
+[(Times-New-Roman) 520.0 -195.0 true true 0 (3)]
+[(Times-New-Roman) 650.0 0.0 true true 0 (\(H)]
+[(Times-New-Roman) 520.0 -195.0 true true 0 (2)]
+[(Times-New-Roman) 650.0 0.0 true true 0 (O\) \(\305\))]
+] -151.7 MCshow
+LTb
+% Begin plot #1
+3.000 UL
+LTb
+0.00 1.00 0.00 C LCb setrgbcolor
+13342 15952 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (MP2 method)]
+] -200.0 MRshow
+3.000 UL
+LTb
+0.00 1.00 0.00 C 13702 15952 M
+1503 0 V
+2220 10279 M
+32 -422 V
+2388 8349 L
+2523 7099 L
+2659 6074 L
+136 -833 V
+136 -667 V
+136 -525 V
+135 -403 V
+136 -300 V
+135 -213 V
+135 -138 V
+135 -76 V
+136 -25 V
+134 19 V
+135 54 V
+135 83 V
+135 105 V
+135 124 V
+134 138 V
+135 148 V
+135 156 V
+134 160 V
+135 163 V
+135 164 V
+135 164 V
+135 162 V
+135 160 V
+135 156 V
+135 151 V
+136 148 V
+135 143 V
+135 137 V
+136 133 V
+136 127 V
+135 122 V
+136 117 V
+136 111 V
+135 107 V
+136 102 V
+136 96 V
+136 93 V
+136 88 V
+135 83 V
+136 80 V
+136 75 V
+136 72 V
+136 69 V
+136 65 V
+136 62 V
+136 59 V
+136 56 V
+136 53 V
+136 50 V
+135 49 V
+136 45 V
+136 44 V
+136 41 V
+136 40 V
+136 38 V
+136 35 V
+136 34 V
+136 33 V
+136 30 V
+136 30 V
+136 28 V
+135 26 V
+136 26 V
+136 24 V
+136 22 V
+136 22 V
+136 21 V
+136 20 V
+136 18 V
+136 18 V
+136 17 V
+136 16 V
+136 16 V
+135 14 V
+136 14 V
+136 14 V
+136 13 V
+136 12 V
+136 11 V
+136 11 V
+136 11 V
+136 10 V
+136 10 V
+136 9 V
+136 9 V
+135 9 V
+136 8 V
+136 8 V
+136 8 V
+136 7 V
+136 7 V
+% End plot #1
+% Begin plot #2
+stroke
+LTb
+0.58 0.00 0.83 C LCb setrgbcolor
+13342 15172 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (DFTB method 1.28 0.14)]
+] -200.0 MRshow
+3.000 UL
+LTb
+0.58 0.00 0.83 C 13702 15172 M
+1503 0 V
+2220 14205 M
+32 -550 V
+136 -2005 V
+2523 9962 L
+2659 8560 L
+2795 7413 L
+136 -922 V
+136 -724 V
+135 -552 V
+136 -405 V
+135 -279 V
+135 -174 V
+135 -87 V
+136 -16 V
+134 40 V
+135 83 V
+135 117 V
+135 140 V
+135 157 V
+134 167 V
+135 172 V
+135 173 V
+134 170 V
+135 166 V
+135 160 V
+135 152 V
+135 143 V
+135 135 V
+135 127 V
+135 117 V
+136 109 V
+135 102 V
+135 93 V
+136 87 V
+136 80 V
+135 73 V
+136 69 V
+136 63 V
+135 59 V
+136 54 V
+136 51 V
+136 47 V
+136 45 V
+135 41 V
+136 39 V
+136 37 V
+136 35 V
+136 33 V
+136 31 V
+136 29 V
+136 28 V
+136 27 V
+136 25 V
+136 25 V
+135 23 V
+136 22 V
+136 21 V
+136 20 V
+136 20 V
+136 18 V
+136 18 V
+136 17 V
+136 16 V
+136 15 V
+136 15 V
+136 14 V
+135 14 V
+136 13 V
+136 12 V
+136 12 V
+136 12 V
+136 11 V
+136 10 V
+136 10 V
+136 10 V
+136 9 V
+136 9 V
+136 9 V
+135 8 V
+136 8 V
+136 8 V
+136 7 V
+136 7 V
+136 7 V
+136 7 V
+136 6 V
+136 6 V
+136 6 V
+136 6 V
+136 5 V
+135 5 V
+136 6 V
+136 4 V
+136 5 V
+136 5 V
+136 4 V
+% End plot #2
+% Begin plot #3
+stroke
+LTb
+0.00 0.00 1.00 C LCb setrgbcolor
+13342 14392 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (DFTB method 1.16 0.12)]
+] -200.0 MRshow
+3.000 UL
+LTb
+0.00 0.00 1.00 C 13702 14392 M
+1503 0 V
+2266 16405 M
+122 -2002 V
+135 -1903 V
+136 -1606 V
+2795 9553 L
+2931 8449 L
+136 -895 V
+135 -714 V
+136 -555 V
+135 -420 V
+135 -304 V
+135 -208 V
+136 -128 V
+134 -62 V
+135 -10 V
+135 32 V
+135 63 V
+135 87 V
+134 104 V
+135 116 V
+135 122 V
+134 126 V
+135 126 V
+135 125 V
+135 121 V
+135 116 V
+135 111 V
+135 106 V
+135 99 V
+136 94 V
+135 87 V
+135 82 V
+136 76 V
+136 71 V
+135 66 V
+136 62 V
+136 58 V
+135 54 V
+136 50 V
+136 47 V
+136 44 V
+136 42 V
+135 39 V
+136 37 V
+136 35 V
+136 33 V
+136 31 V
+136 30 V
+136 29 V
+136 27 V
+136 25 V
+136 25 V
+136 23 V
+135 23 V
+136 21 V
+136 21 V
+136 19 V
+136 19 V
+136 18 V
+136 17 V
+136 16 V
+136 16 V
+136 15 V
+136 15 V
+136 13 V
+135 14 V
+136 12 V
+136 12 V
+136 12 V
+136 11 V
+136 11 V
+136 10 V
+136 10 V
+136 9 V
+136 9 V
+136 9 V
+136 8 V
+135 8 V
+136 8 V
+136 7 V
+136 7 V
+136 7 V
+136 7 V
+136 6 V
+136 6 V
+136 6 V
+136 6 V
+136 5 V
+136 6 V
+135 5 V
+136 5 V
+136 4 V
+136 5 V
+136 4 V
+136 5 V
+% End plot #3
+% Begin plot #4
+stroke
+LTb
+1.00 0.00 0.00 C LCb setrgbcolor
+13342 13612 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (MP2-BSSE)]
+] -200.0 MRshow
+3.000 UL
+LTb
+1.00 0.00 0.00 C 13702 13612 M
+1503 0 V
+2220 12246 M
+32 -455 V
+136 -1641 V
+2523 8776 L
+2659 7633 L
+136 -942 V
+136 -770 V
+136 -620 V
+135 -493 V
+136 -384 V
+135 -291 V
+135 -212 V
+135 -144 V
+136 -88 V
+134 -41 V
+135 -1 V
+135 31 V
+135 58 V
+135 80 V
+134 97 V
+135 110 V
+135 120 V
+134 129 V
+135 133 V
+135 137 V
+135 138 V
+135 139 V
+135 138 V
+135 137 V
+135 133 V
+136 131 V
+135 128 V
+135 124 V
+136 120 V
+136 115 V
+135 112 V
+136 107 V
+136 103 V
+135 98 V
+136 95 V
+136 90 V
+136 86 V
+136 83 V
+135 79 V
+136 75 V
+136 71 V
+136 68 V
+136 65 V
+136 62 V
+136 58 V
+136 56 V
+136 52 V
+136 51 V
+136 47 V
+135 45 V
+136 43 V
+136 40 V
+136 39 V
+136 36 V
+136 35 V
+136 32 V
+136 32 V
+136 29 V
+136 28 V
+136 27 V
+136 26 V
+135 24 V
+136 24 V
+136 22 V
+136 21 V
+136 20 V
+136 20 V
+136 18 V
+136 18 V
+136 17 V
+136 16 V
+136 16 V
+136 15 V
+135 14 V
+136 13 V
+136 13 V
+136 13 V
+136 11 V
+136 12 V
+136 10 V
+136 11 V
+136 10 V
+136 9 V
+136 9 V
+136 9 V
+135 8 V
+136 8 V
+136 8 V
+136 7 V
+136 6 V
+136 6 V
+% End plot #4
+% Begin plot #5
+stroke
+2.000 UL
+LTb
+LT1
+LCb setrgbcolor
+2220 8387 M
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+139 0 V
+138 0 V
+% End plot #5
+stroke
+LTb
+LCb setrgbcolor
+1.000 UL
+LTB
+LCb setrgbcolor
+2220 16405 N
+0 -13885 V
+13705 0 V
+0 13885 V
+-13705 0 V
+Z stroke
+1.000 UP
+1.000 UL
+LTb
+LCb setrgbcolor
+stroke
+grestore
+end
+showpage
+%%Trailer
+%%DocumentFonts: Times-New-Roman
diff --git a/ammoni/figure_ammoni/E-distance-nh4-w-bsse.eps b/ammoni/figure_ammoni/E-distance-nh4-w-bsse.eps
new file mode 100644
index 0000000..bec17b0
--- /dev/null
+++ b/ammoni/figure_ammoni/E-distance-nh4-w-bsse.eps
@@ -0,0 +1,1430 @@
+%!PS-Adobe-2.0 EPSF-2.0
+%%Title: E-distance-NH4-w-bsse.eps
+%%Creator: gnuplot 5.2 patchlevel 2
+%%CreationDate: Fri Oct 18 15:10:26 2019
+%%DocumentFonts: (atend)
+%%BoundingBox: 50 50 900 900
+%%EndComments
+%%BeginProlog
+/gnudict 256 dict def
+gnudict begin
+%
+% The following true/false flags may be edited by hand if desired.
+% The unit line width and grayscale image gamma correction may also be changed.
+%
+/Color true def
+/Blacktext false def
+/Solid false def
+/Dashlength 1 def
+/Landscape false def
+/Level1 false def
+/Level3 false def
+/Rounded false def
+/ClipToBoundingBox false def
+/SuppressPDFMark false def
+/TransparentPatterns false def
+/gnulinewidth 15.000 def
+/userlinewidth gnulinewidth def
+/Gamma 1.0 def
+/BackgroundColor {-1.000 -1.000 -1.000} def
+%
+/vshift -200 def
+/dl1 {
+  10.0 Dashlength userlinewidth gnulinewidth div mul mul mul
+  Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
+} def
+/dl2 {
+  10.0 Dashlength userlinewidth gnulinewidth div mul mul mul
+  Rounded { currentlinewidth 0.75 mul add } if
+} def
+/hpt_ 31.5 def
+/vpt_ 31.5 def
+/hpt hpt_ def
+/vpt vpt_ def
+/doclip {
+  ClipToBoundingBox {
+    newpath 50 50 moveto 900 50 lineto 900 900 lineto 50 900 lineto closepath
+    clip
+  } if
+} def
+%
+% Gnuplot Prolog Version 5.1 (Oct 2015)
+%
+%/SuppressPDFMark true def
+%
+/M {moveto} bind def
+/L {lineto} bind def
+/R {rmoveto} bind def
+/V {rlineto} bind def
+/N {newpath moveto} bind def
+/Z {closepath} bind def
+/C {setrgbcolor} bind def
+/f {rlineto fill} bind def
+/g {setgray} bind def
+/Gshow {show} def   % May be redefined later in the file to support UTF-8
+/vpt2 vpt 2 mul def
+/hpt2 hpt 2 mul def
+/Lshow {currentpoint stroke M 0 vshift R 
+	Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def
+/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
+	Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def
+/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R 
+	Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def
+/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
+  /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
+/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
+ {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
+/BL {stroke userlinewidth 2 mul setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+/AL {stroke userlinewidth 2 div setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+/UL {dup gnulinewidth mul /userlinewidth exch def
+	dup 1 lt {pop 1} if 10 mul /udl exch def} def
+/PL {stroke userlinewidth setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+3.8 setmiterlimit
+% Classic Line colors (version 5.0)
+/LCw {1 1 1} def
+/LCb {0 0 0} def
+/LCa {0 0 0} def
+/LC0 {1 0 0} def
+/LC1 {0 1 0} def
+/LC2 {0 0 1} def
+/LC3 {1 0 1} def
+/LC4 {0 1 1} def
+/LC5 {1 1 0} def
+/LC6 {0 0 0} def
+/LC7 {1 0.3 0} def
+/LC8 {0.5 0.5 0.5} def
+% Default dash patterns (version 5.0)
+/LTB {BL [] LCb DL} def
+/LTw {PL [] 1 setgray} def
+/LTb {PL [] LCb DL} def
+/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
+/LT0 {PL [] LC0 DL} def
+/LT1 {PL [2 dl1 3 dl2] LC1 DL} def
+/LT2 {PL [1 dl1 1.5 dl2] LC2 DL} def
+/LT3 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC3 DL} def
+/LT4 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
+/LT5 {PL [4 dl1 2 dl2] LC5 DL} def
+/LT6 {PL [1.5 dl1 1.5 dl2 1.5 dl1 1.5 dl2 1.5 dl1 6 dl2] LC6 DL} def
+/LT7 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC7 DL} def
+/LT8 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC8 DL} def
+/SL {[] 0 setdash} def
+/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
+/Dia {stroke [] 0 setdash 2 copy vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath stroke
+  Pnt} def
+/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
+  currentpoint stroke M
+  hpt neg vpt neg R hpt2 0 V stroke
+ } def
+/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath stroke
+  Pnt} def
+/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
+  hpt2 vpt2 neg V currentpoint stroke M
+  hpt2 neg 0 R hpt2 vpt2 V stroke} def
+/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath stroke
+  Pnt} def
+/Star {2 copy Pls Crs} def
+/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath fill} def
+/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath fill} def
+/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath stroke
+  Pnt} def
+/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath fill} def
+/DiaF {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath fill} def
+/Pent {stroke [] 0 setdash 2 copy gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath stroke grestore Pnt} def
+/PentF {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath fill grestore} def
+/Circle {stroke [] 0 setdash 2 copy
+  hpt 0 360 arc stroke Pnt} def
+/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
+/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
+/C1 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C2 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C3 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C4 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 180 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C5 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc
+	2 copy moveto
+	2 copy vpt 180 270 arc closepath fill
+	vpt 0 360 arc} bind def
+/C6 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C7 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C8 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 270 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C9 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 270 450 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
+	2 copy moveto
+	2 copy vpt 90 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C11 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 180 arc closepath fill
+	2 copy moveto
+	2 copy vpt 270 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C12 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 180 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C13 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc closepath fill
+	2 copy moveto
+	2 copy vpt 180 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C14 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 360 arc closepath fill
+	vpt 0 360 arc} bind def
+/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
+	neg 0 rlineto closepath} bind def
+/Square {dup Rec} bind def
+/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
+/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
+/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
+/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
+/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
+/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
+/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
+	exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
+/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
+/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
+	2 copy vpt Square fill Bsquare} bind def
+/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
+/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
+/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
+	Bsquare} bind def
+/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
+	Bsquare} bind def
+/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
+/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
+	2 copy vpt Square fill Bsquare} bind def
+/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
+	2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
+/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
+/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
+/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
+/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
+/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
+/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
+/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
+/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
+/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
+/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
+/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
+/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
+/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
+/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
+/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
+/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
+/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
+/DiaE {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath stroke} def
+/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath stroke} def
+/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath stroke} def
+/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath stroke} def
+/PentE {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath stroke grestore} def
+/CircE {stroke [] 0 setdash 
+  hpt 0 360 arc stroke} def
+/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
+/DiaW {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V Opaque stroke} def
+/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V Opaque stroke} def
+/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V Opaque stroke} def
+/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V Opaque stroke} def
+/PentW {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  Opaque stroke grestore} def
+/CircW {stroke [] 0 setdash 
+  hpt 0 360 arc Opaque stroke} def
+/BoxFill {gsave Rec 1 setgray fill grestore} def
+/Density {
+  /Fillden exch def
+  currentrgbcolor
+  /ColB exch def /ColG exch def /ColR exch def
+  /ColR ColR Fillden mul Fillden sub 1 add def
+  /ColG ColG Fillden mul Fillden sub 1 add def
+  /ColB ColB Fillden mul Fillden sub 1 add def
+  ColR ColG ColB setrgbcolor} def
+/BoxColFill {gsave Rec PolyFill} def
+/PolyFill {gsave Density fill grestore grestore} def
+/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
+%
+% PostScript Level 1 Pattern Fill routine for rectangles
+% Usage: x y w h s a XX PatternFill
+%	x,y = lower left corner of box to be filled
+%	w,h = width and height of box
+%	  a = angle in degrees between lines and x-axis
+%	 XX = 0/1 for no/yes cross-hatch
+%
+/PatternFill {gsave /PFa [ 9 2 roll ] def
+  PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
+  PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
+  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
+  clip
+  currentlinewidth 0.5 mul setlinewidth
+  /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
+  0 0 M PFa 5 get rotate PFs -2 div dup translate
+  0 1 PFs PFa 4 get div 1 add floor cvi
+	{PFa 4 get mul 0 M 0 PFs V} for
+  0 PFa 6 get ne {
+	0 1 PFs PFa 4 get div 1 add floor cvi
+	{PFa 4 get mul 0 2 1 roll M PFs 0 V} for
+ } if
+  stroke grestore} def
+%
+/languagelevel where
+ {pop languagelevel} {1} ifelse
+dup 2 lt
+	{/InterpretLevel1 true def
+	 /InterpretLevel3 false def}
+	{/InterpretLevel1 Level1 def
+	 2 gt
+	    {/InterpretLevel3 Level3 def}
+	    {/InterpretLevel3 false def}
+	 ifelse }
+ ifelse
+%
+% PostScript level 2 pattern fill definitions
+%
+/Level2PatternFill {
+/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
+	bind def
+/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} 
+>> matrix makepattern
+/Pat1 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
+	0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
+>> matrix makepattern
+/Pat2 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
+	8 8 L 8 0 L 0 0 L fill}
+>> matrix makepattern
+/Pat3 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
+	0 12 M 12 0 L stroke}
+>> matrix makepattern
+/Pat4 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
+	0 -4 M 12 8 L stroke}
+>> matrix makepattern
+/Pat5 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
+	0 12 M 8 -4 L 4 12 M 10 0 L stroke}
+>> matrix makepattern
+/Pat6 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
+	0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
+>> matrix makepattern
+/Pat7 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
+	12 0 M -4 8 L 12 4 M 0 10 L stroke}
+>> matrix makepattern
+/Pat8 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
+	-4 0 M 12 8 L -4 4 M 8 10 L stroke}
+>> matrix makepattern
+/Pat9 exch def
+/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
+/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
+/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
+/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
+/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
+/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
+/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
+} def
+%
+%
+%End of PostScript Level 2 code
+%
+/PatternBgnd {
+  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
+} def
+%
+% Substitute for Level 2 pattern fill codes with
+% grayscale if Level 2 support is not selected.
+%
+/Level1PatternFill {
+/Pattern1 {0.250 Density} bind def
+/Pattern2 {0.500 Density} bind def
+/Pattern3 {0.750 Density} bind def
+/Pattern4 {0.125 Density} bind def
+/Pattern5 {0.375 Density} bind def
+/Pattern6 {0.625 Density} bind def
+/Pattern7 {0.875 Density} bind def
+} def
+%
+% Now test for support of Level 2 code
+%
+Level1 {Level1PatternFill} {Level2PatternFill} ifelse
+%
+/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
+dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
+currentdict end definefont pop
+%
+%
+% Encoding for ISO-8859-1 (also called Latin1)
+%
+/reencodeISO {
+dup dup findfont dup length dict begin
+{ 1 index /FID ne { def }{ pop pop } ifelse } forall
+currentdict /CharStrings known {
+	CharStrings /Idieresis known {
+		/Encoding ISOLatin1Encoding def } if
+} if
+currentdict end definefont
+} def
+/ISOLatin1Encoding [
+/.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
+/space/exclam/quotedbl/numbersign/dollar/percent/ampersand/quoteright
+/parenleft/parenright/asterisk/plus/comma/minus/period/slash
+/zero/one/two/three/four/five/six/seven/eight/nine/colon/semicolon
+/less/equal/greater/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/backslash/bracketright
+/asciicircum/underscore/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/braceleft/bar/braceright/asciitilde
+/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
+/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
+/.notdef/dotlessi/grave/acute/circumflex/tilde/macron/breve
+/dotaccent/dieresis/.notdef/ring/cedilla/.notdef/hungarumlaut
+/ogonek/caron/space/exclamdown/cent/sterling/currency/yen/brokenbar
+/section/dieresis/copyright/ordfeminine/guillemotleft/logicalnot
+/hyphen/registered/macron/degree/plusminus/twosuperior/threesuperior
+/acute/mu/paragraph/periodcentered/cedilla/onesuperior/ordmasculine
+/guillemotright/onequarter/onehalf/threequarters/questiondown
+/Agrave/Aacute/Acircumflex/Atilde/Adieresis/Aring/AE/Ccedilla
+/Egrave/Eacute/Ecircumflex/Edieresis/Igrave/Iacute/Icircumflex
+/Idieresis/Eth/Ntilde/Ograve/Oacute/Ocircumflex/Otilde/Odieresis
+/multiply/Oslash/Ugrave/Uacute/Ucircumflex/Udieresis/Yacute
+/Thorn/germandbls/agrave/aacute/acircumflex/atilde/adieresis
+/aring/ae/ccedilla/egrave/eacute/ecircumflex/edieresis/igrave
+/iacute/icircumflex/idieresis/eth/ntilde/ograve/oacute/ocircumflex
+/otilde/odieresis/divide/oslash/ugrave/uacute/ucircumflex/udieresis
+/yacute/thorn/ydieresis
+] def
+/Metrics {ExtendTextBox Gswidth} def
+/Lwidth {currentpoint stroke M 0 vshift R Metrics} def
+/Rwidth {currentpoint stroke M dup stringwidth pop neg vshift R Metrics} def
+/Cwidth {currentpoint stroke M dup stringwidth pop -2 div vshift R Metrics} def
+/GLwidth {currentpoint stroke M 0 vshift R {ExtendTextBox} forall} def
+/GRwidth {currentpoint stroke M dup Gwidth vshift R {ExtendTextBox} forall} def
+/GCwidth {currentpoint stroke M dup Gwidth 2 div vshift R {ExtendTextBox} forall} def
+/GLwidth2 {0 Gwidth AddGlyphWidth} def
+/GRwidth2 {Gwidth -1 mul 0 AddGlyphWidth} def
+/GCwidth2 {Gwidth 2 div dup -1 mul AddGlyphWidth} def
+/AddGlyphWidth { dup TBx2 gt {userdict /TBx2 3 -1 roll put} {pop} ifelse
+                 dup TBx1 lt {userdict /TBx1 3 -1 roll put} {pop} ifelse } def
+/MFshow {
+   { dup 5 get 3 ge
+     { 5 get 3 eq {gsave} {grestore} ifelse }
+     {dup dup 0 get findfont exch 1 get scalefont setfont
+     [ currentpoint ] exch dup 2 get 0 exch R dup 5 get 2 ne {dup dup 6
+     get exch 4 get {textshow} {Metrics pop 0 R} ifelse }if dup 5 get 0 eq
+     {dup 3 get {2 get neg 0 exch R pop} {pop aload pop M} ifelse} {dup 5
+     get 1 eq {dup 2 get exch dup 3 get exch 6 get Gswidth pop -2 div
+     dup 0 R} {dup 6 get Gswidth pop -2 div 0 R 6 get
+     textshow 2 index {aload pop M neg 3 -1 roll neg R pop pop} {pop pop pop
+     pop aload pop M} ifelse }ifelse }ifelse }
+     ifelse }
+   forall} def
+/Gswidth {dup type /stringtype eq {stringwidth} {pop (n) stringwidth} ifelse} def
+/MFwidth {0 exch { dup 5 get 3 ge { 5 get 3 eq { 0 } { pop } ifelse }
+ {dup 3 get{dup dup 0 get findfont exch 1 get scalefont setfont
+     6 get Gswidth pop add} {pop} ifelse} ifelse} forall} def
+/MLshow { currentpoint stroke M
+  0 exch R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/MRshow { currentpoint stroke M
+  exch dup MFwidth neg 3 -1 roll R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/MCshow { currentpoint stroke M
+  exch dup MFwidth -2 div 3 -1 roll R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/XYsave    { [( ) 1 2 true false 3 ()] } bind def
+/XYrestore { [( ) 1 2 true false 4 ()] } bind def
+/Times-New-Roman reencodeISO def
+Level1 SuppressPDFMark or 
+{} {
+/SDict 10 dict def
+systemdict /pdfmark known not {
+  userdict /pdfmark systemdict /cleartomark get put
+} if
+SDict begin [
+  /Title (E-distance-NH4-w-bsse.eps)
+  /Subject (gnuplot plot)
+  /Creator (gnuplot 5.2 patchlevel 2)
+%  /Producer (gnuplot)
+%  /Keywords ()
+  /CreationDate (Fri Oct 18 15:10:26 2019)
+  /DOCINFO pdfmark
+end
+} ifelse
+%
+% Support for boxed text - Ethan A Merritt Sep 2016
+%
+/InitTextBox { userdict /TBy2 3 -1 roll put userdict /TBx2 3 -1 roll put
+           userdict /TBy1 3 -1 roll put userdict /TBx1 3 -1 roll put
+	   /Boxing true def } def
+/ExtendTextBox { dup type /stringtype eq
+    { Boxing { gsave dup false charpath pathbbox
+      dup TBy2 gt {userdict /TBy2 3 -1 roll put} {pop} ifelse
+      dup TBx2 gt {userdict /TBx2 3 -1 roll put} {pop} ifelse
+      dup TBy1 lt {userdict /TBy1 3 -1 roll put} {pop} ifelse
+      dup TBx1 lt {userdict /TBx1 3 -1 roll put} {pop} ifelse
+      grestore } if }
+    {} ifelse} def
+/PopTextBox { newpath TBx1 TBxmargin sub TBy1 TBymargin sub M
+               TBx1 TBxmargin sub TBy2 TBymargin add L
+	       TBx2 TBxmargin add TBy2 TBymargin add L
+	       TBx2 TBxmargin add TBy1 TBymargin sub L closepath } def
+/DrawTextBox { PopTextBox stroke /Boxing false def} def
+/FillTextBox { gsave PopTextBox fill grestore /Boxing false def} def
+0 0 0 0 InitTextBox
+/TBxmargin 20 def
+/TBymargin 20 def
+/Boxing false def
+/textshow { ExtendTextBox Gshow } def
+%
+end
+%%EndProlog
+%%Page: 1 1
+gnudict begin
+gsave
+doclip
+50 50 translate
+0.050 0.050 scale
+0 setgray
+newpath
+(Times-New-Roman) findfont 600 scalefont setfont
+BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 17006.00 17006.00 BoxColFill} if
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 2748 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 2748 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-18)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 3658 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 3658 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-16)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 4569 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 4569 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-14)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 5479 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 5479 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-12)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 6390 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 6390 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-10)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 7300 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 7300 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-8)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 8211 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 8211 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-6)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 9121 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 9121 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-4)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 10032 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 10032 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 (-2)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 10942 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 10942 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 0)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 11853 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 11853 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 2)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 12763 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 12763 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 4)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 13674 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 13674 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 6)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 14584 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 14584 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 8)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 15495 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 15495 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 10)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+2580 16405 M
+63 0 V
+13282 0 R
+-63 0 V
+stroke
+2220 16405 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 12)]
+] -233.3 MRshow
+1.000 UL
+LTb
+LCb setrgbcolor
+3793 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+3793 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 3)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+5526 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+5526 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 4)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+7259 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+7259 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 5)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+8993 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+8993 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 6)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+10726 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+10726 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 7)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+12459 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+12459 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 8)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+14192 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+14192 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 9)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+15925 2520 M
+0 63 V
+0 13822 R
+0 -63 V
+stroke
+15925 1920 M
+[ [(Times-New-Roman) 700.0 0.0 true true 0 ( 10)]
+] -233.3 MCshow
+1.000 UL
+LTb
+LCb setrgbcolor
+1.000 UL
+LTB
+LCb setrgbcolor
+2580 16405 N
+0 -13885 V
+13345 0 V
+0 13885 V
+-13345 0 V
+Z stroke
+1.000 UP
+1.000 UL
+LTb
+LCb setrgbcolor
+LCb setrgbcolor
+540 9462 M
+currentpoint gsave translate -270 rotate 0 0 moveto
+[ [(Times-New-Roman) 650.0 0.0 true true 0 (Binding Energy \(kcal.mol)]
+[(Times-New-Roman) 520.0 325.0 true true 0 (-1)]
+[(Times-New-Roman) 650.0 0.0 true true 0 (\))]
+] -281.7 MCshow
+grestore
+LTb
+LCb setrgbcolor
+9252 1020 M
+[ [(Times-New-Roman) 650.0 0.0 true true 0 (Distance of N and O atoms in cluster NH)]
+[(Times-New-Roman) 520.0 -195.0 true true 0 (4)]
+[(Times-New-Roman) 520.0 325.0 true true 0 (+)]
+[(Times-New-Roman) 650.0 0.0 true true 0 (\(H)]
+[(Times-New-Roman) 520.0 -195.0 true true 0 (2)]
+[(Times-New-Roman) 650.0 0.0 true true 0 (O\) \(\305\))]
+] -216.7 MCshow
+LTb
+% Begin plot #1
+4.000 UL
+LTb
+0.00 1.00 0.00 C LCb setrgbcolor
+13342 15952 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (MP2 method)]
+] -200.0 MRshow
+4.000 UL
+LTb
+0.00 1.00 0.00 C 13702 15952 M
+1503 0 V
+2580 11586 M
+2701 9092 L
+2834 7027 L
+2967 5498 L
+3100 4395 L
+133 -766 V
+133 -503 V
+133 -298 V
+133 -139 V
+133 -17 V
+133 74 V
+134 142 V
+133 191 V
+133 225 V
+133 248 V
+133 263 V
+133 269 V
+133 270 V
+133 269 V
+133 263 V
+133 256 V
+133 247 V
+133 238 V
+133 227 V
+133 217 V
+133 206 V
+133 196 V
+133 186 V
+133 176 V
+133 167 V
+133 157 V
+133 149 V
+134 141 V
+133 134 V
+133 126 V
+133 119 V
+133 113 V
+133 108 V
+133 101 V
+133 96 V
+133 92 V
+133 87 V
+133 82 V
+133 78 V
+133 75 V
+133 70 V
+133 68 V
+133 64 V
+133 62 V
+133 58 V
+133 56 V
+133 53 V
+133 51 V
+134 49 V
+133 46 V
+133 45 V
+133 43 V
+133 41 V
+133 40 V
+133 38 V
+133 36 V
+133 35 V
+133 34 V
+133 32 V
+133 32 V
+133 30 V
+133 29 V
+133 28 V
+133 27 V
+133 27 V
+133 25 V
+133 25 V
+133 23 V
+133 24 V
+133 22 V
+134 22 V
+133 21 V
+133 20 V
+133 20 V
+133 19 V
+133 19 V
+133 18 V
+133 17 V
+133 17 V
+133 17 V
+133 16 V
+133 16 V
+133 15 V
+133 15 V
+133 14 V
+133 14 V
+133 14 V
+133 14 V
+133 13 V
+133 13 V
+133 13 V
+% End plot #1
+% Begin plot #2
+stroke
+3.000 UL
+LTb
+0.00 1.00 1.00 C LCb setrgbcolor
+13342 15172 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (DFTB method 1.28 0.14)]
+] -200.0 MRshow
+3.000 UL
+LTb
+0.00 1.00 1.00 C 13702 15172 M
+1503 0 V
+2580 16058 M
+121 -3346 V
+2834 9933 L
+2967 7864 L
+3100 6352 L
+3233 5278 L
+133 -733 V
+133 -469 V
+133 -266 V
+133 -111 V
+133 3 V
+134 88 V
+133 149 V
+133 190 V
+133 217 V
+133 233 V
+133 241 V
+133 242 V
+133 239 V
+133 233 V
+133 225 V
+133 216 V
+133 205 V
+133 195 V
+133 185 V
+133 174 V
+133 165 V
+133 156 V
+133 147 V
+133 139 V
+133 132 V
+133 124 V
+134 117 V
+133 112 V
+133 105 V
+133 101 V
+133 95 V
+133 91 V
+133 86 V
+133 83 V
+133 79 V
+133 75 V
+133 71 V
+133 69 V
+133 65 V
+133 63 V
+133 60 V
+133 58 V
+133 55 V
+133 53 V
+133 51 V
+133 49 V
+133 47 V
+134 45 V
+133 43 V
+133 42 V
+133 40 V
+133 39 V
+133 37 V
+133 36 V
+133 34 V
+133 34 V
+133 32 V
+133 31 V
+133 30 V
+133 30 V
+133 28 V
+133 27 V
+133 26 V
+133 26 V
+133 24 V
+133 24 V
+133 24 V
+133 22 V
+133 22 V
+134 21 V
+133 21 V
+133 19 V
+133 20 V
+133 19 V
+133 18 V
+133 18 V
+133 17 V
+133 16 V
+133 17 V
+133 16 V
+133 15 V
+133 15 V
+133 15 V
+133 14 V
+133 14 V
+133 13 V
+133 13 V
+133 13 V
+133 13 V
+133 12 V
+% End plot #2
+% Begin plot #3
+stroke
+LTb
+1.00 0.00 0.00 C LCb setrgbcolor
+13342 14392 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (DFTB method 1.16 0.12)]
+] -200.0 MRshow
+3.000 UL
+LTb
+1.00 0.00 0.00 C 13702 14392 M
+1503 0 V
+2580 14803 M
+121 -3193 V
+2834 8975 L
+2967 7028 L
+3100 5623 L
+133 -981 V
+133 -651 V
+133 -396 V
+133 -201 V
+133 -54 V
+133 54 V
+134 133 V
+133 188 V
+133 224 V
+133 247 V
+133 259 V
+133 263 V
+133 261 V
+133 255 V
+133 247 V
+133 237 V
+133 225 V
+133 213 V
+133 202 V
+133 190 V
+133 179 V
+133 169 V
+133 158 V
+133 150 V
+133 140 V
+133 133 V
+133 125 V
+134 119 V
+133 112 V
+133 106 V
+133 101 V
+133 96 V
+133 91 V
+133 87 V
+133 82 V
+133 79 V
+133 75 V
+133 72 V
+133 68 V
+133 66 V
+133 63 V
+133 60 V
+133 58 V
+133 55 V
+133 53 V
+133 51 V
+133 49 V
+133 46 V
+134 45 V
+133 44 V
+133 42 V
+133 40 V
+133 38 V
+133 38 V
+133 36 V
+133 34 V
+133 34 V
+133 32 V
+133 31 V
+133 30 V
+133 29 V
+133 29 V
+133 27 V
+133 26 V
+133 26 V
+133 24 V
+133 24 V
+133 24 V
+133 22 V
+133 22 V
+134 21 V
+133 20 V
+133 20 V
+133 20 V
+133 19 V
+133 18 V
+133 17 V
+133 18 V
+133 16 V
+133 17 V
+133 16 V
+133 15 V
+133 15 V
+133 15 V
+133 14 V
+133 14 V
+133 13 V
+133 13 V
+133 13 V
+133 13 V
+133 12 V
+% End plot #3
+% Begin plot #4
+stroke
+LTb
+0.58 0.00 0.83 C LCb setrgbcolor
+13342 13612 M
+[ [(Times-New-Roman) 600.0 0.0 true true 0 (MP2-BSSE)]
+] -200.0 MRshow
+3.000 UL
+LTb
+0.58 0.00 0.83 C 13702 13612 M
+1503 0 V
+2580 12503 M
+2701 9923 L
+2834 7769 L
+2967 6161 L
+3100 4986 L
+133 -832 V
+133 -561 V
+133 -350 V
+133 -186 V
+133 -59 V
+133 38 V
+134 109 V
+133 163 V
+133 200 V
+133 226 V
+133 242 V
+133 252 V
+133 256 V
+133 255 V
+133 252 V
+133 245 V
+133 239 V
+133 229 V
+133 221 V
+133 211 V
+133 200 V
+133 192 V
+133 181 V
+133 172 V
+133 163 V
+133 155 V
+133 146 V
+134 138 V
+133 131 V
+133 124 V
+133 117 V
+133 111 V
+133 105 V
+133 100 V
+133 94 V
+133 90 V
+133 85 V
+133 81 V
+133 76 V
+133 74 V
+133 69 V
+133 66 V
+133 63 V
+133 61 V
+133 57 V
+133 55 V
+133 53 V
+133 50 V
+134 48 V
+133 46 V
+133 44 V
+133 43 V
+133 40 V
+133 39 V
+133 38 V
+133 36 V
+133 35 V
+133 34 V
+133 32 V
+133 31 V
+133 30 V
+133 29 V
+133 29 V
+133 27 V
+133 26 V
+133 26 V
+133 24 V
+133 24 V
+133 23 V
+133 23 V
+134 21 V
+133 21 V
+133 21 V
+133 20 V
+133 19 V
+133 19 V
+133 18 V
+133 17 V
+133 17 V
+133 17 V
+133 16 V
+133 15 V
+133 15 V
+133 15 V
+133 15 V
+133 14 V
+133 14 V
+133 13 V
+133 14 V
+133 13 V
+133 12 V
+% End plot #4
+% Begin plot #5
+stroke
+2.000 UL
+LTb
+LT1
+LCb setrgbcolor
+2580 10942 M
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+135 0 V
+135 0 V
+134 0 V
+135 0 V
+135 0 V
+% End plot #5
+stroke
+LTb
+LCb setrgbcolor
+1.000 UL
+LTB
+LCb setrgbcolor
+2580 16405 N
+0 -13885 V
+13345 0 V
+0 13885 V
+-13345 0 V
+Z stroke
+1.000 UP
+1.000 UL
+LTb
+LCb setrgbcolor
+stroke
+grestore
+end
+showpage
+%%Trailer
+%%DocumentFonts: Times-New-Roman