
HJc        (   @   s  d  Z  d d k Z d d k Z e i d d d d d d d	 d
 d d d d d d d d d d d d d d d d d d d d d d  d! d" d# d$ d% d& d' g%  Z e i d( d) d* d+ d, d- d. d/ d0 d1 d2 d3 d4 d5 d6 d7 d8 d9 d: d; d< d= d> d? d@ dA dB dC dD dE dF dG dH dI dJ dK dL g%  Z e i d dM dN dO dP dQ dR dS dT dU dV dW dX dY dZ g  Z e i d> d? d@ dA dB dC dD dE dF dG dH dI dJ dK d[ g  Z h  d\ d( <d] d+ <d^ d. <d_ d1 <Z h  d` d( <da d+ <db d. <dc d1 <dd d4 <de d7 <df d: <dg d= <dK dA <dh dD <di dJ <dG dj <Z	 dk   Z
 dl   Z dm   Z dn   Z do   Z dp   Z dq   Z dr   Z ds   Z e dt  Z d S(u   s[   
mctscale.py

3He parameter table

26 Jun 2009 Lev - Greywall superfluid 3He phase diagram
iNg|?5^?g`"?gL7A`?gp=
ף?gh|?5?gS?g(\?gV-?gHzG?gM?gsh|??gT㥛 ?gl?gV-?gh|?5 @gK7 @g+ @g/$@g"~j@gy&1@gx@gA`"@gv/@gMbX@gK7@gʡE@gzG@g/$@g)\(@g}?5^I@gffffff@g/$@gS㥛@gNbX9@gn@gS@g!rh@g        g      ?g       @g      @g      @g      @g      @g      @g       @g      "@g      $@g      &@g      (@g      *@g      ,@g      .@g      0@g      1@g      2@g      3@g      4@g      5@gQ85@g      6@g      7@g      8@g      9@g      :@g      ;@g      <@g      =@g      >@g      ?@g      @@g     @@g      A@g%C+A@gjt@gV-@g#~j@gI+@gZd;O@gjt@gS @gDl @g r @gK7 @g+?g/$?guV?gx?gw-A@gQkB@gQ@@gp=
@@g(\µ>@gffffffK@gH@gfffffF@g33333sE@g     @D@g333333C@gYB@gA@gY@@g     >@gGz.A@c           C   s
   t  t  S(   s&   return the melting pressure at AB line(   t   maxt   __P_AB__(    (    (    s   c:\py\lib\he3.pyt   Pmelt8   s    c          C   sP   t  t j }  t t d d d  t  |  g  t t d d d  t |  g  f S(   sw   
    return a bundle (T,P) of arrays defining the outline of the A phase,
    the melting curve is considered flat
    iN(   t   __Tc__t   T_trit   concatenatet   __T_AB__R   t   __Pc__(   t	   above_tri(    (    s   c:\py\lib\he3.pyt   Aphase_contour9   s    c          C   sK   t  t j }  t t  |  t d d g g  t t |  t t   d g g  f S(   sw   
    return a bundle (T,P) of arrays defining the outline of the B phase,
    the melting curve is considered flat
    i    (   R   R   R   R   R   R   R   (   t	   below_tri(    (    s   c:\py\lib\he3.pyt   Bphase_contourA   s    c         C   s%   t  i |  t t d t  i d t  i S(   s6   
    return a Tc [mK] at a given pressure P [bar]
    t   leftt   right(   t   numpyt   interpR   R   t   nan(   t   P(    (    s   c:\py\lib\he3.pyt   TcI   s    c         C   s%   t  i |  t t d t  i d t  i S(   sE   
    return a P [bar] at which the Tc takes a given value T [mK]
    R   R   (   R   R   R   R   R   (   t   T(    (    s   c:\py\lib\he3.pyt   PcO   s    c         C   s%   t  i |  t t d t  i d t  i S(   s8   
    return a T_AB [mK] at a given pressure P [bar]
    R   R   (   R   R   R   R   R   (   R   (    (    s   c:\py\lib\he3.pyt   T_ABU   s    c         C   s%   t  i |  t t d t  i d t  i S(   sG   
    return a P [bar] at which the T_AB takes a given value T [mK]
    R   R   (   R   R   R   R   R   (   R   (    (    s   c:\py\lib\he3.pyt   P_AB[   s    c         C   su   t  i t i d d d t i d  g  i   \ } } t  i |  } t  i |  | | | | d t  i d t  i S(	   s|   
    return Fermi [m / sec] velocity based on Wheatley RMP
    Interpolation is made linearly in coordinates v_F(P^0.8)
    sN   /home/lev/doc/iop_meeting2009/stuff_from_Rob/Data/CriticalThickness-Fit-GL.datt   usecolsi    i   s'   /data/exp/run20/sf/Tc/FermiVelocity.txtR   R   (   i    i   (   R   R   t   flibt	   loadasciit	   transposet   argsortR   R   (   R   t   PPt   vvt   ii(    (    s   c:\py\lib\he3.pyt   v_Fermia   s
    !c         C   s;   d
 d d d t  |   d t i d d t |   d	 S(   sN   
    return zero temperature coherence length [nm] based on Rob's account
    i   gyX5;?i   g      ?gH=NEy8i   gǇ0;gMbP?g    eAgjMS @(   R   R   t   piR   (   R   (    (    s   c:\py\lib\he3.pyt   xi0r   s    c         C   s1   | p | t  |   :} n t |   d | d S(   s6   
    return Ginzburg-Landau coherence length [nm]
    i   g      ?(   R   R!   (   R   R   t   reducedT(    (    s   c:\py\lib\he3.pyt   xi_GLx   s     (   t   __doc__R   R   t   arrayR   R   R   R   t   __Vmol__t   __v_Fermi__R   R	   R   R   R   R   R   R   R!   t   TrueR#   (    (    (    s   c:\py\lib\he3.pys   <module>   sb   																							