|
version 1.40, 2021/05/24 22:17:06
|
version 1.43, 2021/05/26 04:45:48
|
| Line 1293 int computeAbsFlux_los(float *los, int * |
|
| Line 1293 int computeAbsFlux_los(float *los, int * |
|
| *mean_vf_los_ptr = sum*cdelt1*cdelt1*(rsun_ref/rsun_obs)*(rsun_ref/rsun_obs)*100.0*100.0; | *mean_vf_los_ptr = sum*cdelt1*cdelt1*(rsun_ref/rsun_obs)*(rsun_ref/rsun_obs)*100.0*100.0; |
| *count_mask_los_ptr = count_mask_los; | *count_mask_los_ptr = count_mask_los; |
| | |
| printf("USFLUXL=%f\n",*mean_vf_los_ptr); |
//printf("USFLUXL=%f\n",*mean_vf_los_ptr); |
| printf("CMASKL=%f\n",*count_mask_los_ptr); |
//printf("CMASKL=%f\n",*count_mask_los_ptr); |
| | |
| return 0; | return 0; |
| } | } |
| Line 1366 int computeLOSderivative(float *los, int |
|
| Line 1366 int computeLOSderivative(float *los, int |
|
| for (j = 0; j <= ny-1; j++) | for (j = 0; j <= ny-1; j++) |
| { | { |
| if ( bitmask[j * nx + i] < 30 ) continue; | if ( bitmask[j * nx + i] < 30 ) continue; |
| if ( (derx_los[j * nx + i] + dery_los[j * nx + i]) == 0) continue; |
|
| if isnan(los[j * nx + i]) continue; | if isnan(los[j * nx + i]) continue; |
| if isnan(los[(j+1) * nx + i]) continue; | if isnan(los[(j+1) * nx + i]) continue; |
| if isnan(los[(j-1) * nx + i]) continue; | if isnan(los[(j-1) * nx + i]) continue; |
| Line 1381 int computeLOSderivative(float *los, int |
|
| Line 1380 int computeLOSderivative(float *los, int |
|
| | |
| *mean_derivative_los_ptr = (sum)/(count_mask); // would be divided by ((nx-2)*(ny-2)) if shape of count_mask = shape of magnetogram | *mean_derivative_los_ptr = (sum)/(count_mask); // would be divided by ((nx-2)*(ny-2)) if shape of count_mask = shape of magnetogram |
| | |
| printf("MEANGBL=%f\n",*mean_derivative_los_ptr); |
//printf("MEANGBL=%f\n",*mean_derivative_los_ptr); |
| | |
| return 0; | return 0; |
| } | } |
Return to
sw_functions.c
CVS log
Up to [Development] / JSOC / proj / sharp / apps