[ncl-talk] Help Required

Thu Oct 26 13:20:34 MDT 2017

Muhammad,

Is there anyway you can provide the files on a Curtin website?

If not, and since your files are too large to transfer, you are going to
have to provide more information for us to be able to help you.

First, I mentioned in a previous email that you should not be using
"addfiles_GetVar". This function was deprecated a long time ago. Before you
do anything else, replace these two lines:

  LH1 = addfiles_GetVar(a,wrf_file1,"T2")
  LH2 = addfiles_GetVar(b,wrf_file2,"T2")

with:

  LH1 = a[:]->T2
  LH2 = b[:]->T2

However. I am confused about why you are using "addfiles" in the first
place, because you are only opening one file in each case.  Are you
planning to modify this script and later open up multiple WRF output files?

Things should still work as you have them, but it is overkill to use
"addfiles" and the "a[:]" and "b[:]" type of syntax if "a" and "b" are only
pointing to one file.

Anyway, after you change those two lines, add printVarSummary commands for
both of these variables.

  printVarSummary(LH1)
  printVarSummary(LH2)

After this line:

  times = a[:]->Times

add another printVarSummary, and a print:

  printVarSummary(times)
  print(times)

rerun the script and email all the output to ncl-talk.

Also, please do an "ncl_filedump" on both files from the UNIX command line
and attach the output files, so we can see what the individual files look
like:
ncl_filedump
/data/muhdomer/Build_WRF/Martilli/glade/u/home/martilli/wrf_3.8/WRFNCEP_updatedUFRAC/wrfout_d05_2016-04-01_00:00:00
>&! file_ncep_output

ncl_filedump
/data/muhdomer/Build_WRF/Martilli/glade/u/home/martilli/wrf_3.8/WRF_Forest/run/wrfout_d05_2016-04-01_00:00:00
>&! file_forest_output

>From all this output, we might be able to help you figure out what's going
on. The problem is that we don't know what the time array is like for your
data, and whether each wrf output files has multiple timesteps or just one
time step.

--Mary

On Wed, Oct 25, 2017 at 10:23 PM, Muhammad Omer Mughal <
m.mughal1 at postgrad.curtin.edu.au> wrote:

