
    fl                          d dl Zd Zd ZddZy)    Nc                     || z
  }||z
  }||z  ||z  z   }|| z
  |z  ||z
  |z  z   |z  }	|	dkD  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/Eos311202409/getInletX.pyr   r      s    	BB	BB
R%"R%-C7b.BGr>
)S	0A1u	
Q
QVA
QVA	
RB	
RB772b52b5=!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 ]b  \  }	}
|
d-   |z
  t        j                  |t         j                  z  d.z        z  |z  d/z  d+|
d0   |z
  z  |z  d/z  z   ||	<   |
d0   ||	<   |
d-   ||	<   d t        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 ]_  \  }}t        j                  t        j                  ||z
              }||   d4z  |z
  ||<   t        j                  |||   z
        d4z  ||<   a ||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      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1!Q!Q WWR[^Fwwryy!456F12JBuU+IYYy-FYYy-F	9VBGGBK	0B	9VBGGBK	0Bbe8D	B
4xRUU
3''.s22s4x=3G3KKArwwqz\F
lGA, 7B))BFF29-. )$.3svvb4=01$67
 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 ]c  }t         j                  }t        t        |      dz
        D ]7  }t        ||   ||   ||dz      ||dz      |
|   ||         }||k  s1|}|||<   9 e |
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/266$,"56	6	<t	CB
T/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beRUBqsGR!WadAaDABG|A		 VF3q6] )3q6{XQZXQZqTBYqTBYBGnrR"W~-BGGR"WqLBGa<788c
(4.(q	)  i 6)FMr   )r   )numpyr   r   rA   rW    r   r   <module>rZ      s    .;z*r   