function main(args)

* Parse the arguments: date, hour
  if (args = '')
  prompt 'Enter forecast date (example, 20110717) --> ' 
  pull date
  prompt 'Enter forecast hour (example, 00 or 06 or 12 or 18) --> ' 
  pull hour
  else 
  date = subwrd(args,1)
  hour  = subwrd(args,2)
  endif

'reinit'
'open ecmwf_grib.ctl'

* Get info from the descriptor file
'q ctlinfo'
_ctl = result
_undef = getctl(undef)
_tdef = getctl(tdef)
_zdef = getctl(zdef)

maps = 25

* Get the Time axis info
tsize = subwrd(_tdef,2)
_t1 = 1       ;* 2nd half of timeseries
_t2 = 45
tsize = _t2 - _t1 + 1
'set t '_t1' '_t2
'q dims'
times  = sublin(result,5)
_time1 = subwrd(times,6)  
_time2 = subwrd(times,8)
_tdim = _time1' '_time2

tincr = subwrd(_tdef,5)
_tdef = 'tdef 'tsize' linear '_time1' 'tincr

huh = subwrd(_tdef,4)




****************************************
* 500 Height
****************************************
'set annot 0 5'
  i = 1
  while ( i<maps )



* Set up a few preliminary characteristics
setmap(1)
setTeOsPtcols(1)
'set t 'i


'set gxout shaded'
'set csmooth on'
'define geopotdam = HGT500mb(lev=500)/10'
'set clevs 476 480 484 488 492 496 500 504 508 512 516 520 524 528 532 536 540 544 548 552 556 560 564 568 572 576 580 584 588 592 596 600 604 608'
'set ccols 36 36 36 37 38 39 40 41 47 48 49 50 51 52 53 56 58 58 59 60 66 67 69 70 72 73 74 75 77 78 80 81 82 83 84 84'
            
'd geopotdam'

'define PRESmsl = PRESmsl*0.01'

'set gxout contour'
'set ccolor 0'
'set cthick 8'
'set cint 5'
'set clopts 1'
'set clab masked'
'd PRESmsl'

'define HGT500mb=HGT500mb'
'set lev 500'
'set gxout contour'
'set ccolor 1'
'set cint 5'
'set cthick 10'
'set cmax 5520'
'set cmin 5520'
'set clab masked'
'set clopts 1'
'd HGT500mb'
'cbarm'

'set line 1 1 20'


'q dims'
times  = sublin(result,5)

'q time'
timestr=subwrd(result, 3)
weekday1=subwrd(result,6)
day=substr(result, 11, 9)
hr=substr(result, 7, 20)
h=substr(result, 8, 3)
m=substr(result, 13, 3)
d=substr(result, 11, 2)
y=substr(result, 16, 4)
'set strsiz 0.12'
'set string 1 r 4 0' ; 'draw string 9.47 0.6 Valid:  'weekday1' 'd''m''y' 'h



'set strsiz 0.12'
'set string 1 r 4 0' ; 'draw string 5.63 0.6 500hpa & Pressure'
'set strsiz 0.12'
'set string 1 l 4 0' ; 'draw string 0.09 0.6 ECMWF Run: 'huh

* Draw Labels at the top of the page
'set string 1 r 11 0'
'set strsiz 0.12'
'set string 1 l 4 0' ; 'draw string 0.17 1.0 chorleyweather.com (c)'

'printim gfs_NA500hpa'i'.png x1100 y675'

i = i+1
endwhile






* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
* END OF MAIN SCRIPT
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

function setmap(args)
'set mproj nps'
'set mpvals -50 25 30 90'
'set lon -150 100'
'set lat 20 90'
'set mpdset hires'

'set rgb 98 238 238 238'
'set display color white'
'set background 0'
'c'
'set grads off'

'set parea 0.02 9.54 0.8 8.50'
'set rgb 99 1 1 1'
'set line 99 1 6'
*'draw rec 0.02 0.02 9.54 0.81'
*'draw rec 9.54 0.02 10.99 8.48'

return

function setstucols(args)

*light yellow to dark red

'set rgb 21 255 250 170'

'set rgb 22 255 232 120'

'set rgb 23 255 192  60'

'set rgb 24 255 160   0'

'set rgb 25 255  96   0'

'set rgb 26 255  50   0'

'set rgb 27 225  20   0'

'set rgb 28 192   0   0'

'set rgb 29 165   0   0'

*


*light green to dark green


'set rgb 31 230 255 225'

'set rgb 32 200 255 190'

'set rgb 33 180 250 170'

