[gradsusr] Issue with Fwrite

[Davide Sacchetti]

it's a "shift" problem due to 0 - 360 longitude domain

look:
-rw-r--r-- 1 model staff 2304000 Jan 27 14:17
3B42.110423.00z.7.precipitation.bin
2304000 = 4*400*1440 correct!

-rw-r--r-- 1 model staff 2307200 Jan 27 15:50 Obs22apr2011.bin
2307200 = 4*400*1442 (you fwrited two more longitudes)
just to check try to change 1440 into 1442 in your 
Obs22apr2011.ctl and it could work (if it doesn't, just for check,
reduce 0.25 to 0.24 - otherwise it could overlap)

solution may be to force 'set x 1 1440' or 'set x 1 1439' before fwrite


hoping it may help
bye bye

Davide

On Tue, 2015-01-27 at 16:49 +0530, Geeta Geeta wrote:
> Hi  All. 
> 
> 
> I have 8 RF data files (TRMM) and I need to accumulate this Rainfall
> for selected files. I have used following script. Following is the ctl
> file and below is gs file. rfaccum.ctl file::::::::::::::::::::::::
> It looks like, 
> DSET ^3B42.1104%d2.%h2z.7.precipitation.bin
> options template
> TITLE Precipitation estimates
> UNDEF -9999.9
> XDEF 1440 LINEAR -179.875 0.25
> YDEF 400 LINEAR -49.875 0.25
> ZDEF 1 LEVELS 1000
> TDEF 8 LINEAR 06z22apr2011 3hr
> VARS 1
> data 0 99 3B42 TRMM precipitation estimates [mm/3hr]
> ENDVARS
> GS FILE ______________________:::::::::::::::::::::::::::::::::::
> 'reinit''
> open rfaccum.ctl''
> define rf1=data(t=1)'
> 'define rf2=data(t=2)'
> 'define rf3=data(t=3)'
> 'define rf4=data(t=4)'
> 'define rf5=data(t=5)''
> define rf6=data(t=6)'
> 'define rf7=data(t=7)'
> 'define rf8=data(t=8)'
> 'define rf=rf1+rf2+rf3+rf4+rf5+rf6+rf7+rf8'
> 'd rf'
> 'set gxout fwrite'
> 'set fwrite Obs22apr2011.bin'
> 'd rf'
> 'disable fwrite'
> ~The plot of accumulated RF is fine. but when I plot the Obs22apr2011.bin, there is a change in the PLOT can U pls help me. I am attaching the two plots(Figure1.gif and Fig2.gif) these FIGURES are so much different. The file Sizes are also different
> 
> -rw-r--r-- 1 model staff 2304000 Jan 27 14:17 3B42.110423.00z.7.precipitation.bin
> -rw-r--r-- 1 model staff 2304000 Jan 27 14:18 3B42.110423.03z.7.precipitation.bin
> -rw-r--r-- 1 model staff 2304000 Jan 27 14:22 3B42_daily.2011.04.22.7.bin
> -rw-r--r-- 1 model staff 2307200 Jan 27 15:50 Obs22apr2011.bin
> bash-3.2$
> Following ctl file is used for plotting the DAILY data
> dset Obs22apr2011.bin 
> TITLE Precipitation estimates
> UNDEF -9999.9
> XDEF 1440 LINEAR -179.875 0.25
> YDEF 400 LINEAR -49.875 0.25
> ZDEF 1 LEVELS 1000
> TDEF 1 LINEAR 00z22apr2011 1dy
> VARS 1 data 0 99 3B42 TRMM precipitation estimates [mm/3hr]
> ENDVARS
> ~dailyrf.ctl" 11 lines, 305 characters
> Can someone suggest???
> 
> Geeta
> 
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-- 
Davide Sacchetti
Centro Funzionale Meteo Idrologico di Protezione Civile della Regione Liguria
ARPAL Unità Tecnica Complessa di livello Regionale
V.le Brigare Partigiane 2 16121 Genova (I)
tel: +39 010 6437535                    fax: +39 010 6437520
mail: davide.sacchetti at arpal.gov.it     web: www.meteoliguria.it