-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathcodes.html
More file actions
164 lines (126 loc) · 9.25 KB
/
codes.html
File metadata and controls
164 lines (126 loc) · 9.25 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
<!DOCTYPE HTML>
<!--
Spatial by TEMPLATED
templated.co @templatedco
Released for free under the Creative Commons Attribution 3.0 license (templated.co/license)
-->
<html>
<head>
<title>Enrico Camporeale | NOAA Space Weather Prediction Center</title>
<meta charset="utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1" />
<link rel="stylesheet" href="assets/css/main.css" />
</head>
<body class="landing">
<a name="top"></a>
<!-- Header -->
<header id="header" class="alt">
<nav id="nav">
<ul>
<li><a href="index.html">Home</a></li>
<li><a href="research.html">Research</a></li>
<li><a href="mlprojects.html">Machine Learning</a></li>
<li><a href="cv.html">CV</a></li>
<li><a href="codes.html"><font color="red">Software</font></a></li>
<li><a href="publications.html">Publications</a></li>
<li><a href="book.html">Book</a></li>
<li><a href="collaborators.html">Collaborators</a></li>
<li><a href="events.html">Events</a></li>
<li><a href="contact.html">Contacts</a></li>
</ul>
</nav>
</header>
<a href="#menu" class="navPanelToggle"><span class="fa fa-bars"></span></a>
<section id="banner">
<h2>Enrico Camporeale</h2>
</section>
<!-- Three -->
<section id="three" class="wrapper style1">
<div class="container">
<header class="major p">
<h2>Software & Data</h2>
</header>
<h3> Solar wind classification </h3>
We have devised a four-category classification algorithm for the solar wind, based on Gaussian Processes. The four categories are the ones previously adopted in Xu and Borovsky [2015]: ejecta, coronal hole origin plasma, streamer belt origin plasma, and sector reversal origin plasma. The algorithm is trained and tested on a labeled portion of the OMNI2 dataset. </br>
Here, we distribute a database with about 300,000 hours of labeled solar wind data calculated from <a href="ftp://spdf.gsfc.nasa.gov/pub/data/omni/low_res_omni/" target="_blank"> OMNI2</a> for the years 1965-2007.</br>
<b> Please reference the following paper if you publish results based on this database:</b> </br>
E. Camporeale, A. Care', J. Borovsky (2017) <a href="http://onlinelibrary.wiley.com/doi/10.1002/2017JA024383/full" target="_blank"><b>Classification of Solar Wind with Machine Learning</b></a>, <em> J. Geophys. Res.</em>, in press
<p></p>
<ul>
<li>
<a href="./software/OMNI2_classification.dat" > OMNI2_classification.dat </a> (12Mb)</br>
ASCII file with about 300,000 rows and 7 columns. The first three columns are [year, doy, hour]. The last four columns indicate the probability of the event belonging to ejecta, coronal hole origin plasma, streamer belt origin plasma, and sector reversal origin plasma, respectively.
</li>
<li>
<a href="./software/classify_solar_wind.m" > classify_solar_wind.m </a> </br>
MATLAB file with a driver routine to perform Gaussian Process classification on a portion of the OMNI2 dataset. The output is a seven column matrix, as explained above. The routine uses the GPML software, available at <a href="http://www.gaussianprocess.org/gpml/code"> http://www.gaussianprocess.org/gpml/code</a>, and the file below.
</li>
<li>
<a href="./software/parameters_classification.mat" > parameters_classification.mat </a> </br>
This file contains the training set and the optimized hyperparameters to run <a href="./software/classify_solar_wind.m" > classify_solar_wind.m </a>
</li>
<h3> Other codes</h3>
I use a variety of computer codes. I have written some of them myself, while others are part of collaborations. Some are freely available, others are proprietary and not distributable. Here you can find a short description of the codes.
<li><a href="#1">Implicit PIC code</a></li>
<li><a href="#2">Curvilinear PIC code</a></li>
<li><a href="#3">Diffusion code for radiation belt</a></li>
<li><a href="#4">Hot plasma linear dispersion solver</a></li>
<li><a href="#5">Linear kinetic solver for inhomogeneous plasma</a></li>
<li><a href="#6">Mode conversion</a></li>
<li><a href="https://projects.cwi.nl/mlspaceweather/software.html" target="_blank">Software & Data for Machine Learning projects</a></li>
<a name="1"></a>
<b>Implicit PIC code</b><br>
Standard PIC codes use an explicit timestepping scheme that strongly constrains the choice of grid size and time step, to be numerically stable. The implicit PIC code is based on the implicit moment method and, by using a (semi)-implicit time discretization relaxes the stringent stability requirement. Thus one can allow larger grid size, longer time steps, and usually it is possible to resort to physical mass-ratios between species.<br>
I have been involved in the coding of the 2-dimensional PIC code <a href="https://homepages.cwi.nl/~camporea/papers/Markidis_astronum08.pdf" target="_blank"> Parsek2D </a>. More recently a 3D version of the code, named iPIC3D, has been made available on Github. You can download it <a href="https://github.com/CmPA/iPic3D" target="_blank"> here </a>. Contact <a href="https://perswww.kuleuven.be/~u0052182/" target="_blank"> Giovanni Lapenta</a> for more information.
<p></p>
<a name="2"></a>
<b>Curvilinear PIC code</b><br>
I have been involved to the development of a curvilinear PIC code, named CPIC, during my postdoc at the Los Alamos National Laboratory. The code is 3D and electrostatic, and is used for plasma-material interaction studies, allowing fairly arbitrary geometry domains. The details of the code can be found in this <a href="https://homepages.cwi.nl/~camporea/papers/06675865.pdf" target="_blank"> paper </a>.
Unfortunately, the code cannot be openly distributed, as is property of LANL.
<p></p>
<a name="3"></a>
<b>Diffusion code</b><br>
Diffusion codes are extensively used in radiation belt studies. I have coded a 2D diffusion code that was later incorporated as part of the 3D code DREAM3D. I have investigated and compared different schemes for the time discretization, including direct and preconditioned iterative solvers. Details of such study can be found in this <a href="http://onlinelibrary.wiley.com/doi/10.1002/jgra.50293/abstract" target="_blank"> paper</a>. <br>
DREAM3D is LANL property and not distributed, but the 2D part that I wrote can be distributed upon request.
It also includes a routine to calculate the bounce-averaged diffusion coefficients for field-aligned propagating waves.<br>
More recently I have started collaborating with <a href="http://rbm.epss.ucla.edu/person/yuri-shprits/" target="_blank"> Y. Shprits </a> and his group, by using the diffusion code that they have developed in the last several years, VERB, which is also availablae upon request to Dr. Shprits.
<p></p>
<a name="4"></a>
<b>Hot plasma linear dispersion solver</b><br>
The hot plasma linear solver follows the classical derivation of kinetic linear theory from Stix's book. The code has been initially developed by Prof. David Burgess at QMUL (London) and later modified by myself.
In standard mode it uses a Newton-Raphson root finding algorithm, but we have also played with different ideas to overcome the weaknesses of Newton-Raphson (which might actually not converge on a given solution). For instance, we have implemented a bisection method using triangular meshes.
<p></p>
<a name="5"></a>
<b> Mode conversion</b><br>
This code solves a model based on two-fluid linear equations, in 2D. It has been used to study the linear mode conversion of a whistler/lower hybrid wave packet impinging on a density striation. The code is described <a href="http://onlinelibrary.wiley.com/doi/10.1029/2012JA017726/pdf" target="_blank"> here </a>.
<p></p>
<a href="#top">Back to the top</a>
</div>
<!-- Footer -->
<footer id="footer">
<ul class="copyright">
<li>© Enrico Camporeale</li>
<li>Design: <a href="http://templated.co">TEMPLATED</a></li>
<li> The contents of these pages have not been reviewed or approved by CWI.</li></br>
<a rel="license" href="http://creativecommons.org/licenses/by-sa/4.0/"><img alt="Creative Commons License" style="border-width:0" src="https://i.creativecommons.org/l/by-sa/4.0/88x31.png" /></a><br />The databases published here are licensed under a <a rel="license" href="http://creativecommons.org/licenses/by-sa/4.0/">Creative Commons Attribution-ShareAlike 4.0 International License</a>. <br/>
</ul>
</footer>
</div>
</div>
<!-- Default Statcounter code for Homepage
https://ecamporeale.github.io/ -->
<script type="text/javascript">
var sc_project=2920493;
var sc_invisible=1;
var sc_security="0391eb40";
</script>
<script type="text/javascript"
src="https://www.statcounter.com/counter/counter.js"
async></script>
<noscript><div class="statcounter"><a title="Web Analytics"
href="https://statcounter.com/" target="_blank"><img
class="statcounter"
src="https://c.statcounter.com/2920493/0/0391eb40/1/"
alt="Web Analytics"></a></div></noscript>
<!-- End of Statcounter Code -->
</body></html>