'set rgb 34 150 245 140'

'set rgb 35 120 245 115'

'set rgb 36  80 240  80'

'set rgb 37  55 210  60'

'set rgb 38  30 180  30'

'set rgb 39  15 160  15'


*set rgb 39   5 150   5
*


*light blue to dark blue


'set rgb 41 225 255 255'

'set rgb 42 180 240 250'

'set rgb 43 150 210 250'

'set rgb 44 120 185 250'

'set rgb 45  80 165 245'

'set rgb 46  60 150 245'

'set rgb 47  40 130 240'

'set rgb 48  30 110 235'

'set rgb 49  20 100 210'

*


*light purple to dark purple


'set rgb 51 220 220 255'

'set rgb 52 192 180 255'

'set rgb 53 160 140 255'

'set rgb 54 128 112 235'

'set rgb 55 112  96 220'   

'set rgb 56  72  60 200'   

'set rgb 57  60  40 180'

'set rgb 58  45  30 165'

'set rgb 59  40   0 160'

*

*light pink to dark rose  


'set rgb 61 255 230 230'

'set rgb 62 255 200 200'

'set rgb 63 248 160 160'

'set rgb 64 230 140 140'

'set rgb 65 230 112 112'

'set rgb 66 230  80  80'   

'set rgb 67 200  60  60'   

'set rgb 68 180  40  40'

'set rgb 69 164  32  32'

*

*light beige to dark brown   


'set rgb 71 250 240 230'

'set rgb 72 240 220 210'

'set rgb 73 225 190 180'

'set rgb 74 200 160 150'

'set rgb 75 180 140 130'   

'set rgb 76 160 120 110'  

'set rgb 77 140 100  90'

'set rgb 78 120  80  70'

'set rgb 79 100  60  50'

*

return

function setcwcols(args)
'set rgb 33 248 50 60'
'set rgb 34 255 50 89'
'set rgb 35 255 50 185'
'set rgb 36 248 50 255'
'set rgb 37 224 50 255'
'set rgb 38 195 50 255'
'set rgb 39 175 50 255'
'set rgb 40 161 50 255'
'set rgb 41 137 50 255'
'set rgb 42 118 74 255'
'set rgb 43 98 74 255'
'set rgb 44 79 50 255'
'set rgb 45 50 50 255'
'set rgb 46 50 74 255'
'set rgb 47 50 89 255'
'set rgb 48 50 113 255'
'set rgb 49 50 146 255'
'set rgb 50 50 175 255'
'set rgb 51 50 204 255'
'set rgb 52 50 224 255'
'set rgb 53 50 255 253'
'set rgb 54 50 255 228'
'set rgb 55 50 255 200'
'set rgb 56 50 255 161'
'set rgb 57 50 255 132'
'set rgb 58 50 255 103'
'set rgb 59 50 255 79'
'set rgb 60 50 255 60'
'set rgb 61 79 255 50'
'set rgb 62 132 255 50'
'set rgb 63 171 255 50'
'set rgb 64 204 255 50'
'set rgb 65 224 255 50'
'set rgb 66 253 255 50'
'set rgb 67 255 233 50'
'set rgb 68 255 224 50'
'set rgb 69 255 209 50'
'set rgb 70 255 204 50'
'set rgb 71 255 185 50'
'set rgb 72 255 161 50'
'set rgb 73 255 146 50'
'set rgb 74 255 127 50'
'set rgb 75 255 118 50'
'set rgb 76 255 93 50'
'set rgb 77 255 74 50'
'set rgb 78 255 60 50'
'set rgb 79 255 33 33'
'set rgb 80 255 0 0'
'set rgb 81 235 10 0'
'set rgb 82 215 20 0'
'set rgb 83 195 30 0'
'set rgb 84 175 40 0'
'set rgb 85 165 45 0'
'set rgb 86 155 50 0'
'set rgb 87 145 55 0'
'set rgb 88 135 60 0'
'set rgb 89 120 60 0'
'set rgb 90 100 60 0'
'set rgb 91 80 60 0'
'set rgb 20 250 240 230'
'set rgb 21 240 220 210'
'set rgb 22 225 190 180'
'set rgb 23 200 160 150'
'set rgb 24 180 140 130' 
'set rgb 25 160 120 110' 
'set rgb 26 140 100 90'


return

function setNWNcols(args)
'set rgb 99 1 1 1'
'set rgb 16 100 100 100'

