'sdfopen /stornext/online11/eta/ERA5/PSLM_SA/ClimatologiaMensal_1991-2020_ERA5_PSLM.nc'
'sdfopen /stornext/online11/eta/ERA5/PSLM_SA/202101_ERA5_PSLM_MonMean.nc'
'sdfopen /stornext/online11/eta/ERA5/ERA5-OROG_LSMK_SA.nc'

'set mpdset brmap_hires'

'set rgb 37 213 159 213'
'set rgb 38 192 114 192'
'set rgb 16 163  73 164'
'set rgb 17 134  60 136'
'set rgb 18  91  40  91'
*******azul************
'set rgb 39   1  56 101'
'set rgb 19   2  71 129'
'set rgb 20   3  94 151'
'set rgb 21   9 156 219'
'set rgb 22  19 194 233'
'set rgb 23  82 235 254'
'set rgb 24 170 242 254'
'set rgb 25 217 254 254'
*laranja claro ao vermelho
'set rgb 26 255 240 200'
'set rgb 27 255 220 180'
'set rgb 28 255 180 140'
'set rgb 29 255 140 100'
'set rgb 30 255 100  60'
'set rgb 31 255  61  45'
'set rgb 32 255  20   0'
'set rgb 33 236  18   0'
'set rgb 34 166  13   0'
'set rgb 35 121  10   0'

'set display color white'
'c'
*'set lat -30 0'
*'set lon -70 15'
*'set lat -30 0' 
*'set lon -45 -15'
'set grid off'
'set gxout shaded'
'set xlopts 1 4 0.10'
'set ylopts 1 4 0.10'
'set xlint 10'
'set ylint 5'
'set clevs      -15 -10   -5  -4  -3  -2   -1 -0.5 0.5  1    2   3 4  5  10 15'
'set ccols   18      19     20     21      22      23     24    25    0   26   27     28     29     30     31     32      33 '
'set grads off' 
'set stat on'
'd maskout((msl.2(t=1)-msl(t=1))/100,LSM.3(t=1)*-1)'
*'run cbarn 1 1 8.7 4.25'
'run cbarn '
'draw title PSLM Anomaly (hPa) - Jan 2021'
*   'run cbarn'
'printim Anom_PSLM_Jan2021.png'
'quit'
'q pos'
return
say result
'set gxout fwrite'
'set fwrite file.bin' 
'd maskout(Mag1-MagClim,LSM.5(t=1)*-1)'
'disable fwrite'
*'set gxout contour'
*'set clevs 0'
*'d v10.2(t='tempo'))'
*'set clevs 0'
*'d v10.4(t='tempo'))'
***return
*'quit'