ó
‚™xTc           @   sð   d  d l  m Z d „  Z d d „ Z e d k rì d  d l Z d  d l Z e	 e j
 ƒ d k rn e e j ƒ qì e	 e j
 ƒ d k rÆ e e j
 d d ƒ Z e e e j j e j
 d ƒ ƒ e j ƒ  qì e d	 e j j e j
 d
 ƒ f ‚ n  d S(   iÿÿÿÿ(   t   divisionc         C   sl   d d  l  } |  d k rV | j |  ƒ } | |  k  rM | | j |  | ƒ } qb | } n d } d } | | f S(   Niÿÿÿÿi    (   t   matht   sqrt(   t   iobst   isigmaR   t   fobst   fsigma(    (    sM   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/command_line/dtrek.to_cns.pyt   iobs_as_fobs   s    	c         C   s€  g  } g  } d } xü |  j  ƒ  D]î } | d k r5 Pn  | d k r‚ | j d ƒ rd | j | d  ƒ n  | j d ƒ sy q n  d } n  | j ƒ  } t | d ƒ } t | d ƒ } | d k  s | d k  rÌ q n  t | | ƒ \ }	 }
 | j d	 | d | d | d
 |	 |
 f ƒ q W| r'd G| Gd GHn  x | D] } d G| Gd GHq.Wd t | ƒ f GHd GHd GHd GHx | D] } | GHqmWd  S(   Ni    t    t   CRYSTAL_iÿÿÿÿt    i   i   i   s#   INDEX %s %s %s FOBS %.6g SIGMA %.6gi   s   { file:t   }t   {s   NREFlections=%ds   ANOMalous=FALSes2   DECLare NAME=FOBS  DOMAin=RECIprocal TYPE=REAL ENDs2   DECLare NAME=SIGMA DOMAin=RECIprocal TYPE=REAL END(   t
   xreadlinest
   startswitht   appendt   splitt   floatR   t   len(   t   file_objectt	   file_namet   infot   reflectionst   modet   linet   fldsR   R   R   t   sigma(    (    sM   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/command_line/dtrek.to_cns.pyt   dtrek_as_cns_hkl   s<      	 	'   t   __main__Ni   i   t   rs   usage: %s [d*trek_file_name]i    (   t
   __future__R    R   t   NoneR   t   __name__t   syst   os.patht   osR   t   argvt   stdint   opent   ft   patht   abspatht   closet   RuntimeErrort   basename(    (    (    sM   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/command_line/dtrek.to_cns.pyt   <module>   s   	 