Joseph Jensen: Research Log for Summer 2012
Solar Physics REU at Montana
State University, Bozeman, MT |
Home
|
2012_08_09 summary
So I figured that I better give some good information to those that will carry on the work, or help you try to use the code that I did, we'll see how it goes. I wrote up all my instructions for taking the whole process through and it is found in /nfs/hl2/data/jjensen/how_to_calculate_energy_instructions_joes_notes.txt
here is a link to the instructions
instructions
I was able to get 3 active regions done, ar11112, ar11092 at 0500 aug 1 and 1700 july 1, and ar11166. They are all located in the directory
/nfs/hl2/data/jjensen
2012_08_01 ar11112 16 Oct 2011 16:50 UT
Here is all the fancy graphical stuff that I need to put in for ar11112
Here is the results of all the wavelengths, in graphical form

Here are the loops all put together , there are 79 in total

Here are the loops from the four different wavelengths
171 and 193

211 and 335

2012_07_31 ar11112 16 Oct 2011 16:50 UT
For the totals with a combined 79 loops, we get the following energies for ar11112
ERG_FREE = 1.1268797e+31
ERG_FREE_UNC = 4.6774039e+30
ERG_TOT = 8.4537902e+31
ERG_TOT_UNC = 1.2330496e+31
B_POTENTIAL = 8.73035e+31
In the previous active regions that I have done the B_potential has been in very good agreement with the difference of erg_tot and erg_free, but in this case it is off. This active region was much more difficult and the uncertainties that we got in each of the individual wavelengths was much greater than normal. The coronal loops that I could get were very small in length and it made it trickier to fit them. But we did get a number that doesn't seem to inconsistent with what Lucas figured in his paper.
This is the active region ar11112, it was really hard to see any good coronal loops, they were all really low and it was especially difficult to see in the 171 wavelength, so the loop count is really low, but maybe we can get enough when we combine all 4 wavelengths
171 angstrom wavelength with 16 loops, 16 October 2010 16:50 UT
ERG_FREE = 4.7121576e+30
ERG_FREE_UNC = 2.7294846e+30
ERG_TOT = 6.5640676e+31
ERG_TOT_UNC = 1.5541070e+31
193 angstrom wavelength with 20 loops
ERG_FREE = 8.4401558e+30
ERG_FREE_UNC = 3.6710469e+30
ERG_TOT = 7.1039003e+31
ERG_TOT_UNC = 1.4788642e+31
211 angstrom wavelength with 25 loops
ERG_FREE = 4.6339308e+30
ERG_FREE_UNC = 1.7375636e+30
ERG_TOT = 6.2877772e+31
ERG_TOT_UNC = 1.2688882e+31
335 angstrom wavelength with 18 loops
ERG_FREE = 1.7215020e+31
ERG_FREE_UNC = 1.2631532e+31
ERG_TOT = 6.2910905e+31
ERG_TOT_UNC = 2.6109644e+31
2012_07_27 Monte Carlo practice
So I have been experimenting the monte carlo method for computing integrals, and it is cool, I also found out that the monte carlo method was named after a casino that one of the people that used it first's uncle always lost money at, how cool is that(good old wikipedia). any way here is a link to my idl code that has my simple version of the method in IDL. It is accurate to about 3 significant figures.
monte_carlo_practice.pro
2012_07_26 ar11092 at 1700 UT July 31
I put the numbers below, but here is my summary of them. across the board the energy is less 12 hours before our first assesment. Also interesting to note is that for the big calculations (with 80 and 83 loops) the b_potential that is calculated a different way agrees exceptionally well with the total-free energies. The free energy increased by a factor of 2.5 and the total energy increased by about 50%, also the b_potential increased by about 20%.
Here are the graphs of the results, once again the wavelength 193 is a lot lower than the others, with a very small uncertainty

