[gradsusr] Regridding in GrADS

[Davide Sacchetti]

your guess is correct: it may be you did not regrid over time

you could loop over time, eg:
set gxout fwrite
it=1;while(it<NTIME);it=it+1
 set t it
 define newvar ...
 display newvar
endwhile
disable fwrite

bye bye
Davide


On Thu, 2018-07-26 at 09:59 +0800, Joseph Basconcillo wrote:
> Hi, I have NCEP-RA2 data that I want to regrid to N768 resolution,
> which is about 0.17 deg in midlatitude.
> 
> I used:
> 
> 'define newvar=re(var, 0.17)'
> 
> And it worked. 
> 
> However, the file size is smaller than I expected. The original size
> is 63MB (2.5 deg) and when I regrid them to 0.17 deg, the file size
> resulted to about 17.5MB. I have my doubts considering I re-gridded
> the data to higher resolution. File size should be higher. That was
> what I expected.
> 
> Then when I used to tloop and area-average to test the data:
> 
> 'd tloop(aave(var,global))'
> 
> A prompt appeared that "entire grid are set to missing data". And it
> kept looping for few minutes. My guess is that the re-gridding did
> not "regrid" through time and considered only the first timestep. Can
> you help me on this?
> 
> Also, is there an option for regridding method (e.g. bilinear)? 
> 
> Thank you.
> 
> Yours,
> 
> Joseph Basconcillo
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr