-
Notifications
You must be signed in to change notification settings - Fork 1
/
fitj.pro
42 lines (35 loc) · 947 Bytes
/
fitj.pro
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
;+
;
; Name: fitj
;
; Purpose: function to fit HI elongation tracks using Sheeleys formula
;
; Parameters: X
;
; Keywords: -
;
;
; Called by fitelongation.pro
;
;
; History: June 2009/update March 2012
;
; Author: Christian Moestl
; Space Research Institute, Austrian Academy of Sciences
; SSL UC Berkeley
;-
FUNCTION fitj, X
common myfit,xueber,yueber, dst
degtorad=!dpi/180;
;after Sheeley 2008 ApJ
delta=90*degtorad-X[0]; delta wird v. Plane of Sky weg gemessen
rho=X[1]*(xueber-X[2])/dst;
fit=atan(rho*cos(delta)/(1-rho*sin(delta)))/degtorad
sizey=size(yueber); how many data points
residue=double(0); define residue variable
for i=0,sizey(1)-1 do begin
;sum up squared difference between observations and fit for each datapoint
residue=residue+( abs(yueber(i))-abs(fit(i)) )^2;
endfor
return, residue
END