JSOC/proj/sharp/apps/sharp.c - diff

Return to sharp.c CVS log

Up to [Development] / JSOC / proj / sharp / apps

version 1.30, 2014/06/05 21:27:19

version 1.31, 2014/06/16 04:45:18

Line 1049 int findPosition(DRMS_Record_t *inRec, s

float disk_lonc = drms_getkey_float(inRec, "CRLN_OBS", &status);

/* 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)

// 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.;

/* 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;

Legend: