[RegCNET] How to calculate the geopotential height and defintion problem

[rajesh maurya]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20100513/bb98d086/attachment.html>