Changeset 2687 in lmdz_wrf for trunk

Timestamp:

Aug 10, 2019, 11:56:56 AM (5 years ago)

Author:

lfita

Message:

Adding:

• `lonlat_dist': Function to provide the distance in m between two lon,lat coordinates
• `lonlat_Polyarea': Function to provide the area in m of an enclosed polygon in lon,lat coordinates

File:

• trunk/tools/generic_tools.py (modified) (4 diffs)

trunk/tools/generic_tools.py

@@ -16 +16 @@
 from cStringIO import StringIO
 import sys
-import module_ForGen as gen
+import module_ForDef as fdef
+import module_ForGen as fgen
+import module_ForSci as fsci
 
 main = 'nc_var_tools.py'
@@ -168 +170 @@
 # list_norepeatcombos: Function to all possible (Num-1)-combinations of a Num values without repetitions
 # lstring_values: Function to provide a new list-string from a string which is a list of word separated by a character if some values are repeated they are not included
+# lonlat_dist: Function to provide the distance in m between two lon,lat coordinates
+# lonlat_Polyarea: Function to provide the area in m of an enclosed polygon in lon,lat coordinates
 # make_colors: Function to create a palette of colors from a series of high and low values
 # mat_dimreshape: Function to reshape a matrix on a new one according to values along their dimensions
@@ -847 +851 @@
     #sec4 = 3.*((I+4900.+(J-14.)/12.)/100.)/4.
     #juliand = int(sec1 + sec2 + sec3 - sec4)
-    juliand = gen.module_generic.jd(I,J,K)
+    juliand = fgen.module_generic.jd(I,J,K)
 
     juliand = juliand*1. + (hr + mi/60. + ss/3600.)/24.
@@ -16908 +16912 @@
     return finalcol
 
+def lonlat_dist(lonlat1, lonlat2):
+    """ Function to provide the distance in m between two lon,lat coordinates
+        as ds = sqrt(dlon**2+dlat**2)*Earthradii*np.cos(latmean)*DegRad
+      lonlat1: list with longitude, latitude of point 1
+      lonlat2: list with longitude, latitude of point 2
+    >>> lonlat_dist([-70.073759, -38.773790], [-68.414071, -37.284390])
+    2.2299947572458443, 195290.18809
+    """
+    fname = 'lonlat_dist'
+
+    dlon = lonlat2[0] - lonlat1[0]
+    dlat = lonlat2[1] - lonlat1[1]
+
+    latmean = np.mean([lonlat1[1], lonlat2[1]])
+
+    ds = np.sqrt(dlon**2+dlat**2)
+    deg2km = fdef.module_definitions.earthradii*np.cos(latmean*np.pi/180.)*np.pi/180.
+
+    distkm = ds*deg2km
+
+    return ds, distkm
+
+def lonlat_Polyarea(poly):
+    """ Function to provide the area in m of an enclosed polygon in lon,lat 
+          coordinates
+        as ds = ShoelaceArea*Earthradii*np.cos(latmean)*DegRad
+      polygon: list with the consecutive [longitude, latitude] verrtices of the 
+        polygon
+    >>> lonlat_Polyarea([[-70.073759,-38.773790], [-69.993907,-37.264203], 
+      [-68.414071,-37.284390], [-68.437593,-38.867590]])
+    -2.4892578125, -19081159952.039444
+    """
+    fname = 'lonlat_Polyarea'
+
+    apoly = np.array(poly)
+
+    area = fsci.module_scientific.shoelace_area_polygon(nvertex=4, poly=apoly)
+
+    latmean = np.mean(apoly[:,1])
+    deg2km = fdef.module_definitions.earthradii*np.cos(latmean*np.pi/180.)*np.pi/180.
+
+    areakm = area*deg2km*deg2km
+
+    return area, areakm
+
 #quit()