updates

calc_E.py

@@ -26,7 +26,7 @@ import itertools
 # %% Get the data
 
 if argv[0]=='':
-     date_0 = datetime(2016,5,1,6)
+     date_0 = datetime(2016,4,29,6)
      trapath = "../wp21/era5/traj/2016/"
 
 else:
@@ -43,7 +43,7 @@ z_coord_type = "p"
 
 n = 1
 
-path = "./dat_variablekp/" + date_0.strftime("dat_%Y%m%d_%H") + "/" + str(n) + "/"
+path = "./" + date_0.strftime("dat_%Y%m%d_%H") + "/" + str(n) + "/"
 if not exists(path):
     makedirs(path)
 
@@ -62,8 +62,8 @@ T = trajs['T'][::n]
 H = cc.ptoh(P, 1013.25, 8.435) # standard atmosphere
 
 k_p = cc.calc_k(U/1000, V/1000, Omg/100)
-# k_p=15
-k_h = cc.calc_k(U, V, cc.omg2w(Omg/100, T, P))
+k_p=15
+# k_h = cc.calc_k(U, V, cc.omg2w(Omg/100, T, P))
 
 D = list((list(), list(), list()))
 

main.py

@@ -30,16 +30,16 @@ from warnings import catch_warnings, simplefilter
 # only every nth trajectory
 n = 1
 
-force_calc = True
+force_calc = False
 
-date_0 = datetime(2016,5,4,12)
+date_0 = datetime(2016,4,30,12)
 fname = "../wp21/era5/traj/" + date_0.strftime('%Y/traj_%Y%m%d_%H') + ".npy"
 
 B = dd.get_block_dat(date_0)
 coasts = dd.get_coasts(0, 90, -180, 180)
 
 #fname = "traj" + date_0.strftime('_%Y%m%d_%H') + ".npy"
-path = "./dat_variablekp/" + date_0.strftime("dat_%Y%m%d_%H") + "/" + str(n) + "/"
+path = "./" + date_0.strftime("dat_%Y%m%d_%H") + "/" + str(n) + "/"
 if not exists(path):
     makedirs(path)
 
@@ -286,7 +286,7 @@ if True:
 
     boundscs, hullscs = dd.io_bounds(
                 path + "stereop", "$\\alpha$", xcs, ycs, zcs,
-                alpha, force_calc)
+                alpha, False)
 
     bound_meth = "$\\alpha$"