xan/awips2
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions
awips2/cave/com.raytheon.viz.avnconfig/localization/aviation/python/ClimLib.py
Max Schenkelberg 6f60751ec6 Issue #2033 moved avnfps and text workstation files into respective plugins. 
Change-Id: If95cb839ad81ca2a842ff7f6926847ac3928d8f2

Former-commit-id: 77e1a4d8f5237e5fae930c1e00589c752f8b3738
2013-08-15 12:21:43 -05:00

346 lines
12 KiB
Python
Raw Blame History

##
# This software was developed and / or modified by Raytheon Company,
# pursuant to Contract DG133W-05-CQ-1067 with the US Government.
#
# U.S. EXPORT CONTROLLED TECHNICAL DATA
# This software product contains export-restricted data whose
# export/transfer/disclosure is restricted by U.S. law. Dissemination
# to non-U.S. persons whether in the United States or abroad requires
# an export license or other authorization.
#
# Contractor Name: Raytheon Company
# Contractor Address: 6825 Pine Street, Suite 340
# Mail Stop B8
# Omaha, NE 68106
# 402.291.0100
#
# See the AWIPS II Master Rights File ("Master Rights File.pdf") for
# further licensing information.
##
#
# Name:
# ClimLib.py
# GFS1-NHD:A9007.0000-SCRIPT;3
#
# Status:
# DELIVERED
#
# History:
# Revision 3 (DELIVERED)
# Created: 24-MAY-2006 09:05:37 TROJAN
# spr 7144: ficed weather element dictionary, added VLIFR
# category
#
# Revision 2 (DELIVERED)
# Created: 16-FEB-2006 14:35:36 TROJAN
# added check for cache age
#
# Revision 1 (APPROVED)
# Created: 30-JAN-2006 07:51:28 TROJAN
# stdr 945
#
# Change Document History:
# 1:
# Change Document: GFS1-NHD_SPR_7144
# Action Date: 14-FEB-2007 12:27:26
# Relationship Type: In Response to
# Status: CLOSED
# Title: AvnFPS: Incorrect method to specify filtering criteria in Cig/Vis Monthly tool
#
#
# ClimLib.py
# common code for Climate tools
# George Trojan, SAIC/MDL, November 2005
# last update: 05/24/06
import logging, os, cPickle
import Avn
_Logger = logging.getLogger(__name__)
PType = Avn.Bunch(freezing=1<<0, frozen=1<<1, liquid=1<<2)
PBest = Avn.Bunch(drizzle=1<<0, continuous=1<<1, shower=1<<2)
PInt = Avn.Bunch(lgt=1<<0, mod=1<<1, hvy=1<<2)
PTypeStr = { \
(PType.freezing, PBest.drizzle): 'FZDZ', \
(PType.freezing, PBest.continuous): 'FZRA', \
(PType.freezing, PBest.shower): 'SHPL', \
(PType.frozen, PBest.drizzle): 'DZSN', \
(PType.frozen, PBest.continuous): 'SN', \
(PType.frozen, PBest.shower): 'SHSN', \
(PType.liquid, PBest.drizzle): 'DZ', \
(PType.liquid, PBest.continuous): 'RA', \
(PType.liquid, PBest.shower): 'SHRA', \
}
# Codes used in pres_wxm_code
MWxTable = { \
0: {'str': ''}, \
1: {'str': ''}, \
2: {'str': ''}, \
3: {'str': ''}, \
4: {'str': 'FU'}, \
5: {'str': 'HZ'}, \
6: {'str': 'DU'}, \
7: {'str': 'DU'}, \
8: {'str': 'DU'}, \
9: {'str': 'BLSA'}, \
10: {'str': 'BR'}, \
11: {'str': 'BCFG'}, \
12: {'str': 'MIFG'}, \
13: {'str': 'TS'}, \
14: {'str': ''}, \
15: {'str': ''}, \
16: {'str': ''}, \
17: {'str': 'TS'}, \
18: {'str': ''}, \
19: {'str': ''}, \
20: {'str': ''}, \
21: {'str': ''}, \
22: {'str': ''}, \
23: {'str': ''}, \
24: {'str': ''}, \
25: {'str': ''}, \
26: {'str': ''}, \
27: {'str': ''}, \
28: {'str': ''}, \
29: {'str': ''}, \
30: {'str': 'BLSA'}, \
31: {'str': 'BLSA'}, \
32: {'str': 'BLSA'}, \
33: {'str': 'BLSA'}, \
34: {'str': 'BLSA'}, \
35: {'str': 'BLSA'}, \
36: {'str': 'DRSN'}, \
37: {'str': 'DRSN'}, \
38: {'str': 'BLSN'}, \
39: {'str': 'BLSN'}, \
40: {'str': ''}, \
41: {'str': 'BCFG'}, \
42: {'str': 'FG'}, \
43: {'str': 'FG'}, \
44: {'str': 'FG'}, \
45: {'str': 'FG'}, \
46: {'str': 'FG'}, \
47: {'str': 'FG'}, \
48: {'str': 'FG'}, \
49: {'str': 'FG'}, \
50: {'int': PInt.lgt, 'str': 'DZ', \
'best': PBest.drizzle, 'type': PType.liquid}, \
51: {'int': PInt.lgt, 'str': 'DZ', \
'best': PBest.drizzle, 'type': PType.liquid}, \
52: {'int': PInt.mod, 'str': 'DZ', \
'best': PBest.drizzle, 'type': PType.liquid}, \
53: {'int': PInt.mod, 'str': 'DZ', \
'best': PBest.drizzle, 'type': PType.liquid}, \
54: {'int': PInt.hvy, 'str': 'DZ', \
'best': PBest.drizzle, 'type': PType.liquid}, \
55: {'int': PInt.hvy, 'str': 'DZ', \
'best': PBest.drizzle, 'type': PType.liquid}, \
56: {'int': PInt.lgt, 'str': 'FZDZ', \
'best': PBest.drizzle, 'type': PType.freezing}, \
57: {'int': PInt.mod|PInt.hvy, 'str': 'FZDZ', \
'best': PBest.drizzle, 'type': PType.freezing}, \
58: {'int': PInt.lgt, 'str': 'DZRA', \
'best': PBest.drizzle, 'type': PType.liquid}, \
59: {'int': PInt.mod|PInt.hvy, 'str': 'DZRA', \
'best': PBest.drizzle, 'type': PType.liquid}, \
60: {'int': PInt.lgt, 'str': 'RA', \
'best': PBest.continuous, 'type': PType.liquid}, \
61: {'int': PInt.lgt, 'str': 'RA', \
'best': PBest.continuous, 'type': PType.liquid}, \
62: {'int': PInt.mod, 'str': 'RA', \
'best': PBest.continuous, 'type': PType.liquid}, \
63: {'int': PInt.mod, 'str': 'RA', \
'best': PBest.continuous, 'type': PType.liquid}, \
64: {'int': PInt.hvy, 'str': 'RA', \
'best': PBest.continuous, 'type': PType.liquid}, \
65: {'int': PInt.hvy, 'str': 'RA', \
'best': PBest.continuous, 'type': PType.liquid}, \
66: {'int': PInt.lgt, 'str': 'FZRA', \
'best': PBest.continuous, 'type': PType.freezing}, \
67: {'int': PInt.mod|PInt.hvy, 'str': 'FZRA', \
'best': PBest.continuous, 'type': PType.freezing}, \
68: {'int': PInt.lgt, 'str': 'RASN', \
'best': PBest.continuous, 'type': PType.liquid}, \
69: {'int': PInt.mod|PInt.hvy, 'str': 'RASN', \
'best': PBest.continuous, 'type': PType.liquid}, \
70: {'int': PInt.lgt, 'str': 'SN', \
'best': PBest.continuous, 'type': PType.frozen}, \
71: {'int': PInt.lgt, 'str': 'SN', \
'best': PBest.continuous, 'type': PType.frozen}, \
72: {'int': PInt.mod, 'str': 'SN', \
'best': PBest.continuous, 'type': PType.frozen}, \
73: {'int': PInt.mod, 'str': 'SN', \
'best': PBest.continuous, 'type': PType.frozen}, \
74: {'int': PInt.hvy, 'str': 'SN', \
'best': PBest.continuous, 'type': PType.frozen}, \
75: {'int': PInt.hvy, 'str': 'SN', \
'best': PBest.continuous, 'type': PType.frozen}, \
76: {'int': PInt.lgt, 'str': 'IC', \
'best': PBest.continuous, 'type': PType.frozen}, \
77: {'int': PInt.lgt, 'str': 'SG', \
'best': PBest.continuous, 'type': PType.frozen}, \
78: {'int': PInt.lgt, 'str': 'IC', \
'best': PBest.continuous, 'type': PType.frozen}, \
79: {'int': PInt.lgt, 'str': 'PL', \
'best': PBest.continuous, 'type': PType.freezing}, \
80: {'int': PInt.lgt, 'str': 'SHRA', \
'best': PBest.shower, 'type': PType.liquid}, \
81: {'int': PInt.mod|PInt.hvy, 'str': 'SHRA', \
'best': PBest.shower, 'type': PType.liquid}, \
82: {'int': PInt.hvy, 'str': 'SHRA', \
'best': PBest.shower, 'type': PType.liquid}, \
83: {'int': PInt.lgt, 'str': 'SHRASN', \
'best': PBest.shower, 'type': PType.liquid}, \
84: {'int': PInt.mod|PInt.hvy, 'str': 'SHRASN', \
'best': PBest.shower, 'type': PType.liquid}, \
85: {'int': PInt.lgt, 'str': 'SHSN', \
'best': PBest.shower, 'type': PType.frozen}, \
86: {'int': PInt.mod|PInt.hvy, 'str': 'SHSN', \
'best': PBest.shower, 'type': PType.frozen}, \
87: {'int': PInt.lgt, 'str': 'SHGS', \
'best': PBest.shower, 'type': PType.freezing}, \
88: {'int': PInt.mod|PInt.hvy, 'str': 'SHGS', \
'best': PBest.shower, 'type': PType.freezing}, \
89: {'int': PInt.lgt, 'str': 'SHGS', \
'best': PBest.shower, 'type': PType.freezing}, \
90: {'int': PInt.mod|PInt.hvy, 'str': 'SHGS', \
'best': PBest.shower, 'type': PType.freezing}, \
91: {'int': PInt.lgt, 'str': 'TSRA', \
'best': PBest.shower, 'type': PType.liquid}, \
92: {'int': PInt.mod|PInt.hvy, 'str': 'TSRA', \
'best': PBest.shower, 'type': PType.liquid}, \
93: {'int': PInt.lgt, 'str': 'TSGS', \
'best': PBest.shower, 'type': PType.freezing}, \
94: {'int': PInt.mod|PInt.hvy, 'str': 'TSGS', \
'best': PBest.shower, 'type': PType.freezing}, \
95: {'int': PInt.lgt|PInt.mod, 'str': 'TSRA', \
'best': PBest.shower, 'type': PType.liquid}, \
96: {'int': PInt.lgt|PInt.mod, 'str': 'TSGS', \
'best': PBest.shower, 'type': PType.freezing}, \
97: {'int': PInt.hvy, 'str': 'TSRA', \
'best': PBest.shower, 'type': PType.liquid}, \
98: {'int': 0, 'str': 'TS SS', 'best': 0, 'type': 0}, \
99: {'int': PInt.hvy, 'str': 'TSGR', \
'best': PBest.shower, 'type': PType.freezing}, \
}
MOS_obvis = set([5,10]+range(30,36)+[38,39]+range(42,50))
# disjoint
MOS_ptype = {}
for key in PType.values():
MOS_ptype[key] = set([x for x in MWxTable if x>=50 and \
MWxTable[x]['type']&key])
# disjoint
MOS_pbest = {}
for key in PBest.values():
MOS_pbest[key] = set([x for x in MWxTable if x>=50 and \
MWxTable[x]['best']&key])
# not disjoint
MOS_pint = {}
for key in PInt.values():
MOS_pint[key] = set([x for x in MWxTable if x>=50 and \
MWxTable[x]['int']&key])
###############################################################################
# matches GFS MOS
# 1 2 3 4 5 6 7 8
CigCat = [0, 190, 450, 950, 1950, 3050, 6550, 12000, 999999]
# 1 2 3 4 5 6 7
VisCat = [0, 0.45, 0.95, 1.95, 2.9, 5.5, 6.5, 999.0]
UnltdCig = 22000 # as appears in hd5 file
#def get_cat(value, seq):
# for n, (lo, hi) in enumerate(Avn.window(seq)):
# if lo <= value < hi:
# return n
# else:
# raise ValueError('Invalid value: %f' % value)
###############################################################################
# conversion functions
Unlimited = 99999
def hd2us_cig(cig):
'''Converts ceiling from [m] to [ft].
Unlimited ceiling is set to 99999'''
if 0 <= cig <= 12000:
return cig/0.305
elif cig < 30000:
return Unlimited
else:
return -1
def us2hd_cig(cig, unlimited):
'''Converts ceiling from [ft] to [m].
Unlimited (99999) ceiling is set to "unlimited", which should be
set to fh.variables['cig'].Unlimited'''
if cig == Unlimited:
return unlimited
else:
return cig*0.305
def hd2us_vsby(vsby):
'''Converts visibility from [m] to statue miles.'''
if 0 <= vsby <= 100000:
return vsby/1609.0
else:
return -1
def us2hd_vsby(vsby):
'''Converts visibility from statue miles to [m].'''
return vsby*1609.0
def hd2us_wind_speed(speed):
'''Converts speed from [m/s] to [kt].'''
if 0 <= speed <= 100:
return int(speed*3600.0/1852.0+0.5)
else:
return -1
def us2hd_wind_speed(speed):
'''Converts speed from [kt] to [m/s].'''
return speed*1852.0/3600.0
# flight category definition (metric)
FlightCats = { \
0: {'cig': us2hd_cig(190, 22000), 'vis': us2hd_vsby(0.45)}, \
1: {'cig': us2hd_cig(490, 22000), 'vis': us2hd_vsby(0.95)}, \
2: {'cig': us2hd_cig(950, 22000), 'vis': us2hd_vsby(2.9)}, \
3: {'cig': us2hd_cig(3050, 22000), 'vis': us2hd_vsby(5.2)}, \
}
###############################################################################
class Cache:
def __init__(self):
self.home = os.path.join(os.environ['TOP_DIR'], 'data', 'climate',
'cache')
def put(self, id_, key, data):
dir_ = os.path.join(self.home, id_)
try:
if not os.path.isdir(dir_):
os.mkdir(dir_)
path = os.path.join(dir_, key)
cPickle.dump(data, file(path, 'w'), -1)
return True
except IOError, e:
_Logger.exception('put() failed for %s\n%s', path, e)
return False
def get(self, id_, key, cpath=None):
path = os.path.join(self.home, id_, key)
# compare times, return None if climate file newer than cache
try:
if cpath and os.path.getmtime(path) < os.path.getmtime(cpath):
_Logger.info('climate data file newer than cache')
return None
return cPickle.load(file(path))
except (ValueError, IOError, OSError), e:
_Logger.info('get() failed for %s\n%s', path, e)
return None
Reference in a new issue View git blame Copy permalink