[gradsusr] Problem for saving Wind Climatology on GrADS
Sarah Barbier
barbier.sarah at live.fr
Fri Mar 10 07:36:13 EST 2017
Hello GrADS users,
I am trying to attach labels as title and save my wind climatology from ERA INTERIM data from dec 1991 to nov 2007on grads but I do not know why, as soon as I define my climatological variables (u10clim and v10clim) and try for instance to plot climatological winter DJF, the title do not appear and when I try to display my plot after saving, it is blank. The eps file is created but it is blank ..
In addition to that, I cannot overlay none data like seasonal sst .
Here you will find below the code (.gs file) I’m trying to run :
‘reinit’
‘set display color white’
‘c’
‘set grid off’
‘set grads off’
‘set mproj scaled’
‘set mpdset hires’
‘sdfopen wnd10m_erainterim_sst.nc’
‘set lat -20 20
’set lon -80 20’
‘set t 1 12’
‘define u10clim=ave(u10,t+0,t=192,12)’
‘define v10clim=ave(v10,t+0,t=192,12)’
‘modify u10clim seasonal’
'modify v10clim seasonal'
‘ define winteru = (u10clim(t=1),u10clim(t=2),u10clim(t=3))/3'
‘ define winterv = (v10clim(t=1),v10clim(t=2),v10clim(t=3))/3'
‘ define springu = (u10clim(t=4),u10clim(t=5),u10clim(t=6))/3'
‘ define springv = (u10clim(t=4),u10clim(t=5),v10clim(t=6))/3'
‘ define summeru = (v10clim(t=7),v10clim(t=8),v10clim(t=9))/3'
‘ define summerv = (v10clim(t=7),v10clim(t=8),v10clim(t=9))/3'
‘ define fallu = (u10clim(t=10),u10clim(t=11),u10clim(t=12))/3'
‘ define fallv = (v10clim(t=10),v10clim(t=11),v10clim(t=12))/3’
‘ d skip(winter,3,3);winter'
‘draw title DJF Wind climate'
'print essaie.eps’
Did I miss something ?
Regards,
Sarah Barbier
barbier.sarah at live.fr<mailto:barbier.sarah at live.fr>
Student in Master 2 MOCES Paris VI
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170310/8dbcf644/attachment.html
More information about the gradsusr
mailing list