[ncl-talk] Cross sections w/ smoothed terrain.... New issue???

Mon Apr 9 16:02:00 MDT 2018

I couldn't FTP the file over but you should be able to download it from this link:

https://drive.google.com/open?id=1_zNgcHezyGVjA3aLu4xXj70Ho6c7Wwsr

Thank you!
Kacie

---

Kacie N. Shourd, M.S.
드림

Staff Research Scientist
Division of Atmospheric Science
Desert Research Institute, Reno, NV
Kacie.Shourd at dri.edu
Work: (775) 674-7049
________________________________
From: Bill Ladwig [ladwig at ucar.edu]
Sent: Monday, April 09, 2018 1:32 PM
To: Kacie Shourd
Cc: ncl-talk (ncl-talk at ucar.edu)
Subject: Re: [ncl-talk] Cross sections w/ smoothed terrain.... New issue???

Hi Kacie,

Can you upload one of your WRF files by following the instructions here:

https://www.ncl.ucar.edu/ftp_files.shtml

Once you have the file uploaded, can you please send me the name of the file? I'll take a look at your problem.

Thanks,

Bill

On Mon, Apr 9, 2018 at 12:05 PM, Kacie Shourd <Kacie.Shourd at dri.edu<mailto:Kacie.Shourd at dri.edu>> wrote:
Hi everyone,

It's been a while since anyone posted to this thread, but I tried the suggested fix (modifying my WRFUserARW.ncl line 382) to no avail! There is still a good chunk of data missing above the terrain. See the attached output images. The wider image is the default with "dz=0.01*z_max" and the skinnier image is with "dz=0.00001*z_max".

Attempted fix:

"Have you modified line ~382 in WRFUserARW.ncl yet?  dz can be decreased from the default value of 0.01:

...else
        z_max = max(z)
        z_min = 0.
; Changed to:
dz = 0.00001 * z_max
        ; From:
        ;dz = 0.01 * z_max
        nlevels = tointeger( z_max/dz )
        z_var2d = new( (/nlevels/), typeof(z))
        z_var2d(0) = z_min
     end if"

There are quite noticeable differences between the two plots, but no improvement to the terrain. I am close to giving up on this, but thought I would reach out one more time. I have changed other things in the "wrf_intrpd_3d" function also without luck. As mentioned before, the missing data near the terrain was not an issue in a colleague's RIP plot, so I know it is not the data.

Best,
Kacie

---

Kacie N. Shourd, M.S.
드림

Staff Research Scientist
Division of Atmospheric Science
Desert Research Institute, Reno, NV
Kacie.Shourd at dri.edu<mailto:Kacie.Shourd at dri.edu>
Work: (775) 674-7049
________________________________
From: Scott Capps [scapps at atmosdatasolutions.com<mailto:scapps at atmosdatasolutions.com>]
Sent: Tuesday, March 27, 2018 4:47 PM
To: Kacie Shourd
Cc: ncl-talk (ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>)
Subject: Re: [ncl-talk] Cross sections w/ smoothed terrain.... New issue???

Have you modified line ~382 in WRFUserARW.ncl yet?  dz can be decreased from the default value of 0.01:

...else
        z_max = max(z)
        z_min = 0.
; Changed to:
dz = 0.002 * z_max
        ; From:
        ;dz = 0.01 * z_max
        nlevels = tointeger( z_max/dz )
        z_var2d = new( (/nlevels/), typeof(z))
        z_var2d(0) = z_min
     end if

On Tue, Mar 27, 2018 at 4:36 PM, Kacie Shourd <Kacie.Shourd at dri.edu<mailto:Kacie.Shourd at dri.edu>> wrote:
Hi NCL friends,

Attached is an image of the current output I am getting trying to plot a vertical cross section with WRF data in NCL. Also attached is my code, and I can upload the WRF data file for anyone who needs it.

I have used this script tons of times in the past with no issue (i.e. no "missing" data near the ground, or at least the terrain was smoothed to cover that up) and am hoping to figure out what's wrong.

Could it be the data file I have? If yes, is there a simple way I can extrapolate the RH values to below the ground?

Is it the version of NCL I am running?? I used to run this script in v6.4 but am using v6.3 at the moment.

A colleague of mine made a plot with the same WRF file in RIP and did not have this issue. I would prefer to plot this up in NCL if possible as I think the plots come out nicer and this is going into a document for some pretty important folks...

Thanks,
Kacie

---

Kacie N. Shourd, M.S.
드림

Staff Research Scientist
Division of Atmospheric Science
Desert Research Institute, Reno, NV
Kacie.Shourd at dri.edu<mailto:Kacie.Shourd at dri.edu>
Work: (775) 674-7049<tel:(775)%20674-7049>

PUBLIC RECORDS NOTICE: In accordance with NRS Chapter 239, this email and responses, unless otherwise made confidential by law, may be subject to the Nevada Public Records laws and may be disclosed to the public upon request.

_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk<https://linkprotect.cudasvc.com/url?a=http://mailman.ucar.edu/mailman/listinfo/ncl-talk&c=E,1,FvUIVgD7EH9G1hQSxvgiu9p8mf_4k0kEcTABmKT6pVBB-FhMwRwnE4yuShuTwXhNqoezMpeYqatYxdSV83pHVm7uAsGw71owTjHKYqPf7g,,&typo=1>

--
Scott Capps
Principal
Atmospheric Data Solutions<https://linkprotect.cudasvc.com/url?a=http://www.atmosphericdatasolutions.com&c=E,1,ii7FITA5DOfDJKepvgGfOku-AGID3DnAJtS5iWd_oBwb_PIKnfs-1V1XTcN2ZiESPh4EGugxUHZb3Vv0bmuUu2W6sV9LogTv0Nw4PDekVE-r4rFAJI2KDnz5K5wQ&typo=1>
scapps at atmosdatasolutions.com<mailto:scapps at atmosdatasolutions.com>
cell: (949) 910-4385

_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto: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/20180409/217e0692/attachment.html>