| 1 | ; |
|---|
| 2 | ; @coeffhf |
|---|
| 3 | ; |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | ye, r, t |
|---|
| 7 | |
|---|
| 8 | yeyeye = where(t le 0.3) & r0 = r[yeyeye] & t0 = t[yeyeye] |
|---|
| 9 | |
|---|
| 10 | rerr = replicate(1.e-10,n_elements(r0)) |
|---|
| 11 | start = [mean(r0), max(r0)-min(r0), (max(t0)-min(t0))*mean(r0)] |
|---|
| 12 | expr = 'FONC1(X, P)' |
|---|
| 13 | result = MPFITEXPR(expr,t0,r0,rerr,start) |
|---|
| 14 | |
|---|
| 15 | ;plot, r, t, psym=5, yrange=[0,0.6] & oplot, fonc1(t, result), t |
|---|
| 16 | |
|---|
| 17 | yeyeye = where((t gt 0.3) and (t le 1.0)) & r0 = r[yeyeye] & t0 = t[yeyeye] |
|---|
| 18 | ;ye = LINFIT(t0,r0) & oplot, ye(1) * t + ye(0), t |
|---|
| 19 | |
|---|
| 20 | ;; |
|---|
| 21 | ;; verif |
|---|
| 22 | ;; |
|---|
| 23 | ;lim=0.30 & zelim=(-3.85/alog(lim)+0.07*alog(lim))*exp(-4.61*lim) |
|---|
| 24 | ;up=1. & & zeup= - 1.55 * up + 1.27 |
|---|
| 25 | ;slope = (zelim - zeup) / (lim - up) & ord = zeup - slope * up |
|---|
| 26 | ;t1 = findgen(101.)*lim/100. & param1 = (-3.85/alog(t1)+0.07*alog(t1))*exp(-4.61*t1) |
|---|
| 27 | ;t2 = lim + findgen(101.)*(1.0 - lim)/100. & param2 = slope * t2 + ord |
|---|
| 28 | ;plot, param1, t1, yrange=[0.,1.], xrange=[-0.5,1.] & oplot, param2, t2 & oplot, r, t, psym=5 |
|---|
| 29 | |
|---|
| 30 | |
|---|