*Green
'set rgb 17 170 246 130'
'set rgb 18 175 244 116'
'set rgb 19 180 242 102'
'set rgb 20 190 241 88'
'set rgb 21 200 240 74'
'set rgb 22 210 238 60'
'set rgb 23 220 236 46'
'set rgb 24 230 235 32'
'set rgb 25 235 232 22'
'set rgb 26 244 230 11'

'set rgb 27 244 217 11'
'set rgb 28 244 203 11'
'set rgb 29 244 189 11'
'set rgb 30 244 176 11'
'set rgb 31 244 162 11'
'set rgb 32 244 149 11'
'set rgb 33 244 136 11'
'set rgb 34 244 122 11'
'set rgb 35 244 109 11'
'set rgb 36 244 95 11'
'set rgb 37 244 82 11'
'set rgb 38 244 68 11'
'set rgb 39 244 61 11'
'set rgb 40 232 55 9'
'set rgb 41 220 39 8'
'set rgb 42 206 32 7'
'set rgb 43 196 26 10'
'set rgb 44 186 19 15'
'set rgb 45 175 15 20'
'set rgb 46 165 11 25'
'set rgb 47 165 6 30'
'set rgb 48 170 2 35'
'set rgb 49 180 1 40'
'set rgb 50 190 0 50'
'set rgb 51 200 0 60'
'set rgb 52 210 0 70'
'set rgb 53 220 0 80'
'set rgb 54 230 0 90'
*Dark red

*Dark Blue

'set rgb 55 180 20 185'
'set rgb 56 165 16 174'
'set rgb 57 150 12 163'
'set rgb 58 135 8 152'
'set rgb 59 120 4 141'

'set rgb 60 100 0 127'
'set rgb 61 87 0 127'
'set rgb 62 75 0 127'
'set rgb 63 62 0 127'
'set rgb 64 50 0 127'
'set rgb 65 37 0 127'
'set rgb 66 25 0 127'
'set rgb 67 12 0 127'
'set rgb 68 0 0 127'

'set rgb 69 0 24 127'
'set rgb 70 0 30 127'
'set rgb 71 0 40 127'
'set rgb 72 0 50 127'
'set rgb 73 0 60 127'
'set rgb 74 0 70 127'

'set rgb 75 0 82 143'
'set rgb 76 0 98 175'
'set rgb 77 0 114 207'
'set rgb 78 0 130 239'
'set rgb 79 19 146 255'
'set rgb 80 37 154 255'

'set rgb 81 73 172 255'
'set rgb 82 91 180 255'
'set rgb 83 109 188 255'
'set rgb 84 127 196 255'
'set rgb 85 145 204 255'
'set rgb 86 154 208 255'
'set rgb 87 163 212 255'
'set rgb 88 181 220 255'
'set rgb 89 199 228 255'
'set rgb 90 217 236 255'
*Light blue

return

function setBLUEcols(args)

*Dark Blue
'set rgb 56 100 0 127'
'set rgb 57 75 0 127'
'set rgb 58 50 0 127'
'set rgb 59 25 0 127'
'set rgb 60 0 0 127'
'set rgb 61 0 24 127'
'set rgb 62 0 50 127'
'set rgb 63 0 74 127'
'set rgb 64 0 82 143'
'set rgb 65 0 90 159'
'set rgb 66 0 98 175'
'set rgb 67 0 106 191'
'set rgb 68 0 114 207'
'set rgb 69 0 122 223'
'set rgb 70 0 130 239'
'set rgb 71 0 138 255'
'set rgb 72 19 146 255'
'set rgb 73 37 154 255'
'set rgb 74 55 162 255'
'set rgb 75 73 172 255'
'set rgb 76 82 176 255'
'set rgb 77 91 180 255'
'set rgb 78 100 184 255'
'set rgb 79 109 188 255'
'set rgb 80 118 192 255'
'set rgb 81 127 196 255'
'set rgb 82 136 200 255'
'set rgb 83 145 204 255'
'set rgb 84 154 208 255'
'set rgb 85 163 212 255'
'set rgb 86 172 216 255'
'set rgb 87 181 220 255'
'set rgb 88 190 224 255'
'set rgb 89 199 228 255'
'set rgb 90 208 232 255'
'set rgb 91 217 236 255'
'set rgb 92 226 240 255'
'set rgb 93 235 244 255'
'set rgb 94 244 248 255'
'set rgb 95 253 252 255'
*Light blue

return

function setBLACKcols(args)

