<div dir="ltr"><div>The file is about 5GB.<br></div><div><br></div><div>ncl-talk does not need the entire file.</div><div><br></div><div>DO you have the netCDF operators available?</div><div> <br><div>%> ncks -v Times,TAUAER2,TAUAER3,XLAT,XLONG   wrfout_d01_2011-07-01_00:00:12  <a href="http://setareh_chem.nc">setareh_chem.nc</a></div><div><br></div><div>Further, compress  '<a href="http://setareh_chem.nc">setareh_chem.nc</a>'  using gzip</div><div><br></div><div>%>  gzip <a href="http://setareh_chem.nc">setareh_chem.nc</a></div><div><br></div><div><br></div><div>This [setareh_chem.nc.gz] will be substantially smaller.   <br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 7, 2022 at 9:28 AM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Dennis,</div><div>I sent you the file via<span style="background-color:rgb(217,210,233)"> <a href="http://transfer.pcloud.com" target="_blank">transfer.pcloud.com</a></span>. I could not send it here due to the large size of the file. </div><div>Please check the script with the file and let me know your idea.</div><div>Many thanks for your help.</div><div>Best wishes,<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 7, 2022 at 3:18 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>I do not understand 'x' either</div><div><br></div><div><span style="background-color:rgb(213,166,189)"><span style="background-color:rgb(255,255,255)">"I wonder about variable x; I did not define such a variable."</span></span></div><div><br></div><div>ncl 0> f = addfile(".....","r")  ; I do not have your file<br></div><div>ncl 1> print(f)                                             ; like ncdump -h <br></div><div>ncl 2> Times = f->Times                            ; Times(Time, DateStrLen)    (type character)</div><div>ncl 4> Time = wrf_times_c( Times, 0 )      ; "hours since" initial time on file   (double); units recognized by cd_calendar</div><div>ncl 5> print(Time)</div><div><br></div><div>ncl 7> T = f->T    <br>ncl 8> printVarSummary(T)</div><div><br></div><div>ncl 9> T&Time = Time                              ; associate time coordinate<br>ncl 10> printVarSummary(T)<br><br>Variable: T<br>Type: float<br>Total Size: 19683000 bytes<br>            4920750 values<br>Number of Dimensions: 4<br>Dimensions and sizes:    [Time | ??] x [bottom_top | 27] x [south_north | 81] x [west_east | 90]<br>Coordinates: <br>            Time: [   0..  ?]<br>Number Of Attributes: 5<br>  FieldType :  104<br>  MemoryOrder :   XYZ<br>  description :   perturbation potential temperature (theta-t0)<br>  units :       K<br>  stagger : <br>ncl 11> Tavg = calculate_daily_values (T, "avg", 0, True)  <br>ncl 12> printVarSummary(Tavg) <br><br>Variable: Tavg<br>Type: float<br>Total Size: 1574640 bytes<br>            393660 values<br>Number of Dimensions: 4<br>Dimensions and sizes:       [Time | ...] x [bottom_top | 27] x [south_north | 81] x [west_east | 90]<br>Coordinates: <br>            Time: [   0..  12]<br>Number Of Attributes: 8<br>  _FillValue :       9.96921e+36<br>  FieldType :     104<br>  MemoryOrder :   XYZ<br>  description :   perturbation potential temperature (theta-t0)<br>  units :       K<br>  stagger : <br>  Time :        0<br>  NCL_tag :     calculate_daily_values</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 7, 2022 at 5:16 AM Setareh Rahimi <<a href="mailto:setareh.rahimi@gmail.com" target="_blank">setareh.rahimi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Dennis,</div><div>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. <br></div><div><br></div><div>As you can see from the error message, NCL says :<span style="background-color:rgb(255,229,153)">fatal: No coordinate variable exists for dimension (Time) in variable (x)</span></div><div><span style="background-color:rgb(213,166,189)"><span style="background-color:rgb(255,255,255)">I wonder about variable x; I did not define such a variable.</span><br></span></div><div><br></div><div>Please advise me on how to sort this issue out.</div><div>Many thanks in advance, <br></div><div>Best wishes,<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 7, 2022 at 6:02 AM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>As noted in the documentation, <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_daily_values.shtml" target="_blank"><b>calculate_daily_values</b></a> requires that the variable have a **time coordinate** that is recognized by cd_calendar. <br><pre>  Times = f->Times                    ; Times(Time, DateStrLen)    (type character)
  Time = <b>wrf_times_c</b>( Times, 0 )      ; "hours since" initial time on file   (double); units recognized by cd_calendar</pre></div><div>     a2           = f->TAUAER   ; (Time, south_north, west_east)</div><div>     a2&Time = Time   ; associate the 'Time' coordinate with the variable using  standard NCL & syntax<br></div><div>     printVarSummary(a2)</div><div>     print("------------------------------------")     <br><pre>   a2Day = calculate_daily_values (a2, "avg", 0, opt)
   printVarSummary(a2Day)
   printMinMax (a2Day,1)</pre></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Oct 5, 2022 at 9:30 AM Setareh Rahimi via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear NCL users,</div><div><br></div><div>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 ).<br></div><div>I tried the following commands:</div><div> <br></div><div><br></div><div><span style="background-color:rgb(255,229,153)">times = wrf_user_getvar(f,"times",-1)<br>printVarSummary(times)<br>print(times)<br>ymdh = cd_calendar(a2&times, -2)<br>  print(ymdh)</span></div><div><br></div><div>but faced an error.</div><div>Would you please kindly advise me on how can I calculate the daily average of AOD from the hourly output of WRF-chem?</div><div>Many thanks in advance,</div><div>Best wishes,<br></div><div>-- <br><div dir="ltr">S.Rahimi<br><br></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr">S.Rahimi<br><br></div>
</blockquote></div>