version 1.18, 2021/05/26 04:43:52
|
version 1.19, 2021/05/26 19:25:17
|
|
|
DIE("Problem creating sharp cea records."); | DIE("Problem creating sharp cea records."); |
} | } |
| |
|
|
char *keylist = (char *) params_get_str(&cmdparams, "keylist"); | char *keylist = (char *) params_get_str(&cmdparams, "keylist"); |
char *debug = (char *) params_get_str(&cmdparams, "debug"); | char *debug = (char *) params_get_str(&cmdparams, "debug"); |
| |
|
|
int epszflag = (strstr(keylist,"EPSZ") != NULL); | int epszflag = (strstr(keylist,"EPSZ") != NULL); |
int debugflag = (strstr(debug,"debug") != NULL); | int debugflag = (strstr(debug,"debug") != NULL); |
| |
DRMS_Record_t *sharpinrec = sharpinrecset->records[0]; |
for (irec=0;irec<nrecs;irec++) |
DRMS_Record_t *sharpceainrec = sharpceainrecset->records[0]; |
{ |
|
DRMS_Record_t *sharpinrec = sharpinrecset->records[irec]; |
|
DRMS_Record_t *sharpceainrec = sharpceainrecset->records[irec]; |
DRMS_Segment_t *inseg = drms_segment_lookup(sharpceainrec, "Br"); | DRMS_Segment_t *inseg = drms_segment_lookup(sharpceainrec, "Br"); |
int nx = inseg->axis[0]; | int nx = inseg->axis[0]; |
int ny = inseg->axis[1]; | int ny = inseg->axis[1]; |
int nxny = nx * ny; | int nxny = nx * ny; |
int dims[2] = {nx, ny}; | int dims[2] = {nx, ny}; |
|
|
// Temp arrays | // Temp arrays |
float *bh = (float *) (malloc(nxny * sizeof(float))); | float *bh = (float *) (malloc(nxny * sizeof(float))); |
float *bt = (float *) (malloc(nxny * sizeof(float))); | float *bt = (float *) (malloc(nxny * sizeof(float))); |
|
|
float *derx_los = (float *) (malloc(nxny * sizeof(float))); | float *derx_los = (float *) (malloc(nxny * sizeof(float))); |
float *dery_los = (float *) (malloc(nxny * sizeof(float))); | float *dery_los = (float *) (malloc(nxny * sizeof(float))); |
| |
for (irec=0;irec<nrecs;irec++) |
|
//for (irec=5;irec<nrecs;irec++) |
|
{ |
|
|
|
//DRMS_Record_t *sharpinrec = sharpinrecset->records[irec]; |
|
DRMS_Record_t *sharpceainrec = sharpceainrecset->records[irec]; |
|
DRMS_Segment_t *inseg = drms_segment_lookup(sharpceainrec, "Br"); |
|
|
|
// ephemeris variables | // ephemeris variables |
float cdelt1_orig, cdelt1, dsun_obs, imcrpix1, imcrpix2, crpix1, crpix2; | float cdelt1_orig, cdelt1, dsun_obs, imcrpix1, imcrpix2, crpix1, crpix2; |
double rsun_ref, rsun_obs; | double rsun_ref, rsun_obs; |
|
|
free(pmap); | free(pmap); |
free(p1pad); | free(p1pad); |
free(pmapn); | free(pmapn); |
|
|
} //endfor |
|
|
|
free(fx); free(fy); free(fz); | free(fx); free(fy); free(fz); |
free(bh); free(bt); free(jz); | free(bh); free(bt); free(jz); |
free(bpx); free(bpy); free(bpz); | free(bpx); free(bpy); free(bpz); |
|
|
free(err_term2); | free(err_term2); |
free(err_term1); | free(err_term1); |
| |
|
} //endfor |
| |
// Close all the records | // Close all the records |
drms_close_records(sharpinrecset, DRMS_FREE_RECORD); | drms_close_records(sharpinrecset, DRMS_FREE_RECORD); |