'set rgb 50 1 1 1'
'set rgb 51 11 11 11'
'set rgb 52 23 23 23'
'set rgb 53 34 34 34'
'set rgb 54 46 46 46'
'set rgb 55 57 57 57'
'set rgb 56 69 69 69'
'set rgb 57 80 80 80'
'set rgb 58 91 91 91'
'set rgb 59 102 102 102'
'set rgb 60 114 114 114'
'set rgb 61 126 126 126'
'set rgb 62 138 138 138'
'set rgb 63 150 150 150'
'set rgb 64 161 161 161'
'set rgb 65 173 173 173'
'set rgb 66 185 185 185'
'set rgb 67 196 196 196'
'set rgb 68 207 207 207'
'set rgb 69 219 219 219'
'set rgb 70 231 231 231'
'set rgb 71 242 242 242'
'set rgb 72 254 254 254'


return

function setcols(args)


'set rgb 175 73 172 255'
'set rgb 176 82 176 255'
'set rgb 177 91 180 255'
'set rgb 178 100 184 255'
'set rgb 179 109 188 255'
'set rgb 180 118 192 255'
'set rgb 181 127 196 255'
'set rgb 182 136 200 255'
'set rgb 183 145 204 255'
'set rgb 184 154 208 255'
'set rgb 185 163 212 255'
'set rgb 186 172 216 255'
'set rgb 187 181 220 255'
'set rgb 188 190 224 255'
'set rgb 189 199 228 255'
'set rgb 190 208 232 255'
'set rgb 191 217 236 255'
'set rgb 192 226 240 255'
'set rgb 193 235 244 255'
'set rgb 194 244 248 255'
'set rgb 195 253 252 255'
*Light blue


'set rgb 20 234 245 234'
'set rgb 21 200 215 200'
'set rgb 22 160 205 160'
'set rgb 23 120 215 120'
'set rgb 24  80 235  80'
'set rgb 25   0 255   0'
'set rgb 26   0 195   0'
'set rgb 27   0 160   0'
'set rgb 28   0 125   0'

'set rgb 30 255 160 120'
'set rgb 31 160 120 255'
'set rgb 32 160 180 205'

'set rgb 42  32 208  32'
'set rgb 43 208  32 208'
'set rgb 44  64  64 255'
'set rgb 45 255 120  32'
'set rgb 46  32 208 208'
'set rgb 47 240 240   0'

'set rgb 96 139 115  85'
'set rgb 97 100 100 100'
'set rgb 98  64  64  96'
'set rgb 99 254 254 254'

return


function setcols(args)

'set rgb 175 73 172 255'
'set rgb 176 82 176 255'
'set rgb 177 91 180 255'
'set rgb 178 100 184 255'
'set rgb 179 109 188 255'
'set rgb 180 118 192 255'
'set rgb 181 127 196 255'
'set rgb 182 136 200 255'
'set rgb 183 145 204 255'
'set rgb 184 154 208 255'
'set rgb 185 163 212 255'
'set rgb 186 172 216 255'
'set rgb 187 181 220 255'
'set rgb 188 190 224 255'
'set rgb 189 199 228 255'
'set rgb 190 208 232 255'
'set rgb 191 217 236 255'
'set rgb 192 226 240 255'
'set rgb 193 235 244 255'
'set rgb 194 244 248 255'
'set rgb 195 253 252 255'
*Light blue


'set rgb 20 234 245 234'
'set rgb 21 200 215 200'
'set rgb 22 160 205 160'
'set rgb 23 120 215 120'
'set rgb 24  80 235  80'
'set rgb 25   0 255   0'
'set rgb 26   0 195   0'
'set rgb 27   0 160   0'
'set rgb 28   0 125   0'

'set rgb 30 255 160 120'
'set rgb 31 160 120 255'
'set rgb 32 160 180 205'

'set rgb 42  32 208  32'
'set rgb 43 208  32 208'
'set rgb 44  64  64 255'
'set rgb 45 255 120  32'
'set rgb 46  32 208 208'
'set rgb 47 240 240   0'

'set rgb 96 139 115  85'
'set rgb 97 100 100 100'
'set rgb 98  64  64  96'
'set rgb 99 254 254 254'

* farby pre geopotencialny rozdiel
'set rgb 100 0 180 255'
'set rgb 101 0 200 200'
'set rgb 102 0 205 170'
'set rgb 103 0 210 140'
'set rgb 104 0 215 70'
'set rgb 105 0 220 0'
'set rgb 106 80 225 25'
'set rgb 107 160 230 50'
'set rgb 108 195 225 50'
'set rgb 109 230 220 50'
'set rgb 110 230 195 47'
'set rgb 111 230 175 45'
'set rgb 112 235 155 42'
'set rgb 113 240 130 40'
'set rgb 114 248 65 20'
'set rgb 115 255 0 0'
'set rgb 116 255 30 30'
'set rgb 117 250 60 60'
'set rgb 118 245 30 95'
'set rgb 119 240 0 130'
'set rgb 120 235 0 150'
'set rgb 121 230 0 160'
'set rgb 122 225 0 170'

