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

[Dennis Shea]

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170802/76549b55/attachment.html