The final calculation for the energy of the ar11092 at 1700 UT on July 1, 2010 with 83 loops is (B_potential is calculated by a completly different means than just total-free, but it agrees very well with total-free)
ERG_FREE = 1.0015958e+32
ERG_FREE_UNC = 2.8379229e+31
ERG_TOT = 4.7210625e+32
ERG_TOT_UNC = 5.0536270e+31
B_POTENTIAL = 3.77818e+32
the energy for ar11092 at 0500 UT on Aug 1 for 80 loops was
ERG_FREE = 2.5894019e+32
ERG_FREE_UNC = 9.8420592e+31
ERG_TOT = 7.2127592e+32
ERG_TOT_UNC = 1.3928384e+32
B_POTENTIAL = 4.5112687e+32
Here is the energy for the wavelengths of 211 and 335
211 angstrom wavelength with a total of 23 loops 1700 UT july 31
ERG_FREE = 7.3163972e+31
ERG_FREE_UNC = 2.8444000e+31
ERG_TOT = 3.8837488e+32
ERG_TOT_UNC = 5.6001428e+31
211 angstrom wavelength with 22 loops at 0500 UT aug 1
ERG_FREE = 1.1515629e+32
ERG_FREE_UNC = 4.1025398e+31
ERG_TOT = 4.7274411e+32
ERG_TOT_UNC = 8.8767467e+31
335 angstrom wavelength with a total of 10 loops, 1700 UT July 1
ERG_FREE = 1.1682371e+32
ERG_FREE_UNC = 6.0129818e+31
ERG_TOT = 3.9744038e+32
ERG_TOT_UNC = 1.1345828e+32
335 angstrom wavelength with 8 loops, 0500 UT Aug 1
ERG_FREE = 1.9050721e+32
ERG_FREE_UNC = 1.1399854e+32
ERG_TOT = 4.9840895e+32
ERG_TOT_UNC = 1.6048665e+32
2012_07_25 ar11092 at 1700 UT July 31
here is the first results for the separate wavelengths
171 angstrom wavelength with 29 loops at 1700 UT
ERG_FREE = 5.5831548e+31
ERG_FREE_UNC = 3.1909274e+31
ERG_TOT = 3.6561499e+32
ERG_TOT_UNC = 5.8162828e+31
171 angstrom wavelength with 28 loops at 0500 UT the next day (for reference)
ERG_FREE = 1.4118431e+32
ERG_FREE_UNC = 1.0529716e+32
ERG_TOT = 4.2127213e+32
ERG_TOT_UNC = 1.3991548e+32
193 angstrom wavelenght with 21 loops at 1700 UT
ERG_FREE = 1.7252722e+31
ERG_FREE_UNC = 6.4153677e+30
ERG_TOT = 2.6950729e+32
ERG_TOT_UNC = 4.1758684e+31
193 angstrom wavelenght with 22 loops at 0500 UT the next day (for reference)
ERG_FREE = 8.3802320e+31
ERG_FREE_UNC = 3.1845270e+31
ERG_TOT = 5.1409281e+32
ERG_TOT_UNC = 9.3998203e+31
2012_07_23 ar11092 at 0500 UT on August 1, 2010
Here are the graphs of the free and total energies, I have included the ones that still contain the open field lines, as a reference. The numbers are found below the graphs.
Here is the free engery with NO open field lines

Here is the free energy WITH the open field lines

Here is the total energy with NO open field lines

here is the total energy WITH open field lines

Here is the energy calculated with all 80 loops for ar11092
ERG_FREE = 2.5894019e+32
ERG_FREE_UNC = 9.8420592e+31
ERG_TOT = 7.2127592e+32
ERG_TOT_UNC = 1.3928384e+32
just for reference the other calculation with the open loops(total of 130 loops) was
ERG_FREE = 1.9204604e+32
ERG_FREE_UNC = 5.2123262e+31
ERG_TOT = 6.3397548e+32
ERG_TOT_UNC = 8.2420424e+31
171 angstrom wavelength with 28 loops (no open loops)
ERG_FREE = 1.4118431e+32
ERG_FREE_UNC = 1.0529716e+32
ERG_TOT = 4.2127213e+32
ERG_TOT_UNC = 1.3991548e+32
193 angstrom wavelenght with 22 loops (no open loops)
ERG_FREE = 8.3802320e+31
ERG_FREE_UNC = 3.1845270e+31
ERG_TOT = 5.1409281e+32
ERG_TOT_UNC = 9.3998203e+31
211 angstrom wavelength with 22 loops (no open loops)
ERG_FREE = 1.1515629e+32
ERG_FREE_UNC = 4.1025398e+31
ERG_TOT = 4.7274411e+32
ERG_TOT_UNC = 8.8767467e+31
335 angstrom wavelength with 8 loops (no open loops)
ERG_FREE = 1.9050721e+32
ERG_FREE_UNC = 1.1399854e+32
ERG_TOT = 4.9840895e+32
ERG_TOT_UNC = 1.6048665e+32
Here is a picture of all 80 loops that were used for this calculation