* farby pre zmenu teploty
'set rgb 130 80 130 120'
'set rgb 131 120 110 150'
'set rgb 132 180 160 180'
'set rgb 133 160 0 200'
'set rgb 134 135 0 210'
'set rgb 135 110 0 220'
'set rgb 136 70 30 238'
'set rgb 137 30 60 255'
'set rgb 138 15 110 255'
'set rgb 139 0 160 255'
'set rgb 140 54 180 255'
'set rgb 141 60 200 200'
'set rgb 142 144 212 205'
'set rgb 143 185 238 215'
'set rgb 144 224 254 250'
'set rgb 145 255 255 255'
'set rgb 146 246 254 215'
'set rgb 147 240 245 170'
'set rgb 148 232 232 130'
'set rgb 149 230 225 70'
'set rgb 150 230 200 50'
'set rgb 151 230 175 45'
'set rgb 152 235 152 42'
'set rgb 153 240 130 40'
'set rgb 154 245 95 50'
'set rgb 155 250 60 60'
'set rgb 156 245 30 95'
'set rgb 157 240 0 130'
'set rgb 158 235 0 150'
'set rgb 159 220 0 180'
'set rgb 160 200 0 200'

return

function setPrcols(args)

'set rgb 17 240 240 240'
'set rgb 18 230 230 230'
'set rgb 19 220 220 255'
'set rgb 20 190 190 255'
'set rgb 21 150 150 255 '
'set rgb 22 100 100 240 '
'set rgb 23 70 70 200 '
'set rgb 24 40 40 150 '
'set rgb 25 0 120 0 '
'set rgb 26 0 160 0 '
'set rgb 27 0 210 0 '
'set rgb 28 0 240 0'
'set rgb 29 255 250 0 '
'set rgb 30 255 200 0 '
'set rgb 31 255 150 0'
'set rgb 32 255 70 0'
'set rgb 33 180 0 0'
'set rgb 34 90 0 0'
'set rgb 35 40 0 0'

return



function setTeOsPtcols(args)

'set rgb 33 248 50 60'
'set rgb 34 255 50 89'
'set rgb 35 255 50 185'
'set rgb 36 248 50 255'
'set rgb 37 224 50 255'
'set rgb 38 195 50 255'
'set rgb 39 175 50 255'
'set rgb 40 161 50 255'
'set rgb 41 137 50 255'
'set rgb 42 118 74 255'
'set rgb 43 98 74 255'
'set rgb 44 79 50 255'
'set rgb 45 50 50 255'
'set rgb 46 50 74 255'
'set rgb 47 50 89 255'
'set rgb 48 50 113 255'
'set rgb 49 50 146 255'
'set rgb 50 50 175 255'
'set rgb 51 50 204 255'
'set rgb 52 50 224 255'
'set rgb 53 50 255 253'
'set rgb 54 50 255 228'
'set rgb 55 50 255 200'
'set rgb 56 50 255 161'
'set rgb 57 50 255 132'
'set rgb 58 50 255 103'
'set rgb 59 50 255 79'
'set rgb 60 50 255 60'
'set rgb 61 79 255 50'
'set rgb 62 132 255 50'
'set rgb 63 171 255 50'
'set rgb 64 204 255 50'
'set rgb 65 224 255 50'
'set rgb 66 253 255 50'
'set rgb 67 255 233 50'
'set rgb 68 255 224 50'
'set rgb 69 255 209 50'
'set rgb 70 255 204 50'
'set rgb 71 255 185 50'
'set rgb 72 255 161 50'
'set rgb 73 255 146 50'
'set rgb 74 255 127 50'
'set rgb 75 255 118 50'
'set rgb 76 255 93 50'
'set rgb 77 255 74 50'
'set rgb 78 255 60 50'
'set rgb 79 255 33 33'
'set rgb 80 255 0 0'
'set rgb 81 235 10 0'
'set rgb 82 215 20 0'
'set rgb 83 195 30 0'
'set rgb 84 175 40 0'
'set rgb 85 165 45 0'
'set rgb 86 155 50 0'
'set rgb 87 145 55 0'
'set rgb 88 135 60 0'
'set rgb 89 120 60 0'
'set rgb 90 100 60 0'
'set rgb 91 80 60 0'
'set rgb 20 250 240 230'
'set rgb 21 240 220 210'
'set rgb 22 225 190 180'
'set rgb 23 200 160 150'
'set rgb 24 180 140 130'
'set rgb 25 160 120 110'
'set rgb 26 140 100 90'

