From c5b7309c6e07c559700486b425be58652b3f33af Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Mon, 11 Oct 2021 13:40:32 -0500
Subject: [PATCH] start support for Boeing icing data...
---
modules/icing/pireps.py | 14 ++++++++++++++
1 file changed, 14 insertions(+)
diff --git a/modules/icing/pireps.py b/modules/icing/pireps.py
index 07c77f71..3a17d79f 100644
--- a/modules/icing/pireps.py
+++ b/modules/icing/pireps.py
@@ -1,6 +1,8 @@
import datetime
from datetime import timezone
import re
+from pathlib import Path
+import h5py
NO_ICE = '\s*NEG\s*|\s*NONE\s*|\s*NEGATIVE\s*|\s*NO\s*'
ICE_LVL = '\d+-\d+|FL\d+-FL\d+'
@@ -155,3 +157,15 @@ def pirep_icing(filename, lon_range=[-180, 180], lat_range=[-63, 63]):
rpts.append(tup)
return ice_dict, no_ice_dict, neg_ice_dict
+
+
+def pirep_icing_from_boeing(dirname='/ships19/cloud/scratch/ICING/BOEING/icing/2020/', lon_range=[-180, 180], lat_range=[-63, 63]):
+ flist = []
+ for path in Path(dirname).rglob('*.nc'):
+ flist.append(path)
+
+ cnt = 0
+ for fname in flist:
+ h5f = h5py.File(fname, 'r')
+ ice_bool = h5f['ice_bool'][:]
+ cnt += ice_bool.shape[0]
\ No newline at end of file
--
GitLab