ELIASSEN-PALM CROSS SECTIONS

[Eric DeWeaver]

Dear Trevor, Francisco, and Nguyen,

Here's a small grads script that calculates EP fluxes and writes them
to an output file.  The file  'lifecyc.ctl' has winds and temperatures on
standard pressure levels.


'reinit'
'open lifecyc.ctl'

'set gxout fwrite'
'set fwrite data.d/lifecyc_epf.dat'
'set lon 0';'set z 1 13'
'fcor=2*7.292e-5*sin(lat*3.1416/180)'
'ae=6.37e6';'grav=9.81';'pi2=2*3.1416'
'sgk=pow(lev/1000,0.286)'
'cs=cos(lat*3.1416/180)'

tt=1;tmax=360
while(tt<tmax+1)
   'set t 'tt
   'tbar=ave(tmp,lon=0,lon=360,-b)'
   'tsgk=tbar/sgk'
   stab='(sgk*cdiff(tsgk,z)/cdiff(100*lev,z))'
   'vpup=ave(v*(u-ave(u,lon=0,lon=360,-b)),lon=0,lon=360,-b)'
   'vptp=ave(v*(tmp-ave(tmp,lon=0,lon=360,-b)),lon=0,lon=360,-b)'
   'fphi=  -vpup*cs*cs'
   'fprs=  -fcor*ae*cs*cs*vptp'
   'fprs=fprs*(pi2*ae*ae/grav)'
   'fphi=fphi*(pi2*ae*ae/grav)'
   'd fphi';'d fprs/'stab
   say 'wrote record 'tt
   tt=tt+1
endwhile


I haven't used it for awhile, and I don't recall the details,
other than that the output looked ok to me.  It's based on my
reading of the Edmon et al. paper.

I hope this helps, and I welcome any feedback you might have.

Best Wishes,

Eric
-----------------------------------------------
Eric DeWeaver
Assistant Professor
Department of Atmospheric and Oceanic Sciences
University of Wisconsin - Madison
1225 W. Dayton St.
Madison, WI 53706
608-265-5438 (voice) 608-263-4190 (fax)
-----------------------------------------------



On Thu, 15 Mar 2007, Chi Nguyen wrote:

> Hello Alexander,
>
> I am interested in your program as well. It would be great if you could send me that. Thank you very much in advance!
>
> Mai
> -----------------------------------
> Nguyen Chi Mai
> School of Mathematical Sciences
> R. 462, Bld. 28, Clayton Campus
> Monash University
> Work Ph. +61 3 9905 4424
>
>
> Francesco Piani <f.piani at IBIMET.CNR.IT> wrote:>
>> I'm very interested in your program fortran. Can you send me a copy.
>> Thanks in advance.
>> Francesco Piani
>>
>> -------------------------------------------------------------
>> Dott.Ing. Francesco Piani
>> Weather Forecaster
>> Tuscany Region
>> Laboratory for Meteorology and Environmental Modelling - National Council
>> of
>> Research
>> Via Madonna del Piano, 10 - Building D
>> 50019 Sesto Fiorentino (Fi)
>>  Italy
>>
>> Tel.  +39 055 44 830.23
>> Fax  +39 055 44 40 83
>> email: piani at lamma.rete.toscana.it
>>
>> ----- Original Message -----
>> From: "alexander.beck-ratzka" <alexander.beck-ratzka at AEI.MPG.DE>
>> To: <GRADSUSR at LIST.CINECA.IT>
>> Sent: Wednesday, March 14, 2007 2:58 PM
>> Subject: Re: ELIASSEN-PALM CROSS SECTIONS
>>
>>
>>> On Wednesday 03 May 2006 22:20, HALL, Trevor C wrote:
>>>> Hi all,
>>>>  Does anyone have a script to do Eliassen-Palm cross sections.
>>>>
>>>
>>> I have a fortran program which does it, on basis of regular latitude
>>> longitude
>>> grid. If you like, I can provide this to you.
>>>
>>> Cheers
>>>
>>> Alexander
>>>
>>> --
>>> ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
>>> Dr. Alexander Beck-Ratzka - GAT-Developer, group eScience
>>>
>>> MPI for Gravitational Physics (Albert Einstein Institute)
>>> Am Mühlenberg 1
>>> D-14476 Potsdam, Germany
>>>
>>> Tel: 0049-(0)331-567 7192
>>> Email: alexander.beck-ratzka at aei.mpg.de
>>> ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
>