return

function setcols(args)

'set rgb 175 73 172 255'
'set rgb 176 82 176 255'
'set rgb 177 91 180 255'
'set rgb 178 100 184 255'
'set rgb 179 109 188 255'
'set rgb 180 118 192 255'
'set rgb 181 127 196 255'
'set rgb 182 136 200 255'
'set rgb 183 145 204 255'
'set rgb 184 154 208 255'
'set rgb 185 163 212 255'
'set rgb 186 172 216 255'
'set rgb 187 181 220 255'
'set rgb 188 190 224 255'
'set rgb 189 199 228 255'
'set rgb 190 208 232 255'
'set rgb 191 217 236 255'
'set rgb 192 226 240 255'
'set rgb 193 235 244 255'
'set rgb 194 244 248 255'
'set rgb 195 253 252 255'
*Light blue


'set rgb 20 234 245 234'
'set rgb 21 200 215 200'
'set rgb 22 160 205 160'
'set rgb 23 120 215 120'
'set rgb 24  80 235  80'
'set rgb 25   0 255   0'
'set rgb 26   0 195   0'
'set rgb 27   0 160   0'
'set rgb 28   0 125   0'

'set rgb 30 255 160 120'
'set rgb 31 160 120 255'
'set rgb 32 160 180 205'

'set rgb 42  32 208  32'
'set rgb 43 208  32 208'
'set rgb 44  64  64 255'
'set rgb 45 255 120  32'
'set rgb 46  32 208 208'
'set rgb 47 240 240   0'

'set rgb 96 139 115  85'
'set rgb 97 100 100 100'
'set rgb 98  64  64  96'
'set rgb 99 254 254 254'

* farby pre geopotencialny rozdiel
'set rgb 100 0 180 255'
'set rgb 101 0 200 200'
'set rgb 102 0 205 170'
'set rgb 103 0 210 140'
'set rgb 104 0 215 70'
'set rgb 105 0 220 0'
'set rgb 106 80 225 25'
'set rgb 107 160 230 50'
'set rgb 108 195 225 50'
'set rgb 109 230 220 50'
'set rgb 110 230 195 47'
'set rgb 111 230 175 45'
'set rgb 112 235 155 42'
'set rgb 113 240 130 40'
'set rgb 114 248 65 20'
'set rgb 115 255 0 0'
'set rgb 116 255 30 30'
'set rgb 117 250 60 60'
'set rgb 118 245 30 95'
'set rgb 119 240 0 130'
'set rgb 120 235 0 150'
'set rgb 121 230 0 160'
'set rgb 122 225 0 170'

* farby pre zmenu teploty
'set rgb 130 80 130 120'
'set rgb 131 120 110 150'
'set rgb 132 180 160 180'
'set rgb 133 160 0 200'
'set rgb 134 135 0 210'
'set rgb 135 110 0 220'
'set rgb 136 70 30 238'
'set rgb 137 30 60 255'
'set rgb 138 15 110 255'
'set rgb 139 0 160 255'
'set rgb 140 54 180 255'
'set rgb 141 60 200 200'
'set rgb 142 144 212 205'
'set rgb 143 185 238 215'
'set rgb 144 224 254 250'
'set rgb 145 255 255 255'
'set rgb 146 246 254 215'
'set rgb 147 240 245 170'
'set rgb 148 232 232 130'
'set rgb 149 230 225 70'
'set rgb 150 230 200 50'
'set rgb 151 230 175 45'
'set rgb 152 235 152 42'
'set rgb 153 240 130 40'
'set rgb 154 245 95 50'
'set rgb 155 250 60 60'
'set rgb 156 245 30 95'
'set rgb 157 240 0 130'
'set rgb 158 235 0 150'
'set rgb 159 220 0 180'
'set rgb 160 200 0 200'

return



function setTeOsPtcols(args)

