Lifted Index

Hi,

I am having problems with the follow code from Hennaks nordic.gs script:


  • 4 LAYER LIFTED INDEX

i = 3
while ( i<maps )

  • Set up a few preliminary characteristics
    setmap(1)
    setNWNcols(1)

'set t 'i

‘define av4li = (no4lftxsfc+lftxsfc)/2’
‘set gxout shaded’
‘set csmooth on’
‘set clevs -9.25 -9 -8.75 -8.5 -8.25 -8 -7.75 -7.5 -7.25 -7 -6.25 -6.5 -6-25 -6 -5.75 -5.5 -5.25 -5 -4.75 -4.5 -4.25 -4 -3.75 -3.5 -3.25 -3 -2.75 -2.5 -2.25 -2 -1.75 -1.5 -1.25 -1 -0.75 -0.5 -0.25 0’
‘set ccols 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 0’
‘d av4li’

*‘set gxout contour’
*‘set ccolor 0’
*‘set cmax 0’
*‘set cmin -25’
*‘set cint 1’
*‘set clab on’
*‘d no4lftxsfc’

‘cbarm’

‘q dims’
times = sublin(result,5)
hub = subwrd(times,6)

‘set strsiz 0.10’
‘set string 1 r 3 0’ ; ‘draw string 9.35 0.95 Map (c) Finnish WX-server/nordicweather.net’

‘set strsiz 0.12’
‘set string 1 r 3 90’ ; ‘draw string 9.9 4.6 K’
‘set string 1 r 3 0’ ; 'draw string 9.45 0.6 Valid: 'hub
‘set string 1 l 3 0’ ; 'draw string 0.15 0.2 Data: NOAA GFS model, run: 'huh

  • Draw Labels at the top of the page
    ‘set string 1 r 11 0’
    ‘set strsiz 0.14’
    ‘set string 1 l 6 0’ ; ‘draw string 0.15 0.6 Av 4 layer Lifted index’
    ‘printim li’i’.png x900 y675’

i = i+1
endwhile

When the script is run opengrads produces the following error:

Warning: Contour levels are not strictly increasing
This may lead to errors or undesired results
Warning: Contour levels are not strictly increasing
This may lead to errors or undesired results
Cannot plot color bar: No shading information.

The images it produces are blank as seen in the example attached.

Please can anyone help me get this working.

Many Thanks,

Ross