|
version 1.27, 2014/03/19 20:46:26
|
version 1.28, 2014/05/16 21:56:11
|
| Line 1936 void computeSWIndex(struct swIndex *swKe |
|
| Line 1936 void computeSWIndex(struct swIndex *swKe |
|
| | |
| // convert cdelt1_orig from degrees to arcsec | // convert cdelt1_orig from degrees to arcsec |
| float cdelt1 = (atan((rsun_ref*cdelt1_orig*RADSINDEG)/(dsun_obs)))*(1/RADSINDEG)*(3600.); | float cdelt1 = (atan((rsun_ref*cdelt1_orig*RADSINDEG)/(dsun_obs)))*(1/RADSINDEG)*(3600.); |
| int nx1 = nx*cdelt1/2; |
|
| int ny1 = ny*cdelt1/2; |
// define some values for the R calculation |
| |
int scale = round(2.0/cdelt1); |
| |
int nx1 = nx/scale; |
| |
int ny1 = ny/scale; |
| |
|
| |
if (nx1 > floor((nx-1)/scale + 1) ) |
| |
DIE("X-dimension of output array in fsample() is too large."); |
| |
if (ny1 > floor((ny-1)/scale + 1) ) |
| |
DIE("Y-dimension of output array in fsample() is too large."); |
| | |
| // Temp arrays | // Temp arrays |
| float *bh = (float *) (malloc(nxny * sizeof(float))); | float *bh = (float *) (malloc(nxny * sizeof(float))); |
Return to
sharp.c
CVS log
Up to [Development] / JSOC / proj / sharp / apps