[gradsusr] Issue with Fwrite

James T. Potemra jimp at hawaii.edu
Tue Jan 27 14:11:05 EST 2015


Also note you can make use of the GrADS "sum" function to replace
the lines with "define rf1=data(t=1)", "define rf2=data(t=2)" and so
on, e.g., "rf = sum(data,t=1,t=8)"
More details at http://www.iges.org/grads/gadoc/gadocindex.html

Jim

On 1/27/15 2:55 AM, Davide Sacchetti wrote:
> 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



More information about the gradsusr mailing list