Nigel's gosia page

This is not the official gosia web page. This page has various versions of gosia, most of which have not been tested. They are intended to allow users who find bugs, to track down which versions have the bug and which don't.

The official gosia home page is in Rochester at: http://www.pas.rochester.edu/~cline/Gosia/index.html. There you will find the most recent versions of gosia, gosia2, annl, sigma, select, and gremlin. and the manual, which Doug Cline is maintaining.

The Warsaw group also has a gosia page with some talks by Magda and Kasia to introduce beginniers to Gosia: http://www.slcj.uw.edu.pl/en/81.html. It also has a link to the web page for the first gosia workshop (Apr 2008) in Warsaw: http://www.slcj.uw.edu.pl/gosia_workshop" which has the presentations made at that workshop.

Here you can find various versions of gosia. They should (in theory) be similar in the way they work, but if you find a bug in the most recent version, please test some of the older versions to help track down where the bug was introduced. The versions have a date stamp YYYYMMDD and possibly a sub-version. Sub-versions are for bugfixes only (relative to the major version with the date stamp).

Gosia

For target OR projectile excitations, separately.
Date File Comments
02 Apr 2023 gosia_20110524.11.f Bugfix relative to gosia_20110524.10 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
21 Mar 2023 gosia_20110524.10.f Improvement relative to 20110524.9 (Use double precision constants. Fix Leuven efficiency parameterisation).
08 May 2019 gosia_20110524.9.f Improvement relative to 20110524.8 (Asymmetric errors in OP,YIEL; correct finite-distance correction; correct delta; NIST absorption coefficients).
07 Mar 2019 gosia_20110524.8.f Improvement relative to 20110524.7 (Add E0 parameter to Leuven efficiency).
05 Oct 2018 gosia_20110524.7.f Bugfix relative to 20110524.6 (Increase default for CONT option INT from 1 to 1000).
05 Feb 2018 gosia_20110524.6.f Bugfix relative to 20110524.5 (Catch errors in input).
28 Nov 2016 gosia_20110524.5.f Bugfix relative to 20110524.4 (Explicitly convert to avoid warnings).
12 Feb 2016 gosia_20110524.4.f Bugfix relative to 20110524.3 (Allow different efficiency parameters for different experiments).
07 Mar 2014 gosia_20110524.3.f Bugfix relative to 20110524.2 (Fix bug in determination of IKIN in OP,INTI).
07 Nov 2011 gosia_20110524.2.f Bugfix relative to 20110524.1 (Fix bug in lambda=6 mu=3 collision functions).
10 Oct 2011 gosia_20110524.1.f Bugfix relative to 20110524 (Fix bug which causes a failure when a meshpoint in inverse kinematics is at the target angle corresponding to the maximum projectile angle).
24 May 2011 gosia_20110524.f Version including extra options for Rachel.
02 April 2023 gosia_20081208.19.f Bugfix relative to gosia_20081208.18 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
09 May 2019 gosia_20081208.18.f Bug fix relative to 20081208.17 (Better OP,RAND; asymmetric errors in OP,YIEL, NIST absorption coefficients; fix finite-recoil correction; fix multipole-mixing ratio).
07 Mar 2019 gosia_20081208.17.f Improvement relative to 20081208.16 (Add E0 parameter to Leuven efficiency).
05 Oct 2018 gosia_20081208.16.f Bugfix relative to 20081208.15 (Increase default for CONT option INT from 1 to 1000).
05 Feb 2018 gosia_20081208.15.f Bugfix relative to 20081208.14 (Catch errors in input).
07 Mar 2014 gosia_20081208.14.f Bugfix relative to 20081208.13 (Explicitly convert to avoid warnings).
07 Mar 2014 gosia_20081208.13.f Bugfix relative to 20081208.12 (Fix bug in determination of IKIN in OP,INTI).
07 Nov 2011 gosia_20081208.12.f Bugfix relative to 20081208.11 (Fix bug in lambda=6 mu=3 collision functions).
10 Oct 2011 gosia_20081208.11.f Bugfix relative to 20081208.10 (Fix bug which causes a failure when a meshpoint in inverse kinematics is at the target angle corresponding to the maximum projectile angle).
24 Jun 2010 gosia_20081208.10.f Bugfix relative to 20081208.9 (Fix bug in relativistic correction for inverse kinematics).
22 Feb 2010 gosia_20081208.9.f Bugfix relative to 20081208.8 (Fix discontinuity in TRINT).
18 Sep 2009 gosia_20081208.8.f Improvements relative to 20081208.7 (Increased dimensions in ZETA to allow 999 matrix elements correctly this time).
14 Sep 2009 gosia_20081208.7.f Improvements relative to 20081208.6 (Increased dimensions in ZETA to allow 999 matrix elements).
16 Aug 2009 gosia_20081208.6.f Bugfix relative to 20081208.5 (Increased dimensions in VLIN to 101).
20 Jul 2009 gosia_20081208.5.f Bugfix relative to 20081208.4 (integration over PIN diodes was incorrect).
2 Apr 2009 gosia_20081208.4.f Bugfix relative to 20081208.3 (E1 polarization was incorrect).
1 Feb 2009 gosia_20081208.3.f Bugfix relative to 20081208.2 (OP,INTI preserves IKIN flag).
27 Jan 2009 gosia_20081208.2.f Bugfix relative to 20081208.1 (OP,INTI got wrong particle for target excitations).
16 Jan 2009 gosia_20081208.1.f Bugfix relative to 20081208 (give error if too many omega steps).
8 Dec 2008 gosia_20081208.f New OP,INTI option for inverse kinematics (beta version) - see doc here.
02 Apr 2023 gosia_20080630.13.f Bugfix relative to gosia_20080630.12 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
07 Mar 2019 gosia_20080630.12.f Bugfix relative to 20080630.11 (Add E0 parameter to Leuven efficiency).
05 Oct 2018 gosia_20080630.11.f Bugfix relative to 20080630.10 (Increase default INT parameter in CONT option to 1000).
5 Feb 2018 gosia_20080630.10.f Bugfix relative to 20080630.9 (Catch errors in input).
07 Nov 2011 gosia_20080630.9.f Bugfix relative to 20080630.8 (Fix bug in lambda=6 mu=3 collision functions).
24 Jun 2010 gosia_20080630.8.f Bugfix relative to 20080630.7 (Fix bug in relativistic correction for inverse kinematics).
22 Feb 2010 gosia_20080630.7.f Bugfix relative to 20080630.6 (Fix discontinuity in TRINT).
18 Sep 2009 gosia_20080630.6.f Improvements relative to 20080630.5 (Increased dimensions in ZETA to allow 999 matrix elements correctly this time).
14 Sep 2009 gosia_20080630.5.f Improvements relative to 20080630.4 (Increased dimensions in ZETA to allow 999 matrix elements).
16 Aug 2009 gosia_20080630.4.f Bugfix relative to 20080630.3 (Increased dimensions in VLIN to 101).
20 Jul 2009 gosia_20080630.3.f Bugfix relative to 20080630.2 (integration over PIN diodes was incorrect).
02 Apr 2009 gosia_20080630.2.f Bugfix relative to 20080630.1 (E1 polarization was incorrect).
16 Jan 2009 gosia_20080630.1.f Bugfix relative to 20080630 (Bricc files opened 'OLD').
30 Jun 2008 gosia_20080630.f Up to 100 levels, bugfixes, Radware efficiency calibration.
02 April 2023 gosia_20080519.15.f Bugfix relative to 20080519.14 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
07 Mar 2019 gosia_20080519.14.f Bugfix relative to 20080519.13 (Add E0 parameter to Leuven efficiency).
05 Oct 2018 gosia_20080519.13.f Bugfix relative to 20080519.12 (Increase default INT parameter in CONT option to 1000).
05 Feb 2018 gosia_20080519.12.f Bugfix relative to 20080519.11 (Catch errors in input).
07 Nov 2011 gosia_20080519.11.f Bugfix relative to 20080519.10 (Fix bug in lambda=6 mu=3 collision functions).
24 Jun 2010 gosia_20080519.10.f Bugfix relative to 20080519.9 (Fix bug in relativistic correction for inverse kinematics).
22 Feb 2010 gosia_20080519.9.f Bugfix relative to 20080519.8 (Fix dimensions in VLIN - bug introduced in 200080519.5).
22 Feb 2010 gosia_20080519.8.f Bugfix relative to 20080519.7 (Fix discontinuity in TRINT).
18 Sep 2009 gosia_20080519.7.f Improvements relative to 20080519.6 (Increased dimensions in ZETA to allow 999 matrix elements correctly this time).
14 Sep 2009 gosia_20080519.6.f Improvements relative to 20080519.5 (Increased dimensions in ZETA to allow 999 matrix elements).
16 Aug 2009 gosia_20080519.5.f Bugfix relative to 20080519.4 (Increased dimensions in VLIN to 101).
20 Jul 2009 gosia_20080519.4.f Bugfix relative to 20080519.3 (integration over PIN diodes was incorrect).
02 Apr 2009 gosia_20080519.3.f Bugfix relative to 20080519.2 (E1 polarization was incorrect).
16 Jan 2009 gosia_20080519.2.f Bugfix relative to 20080519.1 (give error if too many omega steps).
19 May 2008 gosia_20080519.1.f Bugfix relative to 20080519 (several small bugs).
19 May 2008 gosia_20080519.f Spline, Bricc option for conversion coefficients, OP,SELE.
08 May 2008 gosia_20080508.f Up to 999 matrix elements.
07 May 2008 gosia_20080507.f Corrected constants.
18 Apr 2008 gosia_20080418.f Allow up to 50 meshpoints in energy and theta, Leuven efficiency calibration.
28 Mar 2008 gosia_20080328.f Minor bugfix.
29 Feb 2008 gosia_20080229.f Warn user if NI2 exceeds 50.
26 Jul 2007 gosia_20070726.f Cleanup of code.
21 Jul 2007 gosia_20070721.f Small bugfixes.
17 Jul 2007 gosia_20070717.f Small bugfixes.
09 Jul 2007 gosia_20070709.f 64-bit version.
Tomek's last version (Jun 2006) gosia_06.f 32 bit. Needs special compiler options (e.g. -fno-automatic for g77).

