[ncl-talk] How to define date for WRF-chem to get average daily of a variable ?

Dennis Shea shea at ucar.edu
Fri Oct 7 12:01:40 MDT 2022


  My script follows.. There were no problems. The output is attached.

++++++++++++++++++++++++++++++++++++++++++++++++++++

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

  f = addfile("setareh_chem.nc","r")
  print(f)                              ; same as "nmcdump -h"

  Times = f->Times                      ; [Time | 41] x [DateStrLen | 19]
 (character)

; For fun:   Examine "Times"  in different formats

  Time_0 = wrf_times_c( Times, 0 )      ; "hours since" initial time on
file   (double)
                                                             ;  Time_0 is
recognized by cd_calendar
 ;Time_1 = wrf_times_c( Times, 1 )      ; "hours since 1901-01-01 00:00:00"
   (double)
 ;Time_2 = wrf_times_c( Times, 2 )      ; yyyymmddhhmnss
    (double)
  Time_3 = wrf_times_c( Times, 3 )      ; yyyymmddhh
    (integer)

 ;print(Times)
 ;print(Time_0)
 ;print(Time_1)
 ;print(Time_2)
  print(Time_3)                         ; easy on the eyes
  print("------------------------------------")

  a2           = f->TAUAER2             ; (Time, south_north, west_east)
  a2&Time = Time_0                      ; associate the 'Time' coordinate
with the variable using  standard NCL & syntax
  printVarSummary(a2)
  print("------------------------------------")

  opt_cdv = True             ; option for calculate_daily_values  (cdv)
  opt_cdv at nval_crit = 4      ; default is 1

  a2Day = calculate_daily_values (a2, "avg", 0, opt_cdv)
  printVarSummary(a2Day)
  printMinMax (a2Day,1)

On Fri, Oct 7, 2022 at 9:28 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:

> Dear Dennis,
> I sent you the file via transfer.pcloud.com. I could not send it here due
> to the large size of the file.
> Please check the script with the file and let me know your idea.
> Many thanks for your help.
> Best wishes,
>
> On Fri, Oct 7, 2022 at 3:18 PM Dennis Shea <shea at ucar.edu> wrote:
>
>> I do not understand 'x' either
>>
>> "I wonder about variable x; I did not define such a variable."
>>
>> ncl 0> f = addfile(".....","r")  ; I do not have your file
>> ncl 1> print(f)                                             ; like ncdump
>> -h
>> ncl 2> Times = f->Times                            ; Times(Time,
>> DateStrLen)    (type character)
>> ncl 4> Time = wrf_times_c( Times, 0 )      ; "hours since" initial time
>> on file   (double); units recognized by cd_calendar
>> ncl 5> print(Time)
>>
>> ncl 7> T = f->T
>> ncl 8> printVarSummary(T)
>>
>> ncl 9> T&Time = Time                              ; associate time
>> coordinate
>> ncl 10> printVarSummary(T)
>>
>> Variable: T
>> Type: float
>> Total Size: 19683000 bytes
>>             4920750 values
>> Number of Dimensions: 4
>> Dimensions and sizes: [Time | ??] x [bottom_top | 27] x [south_north |
>> 81] x [west_east | 90]
>> Coordinates:
>>             Time: [   0..  ?]
>> Number Of Attributes: 5
>>   FieldType : 104
>>   MemoryOrder : XYZ
>>   description : perturbation potential temperature (theta-t0)
>>   units : K
>>   stagger :
>> ncl 11> Tavg = calculate_daily_values (T, "avg", 0, True)
>> ncl 12> printVarSummary(Tavg)
>>
>> Variable: Tavg
>> Type: float
>> Total Size: 1574640 bytes
>>             393660 values
>> Number of Dimensions: 4
>> Dimensions and sizes: [Time | ...] x [bottom_top | 27] x [south_north |
>> 81] x [west_east | 90]
>> Coordinates:
>>             Time: [   0..  12]
>> Number Of Attributes: 8
>>   _FillValue : 9.96921e+36
>>   FieldType : 104
>>   MemoryOrder : XYZ
>>   description : perturbation potential temperature (theta-t0)
>>   units : K
>>   stagger :
>>   Time :   0
>>   NCL_tag : calculate_daily_values
>>
>> On Fri, Oct 7, 2022 at 5:16 AM Setareh Rahimi <setareh.rahimi at gmail.com>
>> wrote:
>>
>>> Dear Dennis,
>>> So many thanks for your advice. I made corrections to my script as you
>>> suggested, and attempted to calculate and plot the average AOD for
>>> 2011-06-28 but I still get errors. I attached the script and the errors.
>>>
>>> As you can see from the error message, NCL says :fatal: No coordinate
>>> variable exists for dimension (Time) in variable (x)
>>> I wonder about variable x; I did not define such a variable.
>>>
>>> Please advise me on how to sort this issue out.
>>> Many thanks in advance,
>>> Best wishes,
>>>
>>> On Fri, Oct 7, 2022 at 6:02 AM Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> As noted in the documentation, *calculate_daily_values*
>>>> <https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml>
>>>> requires that the variable have a **time coordinate** that is recognized by
>>>> cd_calendar.
>>>>
>>>>   Times = f->Times                    ; Times(Time, DateStrLen)    (type character)
>>>>   Time = *wrf_times_c*( Times, 0 )      ; "hours since" initial time on file   (double); units recognized by cd_calendar
>>>>
>>>>      a2           = f->TAUAER   ; (Time, south_north, west_east)
>>>>      a2&Time = Time   ; associate the 'Time' coordinate with the
>>>> variable using  standard NCL & syntax
>>>>      printVarSummary(a2)
>>>>      print("------------------------------------")
>>>>
>>>>    a2Day = calculate_daily_values (a2, "avg", 0, opt)
>>>>    printVarSummary(a2Day)
>>>>    printMinMax (a2Day,1)
>>>>
>>>>
>>>>
>>>> On Wed, Oct 5, 2022 at 9:30 AM Setareh Rahimi via ncl-talk <
>>>> ncl-talk at mailman.ucar.edu> wrote:
>>>>
>>>>> Dear NCL users,
>>>>>
>>>>> I am trying to plot the daily average of AOD (aerosol optical depth)
>>>>> from WRF-chem outputs. However, I can not define the time so that NCL plots
>>>>> the daily average of AOD.(the script has been attached, and the model
>>>>> output is too large and I just used ncdump to show the header of some
>>>>> varibles ).
>>>>> I tried the following commands:
>>>>>
>>>>>
>>>>> times = wrf_user_getvar(f,"times",-1)
>>>>> printVarSummary(times)
>>>>> print(times)
>>>>> ymdh = cd_calendar(a2&times, -2)
>>>>>   print(ymdh)
>>>>>
>>>>> but faced an error.
>>>>> Would you please kindly advise me on how can I calculate the daily
>>>>> average of AOD from the hourly output of WRF-chem?
>>>>> Many thanks in advance,
>>>>> Best wishes,
>>>>> --
>>>>> S.Rahimi
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at mailman.ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>>
>
> --
> S.Rahimi
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221007/d7ad6fc5/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: setareh_chem.OUTPUT
Type: application/octet-stream
Size: 8503 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221007/d7ad6fc5/attachment.obj>


More information about the ncl-talk mailing list