version 1.1, 2018/03/28 01:27:06
|
version 1.2, 2020/05/04 21:11:07
|
|
|
MEANGBZ Mean value of the vertical field gradient, in Gauss/Mm | MEANGBZ Mean value of the vertical field gradient, in Gauss/Mm |
R_VALUE Karel Schrijver's R parameter | R_VALUE Karel Schrijver's R parameter |
| |
The indices are calculated on the pixels in which the conf_disambig segment is greater than 70 and |
The indices are calculated on the pixels in which the bitmap segment is greater than 36. |
pixels in which the bitmap segment is greater than 42. These ranges are selected because the CCD |
|
coordinate bitmaps are interpolated for certain data (at the time of this CVS submit, all data |
The SMARP bitmap has 13 unique values because they describe three different characteristics: |
prior to 2013.08.21_17:24:00_TAI contain interpolated bitmaps; data post-2013.08.21_17:24:00_TAI |
the location of the pixel (whether the pixel is off disk or a member of the patch), the |
contain nearest-neighbor bitmaps). |
character of the magnetic field (quiet or active), and the character of the continuum |
|
intensity data (quiet, part of faculae, part of a sunspot). |
In the CCD coordinates, this means that we are selecting the pixels that that equal 33 or 34 in bitmap. Here are the definitions of the pixel values: |
|
|
Here are all the possible values: |
|
|
For bitmap: |
Value Keyword Definition |
1 : weak field outside smooth bounding curve |
0 OFFDISK The pixel is located somewhere off the solar disk. |
2 : strong field outside smooth bounding curve |
1 QUIET The pixel is associated with weak line-of-sight magnetic field. |
33 : weak field inside smooth bounding curve |
2 ACTIVE The pixel is associated with strong line-of-sight magnetic field. |
34 : strong field inside smooth bounding curve |
4 SPTQUIET The pixel is associated with no features in the continuum intensity data. |
|
8 SPTFACUL The pixel is associated with faculae in the continuum intensity data. |
Written by Monica Bobra 15 August 2012 |
16 SPTSPOT The pixel is associated with sunspots in the continuum intensity data. |
Potential Field code (appended) written by Keiji Hayashi |
32 ON_PATCH The pixel is inside the patch. |
Error analysis modification 21 October 2013 |
|
|
Here are all the possible permutations: |
|
|
|
Value Definition |
|
0 The pixel is located somewhere off the solar disk. |
|
5 The pixel is located outside the patch, associated with weak line-of-sight magnetic field, and no features in the continuum intensity data. |
|
9 The pixel is located outside the patch, associated with weak line-of-sight magnetic field, and faculae in the continuum intensity data. |
|
17 The pixel is located outside the patch, associated with weak line-of-sight magnetic field, and sunspots in the continuum intensity data. |
|
6 The pixel is located outside the patch, associated with strong line-of-sight magnetic field, and no features in the continuum intensity data. |
|
10 The pixel is located outside the patch, associated with strong line-of-sight magnetic field, and faculae in the continuum intensity data. |
|
18 The pixel is located outside the patch, associated with strong line-of-sight magnetic field, and sunspots in the continuum intensity data. |
|
37 The pixel is located inside the patch, associated with weak line-of-sight magnetic field, and no features in the continuum intensity data. |
|
41 The pixel is located inside the patch, associated with weak line-of-sight magnetic field, and faculae in the continuum intensity data. |
|
49 The pixel is located inside the patch, associated with weak line-of-sight magnetic field, and sunspots in the continuum intensity data. |
|
38 The pixel is located inside the patch, associated with strong line-of-sight magnetic field, and no features in the continuum intensity data. |
|
42 The pixel is located inside the patch, associated with strong line-of-sight magnetic field, and faculae in the continuum intensity data. |
|
50 The pixel is located inside the patch, associated with strong line-of-sight magnetic field, and sunspots in the continuum intensity data. |
|
|
|
Thus pixels with a value greater than 36 are located inside the patch. |
|
|
|
Written by Monica Bobra |
| |
===========================================*/ | ===========================================*/ |
#include <math.h> | #include <math.h> |
Line 70 int computeAbsFlux(float *bz, int *dims, |
|
Line 90 int computeAbsFlux(float *bz, int *dims, |
|
{ | { |
for (j = 0; j < ny; j++) | for (j = 0; j < ny; j++) |
{ | { |
if ( bitmask[j * nx + i] < 42 ) continue; |
if ( bitmask[j * nx + i] < 36 ) continue; |
if isnan(bz[j * nx + i]) continue; | if isnan(bz[j * nx + i]) continue; |
sum += (fabs(bz[j * nx + i])); | sum += (fabs(bz[j * nx + i])); |
count_mask++; | count_mask++; |
Line 154 int computeBzderivative(float *bz, int * |
|
Line 174 int computeBzderivative(float *bz, int * |
|
{ | { |
for (j = 0; j <= ny-1; j++) | for (j = 0; j <= ny-1; j++) |
{ | { |
if ( bitmask[j * nx + i] < 42 ) continue; |
if ( bitmask[j * nx + i] < 36 ) continue; |
if ( (derx_bz[j * nx + i] + dery_bz[j * nx + i]) == 0) continue; | if ( (derx_bz[j * nx + i] + dery_bz[j * nx + i]) == 0) continue; |
if isnan(bz[j * nx + i]) continue; | if isnan(bz[j * nx + i]) continue; |
if isnan(bz[(j+1) * nx + i]) continue; | if isnan(bz[(j+1) * nx + i]) continue; |