[ncl-talk] Plot vpt-rh

Gurer, Kemal@ARB kemal.gurer at arb.ca.gov
Mon Oct 1 12:18:45 MDT 2018


Hi Bernie,

It is good to hear that the script is working. Good luck for the rest of the averaged plots.

Kemal.

From: Berny Chaimite <bchaimite at gmail.com>
Sent: Monday, October 1, 2018 11:10 AM
To: Gurer, Kemal at ARB <kemal.gurer at arb.ca.gov>
Subject: Re: [ncl-talk] Plot vpt-rh

CAUTION: This email originated from outside of the organization.  Do not click links or open attachments unless you recognize the sender and know the content is safe.

Good evening Kemal!

Thanks, the script already works without any error (in the annex I send an example of the graph that generated). Now I'm working on the script to plot the temporal averages (in a single graphs for the Temp and RH variables).


Thank you.




Best regards,
Berny

 -----
----------------------------------------
Berny A. Chaimite
----------------------------------------
Atmospheric and Climate Sciences Student
Department of Physics
Eduardo Mondlane University

Contact:
Cell: +258 828 957 675(Whatsapp)/846041931/870434233
Tell: +258 2145 6389
Institutional: berny.chaimite at uem.ac.mz<mailto:berny.chaimite at uem.ac.mz>

Maputo, Mozambique


Gurer, Kemal at ARB <kemal.gurer at arb.ca.gov<mailto:kemal.gurer at arb.ca.gov>> escreveu no dia segunda, 1/10/2018 à(s) 19:12:
Hi Bernie,

In my previous email, I only wanted to answer to the first question since the script was not working. Once you see that it is working, then you can tackle plotting the time averages of the variables. For that, you need to comment out the time iteration of reading the data one model output time (hour, day,....) at a time in do loop on line 148 and set “it” time variable to zero when you are reading variables by setting, for example,

u1   = wrf_user_getvar(a,"ua",it)

to

u1   = wrf_user_getvar(a,"ua",0)

The resulting u1 variable here then has the values at all modeling times given in that file, compared to the single time step previously. When your variable has the following dimensions

u1(time,lev,lat,lon),

then you can average your variables at all time steps (first dimension) using:

u1_avg_time = dim_avg_n(u1, 0)

Print the the resulting array using

PrintVarSummary(u1_avg_time)

to confirm that the u1_avg_time variable has only u1(lev,lat,lon) dimensions. Repeat the process for the rest of the variables and now you can plot the vertical cross section at time averaged values of variables at the point of interest in the ncl script.

Kemal.

From: Berny Chaimite <bchaimite at gmail.com<mailto:bchaimite at gmail.com>>
Sent: Sunday, September 30, 2018 11:34 PM
To: Gurer, Kemal at ARB <kemal.gurer at arb.ca.gov<mailto:kemal.gurer at arb.ca.gov>>
Subject: Re: [ncl-talk] Plot vpt-rh

CAUTION: This email originated from outside of the organization.  Do not click links or open attachments unless you recognize the sender and know the content is safe.

Hi Kemal!

1. I made some changes in the script that sent me. I have been trying to run since yesterday, and it gives me the error (which I attached in the form of screenshot). I can not get out of this error.

2. I would like to have a graph that represents the time averages of the differences (ie, for all times of those wrfout_d02 * files) in the vertical profile of temperature and relative humidity. Because my simulations were done for periods of more than a month. So I need to read my variables of interest in all output files and plot in a graphic, do you understand me?

   NB: I was able to make a plot for T2M and RH2M (difference of time averages). now only i need to plot this graph/map of the vertical cross section of Temp. and RH.




Best regards,
Berny Chaimite
 -----
----------------------------------------
Berny A. Chaimite
----------------------------------------
Atmospheric and Climate Sciences Student
Department of Physics
Eduardo Mondlane University

Contact:
Cell: +258 828 957 675(Whatsapp)/846041931/870434233
Tell: +258 2145 6389
Institutional: berny.chaimite at uem.ac.mz<mailto:berny.chaimite at uem.ac.mz>

Maputo, Mozambique


Berny Chaimite <bchaimite at gmail.com<mailto:bchaimite at gmail.com>> escreveu no dia segunda, 1/10/2018 à(s) 02:16:
Thanks a lot.
 -----
----------------------------------------
Berny A. Chaimite
----------------------------------------
Atmospheric and Climate Sciences Student
Department of Physics
Eduardo Mondlane University

Contact:
Cell: +258 828 957 675(Whatsapp)/846041931/870434233
Tell: +258 2145 6389
Institutional: berny.chaimite at uem.ac.mz<mailto:berny.chaimite at uem.ac.mz>

Maputo, Mozambique


Gurer, Kemal at ARB <kemal.gurer at arb.ca.gov<mailto:kemal.gurer at arb.ca.gov>> escreveu no dia domingo, 30/09/2018 à(s) 20:15:
Hi Berry,

The attached ncl script plots the vertical cross section of the difference of a variable obtained from two WRF files over a span of multiple WRF input files. Output of the plot is also attached. The remarks within the script explains the step by step process. You can plot the actual values with a small modification.

Kemal.


From: ncl-talk <ncl-talk-bounces at ucar.edu<mailto:ncl-talk-bounces at ucar.edu>> On Behalf Of Berny Chaimite
Sent: Saturday, September 29, 2018 8:06 PM
To: ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
Subject: [ncl-talk] Plot vpt-rh

CAUTION: This email originated from outside of the organization.  Do not click links or open attachments unless you recognize the sender and know the content is safe.

Good morning, my name is Berny! I want a suggestion of how I can plot the vertical profile of temperature and relative humidity in NCL. My data is output from WRF and is in multiple files (wrfout_d02 *). My final objective is to make a graph of differences in the Vertical Profile of Temperature and Relative Humidity (because I did two simulations with different schemes in wrf) ... etc. I can not see examples of how to do this in ncl, what can you suggest?
the graph that I propose to plot is something similar to the one in annex, which I saw in: http://www.ncl.ucar.edu/Support/talk_archives/2011/1401.html

thank you.


Best regards,
Berny Chaimite
 -----
----------------------------------------
Berny A. Chaimite
----------------------------------------
Atmospheric and Climate Sciences Student
Department of Physics
Eduardo Mondlane University

Contact:
Cell: +258 828 957 675(Whatsapp)/846041931/870434233
Tell: +258 2145 6389
Institutional: berny.chaimite at uem.ac.mz<mailto:berny.chaimite at uem.ac.mz>

Maputo, Mozambique
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181001/e04961ba/attachment.html>


More information about the ncl-talk mailing list