Changeset 2455 in lmdz_wrf

Timestamp:

Apr 22, 2019, 12:12:03 AM (6 years ago)

Author:

lfita

Message:

Adding:

• `zsailing_boat': Function to define an schematic sailing boat from the z-plane with sails
• `z_boat': Function to define an schematic boat from the z-plane

File:

• trunk/tools/geometry_tools.py (modified) (5 diffs)

trunk/tools/geometry_tools.py

@@ -15 +15 @@
 import matplotlib.pyplot as plt
 import os
+import generic_tools as gen
 
 errormsg = 'ERROR -- error -- ERROR -- error'
@@ -47 +48 @@
 # p_triangle: Function to provide the polygon of a triangle from its 3 vertices
 # surface_sphere: Function to provide an sphere as matrix of x,y,z coordinates
+# zsailing_boat: Function to define an schematic sailing boat from the z-plane with sails
+# z_boat: Function to define an schematic boat from the z-plane
 
 ## Plotting
@@ -622 +625 @@
 
 # FROM: http://www.photographers1.com/Sailing/NauticalTerms&Nomenclature.html
-def zsailing_boat(length=10., beam=1., lbeam=0.4, sternbp=0.5):
+def zboat(length=10., beam=1., lbeam=0.4, sternbp=0.5):
     """ Function to define an schematic boat from the z-plane
       length: length of the boat (without stern, default 10)
@@ -629 +632 @@
       sternbp: beam at stern (as percentage of beam, default 0.5)
     """
-    fname = 'zsailing_boat'
+    fname = 'zboat'
 
     bow = np.array([length, 0.])
@@ -671 +674 @@
  
     return boat
+
+def zsailing_boat(length=10., beam=1., lbeam=0.4, sternbp=0.5, lmast=0.6, wmast=0.1, \
+  sd=0., lheads=0.38, lmains=0.55):
+    """ Function to define an schematic sailing boat from the z-plane with sails
+      length: length of the boat (without stern, default 10)
+      beam: beam of the boat (default 1)
+      lbeam: length at beam (as percentage of length, default 0.4)
+      sternbp: beam at stern (as percentage of beam, default 0.5)
+      lmast: position of the mast (as percentage of length, default 0.6)
+      wmast: width of the mast (default 0.1)
+      sd: sails direction respect to center line (default 0.)
+      lheads: length of head sail (as percentage of legnth, defaul 0.38)
+      lmains: length of main sail (as percentage of legnth, defaul 0.55)
+    """
+    import numpy.ma as ma
+    fname = 'zsailing_boat'
+
+    bow = np.array([length, 0.])
+    maxportside = np.array([length*lbeam, -beam])
+    maxstarboardside = np.array([length*lbeam, beam])
+    portside = np.array([0., -beam*sternbp])
+    starboardside = np.array([0., beam*sternbp])
+
+    # forward section
+    fportsaid = circ_sec(bow,maxportside, length*2)
+    fstarboardsaid = circ_sec(maxstarboardside, bow, length*2)
+    dpts = fportsaid.shape[0]
+
+    # aft section
+    aportsaid = circ_sec(maxportside, portside, length*2)
+    astarboardsaid = circ_sec(starboardside, maxstarboardside, length*2)
+    # stern
+    stern = circ_sec(portside, starboardside, length*2)
+    # mast
+    mast = p_circle(wmast,N=dpts)
+    mast = mast + [length*lmast, 0.]
+    # head sails
+    lsail = lheads*length
+    sailsa = np.pi/2. + sd*(1.+0.05)*180./np.pi
+    endsail = np.array([lsail*np.sin(sailsa), lsail*np.cos(sailsa)])
+    endsail[0] = length - endsail[0]
+    hsail = circ_sec(bow, endsail, lsail*1.15)
+    # main sails
+    lsail = lmains*length
+    endsail = np.array([lsail*np.sin(sailsa), lsail*np.cos(sailsa)])
+    endsail[0] = length*lmast - endsail[0]
+    msail = circ_sec(np.array([length*lmast,0.]), endsail, lsail*1.15)
+
+    sailingboat = np.zeros((dpts*8+3,2), dtype=np.float)
+
+    sailingboat[0:dpts,:] = fportsaid
+    sailingboat[dpts:2*dpts,:] = aportsaid
+    sailingboat[2*dpts:3*dpts,:] = stern
+    sailingboat[3*dpts:4*dpts,:] = astarboardsaid
+    sailingboat[4*dpts:5*dpts,:] = fstarboardsaid
+    sailingboat[5*dpts,:] = [gen.fillValueF, gen.fillValueF]
+    sailingboat[5*dpts+1:6*dpts+1,:] = mast
+    sailingboat[6*dpts+1,:] = [gen.fillValueF, gen.fillValueF]
+    sailingboat[6*dpts+2:7*dpts+2,:] = hsail
+    sailingboat[7*dpts+2,:] = [gen.fillValueF, gen.fillValueF]
+    sailingboat[7*dpts+2:8*dpts+2,:] = msail
+
+    sailingboat = ma.masked_equal(sailingboat, gen.fillValueF)
+
+    fname = 'sailboat_L' + str(int(length*100.)) + '_B' + str(int(beam*100.)) +      \
+      '_lb' + str(int(lbeam*100.)) + '_sb' + str(int(sternbp*100.)) +                \
+      '_lm' + str(int(lmast*100.)) + '_wm' + str(int(wmast)) +                       \
+      '_sd' + str(int(sd)) + '_hs' + str(int(lheads*100.)) +                         \
+      '_ms' + str(int(lheads*100.)) + '.dat'
+    if not os.path.isfile(fname):
+        print infmsg
+        print '  ' + fname + ": writting boat coordinates file '" + fname + "' !!"
+        of = open(fname, 'w')
+        of.write('# boat file with Length: ' + str(length) +' max_beam: '+str(beam)+ \
+          'length_at_max_beam:' + str(lbeam) + '% beam at stern: ' + str(sternbp)+   \
+          ' % mast position: '+ str(lmast) + ' % mast width: ' + str(wmast) + ' ' +  \
+          ' sails direction:' + str(sd) + ' head sail length: ' + str(lheads) + ' %'+\
+          ' main sail length' + str(lmains) + '\n')
+        for ip in range(dpts*5):
+            of.write(str(sailingboat[ip,0]) + ' ' + str(sailingboat[ip,1]) + '\n')
+        
+        of.close()
+        print fname + ": Successfull written '" + fname + "' !!"
+ 
+    return sailingboat
 
 def write_join_poly(polys, flname='join_polygons.dat'):