> Hi
>
>
> I repeat my question
>
>
> For the script wrf_ensemble ....... has one problem that the time used for
> the first do loop (time loop) is only 24 hours. This means that it reads
> only the first 24 times from the WRF out file instead of the whole 696
> times(29 days ) in the file. Can you kindly suggest that how the do loop
> can be modified to include the 24 hours for each day ?
>
>
> For the script UHIWind2.ncl I am getting higher temperature difference
> during the afternoon time which is  against the diurnal plots that I have
> produced. This means there is some problem on how the difference of
> temperature is being averaged over 29 days for 24 hours. I will be really
> grateful if you can provide help.
>
>
> Muhammad Omer Mughal
> MSc BSc Mechanical Engineering
> PhD  Research Scholar
> Remote Sensing and Satellite Research Group
> Department of Imaging and Applied Physics
> Curtin University
>
> Curtin University
> Tel | +61 8 9266 7962
> Fax | +61 8 9266 2377
> Mobile | 0470 237 525
>
> Email | m.mughal1 at postgrad.curtin.edu.au <m.lynch at curtin.edu.au>
> Web | http://curtin.edu.au
>
> Curtin University is a trademark of Curtin University of Technology.
> CRICOS Provider Code 00301J (WA), 02637B (NSW)
>
>
> ------------------------------
> *From:* Dennis Shea <shea at ucar.edu>
> *Sent:* Thursday, 26 October 2017 9:48:54 AM
> *To:* Muhammad Omer Mughal
> *Cc:* Mary Haley; ncl-talk at ucar.edu
>
> *Subject:* Re: [ncl-talk] Help Required
>
> How big is the file???
> -------------------
> What does
>
> %> ncdump -h ..
> or, equivalently
> %> ncl_filedump ...
>
> show.
>
> ----
> The divisional ftp server was 100% full. That has been fixed. Please retry.
>
>
>
> On Wed, Oct 25, 2017 at 7:06 PM, Muhammad Omer Mughal <
> m.mughal1 at postgrad.curtin.edu.au> wrote:
>
>> Hi
>>
>>
>> I tried to upload but the server provided time out every time.
>>
>>
>> Any how for the script wrf_ensemble ....... has one problem that the time
>> used for the first do loop (time loop) is only 24 hours. This means that it
>> reads only the first 24 times from the WRF out file instead of the whole
>> 699 times(29 days ) in the file. Can you kindly suggest that how the do
>> loop can be modified to include the 24 hours for each day ?
>>
>>
>> For the script UHIWind2.ncl I am getting higher temperature difference
>> during the afternoon time which is  against the diurnal plots that I have
>> produced. This means there is some problem on how the difference of
>> temperature is being averaged over 29 days for 24 hours. I will be really
>> grateful if you can provide help.
>>
>>
>>
>> Muhammad Omer Mughal
>> MSc BSc Mechanical Engineering
>> PhD  Research Scholar
>> Remote Sensing and Satellite Research Group
>> Department of Imaging and Applied Physics
>> Curtin University
>>
>> Curtin University
>> Tel | +61 8 9266 7962
>> Fax | +61 8 9266 2377
>> Mobile | 0470 237 525
>>
>> Email | m.mughal1 at postgrad.curtin.edu.au <m.lynch at curtin.edu.au>
>> Web | http://curtin.edu.au
>>
>> Curtin University is a trademark of Curtin University of Technology.
>> CRICOS Provider Code 00301J (WA), 02637B (NSW)
>>
>>
>> ------------------------------
>> *From:* Mary Haley <haley at ucar.edu>
>> *Sent:* Thursday, 26 October 2017 1:58:52 AM
>> *To:* Muhammad Omer Mughal
>> *Cc:* ncl-talk at ucar.edu
>> *Subject:* Re: [ncl-talk] Help Required
>>
>> Muhammad,
>>
>> It's hard to debug these kind of scripts without being able to see your
>> data.
>>
>> Can you provide the two data files that are referenced in your script?
>> You can use our ftp:
>>
>> http://www.ncl.ucar.edu/ftp_files.shtml
>>
>> --Mary
>>
>>
>>
>>
>> On Mon, Oct 23, 2017 at 9:11 PM, Muhammad Omer Mughal <
>> m.mughal1 at postgrad.curtin.edu.au> wrote:
>>
>>> Dear NCL team
>>>
>>>
>>> I am attaching two codes. The first code i.e. UHI_Wind2.ncl is run in
>>> conjunction with the attached batch script. The intension is to first
>>> obtain difference between the 2m temperature in the two files and then do
>>> an hourly average of this temperature difference for 29 days. There seems
>>> to be an issue since the difference should appear to be maximum during the
>>> night time. The wrfout files are in UTC and adding 8 hours for Singapore
>>> local time gives me the opposite result i.e the maximum difference is
>>> appearing during day time.
>>>
>>>
>>> In the second code attached, I intend to do obtain the difference in
>>> temperature based on the land use. Here again I am having trouble with
>>> using ntimes since the plot appearing is inconsistent with the definition
>>> of UHI i.e a lowest values appearing during the noon time in the diurnal
>>> cycle.
>>>
>>>
>>> I would be really grateful to seek help from the  NCL community
>>>
>>>
>>>
>>> Muhammad Omer Mughal
>>> MSc BSc Mechanical Engineering
>>> PhD  Research Scholar
>>> Remote Sensing and Satellite Research Group
>>> Department of Imaging and Applied Physics
>>> Curtin University
>>>
>>> Curtin University
>>> Tel | +61 8 9266 7962
>>> Fax | +61 8 9266 2377
>>> Mobile | 0470 237 525
>>>
>>> Email | m.mughal1 at postgrad.curtin.edu.au <m.lynch at curtin.edu.au>
>>> Web | http://curtin.edu.au
>>>
>>> Curtin University is a trademark of Curtin University of Technology.
>>> CRICOS Provider Code 00301J (WA), 02637B (NSW)
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>> _______________________________________________
>> 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/20171026/4489ab5f/attachment.html>