Changeset 2636 in lmdz_wrf

Timestamp:

Jun 25, 2019, 5:01:01 AM (5 years ago)

Author:

lfita

Message:

Adding:

• `EstuarioRioPlata?': Function to plot an eschematic representation of the Estuario of Rio de la Plata
• `NotablePostitions?': Notable points

File:

• trunk/tools/nautical.py (modified) (3 diffs)

trunk/tools/nautical.py

@@ -19 +19 @@
 # buoy1: Function to draw a buoy as superposition of prism and section of ball
 # band_lighthouse: Function to plot a lighthouse with spiral bands
+# EstuarioRioPlata: Function to plot an eschematic representation of the Estuario of Rio de la Plata
 # green_buoy1: Function to draw a green mark buoy using buoy1
 # isolateddanger_buoy1: Function to draw an isolated danger buoy using buoy1
@@ -33 +34 @@
 #   circular segments
 # [north/east/south/west_buoy1: Function to draw a [North/East/South/West] danger buoy using buoy1
+
+# Definitions [Name, lat, lon]
+NotablePostitions = {                                                                \
+  'PuntaMedanos': ['Punta Medanos', np.array([-36.8494, -56.6395])],                 \
+  'PuntaRaza': ['Punta Raza', np.array([-36.2929, -56.7474])],                       \
+  'RioSalado': ['Rio Salado', np.array([-35.7423, -57.3635])],                       \
+  'PuntaIndio': ['Punt Indio', np.array([-35.4179, -57.0959])],                      \
+  'PuntaAtalaya': ['Punta Atalaya', np.array([-35.01868, -57.5181])],                \
+  'Tigre': ['Tigre', np.array([-34.4486, -58.4989])],                                \
+  'MartinChico': ['Martin Chico', np.array([-34.1681, -58.2118])],                   \
+  'Colonia': ['Colonia del Sacramento', np.array([-34.4724, -57.8556])],             \
+  'ArroyoRosario': ['Arroyo Rosario', np.array([-34.4331, -57.3504])],               \
+  'Montevideo': ['Montevideo', np.array([-34.9216, -56.1574])],                      \
+  'PuntaEste': ['Punta del Este', np.array([-34.9830, -54.9533])],                   \
+  'CaboPolonio': ['Cabo Polonio', np.array([-34.4083, -53.7782])]}
 
 # FROM: http://www.photographers1.com/Sailing/NauticalTerms&Nomenclature.html
@@ -1174 +1190 @@
     return buoy, buoysecs, buoydic
 
+def EstuarioRioPlata(N=300):
+    """ Function to plot an eschematic representation of the Estuario of Rio de la Plata
+      N: total number of vertices to use
+    """
+    fname = 'EstuarioRioPlata'
+
+    secs0 = ['PuntaMedanos', 'PuntaRaza', 'RioSalado', 'PuntaIndio', 'PuntaAtalaya', \
+      'Tigre', 'MartinChico', 'Colonia', 'ArroyoRosario', 'Montevideo', 'PuntaEste', \
+      'CaboPolonio']
+    secs = []
+    dic = {}
+    rads = [5., 1.0, 5., 5., 5., 5., 5., 5., 5., 5., 5.]
+    lengths = ['short', 'short', 'short', 'short', 'short', 'short', 'short',        \
+      'short', 'short', 'short', 'short']
+    sides = ['right', 'left', 'left', 'right', 'left', 'left', 'left', 'left',       \
+      'right', 'left', 'right']
+    Nsecs = len(secs0)
+    Nn = N/Nsecs
+    estuario = np.zeros((N,2), dtype=np.float)
+
+    iip = 0
+    # Atlantic_PuntaRaza
+    prevn = 'PuntaMedanos'
+    pv = NotablePostitions[prevn]
+    ip = pv[1]
+    for isec in range(1,Nsecs-1):
+        iisec = isec - 1
+        aname = secs0[isec]
+        pv = NotablePostitions[aname]
+        ep = pv[1]
+        dps = geo.dist_points(ip,ep)
+        estuario[iip:iip+Nn,:] = geo.circ_sec(ip,ep, dps*rads[iisec], lengths[iisec],\
+          sides[iisec], Nn)
+        secs.append(prevn+'_'+aname)
+        dic[prevn+'_'+aname] = [estuario[iip:iip+Nn,:], ['-', 'k', 1.]]
+        ip = ep + 0.
+        prevn = aname + ''
+        iip = iip + Nn
+
+    Nn2 =  N - (Nsecs-2)*Nn
+    isec = Nsecs-1
+    iisec = isec - 1
+    aname = secs0[isec]
+    pv = NotablePostitions[aname]
+    ep = pv[1]
+    dps = geo.dist_points(ip,ep)
+    isec = Nsecs - 1
+    estuario[iip:N,:] = geo.circ_sec(ip, ep, dps*rads[iisec], lengths[iisec],        \
+      sides[iisec], Nn2)
+    secs.append(prevn+'_'+aname)
+    dic[prevn+'_'+aname] = [estuario[iip:N,:], ['-', 'k', 1.]]
+
+    return estuario, secs, dic
+