version 1.5, 2012/10/23 19:35:37
|
version 1.6, 2012/12/18 22:17:25
|
|
|
* v0.0 Jul 02 2012 | * v0.0 Jul 02 2012 |
* v0.1 Jul 23 2012 | * v0.1 Jul 23 2012 |
* v0.2 Sep 04 2012 | * v0.2 Sep 04 2012 |
|
* v0.3 Dec 18 2012 |
* | * |
* Notes: | * Notes: |
* v0.0 | * v0.0 |
|
|
* SW indices fixed | * SW indices fixed |
* Added doppler and continuum | * Added doppler and continuum |
* Added other keywords: HEADER (populated by cvs build version), DATE_B | * Added other keywords: HEADER (populated by cvs build version), DATE_B |
|
* v0.3 |
|
* Fixed memory leakage of 0.15G per rec; denoted with "Dec 18" |
* | * |
* Example: | * Example: |
* sharp "mharp=hmi.Mharp_720s[1404][2012.02.20_10:00]" \ | * sharp "mharp=hmi.Mharp_720s[1404][2012.02.20_10:00]" \ |
Line 294 char *CEASegs[] = {"magnetogram", "bitma |
|
Line 297 char *CEASegs[] = {"magnetogram", "bitma |
|
| |
| |
char *module_name = "sharp"; | char *module_name = "sharp"; |
char *version_id = "2012 Jul 02"; /* Version number */ |
char *version_id = "2012 Dec 18"; /* Version number */ |
| |
ModuleArgs_t module_args[] = | ModuleArgs_t module_args[] = |
{ | { |
|
|
| |
} // irec | } // irec |
| |
|
|
drms_close_records(mharpRS, DRMS_FREE_RECORD); | drms_close_records(mharpRS, DRMS_FREE_RECORD); |
drms_close_records(bharpRS, DRMS_FREE_RECORD); | drms_close_records(bharpRS, DRMS_FREE_RECORD); |
|
drms_close_records(dopRS, DRMS_FREE_RECORD); // Dec 18 2012 |
|
drms_close_records(contRS, DRMS_FREE_RECORD); // Dec 18 2012 |
| |
return 0; | return 0; |
| |
Line 748 int mapScaler(DRMS_Record_t *sharpRec, D |
|
Line 754 int mapScaler(DRMS_Record_t *sharpRec, D |
|
if (status) return 0; | if (status) return 0; |
| |
if (inArray) drms_free_array(inArray); | if (inArray) drms_free_array(inArray); |
|
if ((xsz != FOURK) || (ysz != FOURK)) free(inData); // Dec 18 2012 |
if (outArray) drms_free_array(outArray); | if (outArray) drms_free_array(outArray); |
return 0; | return 0; |
| |
Line 1107 int performSampling(float *outData, floa |
|
Line 1114 int performSampling(float *outData, floa |
|
// Rebinning, smoothing | // Rebinning, smoothing |
| |
frebin(outData0, outData, ncol0, nrow0, mInfo->nbin, 1); // Gaussian | frebin(outData0, outData, ncol0, nrow0, mInfo->nbin, 1); // Gaussian |
|
free(outData0); // Dec 18 2012 |
| |
// | // |
| |
Line 1589 int writeCutout(DRMS_Record_t *outRec, D |
|
Line 1597 int writeCutout(DRMS_Record_t *outRec, D |
|
| |
#if DISAMB_AZI | #if DISAMB_AZI |
if (!strcmp(SegName, "azimuth")) { | if (!strcmp(SegName, "azimuth")) { |
DRMS_Segment_t *disambSeg = drms_segment_lookup(inRec, "disambig"); |
DRMS_Segment_t *disambSeg = NULL; |
|
disambSeg = drms_segment_lookup(inRec, "disambig"); |
if (!disambSeg) {drms_free_array(cutoutArray); return 1;} | if (!disambSeg) {drms_free_array(cutoutArray); return 1;} |
DRMS_Array_t *disambArray; | DRMS_Array_t *disambArray; |
if (disambSeg->axis[0] == nx && disambSeg->axis[1] == ny) { | if (disambSeg->axis[0] == nx && disambSeg->axis[1] == ny) { |
disambArray = drms_segment_read(disambSeg, DRMS_TYPE_CHAR, &status); | disambArray = drms_segment_read(disambSeg, DRMS_TYPE_CHAR, &status); |
if (status) {drms_free_array(cutoutArray); return 1;} | if (status) {drms_free_array(cutoutArray); return 1;} |
} else { | } else { |
|
drms_free_array(cutoutArray); |
return 1; | return 1; |
} | } |
double *azimuth = (double *) cutoutArray->data; | double *azimuth = (double *) cutoutArray->data; |
Line 1774 void computeSWIndex(struct swIndex *swKe |
|
Line 1784 void computeSWIndex(struct swIndex *swKe |
|
| |
// Clean up | // Clean up |
| |
|
drms_free_array(bitmaskArray); // Dec 18 2012 Xudong |
drms_free_array(maskArray); | drms_free_array(maskArray); |
drms_free_array(bxArray); | drms_free_array(bxArray); |
drms_free_array(byArray); | drms_free_array(byArray); |