2012_07_23 ar11166
Here are the loops for the separate wavelengths for ar11166
171 angstrom 46 loops used for computation 82 identified________________ 193 angstrom 49 loops used for computation 65 identified

211 angstrom 38 loops used for computation 45 identified________________ 335 angstrom 29 loops used for computation 36 identified

I have gone through and I have taken out the loops that were open in active region 11166, and these are the new energies that I have calculated, I put the graphs first, and then the numbers are below.


This is using all the loops for ar11166, there are 164 loops total, and none of them are open loops
ERG_FREE = 2.5109260e+32
ERG_FREE_UNC = 1.0641941e+32
ERG_TOT = 9.6218660e+32
ERG_TOT_UNC = 1.7317870e+32
171 angstrom energy calculations with 46 loops (with no open loops)
ERG_FREE = 2.4201141e+32
ERG_FREE_UNC = 1.5912902e+32
ERG_TOT = 9.7237749e+32
ERG_TOT_UNC = 2.7254838e+32
193 angstrom energy calculations with 49 loops(with no open loops)
ERG_FREE = 3.6577528e+31
ERG_FREE_UNC = 1.8569837e+31
ERG_TOT = 7.1546491e+32
ERG_TOT_UNC = 1.7992224e+32
211 angstrom energy calculations with 38 loops (with no open loops)
ERG_FREE = 3.1796010e+32
ERG_FREE_UNC = 1.6965099e+32
ERG_TOT = 9.7700127e+32
ERG_TOT_UNC = 2.5826349e+32
335 angstrom energy calculations with 29 loops(with no open loops)
ERG_FREE = 6.8843366e+31
ERG_FREE_UNC = 3.2649502e+31
ERG_TOT = 6.3550476e+32
ERG_TOT_UNC = 1.4238922e+32
Here is a picture of all 164 loops that were drawn

2012_07_17
So here are the first of the energies for the active region 11166
171 angstrom energy calculations with 67 loops
ERG_FREE = 3.0284701e+32
ERG_FREE_UNC = 1.6882441e+32
ERG_TOT = 1.0667411e+33
ERG_TOT_UNC = 2.7894277e+32
193 angstrom energy calculations with 60 loops
ERG_FREE = 4.3170972e+31
ERG_FREE_UNC = 1.6581506e+31
ERG_TOT = 7.9505390e+32
ERG_TOT_UNC = 1.8053155e+32
211 angstrom energy calculations with 43 loops
ERG_FREE = 2.6178505e+32
ERG_FREE_UNC = 1.3706560e+32
ERG_TOT = 9.3362093e+32
ERG_TOT_UNC = 2.2209015e+32
335 angstrom energy calculations with 29 loops(I made sure there were no open loops, I will fix the others later)
ERG_FREE = 6.8843366e+31
ERG_FREE_UNC = 3.2649502e+31
ERG_TOT = 6.3550476e+32
ERG_TOT_UNC = 1.4238922e+32
2012_07_09
Here is the grand total with all the loops put together for ar11092 with a total of 130 loops
ERG_FREE = 1.9204604e+32
ERG_FREE_UNC = 5.2123262e+31
ERG_TOT = 6.3397548e+32
ERG_TOT_UNC = 8.2420424e+31
These are the calculations of the four separate wavelengths for the active region ar11092 on august 1.
211 angstrom energy calculations with 30 loops
ERG_FREE = 1.3389243e+32
ERG_FREE_UNC = 4.2251313e+31
ERG_TOT = 5.1831448e+32
ERG_TOT_UNC = 9.1464023e+31
171 angstrom energy calculations with 50 loops
ERG_FREE = 1.7747051e+32
ERG_FREE_UNC = 9.6831359e+31
ERG_TOT = 6.1052848e+32
ERG_TOT_UNC = 1.4735266e+32
193 angstrom energay calculations with 35 loops
ERG_FREE = 1.0359039e+32
ERG_FREE_UNC = 3.8393463e+31
ERG_TOT = 5.6236923e+32
ERG_TOT_UNC = 1.1393917e+32
335 angstrom energy calculations with 15 loops
ERG_FREE = 1.8752290e+32
ERG_FREE_UNC = 8.7563738e+31
ERG_TOT = 5.2983089e+32
ERG_TOT_UNC = 1.2035350e+32
Here are all 130 loops that went into the calculation for ar11092

