You need to set dfile 1 or 2 ... before defining your variable in the corresponding file.
Hope this will help
Regards !
Mouhamadou Bamba SYLLA, PhD
Post-Doc Fellow, PWC/ESP
International Centre for Theoretical Physics
Strada Costiera 11, 34100 Trieste, Italy
Office tel: 0039 040 2240 225
Cell: 0039 34 71 94 06 55
emails: msylla@ictp.it, syllabamba@yahoo.fr
***************************************************
De : rajesh maurya <rksmaurya@gmail.com>
À : regcnet@lists.ictp.it
Envoyé le : Jeu 13 mai 2010, 9h 54min 50s
Objet : [RegCNET] Fwd: How to calculate the geopotential height and defintion problem
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
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@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@gmail.com
Cc: regcnet@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@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@cdac.in
Tel (O): +91-020-2570 4226, (Cell): +91-094223 56326
--------------------------------------------------------------------------------
From: Satyaban Bishoyi Ratna <satyaban_2001@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@gmail.com
Cc: regcnet@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 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@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@cdac.in
Tel (O): +91-020-2570 4226, (Cell): +91-094223 56326
--------------------------------------------------------------------------------
--
Rajesh Kumar Maurya
M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research Scholar
Department of Geophysics
Banaras Hindu University
Varanasi 221005