[gradsusr] smoothing large resolution data
wendi harjupa
wendiharjupa at gmail.com
Wed Jan 25 20:38:33 EST 2017
Dear Mr. Andrew
Thank you very much for your instructions.
I am very sorry for asking you many questions.
But, I got error message as below :
Data Request Warning: Request is completely outside file limits
Entire grid contents are set to missing data
Grid limits of file: X = 1 280 Y = 1 240 Z = 1 1 T = 1 1 E = 1 1
Grid limits of request: X = 134 160 Y = 107 128 Z = 1 1 T = -156 -156
E =1 1
Would you mind to explain me the error meaning?
Thank you very much.
Best regards,
wendi
On Tue, Jan 24, 2017 at 12:33 AM, Andrew Friedman <andfried at gmail.com>
wrote:
> Hi Wendy,
>
> I would place the command at the end of your script.
>
> 'define try=(rr, albedo.2, aave)’
>
> Do you get an error? Based on your script, I’m wondering if you might be
> closing your files too soon.
> Assuming no error message, you can go to a specific lat/lon point where
> you know there is data in the original file and type
> 'd try’
>
> You can also try ‘set gxout grid’ to look at the values in specific
> gridboxes.
>
> Hope this helps.
> Andrew
>
> > On Jan 23, 2017, at 1:20 PM, wendi harjupa <wendiharjupa at gmail.com>
> wrote:
> >
> > Dear Mr. Andrew
> >
> > Thank you very much for your excellent supervision
> >
> > Below is my simple script
> >
> > 'open ./Rain.ctl'
> > 'set lat 33.5 35.5'
> > 'set lon 134.3 136.6'
> > # 'd rr' ! variable
> > 'xcbar 8.85 9.05 0.8 4 -fw 0.11 -fh 0.15 -edge triangle -fstep 1
> -line on -lc 1'
> > 'close 1'
> >
> > ##second data/resolution 0.1x0.1 (lat, lon)
> >
> > 'open ./Him.ctl'
> > 'set lat 33.5 35.5'
> > 'set lon 134.3 136.6'
> > # 'd albedo' ! variable
> > 'xcbar 6.3 9 0.23 0.37 -fw 0.11 -fh 0.15 -edge triangle -fs 2 -fo 1
> -line on'
> > 'close 1'
> >
> > 'd lterp(rr, albedo.2, aave)'
> >
> > my control file as below
> > 1. Rain.ctl
> >
> > DSET ~/data1
> > OPTIONS yrev
> > UNDEF -1.E35
> > XDEF 800 LINEAR 134.375 0.003125
> > YDEF 960 LINEAR 33.5833 0.00208333
> > ZDEF 1 LINEAR 1 1
> > TDEF 120 LINEAR 08Aug2015 1hr
> > VARS 1
> > rr 0 98 rairate [mm/hr]
> > ENDVARS
> >
> > 2. Him.ctl
> > DSET ${file_in_vor2}
> > UNDEF -1.0
> > OPTIONS yrev
> > XDEF 280 LINEAR 121.05 0.1
> > YDEF 240 LINEAR 22.95 0.1
> > ZDEF 1 LINEAR 1 1
> > TDEF 1 LINEAR 13:05Z08AUG2015 5mn
> > VARS 1
> > ALBEDO 0 1 reflectivity
> > ENDVARS
> >
> > My next question where should I locate my lterp command??
> > and how it would be like?
> > I have tried this command below, it doesnor work
> >
> > 'define try=(rr, albedo.2, aave)'
> > 'd try'
> >
> > Thank you very much for your help. I am very sorry for ask many
> questions.
> >
> > Best regards,
> > wendi
> >
> >
> >
> >
> >
> >
> > On Mon, Jan 23, 2017 at 6:48 PM, Andrew Friedman <andfried at gmail.com>
> wrote:
> > Hi Wendi,
> >
> > You can use the command as in examples 4 & 5 for the lterp documentation:
> > 'define newvar = (rr, albedo.2, aave)’
> > What does newvar look like over specific gridboxes where you know there
> is data?
> >
> > Again, I would try first running a simplified version of your script
> that doesn’t call other scripts in order to simplify the problem as much as
> possible.
> >
> > Best, Andrew
> >
> > > On Jan 23, 2017, at 3:08 AM, wendi harjupa <wendiharjupa at gmail.com>
> wrote:
> > >
> > > Dear Mr. Andrew
> > >
> > > Many thanks for your great support.
> > >
> > > I think my problem is, how to use the lterp function. where should I
> write the lterp funtion
> > > and how to write is. I just found some sample in using lterp command,
> but It did work.
> > >
> > > Based on the sample script that I have sent in the previous email,
> where should I locate my lterp
> > > command?
> > >
> > > Thank you very much.
> > >
> > > Best regards,
> > > wendi
> > >
> > > On Sun, Jan 22, 2017 at 11:43 PM, Andrew Friedman <andfried at gmail.com>
> wrote:
> > > Hi Wendi,
> > >
> > > I don’t see an error offhand with your dimensions.
> > > What does the numerical output of lterp look like in the regions where
> you have data? Do you get any values?
> > >
> > > Andrew
> > >
> > > > On Jan 21, 2017, at 4:11 AM, wendi harjupa <wendiharjupa at gmail.com>
> wrote:
> > > >
> > > > Dear Mr. Andrew,
> > > >
> > > > Many thanks for your great support,
> > > >
> > > > Below is q dims and q file of my data,
> > > >
> > > > ## Data 1 (large scale : 0.1x0.1 deg)
> > > >
> > > > q dims 1
> > > > X is varying Lon = 121.05 to 148.95 X = 1 to 280
> > > > Y is varying Lat = 22.95 to 46.85 Y = 1 to 240
> > > > Z is fixed Lev = 1 Z = 1
> > > > T is fixed Time = 13:05Z08AUG2015 T = 1
> > > > E is fixed Ens = 1 E = 1
> > > >
> > > > q file 1
> > > >
> > > > Descriptor: Him.ctl
> > > > Binary: /TAKAYASU/wendi/RDCA/rdca_201508080520.dat
> > > > Type = Gridded
> > > > Xsize = 280 Ysize = 240 Zsize = 1 Tsize = 1 Esize = 1
> > > > Number of Variables = 1
> > > > albedo 0 1 reflectivity
> > > >
> > > >
> > > > ## data 2 (small scale : 0.002x0.003 deg)
> > > >
> > > > q dims 2
> > > >
> > > > X is varying Lon = 134.375 to 136.872 X = 1 to 800
> > > > Y is varying Lat = 33.5833 to 35.5812 Y = 1 to 960
> > > > Z is fixed Lev = 1 Z = 1
> > > > T is fixed Time = 00Z08AUG2015 T = 1
> > > > E is fixed Ens = 1 E = 1
> > > >
> > > > q file 2
> > > >
> > > > Binary: /TAKAYASU/wendi/xband/G000.grads/KINKI/20150808_0900-1700/data/KINKI00001-
> 20150808-1420-G000-EL000000.bin
> > > > Type = Gridded
> > > > Xsize = 800 Ysize = 960 Zsize = 1 Tsize = 120 Esize = 1
> > > > Number of Variables = 1
> > > > rr 0 98 rairate [mm/hr]
> > > >
> > > > is there any mistakes on that information?
> > > >
> > > > Best regards,
> > > > wendi
> > > >
> > > >
> > > >
> > > >
> > > > On Fri, Jan 20, 2017 at 6:25 PM, Andrew Friedman <andfried at gmail.com>
> wrote:
> > > > Hi Wendi,
> > > > What is your dimension environment?
> > > > Can you post the output of ‘q dims’, ‘q file 1’ and ‘q file 2'
> > > > Andrew
> > > >
> > > > > On Jan 20, 2017, at 1:28 AM, wendi harjupa <wendiharjupa at gmail.com>
> wrote:
> > > > >
> > > > > Dear Mr. Andrew
> > > > >
> > > > > Thank you very much for your explanation,
> > > > > I am still struggling to solve it.
> > > > >
> > > > > Best regards,
> > > > > wendi
> > > > >
> > > > > On Wed, Jan 18, 2017 at 10:16 PM, Andrew Friedman <
> andfried at gmail.com> wrote:
> > > > > Hi Wendi,
> > > > >
> > > > > You can check ‘q dims’ to see what your time dimensions are.
> > > > >
> > > > > I recommend going through the website examples for lterp.gs,
> particularly #3-5 - these show situations that often come up about
> mismatched dimensions.
> > > > > As they show, you often need to specify time or dfile, or some
> combination. Another possibility could be wrapped longitude coordinates.
> > > > >
> > > > > If none of these steps work, can you post a simplified version of
> your problem, without the anciallary scripts in your code?
> > > > >
> > > > > Andrew
> > > > >
> > > > > > On Jan 18, 2017, at 2:56 AM, wendi harjupa <
> wendiharjupa at gmail.com> wrote:
> > > > > >
> > > > > > Dear Mr Andrew,
> > > > > >
> > > > > > Many thanks for your response
> > > > > >
> > > > > > There is no specific error, but the is no change of result when
> I use
> > > > > > 'd lterp command'.
> > > > > >
> > > > > > I set my dimension environment as below :
> > > > > >
> > > > > > 'set lat 33.000 36.000'
> > > > > > 'set lon 134.000 137.000'
> > > > > >
> > > > > > I am little confused how to use the 'd lterp' function in my
> script program.
> > > > > >
> > > > > > Would you mind to give me more clue?
> > > > > >
> > > > > > Best regards,
> > > > > > wendi
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > On Wed, Jan 18, 2017 at 5:55 AM, Andrew Friedman <
> andfried at gmail.com> wrote:
> > > > > > Hi Wendi,
> > > > > > What specific error message are you getting?
> > > > > > A typical problem is that the time axes of the two files have to
> be the same - you may want to check your dimension environment.
> > > > > > Andrew
> > > > > >
> > > > > > > On Jan 17, 2017, at 1:53 PM, wendi harjupa <
> wendiharjupa at gmail.com> wrote:
> > > > > > >
> > > > > > > Dear Mr. Andrew and colleagues
> > > > > > >
> > > > > > > Please tell me what wrong with this script. It does not work
> to interpolate the coarse grid to fine grid
> > > > > > >
> > > > > > > ##first data/resolution 0.002x0.003 (lat, lon)
> > > > > > >
> > > > > > > 'open ./Rain.ctl'
> > > > > > > 'setXrain.gs'
> > > > > > > 'setCOMMON.gs' ! set parameter (lon, lat, parea, grid
> off, etc)
> > > > > > > # 'd rr' ! variable
> > > > > > > 'xcbar 8.85 9.05 0.8 4 -fw 0.11 -fh 0.15 -edge triangle
> -fstep 1 -line on -lc 1'
> > > > > > > 'close 1'
> > > > > > >
> > > > > > > ##second data/resolution 0.1x0.1 (lat, lon)
> > > > > > >
> > > > > > > 'open ./Him.ctl'
> > > > > > > 'setCOMMON.gs'
> > > > > > > 'setHimawari.gs'
> > > > > > > # 'd albedo' ! variable
> > > > > > > 'xcbar 6.3 9 0.23 0.37 -fw 0.11 -fh 0.15 -edge triangle
> -fs 2 -fo 1 -line on'
> > > > > > > 'close 1'
> > > > > > >
> > > > > > > 'd lterp(rr, albedo.2, aave)'
> > > > > > >
> > > > > > > Any help will be very appreciated.
> > > > > > > Thank you very much.
> > > > > > >
> > > > > > > Best regards,
> > > > > > > wendi
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > On Tue, Jan 17, 2017 at 7:12 PM, wendi harjupa <
> wendiharjupa at gmail.com> wrote:
> > > > > > > Dear Mr. Andrew..
> > > > > > >
> > > > > > > Many thanks for your valuable information.
> > > > > > >
> > > > > > > Best regards,
> > > > > > > wendi
> > > > > > >
> > > > > > > On Tue, Jan 17, 2017 at 6:45 PM, Andrew Friedman <
> andfried at gmail.com> wrote:
> > > > > > > Hi Wendi,
> > > > > > > You can use the function lterp for spatial interpolation:
> http://cola.gmu.edu/grads/gadoc/gradfunclterp.html
> > > > > > > -Andrew
> > > > > > >
> > > > > > > > On Jan 17, 2017, at 10:29 AM, wendi harjupa <
> wendiharjupa at gmail.com> wrote:
> > > > > > > >
> > > > > > > > Dear All
> > > > > > > >
> > > > > > > > I got problem in overlaying two type of data.
> > > > > > > > As I have attached the plot of data in this email.
> > > > > > > >
> > > > > > > > First data is rain data with resolution is about 0.002x0.003
> deg (lon ,lat).
> > > > > > > > The rain data is indicated by rainbow color.
> > > > > > > > and second data is cloud data with resolution is 0.1x0.1 deg
> (lon, lat)
> > > > > > > > The cloud data is indicated by gray color
> > > > > > > >
> > > > > > > > The plot looks not good. Is there any tools in grads to
> smooth the cloud data?
> > > > > > > > Then, we can interpolate cloud data to become 0.005 or
> 0.001, and plot data will be more nice.
> > > > > > > >
> > > > > > > > Thank you very much for your attention and help.
> > > > > > > >
> > > > > > > > Best regards,
> > > > > > > > wendi
> > > > > > > >
> > > > > > > > <large skale.png>____________________
> ___________________________
> > > > > > > > gradsusr mailing list
> > > > > > > > gradsusr at gradsusr.org
> > > > > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > > > > >
> > > > > > >
> > > > > > > _______________________________________________
> > > > > > > gradsusr mailing list
> > > > > > > gradsusr at gradsusr.org
> > > > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > > > > >
> > > > > > >
> > > > > > > _______________________________________________
> > > > > > > gradsusr mailing list
> > > > > > > gradsusr at gradsusr.org
> > > > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > > > >
> > > > > >
> > > > > > _______________________________________________
> > > > > > gradsusr mailing list
> > > > > > gradsusr at gradsusr.org
> > > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > > > >
> > > > > > _______________________________________________
> > > > > > gradsusr mailing list
> > > > > > gradsusr at gradsusr.org
> > > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > > >
> > > > >
> > > > > _______________________________________________
> > > > > gradsusr mailing list
> > > > > gradsusr at gradsusr.org
> > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > > >
> > > > > _______________________________________________
> > > > > gradsusr mailing list
> > > > > gradsusr at gradsusr.org
> > > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > >
> > > >
> > > > _______________________________________________
> > > > gradsusr mailing list
> > > > gradsusr at gradsusr.org
> > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > > >
> > > > _______________________________________________
> > > > gradsusr mailing list
> > > > gradsusr at gradsusr.org
> > > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170126/3bb3ae2c/attachment-0001.html
More information about the gradsusr
mailing list