[RegCNET] Fwd: How to calculate the geopotential height and defintion problem
Marcelino II Villafuerte
mqvillafuerte at gmail.com
Thu May 13 10:25:40 CEST 2010
hi rajesh,
I attached a grads script I used to produce a monthly average rainfall data
from daily values netcdf files separated
every year. It works at grads 2.0 version or later. you can modify it for
your need. for example you need to replace
sdfopen with simply open because your opening a ctl file. Notice that I used
undefine to allow other variables be defined.
hope it helps! good luck!
On Thu, May 13, 2010 at 3:54 PM, rajesh maurya <rksmaurya at gmail.com> wrote:
> Hi Satyaban Bishoyi Ratna
>
> I did the same as you suggested but i am getting the warning messages.
> Please have a look on this and give your suggestion.
>
> ga-> define r1=ave(tpr.1,t=1,t=120)*30/10
> Averaging. dim = 3, start = 1, end = 120
> Define memory allocation size = 133584 bytes
> ga-> define r2=ave(tpr.2,t=1,t=124)*31/10
> Averaging. dim = 3, start = 1, end = 124
>
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -119 -119 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -118 -118 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -117 -117 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -116 -116 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -115 -115 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -114 -114 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -113 -113 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -112 -112 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = -111 -111 E
> = 1 1
>
>
> Warning issued for variable = tpr.2
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
>
>
> ---------- Forwarded message ----------
> From: Satyaban Bishoyi Ratna <satyaban_2001 at rediffmail.com>
> Date: Thu, May 13, 2010 at 12:50 PM
> Subject: Re: [RegCNET] How to calculate the geopotential height and
> defintion problem
> To: rksmaurya at gmail.com
> Cc: regcnet at lists.ictp.it
>
>
>
> Hi Rajesh,
> Why don't you use the Post Processing tool of RegCM3.
> This post processing tool will solve all your query in this mail.
> This tool can help you to derive monthly mean of your desired variables and
> also it can get geopotential height.
>
> Anyway, there is some wrong with defining time with your grads script.
> You see each CTL file and average the variable as per the number of time
> stamp available. You can use this example script with slight modification of
> your script.
>
>
> define r1=ave(tpr.1,t=1,t=120)*30/10
> define r2=ave(tpr.2,t=1,t=124)*31/10
> define r3=ave(tpr.3,t=1,t=124)*31/10 ....and so on.
>
> The same method you can apply for wind. Hope this helps you.
>
> Thanks.
>
> Regards
> Satyaban
>
>
>
>
> On Thu, 13 May 2010 12:01:23 +0530 wrote
>
> >Hello All,
>
> I am plotting monthly and seasonal precipitaion , circulation
> field(200,850).
> These are the step(Given below) i am following for plot like precipitation
> and vector wind.
> Let me correct , if i am missing something for correct plot.
> I need help on how to plot geopotential height(monthly and
> seasonal) and give the idea(with example) .
> Please help me regarding the above problem.
>
> ga-> open ATM.2005060100.ctl
> Scanning description file: ATM.2005060100.ctl
> Data file ATM.2005060100 is open as file 1
> LON set to 49.4242 111.019
> LAT set to -10.48 39.77
> LEV set to 1008.43 1008.43
>
> Time values set: 2005:6:1:6 2005:6:1:6
> E set to 1 1
> ga-> open ATM.2005070100.ctl
> Scanning description file: ATM.2005070100.ctl
> Data file ATM.2005070100 is open as file 2
> ga-> open ATM.2005080100.ctl
>
> Scanning description file: ATM.2005080100.ctl
> Data file ATM.2005080100 is open as file 3
> ga-> open ATM.2005090100.ctl
> Scanning description file: ATM.2005090100.ctl
> Data file ATM.2005090100 is open as file 4
>
>
> Define function
>
> 1)- precipitation
>
> ga-> define r1=ave(tpr.1,t=1,t=120)*30/10
> Averaging. dim = 3, start = 1, end = 120
> Define memory allocation size = 133584 bytes
> ga-> define r2=ave(tpr.2,t=120,t=244)*31/10
> Averaging. dim = 3, start = 120, end = 244
>
> Data Request Warning: Request is completely outside file limits
> Entire grid contents are set to missing data
> Grid limits of file: X = 1 138 Y = 1 121 Z = 1 23 T = 1 124 E = 1
> 1
> Grid limits of request: X = 1 138 Y = 1 121 Z = 1 1 T = 0 0 E = 1 1
>
> Warning issued for variable = tpr.2
> Define memory allocation size = 133584 bytes
>
> ga-> define r3=ave(tpr.3,t=244,t=368)*31/10
> ...
> ..
> ga-> define r4=ave(tpr.4,t=368,t=488)*30/10
> ..
> ..
> ga-> define r=(r1+r2+r3+r4)
>
>
> 2) vector wind
>
> ga-> define u1=ave(u.1,t=1,t=120)
>
>
>
> ga-> define v1=ave(v.1,t=1,t=120)
>
>
> ga-> define u2=ave(u.2,t=120,t=244)
>
>
>
> ga-> define v2=ave(v.2,t=120,t=244)
>
> ga-> d skip(u1,5,5);v1
> ..
> ...
> ..
> ga-> define u=(u1+u2+u3+u4)
>
> ga-> define v=(v1+v2+v3+v4)
>
>
> --
> Waiting for your necessary suggestion on it.
>
> Thanks
> Rajesh Kumar Maurya
> M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research
> Scholar
> Department of Geophysics
> Banaras Hindu University
> Varanasi 221005
>
> _______________________________________________
>
> RegCNET mailing list
>
> RegCNET at lists.ictp.it
>
> https://lists.ictp.it/mailman/listinfo/regcnet
>
>
>
> -------------------------------------------------------------------------------
> Dr. Satyaban Bishoyi Ratna
> Computational Atmospheric Sciences Team
> Scientific & Engineering Computing Group
> Centre for Development of Advanced Computing (C-DAC)
> Pune University Campus, Ganeshkhind
> Pune - 411 007, Maharashtra, India
> E-mail: satyabanb at cdac.in
> Tel (O): +91-020-2570 4226, (Cell): +91-094223 56326
>
> --------------------------------------------------------------------------------
>
> <http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?>
>
>
>
> --
> Rajesh Kumar Maurya
> M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research
> Scholar
> Department of Geophysics
> Banaras Hindu University
> Varanasi 221005
>
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20100513/373174bd/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aphrotomonthly.gs
Type: application/octet-stream
Size: 6025 bytes
Desc: not available
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20100513/373174bd/attachment-0002.obj>
More information about the RegCNET
mailing list