[ncl-talk] Help Required

Barry Lynn barry.h.lynn at gmail.com
Thu Oct 26 15:12:43 MDT 2017


Muhammed:

Are you reading a single file or multiple files?

I've never heard of a WRF file whose time index ends at one day, but the
run goes for 29 days.

Barry

On Thu, Oct 26, 2017 at 10:20 PM, Mary Haley <haley at ucar.edu> wrote:

> 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
>>>
>>>
>>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>


-- 
Barry H. Lynn, Ph.D
Senior Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
US 914 432 3108
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171027/a262a8b7/attachment.html>


More information about the ncl-talk mailing list