Jun 10, 2018

Bering Sea Rule Verification

For the Bering Sea Rule (Renken et al.) fans out there, I am slowly automating verification results. Not sure what metrics to correlate or at what time frequency, but it is a start. Attached is an example of a 5-day H5 anomaly for 6/2-6.



A popular way to visualize the Bering Sea Rule is to utilize the maps @crankywxguy generates (http://www.stormhamster.com/bsr/bsr.php) which showcases daily surface and 500MB. I may just run with those metrics for the time being and concentrate on the "pattern" scoring 5-day anomalies.

For those interested in the verification method; the verification data are generated from ESRL NetCDF files utilizing the scorr function within GrADS. A snippet of the GrADS script is below. Beware, my GrADS scripting is amatuer at best. Enjoy!

function script(args)
myVarA=subwrd(args,1)
myVarB=subwrd(args,2)
myVarC=subwrd(args,3)
myVarD=subwrd(args,4)

'reinit'
rc = gsfallow("on")
nrows=2
ncols=1
args=nrows%' '%ncols
'sdfopen Z:\VMSharedFolder\Projects\Verification\esrl-ncdeflect\'myVarC
'set mproj latlon'
panels(args)

p=1
_vpg.1
'set parea 1 9.8 1 7.5'
'set t 1'
'set grads off'
'set gxout shaded'
'set rgb 16 0 0 255'
'set rgb 17 55 55 255'
'set rgb 18 110 110 255'
'set rgb 19 165 165 255'
'set rgb 20 220 220 255'
'set rgb 21 255 220 220'
'set rgb 22 255 165 165'
'set rgb 23 255 110 110'
'set rgb 24 255 55 55'
'set rgb 25 255 0 0'
'set rgb 26 255 250 250'
'set rgb 27 125 125 125'
'set clevs -100 -80 -60 -40 -20 20 40 60 80 100'
'set ccols 16 17 18 19 20 27 21 22 23 24 25'
'set lon 140 210'
'set lat 39 65'
'set frame off'
'set mpdraw on'
'set gxout shaded'
'set csmooth on'
'd hgt1'
'd scorr(hgt1,hgt2,lon=0,lon=75,lat=57.5,lat=90)'
myVar=result
'run cbarc.gs'
'draw string 1.6 .35 Bering Sea Rule Intialized: 'myVarB
'draw string 5.8 .35 Bering Sea Rule Forecast Dates: 'myVarA
'draw string 2.8 4.12 500mb GPH (m) Anomaly Spatial Correlation (Bering Sea to CONUS Teleconnection)'
'draw title 'myVar

p=2
_vpg.2
'set parea 1 9.8 1 7.5'
'set t 1'
'set grads on'
'set gxout shaded'
'set rgb 16 0 0 255'
'set rgb 17 55 55 255'
'set rgb 18 110 110 255'
'set rgb 19 165 165 255'
'set rgb 20 220 220 255'
'set rgb 21 255 220 220'
'set rgb 22 255 165 165'
'set rgb 23 255 110 110'
'set rgb 24 255 55 55'
'set rgb 25 255 0 0'
'set rgb 26 255 250 250'
'set rgb 27 125 125 125'
'set clevs -100 -80 -60 -40 -20 20 40 60 80 100'
'set ccols 16 17 18 19 20 27 21 22 23 24 25'
'set lon 232.5 302.5'
'set lat 24.5 50'
'set frame off'
'set mpdraw on'
'set gxout shaded'
'set csmooth on'
'd hgt2'
'draw string 3.9 .35 CONUS Observed Dates: 'myVarA
'draw string 0.05 .27 @OSNW3'

'set vpage off'
'printim Z:\VMSharedFolder\Projects\Verification\result-images\'myVarD
'set gxout print'
'set prnopts %6.2f'
write('bsr-hgt.txt', myVarB' 'myVar, append)
;


Also, there is some pre-GrADS NetCDF magic to line up the two domains for use with the scorr function. If there are any questions, comments, or suggestions on the material presented please let me know. Thank you for reading.

Jun 7, 2018

RRWT Forecast Verification

Slowly automating RRWT forecast verification. Utilizing GrADS scorr function at the moment. Unsure hemispheric spatial correlation is the best path, but I find it entertaining that generating a forecast verification is considerably more intensive than generating a forecast.



Experimenting with the verification strategy and metrics. Attempting to verify each of the current outlooks.



Scores located here: http://www.consonantchaos.com/a-all.html

For those interested in the verification method, it is a quite simple. The verification data are generated from ESRL NetCDF files utilizing the scorr function within GrADS. If there are any questions, comments, or suggestions on the material presented please let me know. Thank you for reading.

Jun 1, 2018

2018 Indian Summer Monsoon

14 of the 31 May RRWT precipitation rate outlooks for June, July, and August are suggesting a potential deficit in Indian rainfall compared to the average of the past 10 years.




8 of the 31 are suggesting similar accumulations of the past 10 years.





If there are any questions, comments, or suggestions on the material presented please let me know. Thank you for reading.