PDE [message #93847] |
Sat, 05 November 2016 02:19  |
Ali Gamal
Messages: 98 Registered: June 2013
|
Member |
|
|
Hi,
How can I write this equation in idl
∑_i〖V_i*〖(∂I/∂λ)〗_i 〗
|
|
|
|
|
Re: PDE [message #93850 is a reply to message #93847] |
Sat, 05 November 2016 06:58   |
Ali Gamal
Messages: 98 Registered: June 2013
|
Member |
|
|
On Saturday, November 5, 2016 at 11:19:57 AM UTC+2, AGW wrote:
> Hi,
> How can I write this equation in idl
>
> ∑_i〖V_i*〖(∂I/∂λ)〗_i 〗
where
vint=[4,18]
λ=findgen(vint(1)-vint(0)+1)-(vint(1)-vint(0))/2.
for i=0,nx-1 do begin
for j=0,ny-1 do begin
SI=stokesout[i,j,vint(0):vint(1),0]
SV=stokesout[i,j,vint(0):vint(1),3]
|
|
|
Re: PDE [message #93853 is a reply to message #93850] |
Mon, 07 November 2016 04:29  |
Markus Schmassmann
Messages: 129 Registered: April 2016
|
Senior Member |
|
|
On 11/05/2016 02:58 PM, AGW wrote:
> On Saturday, November 5, 2016 at 11:19:57 AM UTC+2, AGW wrote:
>> Hi,
>> How can I write this equation in idl
>>
>> ∑_i〖V_i*〖(∂I/∂λ)〗_i 〗
>
> where
> vint=[4,18]
> λ=findgen(vint(1)-vint(0)+1)-(vint(1)-vint(0))/2.
>
> for i=0,nx-1 do begin
> for j=0,ny-1 do begin
>
> SI=stokesout[i,j,vint(0):vint(1),0]
> SV=stokesout[i,j,vint(0):vint(1),3]
IdL=(stokesout[*,*,vint[0]+1:vint(1)+1,0]-stokesout[*,*,vint [0]+1:vint(1)+1,0])/2.
out=total(sv*IdL,1)
if the wavelength positions form [3:19] are not uniform, your have to
divide out an array instead of 2.
I hope this is what you are looking for, Markus
|
|
|