'set rgb 33 248 50 60'
'set rgb 34 255 50 89'
'set rgb 35 255 50 185'
'set rgb 36 248 50 255'
'set rgb 37 224 50 255'
'set rgb 38 195 50 255'
'set rgb 39 175 50 255'
'set rgb 40 161 50 255'
'set rgb 41 137 50 255'
'set rgb 42 118 74 255'
'set rgb 43 98 74 255'
'set rgb 44 79 50 255'
'set rgb 45 50 50 255'
'set rgb 46 50 74 255'
'set rgb 47 50 89 255'
'set rgb 48 50 113 255'
'set rgb 49 50 146 255'
'set rgb 50 50 175 255'
'set rgb 51 50 204 255'
'set rgb 52 50 224 255'
'set rgb 53 50 255 253'
'set rgb 54 50 255 228'
'set rgb 55 50 255 200'
'set rgb 56 50 255 161'
'set rgb 57 50 255 132'
'set rgb 58 50 255 103'
'set rgb 59 50 255 79'
'set rgb 60 50 255 60'
'set rgb 61 79 255 50'
'set rgb 62 132 255 50'
'set rgb 63 171 255 50'
'set rgb 64 204 255 50'
'set rgb 65 224 255 50'
'set rgb 66 253 255 50'
'set rgb 67 255 233 50'
'set rgb 68 255 224 50'
'set rgb 69 255 209 50'
'set rgb 70 255 204 50'
'set rgb 71 255 185 50'
'set rgb 72 255 161 50'
'set rgb 73 255 146 50'
'set rgb 74 255 127 50'
'set rgb 75 255 118 50'
'set rgb 76 255 93 50'
'set rgb 77 255 74 50'
'set rgb 78 255 60 50'
'set rgb 79 255 33 33'
'set rgb 80 255 0 0'
'set rgb 81 235 10 0'
'set rgb 82 215 20 0'
'set rgb 83 195 30 0'
'set rgb 84 175 40 0'
'set rgb 85 165 45 0'
'set rgb 86 155 50 0'
'set rgb 87 145 55 0'
'set rgb 88 135 60 0'
'set rgb 89 120 60 0'
'set rgb 90 100 60 0'
'set rgb 91 80 60 0'
'set rgb 20 250 240 230'
'set rgb 21 240 220 210'
'set rgb 22 225 190 180'
'set rgb 23 200 160 150'
'set rgb 24 180 140 130'
'set rgb 25 160 120 110'
'set rgb 26 140 100 90'

return

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
function vrng(f1,f2)
'set gxout stat'
'd 'f1
data = sublin(result,8)
ymx = subwrd(data,5)
ymn = subwrd(data,4)
'd 'f2
data = sublin(result,8)
zmx = subwrd(data,5)
zmn = subwrd(data,4)
if (zmx > ymx) ; ymx = zmx ; endif
if (zmn < ymn) ; ymn = zmn ; endif
dy = ymx-ymn
ymx = ymx + 0.08 * dy
ymn = ymn - 0.08 * dy
if ((ymx-ymn)/2.2 < 1)
  incr = (ymx-ymn)/4
  incr = 0.01 * (math_nint(100*incr))
else
  incr = math_nint((ymx-ymn)/4)
endif
'set vrange 'ymn' 'ymx
'set ylint 'incr
if (ymn=0 & ymx=0 & incr=0)
  'set vrange -.9 .9'
  'set ylint 1'
endif
'set gxout line'
return

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
function rh2vrng(f1)
'set gxout stat'
'd 'f1
data = sublin(result,8)
ymn = subwrd(data,4)
ymx = subwrd(data,5)
if (ymn < 20) 
  miny = 0 
  'set ylevs 20 40 60 80'
endif
if (ymn >= 20 & ymn < 30) 
  miny = 20 
  'set ylevs 30 50 70 90'
endif
if (ymn >= 30 & ymn < 40) 
  miny = 30 
  'set ylevs 40 50 60 70 80 90'
endif
if (ymn >= 40 & ymn < 50) 
  miny = 40 
  'set ylevs 50 60 70 80 90'
endif
if (ymn >= 50 & ymn < 60) 
  miny = 50
  'set ylevs 60 70 80 90'
endif
if (ymn >= 60) 
  miny = 60
  'set ylevs 70 80 90'
endif
'set vrange 'miny' 'ymx+3
return


* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
function getctl(handle)
line = 1
found = 0
while (!found)
  info = sublin(_ctl,line)
  if (subwrd(info,1)=handle)
    _handle = info
    found = 1
  endif
  line = line + 1
endwhile
return _handle

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
function getgrid(dodsvar,myvar)

'set lon '_xdim
'set lat '_ydim
'set lev '_zgrd
'set time '_tdim

* Write the variable to a file
'set gxout fwrite'
'set fwrite dummy.dat'
'd 'dodsvar
'disable fwrite'

