[ncl-talk] hovmoller plotting issue using high-resolution data

[Jeremy Berman]

Hi Jack,

Thanks for the response. I double-checked, but the variables read in my script are correct: var1_all (ensemble-mean values) is stored correctly in cntfld (for contour lines), var2_all (sensitivity values) is stored correctly in shdfld (for the color fill or shading), and for that matter tscore_all is stored correctly in stsfld (for the stippling). I should mention that the color fill and contouring in the previously attached image for the most part is what I expected --- with the exception being the issue I brought up.

To reiterate my issue: I don't understand how or why NCL is doing the color fill of low positive sensitivity between the "boxes" of high positive sensitivity (the color fill) (e.g., at 1400 km at 21 h, 1600 km at 24 h, 1800 at 27 h, etc.), rather than just having the signal be persistent between those boxes. Again, there is only data every 3 hours (as labeled on the y-axis) so there are no values in-between those "boxes" but NCL appears to be treating those regions as some other value.

Thus my problem still persists. If you or anyone have any ideas to address this issue, that would be much appreciated. 

Thanks,
Jeremy

On Feb 23, 2016, at 6:13 AM, Jack Ogaja <jack.ogaja at gmail.com>
 wrote:

> It looks like you are just mixing up the variables. If you switch var1_all and var2_all and re-evaluate your level limits you will probably see what you expect.
> 
> Jack
> 
>> On 22 Feb 2016, at 21:51, Jeremy Berman <jdberman at albany.edu> wrote:
>> 
>> Hello,
>> 
>> I am creating a Hovmoller using data that was computed based off of a WRF simulation (15 km grid, output every 3 hours). The data consists of an ensemble-mean isentropic level of potential vorticity, sensitivity values, and a student's t-test; all as a function of time (12 different forecast hours) and distance along the hovmoller (~2400 km). My NCL script reads in this data, and overlays these three variables on a single plot, as a function of forecast hour and distance --- the resulting plot is attached.
>> 
>> The issue is focused on the shading and how NCL is plotting it. The shading corresponds to the sensitivity values, with positive (negative) values in the warm (cool) colors. However, NCL seems to be doing something weird with interpolating between the data points. For example: there are boxes of high positive sensitivity at 1400 km at 21 h, 1600 km at 24 h, 1800 at 27 h, etc.
>> 
>> Thus my question is: Is there something wrong with how I am plotting the data in NCL? Such as with my resources? I am using Rasterfill for example, which has worked in previous applications, but have tried other ones with the same result. Or is there something inherently wrong with using NCL for this application?
>> 
>> Has anyone experienced this plotting problem before, and have any suggestions? Attached is the .nc file (combined_fields.nc), the main NCL script (hovmoller_clean_pv.ncl), and the outputted image (NCL_hovmoller_335K_PV.png). I am using NCL version 6.1.2.
>> 
>> Thank you for your time and consideration.
>> 
>> Many Thanks,
>> 
>> Jeremy Berman
>> PhD Graduate Student
>> University at Albany/SUNY
>> 
>> <NCL_hovmoller_335K_PV.png><hovmoller_clean_pv.ncl><combined_fields.nc>
>> 
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>