Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- FUNCTION maketrace_defringe, prm, lm, idx, period, fitname, noval
- COMPILE_OPT idl2, hidden
- ;Rebuild trace
- nwa = DBLARR(N_ELEMENTS(idx), N_ELEMENTS(period))
- cnt = prm[*,1]
- nx = (SIZE(prm, /DIM))[0]
- nprm = (SIZE(prm, /DIM))[1]
- FOR xi = 0, nx-1 DO IF (cnt[xi] NE noval) THEN BREAK
- SWITCH fitname OF
- 'GFCnb':
- 'GVCnb': BEGIN
- FOR i = xi, nx-1 DO BEGIN
- IF (cnt[i] EQ noval) THEN CONTINUE
- l1 = lm[0,i]
- l2 = lm[1,i]
- gauss = prm[i,0] * EXP(-0.5*((period-prm[i,1])/prm[i,2])^2)
- nwa[i,l1:l2] = gauss[l1:l2]
- ENDFOR
- BREAK
- END
- 'GFC':
- 'GVC':
- 'GFCest': BEGIN
- FOR i = xi, nx-1 DO BEGIN
- IF (cnt[i] EQ noval) THEN CONTINUE
- l1 = lm[0,i]
- l2 = lm[1,i]
- gauss = prm[i,0] * EXP(-0.5*((period-prm[i,1])/prm[i,2])^2) + $
- prm[i,3]
- nwa[i,l1:l2] = gauss[l1:l2]
- ENDFOR
- BREAK
- END
- ELSE: message, "Method '" + fitname + $
- "' not enabled in maketrace_defringe"
- ENDSWITCH
- RETURN, nwa
- END
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement