Re: Plotting vertical profiles from model output and radio sondes on the same plot

From: BasitAli Khan <BasitAli.Khan_at_nyahnyahspammersnyahnyah>
Date: Thu Mar 21 2013 - 05:52:23 MDT

Hi Mary,
Thanks for digging this out for me. Unfortunately this is not what i need. I want to do a simple xy plot with height/pressure on y axis and magnitude of the desired variable on x-axis. Please see the attached figure. It remains simple as long as I plot data only from one source that is either model output (as the case in the attached figure) or soundings (observed data from radar/lidar etc.). But, it becomes complex when i want to plot data from both sources i.e. model and sounding on the same plot. Why …… because both model and observed values are recorded at different height/pressure levels. For example data point in model output in z dimension are 1000, 900, 850, 700, 500 mb while from observed data are recorded at 1010, 950, 870, 800, 550 mb. Now, there are two options that i can think of handling this problem, one is to interpolate the variable from soundings to the model levels. And second option is to attach pressure/height from sounding as a vertical coordinate variable to the desired sounding variable e.g. potential temperature, and then plot theta from sounding without interpolation to model levels but on the same y-axis and hope that sounding variable will be plotted at the right pressure on the plot. This simply means that i want to add curves for theta, wd, ws, and specific humidity from sounding data to the attached plot. I want to go for the second option(without interpolation) but don't know how should i go about this. To use 2nd option, i even can't use overlay function because the data from model and observed are on different levels. I will greatly appreciate if you could advise some solution, share some script and/or give some idea as how should I handle this.

Thanks and best regards,

----
Basit A. Khan, Ph.D.
Postdoctoral Research Fellow
Division of Physical Sciences & Engineering
Office# 3204, Level 3, Building 1,
King Abdullah University of Science & Technology
4700 King Abdullah Blvd, Box 2753, Thuwal 23955 –6900,
Kingdom of Saudi Arabia.
Office: +966(0)2 808 0276,  Mobile: +966(0)5 0860 3617
E-mail: basitali.khan@kaust.edu.sa<mailto:basitali.khan@kaust.edu.sa>
Skype name: basit.a.khan
From: Mary Haley <haley@ucar.edu<mailto:haley@ucar.edu>>
Date: Wed, 20 Mar 2013 19:00:21 +0300
To: KAUST <basitali.khan@kaust.edu.sa<mailto:basitali.khan@kaust.edu.sa>>
Cc: "ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>>
Subject: Re: Plotting vertical profiles from model output and radio sondes on the same plot
Basit,
See the attached plot and let me know if this will work for you. Note that just for "fun", I reversed the Y axis of the first plot. You may not want this, so just remove the "trYReverse" lines.
I basically created two individual contour plots, and then "added"' the line contour plot to the filled contour plot using gsn_add_annotation.
Doing it this way allows me to retain the Y axis titles and labels.
I believe if you try to use "overlay", then the Y axis titles and labels will be removed from the second plot.
--Mary
[cid:a74907c0-7ddc-4608-9258-eff77c1bcfa9@KAUST.EDU.SA]
On Mar 18, 2013, at 10:27 AM, BasitAli Khan wrote:
Dear all,
I want to plot vertical profiles of theta, wind speed etc. of model output and radio sondes on the same plot. However, model vertical layers  and therefore height levels are different from the radio sondes levels.  Just wondering how would i do it in ncl?
I will greatly appreciate any help.
Best regards,
----
Basit A. Khan, Ph.D.
Postdoctoral Research Fellow
Division of Physical Sciences & Engineering
Office# 3204, Level 3, Building 1,
King Abdullah University of Science & Technology
4700 King Abdullah Blvd, Box 2753, Thuwal 23955 –6900,
Kingdom of Saudi Arabia.
Office: +966(0)2 808 0276,  Mobile: +966(0)5 0860 3617
E-mail: basitali.khan@kaust.edu.sa<mailto:basitali.khan@kaust.edu.sa>
Skype name: basit.a.khan
________________________________
This message and its contents including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

overlay.png
Received on Thu Mar 21 05:52:52 2013

This archive was generated by hypermail 2.1.8 : Fri Mar 22 2013 - 16:45:18 MDT