[ncl-talk] Fw: script for daily average

Muhammad Omer Mughal m.mughal1 at graduate.curtin.edu.au
Thu Feb 1 18:23:09 MST 2018


Hi Dennis


Thanks for the reply. I think its working now since I receive the following output.


Variable: tc2DayAvg
Type: float
Total Size: 3250800 bytes
            812700 values
Number of Dimensions: 3
Dimensions and sizes:    [Time | 30] x [south_north | 129] x [west_east | 210]
Coordinates:
            Time: [   0.. 696]
Number Of Attributes: 9
  _FillValue :    9.96921e+36
  FieldType :    104
  MemoryOrder :    XY
  description :    Daily Mean: TEMP at 2 M
  units :    K
  stagger :
  coordinates :    XLONG XLAT XTIME
  Time :       0
  NCL_tag :    calculate_daily_values: arith=avg
(0)    Daily Mean: TEMP at 2 M (K) : min=298.312   max=306.723
(0)    -----



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<mailto:m.lynch at curtin.edu.au>
Web | http://curtin.edu.au<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: Friday, 2 February 2018 3:16:17 AM
To: Mary Haley
Cc: Muhammad Omer Mughal; ncl-talk at ucar.edu
Subject: Re: [ncl-talk] Fw: script for daily average

You had the wrong syntax.

@ is for attributes
&  is for coordinates

---
So I try to do the following to correct this error

   tc2 at Time=wrf_times_c(a->Times,0)

and repeat the command but I still receive the same error as above.
---

Please try:
   tc2&Time=wrf_times_c(a->Times,0)

On Thu, Feb 1, 2018 at 9:43 AM, Mary Haley <haley at ucar.edu<mailto:haley at ucar.edu>> wrote:
The calculate_daily_values documentation states that the variable must have a "time" coordinate attached to the leftmost dimension, with units recognized by cd_calendar, i.e., units like "hours since 1-1-1 00:00:0.0"

You are working with WRF data, which doesn't have these kind of units
​for time by default. You may need to use "wrf_times_c" to first convert time to have the correct units, and then attach this new time array to your variable before you call calculate_daily_values:

     Time_0
​  ​
= wrf_times_c(
​
a->Times, 0 )
​        ​
;
​ ​
opt=0 ==>
​ ​
"hours since" initial time on file
​
tc2
​&Time = Time_0​
​

​                        ​
; a
​ttach Time_0 to the "Time" dimension of tc2​


​     ​
tc2DayAvg
​ ​
= calculate_daily_values(tc2, "avg", nDim, False)

Please read the documentation on wrf_times_c to see what other options are available, as the above may not be appropriate for your data.​

http://www.ncl.ucar.edu/Document/Functions/WRF_contributed/wrf_times_c.shtml

--Mary

​

On Wed, Jan 31, 2018 at 6:31 PM, Muhammad Omer Mughal via ncl-talk <ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>> wrote:

Hi


I tried to use the suggestion provided by Dennis below and it shows the error.I am using NCL version 6.4.0. I am repeating the steps here for convenience.

load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

a=addfile("wrfout_d05_2016-04-01_00:00:00","r")

 tc2= wrf_user_getvar(a,"T2",-1)

nDim    = 0

tc2DayAvg= calculate_daily_values(tc2, "avg", nDim, False)


At this stage I get the following error

fatal:No coordinate variable exists for dimension (Time) in variable (x)
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 18125 in file /usr/local/ncl-6.4.0/lib/ncarg/nclscripts/csm/contributed.ncl

So I try to do the following to correct this error

tc2 at Time=wrf_times_c(a->Times,0)

and repeat the command but I still receive the same error as above.


I would appreciate the help depending upon your convenience.

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<mailto:m.lynch at curtin.edu.au>
Web | http://curtin.edu.au<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<mailto:shea at ucar.edu>>
Sent: Thursday, 1 February 2018 5:51:50 AM
To: Muhammad Omer Mughal
Cc: ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
Subject: Re: [ncl-talk] Fw: script for daily average

I am not sure what you were doing but *read the following documentation":

https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml
requires 6.4.0

========
  nDim    = 0
  tc2DayAvg= calculate_daily_values(tc2, "avg", nDim, False)
  tc2DayAvg at description = "Daily Mean: " + tc2 at description

  printVarSummary(tc2DayAvg)
  printMinMax(tc2DayAvg,0)
  print("-----")

On Wed, Jan 31, 2018 at 2:24 PM, Muhammad Omer Mughal via ncl-talk <ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>> wrote:



​
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<tel:+61%208%209266%207962>
Fax | +61 8 9266 2377<tel:+61%208%209266%202377>
Mobile | 0470 237 525

Email | m.mughal1 at postgrad.curtin.edu.au<mailto:m.mughal1 at postgrad.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: Muhammad Omer Mughal
Sent: Wednesday, 31 January 2018 12:17:57 PM
To: ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
Subject: script for daily average


Hi


I wrote the following script to calculate the daily averages


a=addfile("./wrfout_d05_2016-04-01_00:00:00","r")
tc2= wrf_user_getvar(a,"T2",-1)
tc2_dims=dimsizes(tc2)
ntim=tc2_dims(0)
nlat=tc2_dims(1)
nlon=tc2_dims(2)

ndays=ntim/24

tc2D=reshape(tc2,(/ndays,nlat,nlon/))
copy_VarAtts(tc2,tc2D)
tc2_hr_avg=dim_avg_n_Wrap(tc2D,0)

printVarSummary(tc2_hr)

but I am seeing the following error


input dimension sizes cannot conform to dimensions of input array

The Var summaries are as follows
Variable: tc2
Type: float
Total Size: 75526920 bytes
            18881730 values
Number of Dimensions: 3
Dimensions and sizes:    [Time | 697] x [south_north | 129] x [west_east | 210]
Coordinates:
Number Of Attributes: 6
  coordinates :    XLONG XLAT XTIME
  stagger :
  units :    K
  description :    TEMP at 2 M
  MemoryOrder :    XY
  FieldType :    104

Variable: nlat
Type: integer
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes:    [1]
Coordinates:
(0)    129







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<mailto:m.lynch at curtin.edu.au>
Web | http://curtin.edu.au<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<mailto: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<mailto: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/20180202/72c10253/attachment.html>


More information about the ncl-talk mailing list