2012_07_02
So the majority of last week was spent trying to get the magnetic field ready to go. We finally got it going and now I am running the 171 wavelength on the active region 11092. Things are moving along nicely and as I get the calculations for the energies I will post them here.
193 angstrom

211 angstrom

2012_06_26
Lucas and I spent a fair amount of the morning, getting the AIA images aligned so I could start drawing lines. We finally got it going and I started drawing lines, I am about up to 50. I already started the Lothian_browning program because it takes forever, so we'll see how that region comes. After 171 I'll move on and draw lines on the other wavelengths
The rest of the day was spent finishing up the calculation for the ar11158 region, and here is is
ERG_FREE = 2.09205e+32
ERG_TOT = 8.82972e+32
This is all the loops that I have, so this wavelength is done. Here is the picture with all the loops,

2012_06_25
I started the program for the other active region that we want to do, ar11092. so the alpha fields are getting generated, that should be done today, so I will start on the loops tomorrow, and will start in the 171 angstrom wavelength to draw loops.
This is the most recent computation of the ar11158 region, we have about 105 loops gone through, and this is the energies that we get.
ERG_FREE FLOAT = 1.21176e+32
ERG_TOT FLOAT = 7.71093e+32
and here is the picture of the area with the loops and footpoints on,

2012_06_22
So I finished the calculation with 46 loops and for the energy I got
ERG_FREE = 1.42730e+32
ERG_TOT = 5.06741e+32
and here is the picture of the loops and footpoints, it is cool

2012_06_21
So I calculated the energy for this active region, but only with 10 loops, The picture looks great,

The vallues that I got for the energy are
erg_free = 1.36248e+31
erg_tot = 9.23385e+31
I will put in the other loops that I have calculated, the program has finished about 30 curves, so I'll let it run overnight agian(for the second night!!), and then I will input the next curves and get a better result.
2012_06_20
Today I spent most of the time getting the next part of the program to work, we were using a two programs, one called lothian_n_browning_visual_fit_v2.pro, and the other one was select_minima.pro We had to modify each due to the in_model/outfolder problem, but it was fairly simple so we got both running, I let the lothian program run all night and finished about 20 curves, of 123 so now I will go through and look at the best fit.
This is the most up to date curve picture

2012_06_19
Here we have a picture of the bezier curves that I have drawn for active region 11158 in the 211 angstrom wavelength. There are a lot of curves, and I am probably not done.

2012_06_18
So I finally got some lines drawn we got the programs all set, and now I am drawing lines on ar11158 in the 211 angdtrom wavelength, we finally got it done at the end of the day, so I'll get some pictures tomorrow.
2012_06_13
The day was spent cheking the code making sure the conversions were correct, and also trying to get the pictures on the webpage a little better quality. I figured out how to do it, so from now on they should look a little better, here is an example,

2012_06_12
here is a better image of the program and the lines that have been drawn,

this is an output of the program that gives the length of the lines in pixals length is; 591.320
length is; 1004.40
length is; 444.758
length is; 603.810
These are the lengths in pixals so they look big. I could eventually work on labeling the lines to make more sense, but that is a task for another day.

So this is the final one,and here are the lengths of the loops in Megameters,
length is; 191.942
length is; 104.666
length is; 192.763
length is; 275.656
I think that I have the program down, I can set the number of loops that I draw, and the number of points on each line, also I sized the axis so that it is in Mega Meters, here is a link to my code as of right now,
2012_06_12_tracer_j.pro
I am sure that it is not the most optimal code, but it gets the job done
2012_06_08
Today I met with Dr. Longcope, and he gave me an assignment for IDL. I was to print out a .fits file and make a program that would trace a line that followed the clicks of my mouse pointer and then superimpose the drawn loops on the sun's loops I was able to get it done and also able to calculate the length of the drawn loops. I am still playing around with saving it to a file with the loops drawn, so that will be for next week. Also I need to figure out a way to make the x and y axis in arcsecs.

you can't see the lines very well but they are there, I'll play with it some more tuesday, there are four of them drawn on there.
2012_06_07
This week has been spent learning IDL, the logic is more or less the same throughout the programming languages, it is just learning the sintax that is taking the time.
2012_06_02 through 06
We made it to Montana safe and sound and hungry, it was a beautiful drive up and we had a good trip. The first day was good we had a the classic orientation stuff and it made me excited for the REU.
| |