Q Evaluation routines

From magneticfields
Revision as of 09:09, 20 May 2014 by Hiroyuki (Talk | contribs)

Jump to: navigation, search

We have evaluated Q using two independent routines, one by Tashiro, second by Ferrer, and performed a triple check using Vachaspati's Mathematica notebook.


Tashiro's C routines

Here are Tashiro's Fortran routines.File:Gamma codes.tar

Ferrer's Python routines

The following Python script takes the files in nosourceoutput.tar.gz as input and calculates Q for each E1 and E2 combination Media:qcalc.py. A small modification of the script that calculates Q separately for the northern and southern hemispheres is given here Media:qcalcnorthsouth.py. Both scripts generate a figure of the output and write the value of Q and the standard error to file. For example, the first script produces the figure Media:jan14z100bgt80.pdf and the file arrays Media:qarrays.tar.gz.

Vachaspati's Mathematica routine

Intended as a check of some Q values generated by other routines: File:Cpgamma-northsouth.nb

Retrieved from "https://sites.physics.wustl.edu/magneticfields/wiki/index.php?title=Q_Evaluation_routines&oldid=117"