
    ۨhl                     "    S SK rS rS rSS jrg)    Nc                     X -
  nX1-
  nXf-  Xw-  -   nX@-
  U-  XQ-
  U-  -   U-  n	U	S:  a  Sn	OU	S:  a  Sn	X	U-  -   n
XU-  -   nX-
  nX-
  n[         R                  " X-  X-  -   5      nU$ )N   r   )npsqrt)x1y1x2y2x3y3pxpydsquxydxdydists                  G/Users/jklymak/Dropbox/Teaching/OSMCruiseData/CruiseData25/getInletX.pyr   r      s    	B	B
%"%-C7b.BGr>
)S	0A1u	
Q
VA
VA	
B	
B772525=!DK    c                    SS/SS/SS/SS/S	S
/SS/SS/SS/SS/SS/SS/SS/SS/SS/SS/SS/S S!/S"S#/S$S%/S&S'//n[         R                  " U 5      n [         R                  " U5      nS(nS)nS*n[         R                  " [        U5      5      S+[         R                  " [        U5      5      -  -   nS,[         R                  " U5      -  nS,[         R                  " U5      -  n[        U5       H`  u  pU
S-   U-
  [         R                  " U[         R                  -  S.-  5      -  U-  S/-  S+U
S0   U-
  -  U-  S/-  -   Xi'   U
S0   Xy'   U
S-   X'   Mb     [         R                  " U5      S,-  n[         R                  " [         R                  " [         R                  " [         R                  " U5      5      5      5      US0S & [         R                  " S,S1S25      n[         R                  " XU5      n[         R                  " XU5      n[         R                  " X[         R                  " U5      5      n[         R                  " X[         R                  " U5      5      nUS+U-  -   nS3nX-
  [         R                  " U[         R                  -  S.-  5      -  U-  S/-  S+X-
  -  U-  S/-  -   nS-[         R                  " U5      -  nS-[         R                  " U5      -  n[        U5       Hc  u  nn[         R                  " [         R                  " UU-
  5      5      nUU   S4-  U-
  UU'   [         R                  " UUU   -
  5      S4-  UU'   Me     UU4$ )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      s    
y	!)))))))))))))))))))'	C* --
C
--
C	D	DF 
#c(	BrxxC11	1Bbggbk>Dbggbk>D~tDy"&&beeC007;B!T	N6<QRU<UU!!  WWR[^Fwwryy!456F12JBuU+IYYy-FYYy-F	9BGGBK	0B	9BGGBK	0Bbe8D	B
xRUU
3''.s22sx=3G3KKArwwqz\F
lGA,B))BFF29-. )$.3svvb4=01$6 
 7?r   c                    X-
  [         R                  " U[         R                  -  S-  5      -  S-  S-  nX-
  S-  S-  n[         R                  " [         R                  " [         R
                  " U5      S-  [         R
                  " U5      S-  -   5      5      n	[         R                  " S/U	5      n	[         R                  " U5      U-
  [         R                  " U[         R                  -  S-  5      -  S-  S-  n
[         R                  " U5      U-
  S-  S-  n[         R                  " [        U
5      5      n[        [        U
5      5       Ha  n[         R                  n[        [        U5      S-
  5       H3  n[        X   X   XS-      XS-      X   X   5      nUU:  d  M-  UnXU'   M5     Mc     U
S-  n[        [        U
5      5       H  n[        X   5      nUU   nUUS-      nUU   nUUS-      nX   U-
  nX   U-
  nUUU-
  -  UUU-
  -  -   nU[         R                  " UU-
  S-  UU-
  S-  -   5      -  n[        U5        UU	U   -   UU'   M     X   U-
  nU$ )z;
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   s4    /266$,"56	6	<t	CB
/3	t	#Byya"''"+q.!@ABHyy"x(H	s	d	"BFF4"%%<$+>$??#ELA	s	d	"C',A
((3q6
C3q6]&&s2wqy!AbeRUBsGR!WadADABG|A	 "  VF3q6]36{XQZXQZTBYTBYBGnrR"W~-BGGR"WqLBGa<788c
(4.(q	    6)FMr   )r   )numpyr   r   rA   rW    r   r   <module>rZ      s    .;z*r   