'reinit' * para deixar o fundo branco 'set display color white' 'clear' 'sdfopen output.nc' 'set mpdset mres' 'set map 1 1 0' 'set lon -85 -30' 'set lat -55 10' 'set gxout contour' 'set cint 3' 'set grads off' 'set parea 0.5 5.25 0 8.5' 'set digsiz 0.105' 'set t 1' *transforma latitude e longitude de graus para radianos 'define latr=lat*3.1415/180' 'define lonr= lon*3.1415/180' *raio da Terra 'define r=6.371e6' *variacao em relacao a latitude 'define dy=r*cdiff(latr,y)' *variacao em relacao a longitude 'define dx=r*cos(latr)*cdiff(lonr,x)' *variacao de temperatura com long 'define dtx=cdiff(t2m,x)' *variacao da temperatura com relacao a lat 'define dty=cdiff(t2m,y)' *unidade: oC/km 'define dtdx=dtx/dx*1000' 'define dtdy=dty/dy*1000' 'set cint 3' 'd t2m-273' 'set arrscl 0.3 0.01' 'd skip(dtdx,5,5);dtdy;mag(dtdx,dtdy)' 'draw string 0.9 7.9 Temperatura 2m (oC) e Grad T (oC/km) - JAN 2011' 'printim gradt2m.gif'