
    fl                     "    d dl Zd Zd ZddZdS )    Nc                     || z
  }||z
  }||z  ||z  z   }|| z
  |z  ||z
  |z  z   |z  }	|	dk    rd}	n|	dk     rd}	| |	|z  z   }
||	|z  z   }|
|z
  }||z
  }t          j        ||z  ||z  z             }|S )N   r   )npsqrt)x1y1x2y2x3y3pxpydsquxydxdydists                  G/Users/jklymak/Dropbox/Teaching/CruiseData314/CruiseData24/getInletX.pyr   r      s    	BB	BB
R%"R%-C7b.BGr>
)S	0A1uu	
Q
QVA
QVA	
RB	
RB72b52b5=!!DK    c                 ~   ddgddgddgddgd	d
gddgddgddgddgddgddgddgddgddgddgddgd d!gd"d#gd$d%gd&d'gg}t          j        |           } t          j        |          }d(}d)}d*}t          j        t          |                    d+t          j        t          |                    z  z   }d,t          j        |          z  }d,t          j        |          z  }t          |          D ]i\  }	}
|
d-         |z
  t          j        |t           j        z  d.z            z  |z  d/z  d+|
d0         |z
  z  |z  d/z  z   ||	<   |
d0         ||	<   |
d-         ||	<   jt          j        |          d,z  }t          j        t          j        t          j	        t          j
        |                                        |d0d <   t          j        d,d1d2          }t          j        |||          }t          j        |||          }t          j        ||t          j        |                    }t          j        ||t          j        |                    }|d+|z  z   }d3}| |z
  t          j        |t           j        z  d.z            z  |z  d/z  d+||z
  z  |z  d/z  z   }d-t          j        |          z  }d-t          j        |          z  }t          |          D ]b\  }}t          j        t          j	        ||z
                      }||         d4z  |z
  ||<   t          j	        |||         z
            d4z  ||<   c||fS )5Ng 4^g[ CH@g^gCͩDH@gzj^gO9FH@gI,)w^gVGH@goc#^g4)HH@g\<^gR臭JH@g^grQH@g_^g .VH@gƉvYH@gQy9^gpi[H@ge^g*ʥ]H@g"1^gτ&%_H@gjO9^gI2`H@gV^gTb.bH@gc^gQUaH@gݰmQ^gC6^H@gN^gl%[H@g˛õ^g<.[H@gKƼ^g>;XH@gI^gTJH@g     ^g33333SH@i<  y              ?        r            N@r   g     L@i'  goC.@g     @@)r   
atleast_1dzeroslenreal	enumeratecospicumsumabsdifflinspaceinterpimagargmin)lonlatdatlon0lat0mpernmxxlattlonnndalongx	alongxnewlatnewlonnewxrxixnewx0r   xalongxacrossindclosests                           r   	getInletXrA      sU   
y	!)))))))))))))))))))'	C* -

C
-

C	D	DF 
#c((		BrxC111	1Bbgbkk>Dbgbkk>D~~  !tDy"&beC0007;B!T	N6<QRU<UU1!Q!QWR[[^Fwry!4!45566F122JBuU++IYy--FYy--F	9VBGBKK	0	0B	9VBGBKK	0	0Bbe8D	B
4xRU
3'''.s22s4x=3G3KKArwqzz\F

lGA,, 7 7B)BF29--.. )$.3svb4=011$67?r   c                 &   | |z
  t          j        |t           j        z  dz            z  dz  dz  }||z
  dz  dz  }t          j        t          j        t          j        |          dz  t          j        |          dz  z                       }	t          j        dg|	          }	t          j        |          |z
  t          j        |t           j        z  dz            z  dz  dz  }
t          j        |          |z
  dz  dz  }t          j        t          |
                    }t          t          |
                    D ]{}t           j        }t          t          |          dz
            D ]M}t          ||         ||         ||dz            ||dz            |
|         ||                   }||k     r|}|||<   N||
dz  }t          t          |
                    D ]}t          ||                   }||         }||dz            }||         }||dz            }|
|         |z
  }||         |z
  }|||z
  z  |||z
  z  z   }|t          j        ||z
  dz  ||z
  dz  z             z  }t          |           ||	|         z   ||<   |	|         |z
  }|S )zC
    get potison along an inlet line   Anchor Ind should be S4
    g     f@r   g?   r   r   )r   r!   r"   r#   r   r%   appendr   r   r   rangeInfr   intprint)linelonslinelatsr*   r+   r-   r.   	anchorindr0   yydistliner   r   r?   jthedistiddr5   inddr<   r   y0r   xpypdots                             r   getInletXNewrW   W   s    T/26$,"566	6	<t	CB
T/3	t	#Bya"'"++q.!@AABBHy"x((H	s		d	"BF4"%<$+>$?$??#ELA	s		d	"C',A
(3q66

C3q66]]  &s2wwqy!! 	 	AbeRUBqsGR!WadAaDAABG||A		
 VF3q66]] ) )3q6{{XQZXQZqTBYqTBYBGnrR"W~-BGR"WqLBGa<7888c


(4.(q		i 6)FMr   )r   )numpyr   r   rA   rW    r   r   <module>rZ      sN         .; ; ;z* * * * * *r   