|
version 1.23, 2014/03/05 19:39:32
|
version 1.24, 2014/03/07 21:15:54
|
| Line 1432 int getBErr(float *bx_err, float *by_err |
|
| Line 1432 int getBErr(float *bx_err, float *by_err |
|
| if (sphere2img (lat, lon, disk_latc, disk_lonc, &xi, &zeta, | if (sphere2img (lat, lon, disk_latc, disk_lonc, &xi, &zeta, |
| disk_xc, disk_yc, 1.0, pa, 0., 0., 0., 0.)) { | disk_xc, disk_yc, 1.0, pa, 0., 0., 0., 0.)) { |
| bx_err[ind_map] = DRMS_MISSING_FLOAT; | bx_err[ind_map] = DRMS_MISSING_FLOAT; |
| bx_err[ind_map] = DRMS_MISSING_FLOAT; |
by_err[ind_map] = DRMS_MISSING_FLOAT; |
| bx_err[ind_map] = DRMS_MISSING_FLOAT; |
bz_err[ind_map] = DRMS_MISSING_FLOAT; // Mar 7 |
| continue; | continue; |
| } | } |
| | |
| Line 2060 void computeSWIndex(struct swIndex *swKe |
|
| Line 2060 void computeSWIndex(struct swIndex *swKe |
|
| swKeys_ptr->totpot_err = DRMS_MISSING_FLOAT; | swKeys_ptr->totpot_err = DRMS_MISSING_FLOAT; |
| } | } |
| | |
| |
|
| if (computeShearAngle(bx_err, by_err, bz_err, bx, by, bz, bpx, bpy, bpz, dims, | if (computeShearAngle(bx_err, by_err, bz_err, bx, by, bz, bpx, bpy, bpz, dims, |
| &(swKeys_ptr->meanshear_angle), &(swKeys_ptr->meanshear_angle_err), &(swKeys_ptr->area_w_shear_gt_45), | &(swKeys_ptr->meanshear_angle), &(swKeys_ptr->meanshear_angle_err), &(swKeys_ptr->area_w_shear_gt_45), |
| mask, bitmask)) { | mask, bitmask)) { |
| Line 2068 void computeSWIndex(struct swIndex *swKe |
|
| Line 2069 void computeSWIndex(struct swIndex *swKe |
|
| swKeys_ptr->meanshear_angle_err= DRMS_MISSING_FLOAT; | swKeys_ptr->meanshear_angle_err= DRMS_MISSING_FLOAT; |
| } | } |
| | |
| |
/* |
| if (computeR(bz_err, los , dims, &(swKeys_ptr->Rparam), cdelt1, rim, p1p0, p1n0, | if (computeR(bz_err, los , dims, &(swKeys_ptr->Rparam), cdelt1, rim, p1p0, p1n0, |
| p1p, p1n, p1, pmap, nx1, ny1)) | p1p, p1n, p1, pmap, nx1, ny1)) |
| { | { |
| swKeys_ptr->Rparam = DRMS_MISSING_FLOAT; // If fail, fill in NaN | swKeys_ptr->Rparam = DRMS_MISSING_FLOAT; // If fail, fill in NaN |
| } | } |
| |
*/ |
| | |
| // Clean up the arrays | // Clean up the arrays |
| | |
| Line 2082 void computeSWIndex(struct swIndex *swKe |
|
| Line 2084 void computeSWIndex(struct swIndex *swKe |
|
| drms_free_array(bxArray); | drms_free_array(bxArray); |
| drms_free_array(byArray); | drms_free_array(byArray); |
| drms_free_array(bzArray); | drms_free_array(bzArray); |
| |
drms_free_array(losArray); // Mar 7 |
| |
drms_free_array(bx_errArray); |
| |
drms_free_array(by_errArray); |
| |
drms_free_array(bz_errArray); |
| | |
| free(bh); free(bt); free(jz); free(jz_smooth); | free(bh); free(bt); free(jz); free(jz_smooth); |
| free(bpx); free(bpy); free(bpz); | free(bpx); free(bpy); free(bpz); |
Return to
sharp.c
CVS log
Up to [Development] / JSOC / proj / sharp / apps