xan/awips2
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

ASM #645 - GFE: CarSnowAmt smart tool can return incorrect amount

Browse source 
Change-Id: I8aea3ab2ccfec7e411e605ddfdfeef0a823687c2

Former-commit-id: 25f4c8acdd8975d3001528212cbc7617cfd39500
This commit is contained in:
Bruce.Hundermark 2014-12-04 15:56:01 +00:00
parent 52d59a2bc9
commit 6f73b491f8
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
cave/com.raytheon.viz.gfe/localization/gfe/userPython/smartTools/CarSnowAmt.py
View file

@ -167,8 +167,8 @@ class Tool (SmartScript.SmartScript):
cubeShape = (len(tCube), tCube.shape[1], tCube.shape[2])
gridShape = (tCube.shape[1], tCube.shape[2])
layerSR = zeros(cubeShape, 'int8')
pvvAvg = zeros(cubeShape, 'int8')
pvvSum = zeros(gridShape, 'int8')
pvvAvg = zeros(cubeShape, 'float32')
pvvSum = zeros(gridShape, 'float32')
for i in range(len(ghCube) - 1):
#for i in range(len(ghCube)):
@ -235,10 +235,10 @@ class Tool (SmartScript.SmartScript):
totalSnowRatio = where(less(pvvSum, 1.0), thicknessSnowRatio, totalSnowRatio)
# If there's any layer above 0.0C, snowRatio gets 0
mask = greater(tCube, 272.65)
mask = any(less_equal(tCube, 272.65), axis = 0)
mask = sum(mask) # reduce to single level by adding bits vertically
# if mask == 0, nowhere in the column is temp < 0.5C
totalSnowRatio = where(equal(mask, 0), totalSnowRatio, 0.0)
totalSnowRatio = where(mask, totalSnowRatio, 0.0)
#Calculate Snowfall - taper to zero from 31 to 34 F.
snowfall = QPF * totalSnowRatio
snowfall = where(greater(T, 31.0), pow(35.0 - T,2)/16.0 * snowfall , snowfall)