* Write a descriptor file 
rc = write(dummy.ctl,'dset ^dummy.dat')
rc = write(dummy.ctl,_undef,append)
rc = write(dummy.ctl,'xdef 1 linear 1 1',append)
rc = write(dummy.ctl,'ydef 1 linear 1 1',append)
rc = write(dummy.ctl,_zdef,append)
rc = write(dummy.ctl,_tdef,append)
rc = write(dummy.ctl,'vars 1',append)
rc = write(dummy.ctl,'dummy '_newzsize' -999 dummy',append)
rc = write(dummy.ctl,'endvars',append)
rc = close (dummy.ctl)

* Open the dummy file, define variable, close dummy file
'open dummy.ctl'
line = sublin(result,2)
dummyfile = subwrd(line,8)
'set dfile 'dummyfile
'set lon 1'
'set lat 1'
'set lev '_zbot' '_ztop
'set time '_time1' '_time2
'define 'myvar' = dummy.'dummyfile
'close 'dummyfile
'set dfile 1'
return

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
function getetarh(dodsvar,myvar)

* swap out original pressure vars
tmpzgrd = _zgrd
tmpzdef = _zdef
tmpzbot = _zbot
tmpztop = _ztop
tmpzsize = _newzsize

* retrieve rh data over the rh pressure range
_zgrd = _rhzgrd
_zdef = _trhzdef
_ztop = _rhztop
_zbot = _rhzbot
_newzsize = _trhzsize
getgrid(dodsvar,tmprh)

* swap in original pressure vars
_zgrd = tmpzgrd
_zdef = tmpzdef
_zbot = tmpzbot
_ztop = tmpztop
_newzsize = tmpzsize

'set lon '_xdim
'set lat '_ydim
'set lev '_rhzgrd
'set time '_tdim

* Write the variable to a file
'set gxout fwrite'
'set fwrite dummy.dat'
t = _t1
while (t <= _t2)
  'set t 't
  z = 1
  while (z <= _newrhzsize)
    level = subwrd(_rhlevs,z) 
    'set lev 'level
    'd tmprh'
    z = z + 1
  endwhile
  t = t + 1
endwhile
'disable fwrite'

* Write a descriptor file 
rc = write(dummy.ctl,'dset ^dummy.dat')
rc = write(dummy.ctl,_undef,append)
rc = write(dummy.ctl,'xdef 1 linear 1 1',append)
rc = write(dummy.ctl,'ydef 1 linear 1 1',append)
rc = write(dummy.ctl,_rhzdef,append)
rc = write(dummy.ctl,_tdef,append)
rc = write(dummy.ctl,'vars 1',append)
rc = write(dummy.ctl,'dummy '_newrhzsize' -999 dummy',append)
rc = write(dummy.ctl,'endvars',append)
rc = close (dummy.ctl)

* Open the dummy file, define variable, close dummy file
'open dummy.ctl'
line = sublin(result,2)
dummyfile = subwrd(line,8)
'set dfile 'dummyfile
'set lon 1'
'set lat 1'
'set lev '_rhzbot' '_rhztop
'set time '_time1' '_time2
'q dims'
'define 'myvar' = dummy.'dummyfile
'close 'dummyfile
'set dfile 1'

return

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
function getseries(dodsvar,myvar,level)

'set lon '_xdim
'set lat '_ydim
'set lev 'level' 'level
'set time '_tdim

* Write the variable to a file
'set fwrite dummy.dat'
'set gxout fwrite'
'd 'dodsvar
'disable fwrite'

* Write a descriptor file 
rc = write(dummy.ctl,'dset ^dummy.dat')
rc = write(dummy.ctl,_undef,append)
rc = write(dummy.ctl,'xdef 1 linear 1 1',append)
rc = write(dummy.ctl,'ydef 1 linear 1 1',append)
rc = write(dummy.ctl,'zdef 1 linear 1 1',append)
rc = write(dummy.ctl,_tdef,append)
rc = write(dummy.ctl,'vars 1',append)
rc = write(dummy.ctl,'dummy 0 -999 dummy',append)
rc = write(dummy.ctl,'endvars',append)
rc = close(dummy.ctl)

* Open the dummy file, define variable, close dummy file
'open dummy.ctl'
line = sublin(result,2)
dummyfile = subwrd(line,8)
'set dfile 'dummyfile
'set lon 1'
'set lat 1'
'set lev 'level
'set time '_time1' '_time2
'define 'myvar' = dummy.'dummyfile
'close 'dummyfile
'set dfile 1'
'set gxout contour'

return