iemand een idee waarom de gfs0.25 data niet werkt bij mijn oude scripts die op 0.5° perfect werken? Als test een klein script hierna, hij geeft fout bij
'define t825=tmpprs(lev=825)'. Fout luidt 'World coordinates convert to non-integer grid coordinates....
thx
'set lat 49 55'
'set lon 2 7'
'define t825=tmpprs(lev=825)'
'define t800=tmpprs(lev=800)'
n=1
while (n<2)
timer=((n*3)-3)
'set t 'n''
'q time'
res = subwrd(result,3)
datum2 = substr(res,4,9)
uur = substr(res,1,2)
dag = subwrd(result,6)
fff = substr(dag,1,2)
clear
'set display color white'
clear
'define pmsl= prmslmsl/100'
'define t2m = tmp2m-273.15'
********************TWmsl*********************
'define tw = ((-5.806+(0.672*t2m))-(0.006*(t2m*t2m))+((0.061+(0.004*t2m)+(0.000099*(t2m*t2m)))*rh2m)+(((-0.000033-(0.000005*t2m)-(0.0000001*(t2m*t2m))))*(rh2m*rh2m)))'
********************TWmsl*********************
'set gxout grid'
'set gridln off'
'set ccolor 1'
'set digsize 0.10'
'set dignum 1'
'display snodsfc/100'
'draw title T2m (C) 'dag' 'datum2' 'uur'z +'timer'h (run: GFS'run'z)'
n=n+1
endwhile
Quote selectie