The changes from one major version to the next are given here: gosia_changelog.txt. The differences between minor versions are bug fixes, correction of typos etc.

Note that in addition to the changes listed in the change log, there are lots of changes, which involve the addition of comments.

Gosia2

For handling target and projectile excitations simultaneously.
Date File Comments
03 Apr 2023 gosia2_20081208.26.f Bugfix relative to gosia2_20081208.25 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
22 Mar 2021 gosia2_20081208.25.f Fix bug in getting conversion coefficients with BrIcc during minimisation
09 May 2019 gosia2_20081208.24.f Asymmetric errors in OP,YIEL; fix multipole-mixing ratios; fix finite-recoil correction; use NIST absorption coefficients; better OP,RAND
07 May 2019 gosia2_20081208.23.f Liam's improvements for convergence of OP,MINI
07 Mar 2019 gosia2_20081208.22.f Add E0 parameter to Leuven efficiency.
28 May 2018 gosia2_20081208.21.f Rewind matrix element file before writing to it.
13 Apr 2018 gosia2_20081208.20.f Use different files for conversion coefficients of projectile/target.
06 Feb 2018 gosia2_20081208.19.f Trap some common errors in inputs
28 Nov 2016 gosia2_20081208.18.f Clean up some warnings by explicitly converting types.
13 Oct 2016 gosia2_20081208.17.f Increase the maximum number of iterations.
12 Feb 2016 gosia2_20081208.16.f Allow different efficiency parameters for different experiments.
7 Mar 2014 gosia2_20081208.15.f Bugfix relative to 20081208.14 (Fix bug in determination of IKIN in OP,INTI).
07 Nov 2011 gosia2_20081208.14.f Bugfix relative to 20081208.13 (Fix bug in lambda=6 mu=3 collision functions).
20 May 2011 gosia2_20081208.13.f Bugfix relative to 20081208.12 (Fix bug which causes a failure when a meshpoint in inverse kinematics is at the target angle corresponding to the maximum projectile angle).
20 May 2011 gosia2_20081208.12.f Bugfix relative to 20081208.11 (Fix bug which overwrites files).
24 Jun 2010 gosia2_20081208.11.f Bugfix relative to 20081208.10 (Fix bug in relativistic correction for inverse kinematics).
23 Apr 2010 gosia2_20081208.10.f Bugfix relative to 20081208.9 (Fix LCK, control option bug).
22 Feb 2010 gosia2_20081208.9.f Bugfix relative to 20081208.8 (Fix TRINT discontinuity bug).
14 Sep 2009 gosia2_20081208.8.f Improvement relative to 20081208.7 (Increased dimensions in ZETA to allow 999 matrix elements).
16 Aug 2009 gosia2_20081208.7.f Bugfix relative to 20081208.6 (Increased dimensions in VLIN to 101).
20 Jul 2009 gosia2_20081208.6.f Bugfix relative to 20081208.5 (integration over PIN diodes was incorrect).
02 Apr 2009 gosia2_20081208.5.f Bugfix relative to 20081208.4 (E1 polarization was incorrect).
01 Feb 2009 gosia2_20081208.4.f Bugfix relative to 20081208.3 (OP,INTI preserves IKIN flag).
27 Jan 2009 gosia2_20081208.3.f Bugfix relative to 20081208.2 (OP,INTI got wrong particle for target excitations).
23 Jan 2009 gosia2_20081208.2.f Bugfix relative to 20081208.1 (multiply by dsigma*SIN(theta) for OP,INTI as well as OP,INTG).
19 Jan 2009 gosia2_20081208.1.f Bugfix relative to 20081208 (give error if too many omega steps).
08 Dec 2008 gosia2_20081208.f New OP,INTI option for inverse kinematics (beta version) - see doc here.
03 Apr 2023 gosia2_20080630.14.f Improvement relative to gosia2_20080630.13 Liam's improvements for convergence of OP,MINI
03 Apr 2023 gosia2_20080630.13.f Bugfix relative to gosia2_20080630.12 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
07 Mar 2019 gosia2_20080630.12.f Improvement relative to 20080630.11 (Add E0 parameter to Leuven efficiency).
29 May 2018 gosia2_20080630.11.f Bugfix relative to 20080630.10 (Fix precision of adiabatic constant and sinh/cosh and rewind matrix element file before writing to it).
07 Nov 2011 gosia2_20080630.10.f Bugfix relative to 20080630.9 (Fix bug in lambda=6 mu=3 collision functions).
24 Jun 2010 gosia2_20080630.9.f Bugfix relative to 20080630.8 (Fix bug in relativistic correction for inverse kinematics).
23 Apr 2010 gosia2_20080630.8.f Bugfix relative to 20080630.7 (Fix LCK, control option bug).
22 Feb 2010 gosia2_20080630.7.f Bugfix relative to 20080630.6 (Fix TRINT discontinuity bug).
14 Sep 2009 gosia2_20080630.6.f Improvement relative to 20080630.5 (Increased dimensions in ZETA to allow 999 matrix elements).
16 Aug 2009 gosia2_20080630.5.f Bugfix relative to 20080630.4 (Increased dimensions in VLIN to 101).
20 Jul 2009 gosia2_20080630.4.f Bugfix relative to 20080630.3 (integration over PIN diodes was incorrect).
02 Apr 2009 gosia2_20080630.3.f Bugfix relative to 20080630.2 (E1 polarization was incorrect).
19 Jan 2009 gosia2_20080630.2.f Bugfix relative to 20080630.1 (give error if too many omega steps).
30 Jun 2008 gosia2_20080630.1.f Bugfix relative to 20080630 (Bricc files opened 'OLD').
30 Jun 2008 gosia2_20080630.f Up to 100 levels and 999 matrix elements, bugfixes, Radware efficiency calibration. KNOWN BUG: OP,BRIC is broken (use 20080630-1 instead).
03 Apr 2023 gosia2_20080519.13.f Improvement relative to gosia2_20080519.12 Liam's improvements for convergence of OP,MINI
03 Apr 2023 gosia2_20080519.12.f Bugfix relative to gosia2_20080519.11 (Explicitly initialise idr. Correct errorbars on B(E2) values. Write absolute value of lab angle. Fix format for long lifetimes. Change default value of DIPOL.)
07 Mar 2019 gosia2_20080519.11.f Improvement relative to 20080519.10 (Add E0 parameter to Leuven efficiency).
29 May 2018 gosia2_20080519.10.f Bugfix relative to 20080519.9 (Fix precision of adiabatic constant and sinh/cosh and rewind matrix element file before writing to it).
07 Nov 2011 gosia2_20080519.9.f Bugfix relative to 20080519.8 (Fix bug in lambda=6 mu=3 collision functions).
24 Jun 2010 gosia2_20080519.8.f Bugfix relative to 20080519.7 (Fix bug in relativistic correction for inverse kinematics).
23 Apr 2010 gosia2_20080519.7.f Bugfix relative to 20080519.6 (Fix LCK, control option bug).
22 Feb 2010 gosia2_20080519.6.f Bugfix relative to 20080519.5 (Use FILE= not NAME= in open statement).
22 Feb 2010 gosia2_20080519.5.f Bugfix relative to 20080519.4 (Fix TRINT discontinuity bug).
16 Aug 2009 gosia2_20080519.4.f Bugfix relative to 20080519.3 (Increased dimensions in VLIN to 101).
20 Jul 2009 gosia2_20080519.3.f Bugfix relative to 20080519.2 (integration over PIN diodes was incorrect).
02 Apr 2009 gosia2_20080519.2.f Bugfix relative to 20080519.1 (E1 polarization was incorrect).
19 Jan 2009 gosia2_20080519.1.f Bugfix relative to 20080519 (give error if too many omega steps).
19 May 2008 gosia2_20080519.f Spline, Bricc option for conversion coefficients, OP,SELE, corrected constants.
18 Apr 2008 gosia2_20080418.f Cleanup a bit.
17 Apr 2008 gosia2_20080417.f Allow up to 50 meshpoints in energy and theta, give error if NI2 exceeds 50.
14 Apr 2008 gosia2_20080414.f 64-bit version.

The changes from one major version to the next are given here: gosia2_changelog.txt. The differences between minor versions are bug fixes, correction of typos etc.

Note that in addition to the changes listed in the change log, there are lots of changes, which involve the addition of comments.