function plt(args)
caso=subwrd(args,1)
if (caso='') ; caso='ERA_5' ; endif
date=subwrd(args,2)
if (date='') ; date='2023021312' ; endif
casestudy='SaoSebastiao'

'reinit'
dir='/pesq/dados/bamc/paulo.bonatti/Modal_Energetics/Version_3.0/Horizontal_Decomposition/dataout/'
filea='ENCL'caso''date''date'P.icn.RB0125L037.ctl'
fileb='ENCL'caso'20230201002023020100P.icn.RB0125L037.ctl'
'open 'dir''filea
'open 'dir''fileb

'query time'
dg=subwrd(result,3)
dateg=substr(dg,4,9)

'set display color white'
'clear'

'set grads off'
'set mpdset brmap_mres'
'set map 1 1 3'
'set xlopts 1 8 0.15'
'set ylopts 1 8 0.15'
'set clopts 1 8 0.12'
'set lat -45 45'
'set cthick 8'
'display (etkvkv1.1-etkvkv1.2(t=1))/1000'
pa=subwrd(result,2)
pb=subwrd(result,4)
dp=subwrd(result,6)
ct='Contours: 'pa' to 'pb' by 'dp' [kJ/kg]'
ha=434.5
hb=125.1
za=05
zb=10
tr='Case Study : 'casestudy'  RB0125L037'
'draw title `1'caso' Kelvin Energy Class 1 : H'za'='ha' to H'zb'='hb'\`1 Daily Mean 'dateg' Minus February Monthly Mean'
'draw xlab `1'ct'\`1'tr

fig='../png/Energy_Class1_Anom_Kelvin_'caso'_'casestudy'_'dateg'_TR_RB0125L037.png'
say ' '
say ' fig = 'fig
say ' '
'printim 'fig' png white'

return