[ncl-talk] Scripts for 1000-500mb thickness

[Dennis Shea]

I think that someone who is more familiar with WRF (? wrfhelp at ucar.edu ?)
will have to help. I am not aware of any function that would extrapolate
WRF output to 1000 hPa.

--
A post to the wrf forum in 2014 asked your question:

        http://forum.wrfforum.com/viewtopic.php?f=32&t=7610

It was never answered.

===
The poster tried

z = wrf_user_getvar(a, "z",it)             ; grid point height
p = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical coordinat

z_plane_500   = wrf_user_intrp3d( z,p,"h",500.,0.,False)
z_plane_1000 = wrf_user_intrp3d( z,p,"h",1000.,0.,False)

thick_500_1000 = z_plane_500-z_plane_1000 ;

The poster said this did not work. However, there may have been updates
since that time.

Good luck

On Wed, Aug 2, 2017 at 8:04 AM, Daniel Vecellio <djvecellio at tamu.edu> wrote:

> Oh, sorry. Yeah, it was a general question, so I didn't know how much data
> to put in my email. It is the normal output data from a WRF run on model
> levels. I imagined something like you had above, but was wondering how
> places above 1000 mb might need to be handled.
>
> Thanks,
> Dan
>
> On Wed, Aug 2, 2017 at 8:42 AM, Dennis Shea <shea at ucar.edu> wrote:
>
>> What data do you have available?
>> People are reluctant to get involved unless some information is provided.
>>
>> Are the data on isobaric levels?
>> Are they gridded?
>> Station rawindsonde data?
>> What format?
>> etc
>>
>> For (say) gridded NCEP Reanalysis on isobaric levels, the thickness would
>> simple:
>>
>>    f = addfile("hgt....nc","r")   ; (time,lev,lat,lon)
>>    z = f->hgt
>>    thick = z(:,{500},:,:) - z(:,{1000},:,:)   ; (ntim,nlat,mlon)
>>    copy_VarCoords(z(:,0,:,:), thick)
>>    thick at long_name = "Thickness: 500-1000"
>>    thick at units            = z at units
>>    printVarSummary(thick)
>>
>>
>>
>>
>>
>>
>>
>> On Tue, Aug 1, 2017 at 1:30 PM, Daniel Vecellio <djvecellio at tamu.edu>
>> wrote:
>>
>>> Hello all,
>>>
>>> I know there are NCL scripts already out on the Internet for nearly
>>> everything, but I'm having a ton of trouble finding one for 1000-500mb
>>> thickness. Does anyone have one that they could point me to?
>>>
>>> Thanks,
>>> Dan
>>>
>>> --
>>>
>>> *Daniel J. Vecellio*
>>>
>>> *Ph.D. Student - Climate Science Lab, **Department of Geography, **Texas
>>> A&M University*
>>> *Head Teaching Assistant - Geography 213, Texas A&M University*
>>> *Education Chair - United States Permafrost Association*
>>> *ExCom and United States National Representative - Permafrost Young
>>> Researchers Network*
>>> *Communications Director - Students and New Professionals Group,
>>> International Society of Biometeorology *
>>>
>>> *M.S. Atmospheric Science - Texas Tech University*
>>> *B.S. Meteorology - Pennsylvania State University*
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=pHuv8gN_LxmdVzqdE0ozhLEwZ41Lw9GQB_n3CEBvO0U&m=bR3jpg-IUuDwbCwoBH1CyJp--PeAuHLKl5WGMdWZtwg&s=i7v223VQTTJLxHzJBuSC4FXEW48XRqosYGVl97poNKo&e=>
>>>
>>>
>>
>
>
> --
>
> *Daniel J. Vecellio*
>
> *Ph.D. Student - Climate Science Lab, **Department of Geography, **Texas
> A&M University*
> *Head Teaching Assistant - Geography 213, Texas A&M University*
> *Education Chair - United States Permafrost Association*
> *ExCom and United States National Representative - Permafrost Young
> Researchers Network*
> *Communications Director - Students and New Professionals Group,
> International Society of Biometeorology *
>
> *M.S. Atmospheric Science - Texas Tech University*
> *B.S. Meteorology - Pennsylvania State University*
>
> Office: 814-1 O&M
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170802/235c0d21/attachment.html