|
version 1.30, 2014/06/05 21:27:19
|
version 1.33, 2015/01/25 18:28:41
|
| Line 1049 int findPosition(DRMS_Record_t *inRec, s |
|
| Line 1049 int findPosition(DRMS_Record_t *inRec, s |
|
| float disk_lonc = drms_getkey_float(inRec, "CRLN_OBS", &status); | float disk_lonc = drms_getkey_float(inRec, "CRLN_OBS", &status); |
| | |
| /* Center coord */ | /* Center coord */ |
| |
// Changed into double Jun 16 2014 XS |
| | |
| float minlon = drms_getkey_float(inRec, "LONDTMIN", &status); if (status) return 1; // Stonyhurst lon |
double minlon = drms_getkey_double(inRec, "LONDTMIN", &status); if (status) return 1; // Stonyhurst lon |
| float maxlon = drms_getkey_float(inRec, "LONDTMAX", &status); if (status) return 1; |
double maxlon = drms_getkey_double(inRec, "LONDTMAX", &status); if (status) return 1; |
| float minlat = drms_getkey_float(inRec, "LATDTMIN", &status); if (status) return 1; |
double minlat = drms_getkey_double(inRec, "LATDTMIN", &status); if (status) return 1; |
| float maxlat = drms_getkey_float(inRec, "LATDTMAX", &status); if (status) return 1; |
double maxlat = drms_getkey_double(inRec, "LATDTMAX", &status); if (status) return 1; |
| | |
| // A bug fixer for HARP (per M. Turmon) | // A bug fixer for HARP (per M. Turmon) |
| // When AR is below threshold, "LONDTMIN", "LONDTMAX" will be wrong | // When AR is below threshold, "LONDTMIN", "LONDTMAX" will be wrong |
| Line 1084 int findPosition(DRMS_Record_t *inRec, s |
|
| Line 1085 int findPosition(DRMS_Record_t *inRec, s |
|
| mInfo->yc = (maxlat + minlat) / 2.; | mInfo->yc = (maxlat + minlat) / 2.; |
| | |
| /* Size */ | /* Size */ |
| |
// Rounded to 1.d3 precision first. Jun 16 2014 XS |
| | |
| mInfo->ncol = round((maxlon - minlon) / mInfo->xscale); |
mInfo->ncol = round(round((maxlon - minlon) * 1.e3) / 1.e3 / mInfo->xscale); |
| mInfo->nrow = round((maxlat - minlat) / mInfo->yscale); |
mInfo->nrow = round(round((maxlat - minlat) * 1.e3) / 1.e3 / mInfo->yscale); |
| | |
| return 0; | return 0; |
| | |
| Line 1971 void computeSWIndex(struct swIndex *swKe |
|
| Line 1973 void computeSWIndex(struct swIndex *swKe |
|
| float *jz_err_squared = (float *) (malloc(nxny * sizeof(float))); | float *jz_err_squared = (float *) (malloc(nxny * sizeof(float))); |
| float *jz_err_squared_smooth = (float *) (malloc(nxny * sizeof(float))); | float *jz_err_squared_smooth = (float *) (malloc(nxny * sizeof(float))); |
| float *jz_rms_err = (float *) (malloc(nxny * sizeof(float))); | float *jz_rms_err = (float *) (malloc(nxny * sizeof(float))); |
| |
float *err_term1 = (float *) (malloc(nxny * sizeof(float))); |
| |
float *err_term2 = (float *) (malloc(nxny * sizeof(float))); |
| | |
| // define some values for the R calculation | // define some values for the R calculation |
| int scale = round(2.0/cdelt1); | int scale = round(2.0/cdelt1); |
| Line 2039 void computeSWIndex(struct swIndex *swKe |
|
| Line 2043 void computeSWIndex(struct swIndex *swKe |
|
| } | } |
| | |
| computeJz(bx_err, by_err, bx, by, dims, jz, jz_err, jz_err_squared, mask, bitmask, cdelt1, rsun_ref, rsun_obs, | computeJz(bx_err, by_err, bx, by, dims, jz, jz_err, jz_err_squared, mask, bitmask, cdelt1, rsun_ref, rsun_obs, |
| derx, dery); |
derx, dery, err_term1, err_term2); |
| | |
| | |
| if(computeJzsmooth(bx, by, dims, jz, jz_smooth, jz_err, jz_rms_err, jz_err_squared_smooth, &(swKeys_ptr->mean_jz), | if(computeJzsmooth(bx, by, dims, jz, jz_smooth, jz_err, jz_rms_err, jz_err_squared_smooth, &(swKeys_ptr->mean_jz), |
| Line 2140 void computeSWIndex(struct swIndex *swKe |
|
| Line 2144 void computeSWIndex(struct swIndex *swKe |
|
| free(bt_err); free(bh_err); free(jz_err); | free(bt_err); free(bh_err); free(jz_err); |
| free(jz_err_squared); free(jz_rms_err); | free(jz_err_squared); free(jz_rms_err); |
| free(jz_err_squared_smooth); | free(jz_err_squared_smooth); |
| |
free(err_term2); |
| |
free(err_term1); |
| | |
| // free the arrays that are related to the r calculation | // free the arrays that are related to the r calculation |
| free(rim); | free(rim); |
Return to
sharp.c
CVS log
Up to [Development] / JSOC / proj / sharp / apps