[ncl-talk] Daily averages from 6 hourly data
Dennis Shea
shea at ucar.edu
Thu Apr 18 21:55:02 MDT 2019
This is a subtle error which I can not seem to isolate. The issue is tha
dimensions:
lon = 192 ;
lat = 94 ;
*level = 1* ;
time = UNLIMITED ; // (1460 currently)
variables:
short uwnd(time, *level, *lat, lon) ;
===
It works fine on variable with multiple level.
I speculate that NCL's dimension reduction that leads to the error.
===
Fortunately, there is a work-around:
Change:
uwnd_practice = short2flt(fil1->uwnd) ; [time | 1460]
x *[level
| 1]* x [lat | 94] x [lon | 192]
to
uwnd_practice = short2flt(fil1->uwnd*(:,0,:,:)) *; [time | 1460] x
[lat | 94] x [lon | 192]
The result will be:
uwnd_daily: [time | 365] x [lat | 94] x [lon | 192]
On Thu, Apr 18, 2019 at 6:42 AM Barry Lynn <barry.h.lynn at gmail.com> wrote:
> Hi:
>
> Please see here -- the explanation of the error message might be of help.
>
>
> https://www.ncl.ucar.edu/Document/Language/error_messages.shtml#AssignMismatch
>
> On Thu, Apr 18, 2019 at 3:18 PM Nkese Mc Shine <Nkese.McShine at sta.uwi.edu>
> wrote:
>
>> Hi,
>>
>> This line causes the error.
>>
>>
>>
>> uwnd_daily = calculate_daily_values(uwnd_practice, "avg", 0, opt)
>>
>>
>>
>> Regards,
>>
>> Nkese.
>>
>>
>>
>> *From:* Barry Lynn <barry.h.lynn at gmail.com>
>> *Sent:* Thursday, April 18, 2019 8:16 AM
>> *To:* Nkese Mc Shine <Nkese.McShine at sta.uwi.edu>
>> *Cc:* ncl-talk at ucar.edu
>> *Subject:* Re: [ncl-talk] Daily averages from 6 hourly data
>>
>>
>>
>> Hello:
>>
>>
>>
>> What line causes your error?
>>
>>
>>
>> On Thu, Apr 18, 2019 at 2:52 PM Nkese Mc Shine <Nkese.McShine at sta.uwi.edu>
>> wrote:
>>
>> Good morning,
>>
>> I have 6 hourly wind data and I am trying to get the daily averages
>> however I am getting and error when I use calculate_daily_values. I have
>> ncl 6.6.2 installed, can anyone help me? Below is the script used and the
>> output.
>>
>> Script:
>> fil1 = addfile("/media/wind/NCEP-DOE_Reanalysis_2/
>> uwnd.10m.gauss.1979.nc","r")
>> uwnd_practice = short2flt(fil1->uwnd)
>>
>>
>> printVarSummary(uwnd_practice)
>> printMinMax(uwnd_practice, True)
>>
>> opt = True
>> opt at nval_crit = 4
>>
>> uwnd_daily = calculate_daily_values(uwnd_practice, "avg", 0, opt)
>> printVarSummary(uwnd_daily)
>> printMinMax(uwnd_daily, True)
>>
>>
>> Output:
>> Variable: uwnd_practice
>> Type: float
>> Total Size: 105400320 bytes
>> 26350080 values
>> Number of Dimensions: 4
>> Dimensions and sizes: [time | 1460] x [level | 1] x [lat | 94] x [lon
>> | 192]
>> Coordinates:
>> time: [1569072..1577826]
>> level: [10..10]
>> lat: [88.542..-88.542]
>> lon: [ 0..358.125]
>> Number Of Attributes: 19
>> long_name : 6-Hourly Forecast of U-wind at 10 m
>> unpacked_valid_range : ( -120, 120 )
>> actual_range : ( -38.15, 46.84 )
>> units : m/s
>> precision : 2
>> least_significant_digit : 1
>> GRIB_id : 33
>> GRIB_name : U GRD
>> var_desc : u-wind
>> dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)
>> level_desc : 10 m
>> statistic : Individual Obs
>> parent_stat : Other
>> standard_name : eastward_wind
>> _FillValue_original : -32767
>> _FillValue : -32767
>> missing_value_original : 32766
>> valid_range : ( -120, 120 )
>> missing_value : -32767
>> (0)
>> (0) 6-Hourly Forecast of U-wind at 10 m (m/s) : min=-38.15 max=46.84
>>
>> fatal:["NclVar.c":1390]:Assignment type mismatch, right hand side can't
>> be coerced to type of left hand side
>> fatal:["NclAtt.c":262]:Attribute assignment type mismatch.
>>
>> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 19034 in
>> file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
>>
>> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 35 in
>> file dailyaverages2.ncl
>>
>>
>> Regards,
>> Nkese.
>>
>> CONFIDENTIALITY: This email (including any attachments) may contain
>> confidential, proprietary and/or privileged information. Any duplication,
>> copying, distribution, dissemination, transmission, disclosure or use in
>> any manner of this email (including any attachments) without the
>> authorisation of the sender is strictly prohibited. If you receive this
>> email (including any attachments) in error, please notify the sender and
>> delete this email (including any attachments) from your system. Thank you.
>>
>> _______________________________________________
>> 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 Associate Scientist, 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
>> CONFIDENTIALITY: This email (including any attachments) may contain
>> confidential, proprietary and/or privileged information. Any duplication,
>> copying, distribution, dissemination, transmission, disclosure or use in
>> any manner of this email (including any attachments) without the
>> authorisation of the sender is strictly prohibited. If you receive this
>> email (including any attachments) in error, please notify the sender and
>> delete this email (including any attachments) from your system. Thank you.
>>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Associate Scientist, 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
>
> _______________________________________________
> 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/20190418/ddae9639/attachment.html>
More information about the ncl-talk
mailing list