1c1
<       program geh_main
---
>       program ph2_main
7,10c7
< C     Modified to keep beta constant, second component peak dropped,
< C     New intrinsic profile used, adb bugs removed etc.
< C
< C      ****************************************************************
---
> 
13,14c10,11
<       real phi,m1,m2,E0,W1,W2,beta,corner(3),minchi2,labmin
<       real vertix(4,3),labprob(1000),labenergy(1000),labmax
---
>       real phi,m1,m2,E0,W1,W2,beta,corner(4),minchi2,labmin
>       real vertix(5,4),labprob(1000),labenergy(1000),labmax
119c116
< 	call pgenv((7.22e-21),(7.28e-21),(2e-23),(labmax + (labmax/10)),
---
> 	call pgenv((7.22e-21),(7.28e-21),(2e-23),(labmax + (labmax/100)),
134c131
< C       betac = (betac / 100) * beta
---
>        betac = (betac / 100) * beta
138c135
< C       betac = beta + betac
---
>        betac = beta + betac
144,145c141,142
< c     * corner(3) = E0
< c     * corner(4) = beta removed
---
> c     * corner(3) = beta
> c     * corner(4) = E0
147c144
<        do i=1,4
---
>        do i=1,5
150,151c147,148
<           vertix(i,3) = E0
< C         vertix(i,4) = E0
---
>           vertix(i,3) = beta
>           vertix(i,4) = E0
156,157c153,154
<        vertix(4,3) = E0c
< C      vertix(5,4) = E0c
---
>        vertix(4,3) = betac
>        vertix(5,4) = E0c
159c156
<        do i=1,4
---
>        do i=1,5
170c167
<        call amoeba(vertix,theresult,4,3,3,ftol,phi2,iter)
---
>        call amoeba(vertix,theresult,5,4,4,ftol,phi2,iter)
172c169
<        do i=1,4
---
>        do i=1,5
175,176c172,173
<           vertix(i,3) = bestE0
< C          vertix(i,4) = bestE0
---
>           vertix(i,3) = bestbeta
>           vertix(i,4) = bestE0
181,182c178,179
< C       vertix(4,3) = betac - (betac / 15)
<        vertix(4,3) = E0c - (E0c /500)
---
>        vertix(4,3) = betac - (betac / 15)
>        vertix(5,4) = E0c - (E0c /500)
188,189c185,186
<        do i=1,4
<           do k=1,3
---
>        do i=1,5
>           do k=1,4
199c196
<        call amoeba(vertix,theresult,4,3,3,ftol,phi2,iter)
---
>        call amoeba(vertix,theresult,5,4,4,ftol,phi2,iter)
206c203
<        write (6,*)'Best beta is ',beta
---
>        write (6,*)'Best beta is ',bestbeta
222c219
< C        call pgtext(7.27e-21,labmax-200,'bestW1(cm^-1)=')
---
>         call pgtext(7.27e-21,labmax-200,'bestW1(cm^-1)=')
225c222
< C        call pgtext(7.22e-21,labmax-400,'bestbeta=')
---
>         call pgtext(7.22e-21,labmax-400,'bestbeta=')
230c227
< C        call pgtext(7.275e-21,labmax-200,charw1cm)
---
>         call pgtext(7.275e-21,labmax-200,charw1cm)
234c231
< C        call pgtext(7.225e-21,labmax-400,charbeta)
---
>         call pgtext(7.225e-21,labmax-400,charbeta)
256c253
< 	real phi,m1,m2,corner(3),calcenergy1(1000),calcprob(1000),beta
---
> 	real phi,m1,m2,corner(4),calcenergy1(1000),calcprob(1000),beta
278c275
< 	  calcprob(iphi)=(1.0/(2*3.14159))*(1.0+beta*0.5*
---
> 	  calcprob(iphi)=(1.0/(2*3.14159))*(1.0+corner(3)*0.5*
280c277
< C	  calcenergy1(iphi)=(m1*m2/(m1+m2))*(corner(3)/m2+2.0*sqrt(E0*
---
> C	  calcenergy1(iphi)=(m1*m2/(m1+m2))*(corner(4)/m2+2.0*sqrt(E0*
284c281
<      &     (corner(3)/m2+(2.0*(sqrt(corner(3)*
---
>      &     (corner(4)/m2+(2.0*(sqrt(E0*
288c285
<              calcenergy2(iphi)=(m1*m2/(m1+m2))*(corner(3)/m2+2.0*sqrt
---
>              calcenergy2(iphi)=(m1*m2/(m1+m2))*(corner(4)/m2+2.0*sqrt
291c288
<         write(9,*)calcprob(iphi),calcenergy1(iphi)
---
> C        write(9,*)calcprob(iphi),calcenergy1(iphi)
337c334
<                  k=min(max(j-3/2,1),177)
---
>               k=min(max(j-3/2,1),177)
339c336
<            call polint(calcenergy(k),calcprob(k),3,labenergy(i),
---
>            call polint(calcenergy(k),calcprob(k),4,labenergy(i),
446,447c443,444
<          bestE0=corner(3)
< C         bestE0=corner(4)
---
>          bestbeta=corner(3)
>          bestE0=corner(4)
453,469d449
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
< 
