;
; WRITTEN: 09Dec2002  Xuepu Zhao
;
function fte2pv,fte,model 

   absfte=abs(fte)
   CASE 1 OF
    model eq 'wang90': pv = 415 * 0.938385^absfte + 316
    model eq 'wang902': begin
                        nd=n_elements(fte)
                        pv=fte
                        for i=0,nd-1 do begin
                            ftei=absfte(i)
                            if finite(ftei) eq 1 then begin
                               CASE 1 OF
                                ftei le 3.5:                pv(i)=700
                                ftei gt 3.5 and ftei le 9:  pv(i)=600
                                ftei gt 9   and ftei le 18: pv(i)=500
                                ftei gt 18  and ftei le 54: pv(i)=400
                                ftei gt 54:                 pv(i)=330
                                else:
                               ENDCASE
                            endif else pv(i)=!values.f_nan
                        endfor
                        end       
    model eq 'zhao03' : begin
                        nd=n_elements(fte)
                        pv=fte
                        for i=0,nd-1 do begin
                            ftei=absfte(i)
                            if finite(ftei) eq 1 then begin
                               CASE 1 OF
                                ftei le 3.5:                pv(i)=700
                                ftei gt 3.5 and ftei le 12:  pv(i)=600
                                ftei gt 12   and ftei le 16: pv(i)=500
                                ftei gt 16  and ftei le 54: pv(i)=400
                                ftei gt 54:                 pv(i)=330
                                else:
                               ENDCASE
                            endif else pv(i)=!values.f_nan
                        endfor
                        end
    model eq 'riley' : begin
                        nd=n_elements(fte)
                        pv=fte
                        for i=0,nd-1 do begin
                            ftei=absfte(i)
                            if finite(ftei) eq 1 then begin
                               CASE 1 OF
                                ftei le 3.5:                pv(i)=650
                            ;   ftei gt 3.5 and ftei le 12:  pv(i)=625
                                ftei gt 3.5 and ftei le 10:  pv(i)=625
                                ftei gt 10 and ftei le 12:  pv(i)=600
                                ftei gt 12   and ftei le 30: pv(i)=300
                            ;   ftei gt 12   and ftei le 16: pv(i)=600
                            ;   ftei gt 16  and ftei le 54: pv(i)=300
                                ftei gt 30:                 pv(i)=265
                                else:
                               ENDCASE
                            endif else pv(i)=!values.f_nan
                        endfor
                        end
    model eq 'wang97' : pv = 639 * 0.877432^absfte + 400
    model eq 'bala00' : pv = 267 + 410.0/absfte^(1.0/2.5)
    model eq 'arg00'  : pv = 285 + 650.0/absfte^(5./9.)
    model eq 'zhao942': begin
                        xx=alog10(absfte)
                        pv = 868.2 - 264.7*xx + 9.4*xx*xx
                        end
    model eq 'zhao943': begin
                        xx=alog(absfte)
                        pv = 868.1 - 114.7*xx + 1.7*xx*xx
                        end
    model eq 'zhao951': pv = 683 * 0.971158^absfte + 147
    else:
   ENDCASE
   return,pv

end