[ncl-talk] Hourly mean csv file

[Toni Klemm]

Sorry, there was a typo in my code, in the ListSetType command. Here it is again:


do year = [from],[to]			; years always have four digits

  do month = 1,12

    if (month .lt. 10) then
      month_zero = "0" + month		; turns "1"..."9" into "01"..."09" to match filename structure
    end if

    if(month .gt. 9) then
      month_zero = month		; leave “10”..."12" unchanged
    end if

    do day = 1,31

      if (day .lt. 10) then
        day_zero = "0" + day		; turns "1"..."9" into "01"..."09" to match filename structure
      end if

      if(hour .gt. 9) then
        day_zero = day			; leave "10…"31" unchanged
      end if

      do hour = 0,23

        if (hour .lt. 10) then
          hour_zero = "0" + hour	; turns "1"..."9" into "01"..."09" to match filename structure
        end if

        if(hour .gt. 9) then
          hour_zero = hour		; leave "10…"23" unchanged
        end if
     

       AOD_files = systemfunc("ls ~/your/file/directory/" + year + "/" + month_zero + "/" + day_zero + "/" + hour_zero + “_AOD.csv")
       AOD_f = addfiles(AOD_files,"r")
       ListSetType(AOD_f,"cat”)


...
proceed with your analysis, e.g., dim_avg_n
...


      end do
    end do
  end do
end do



Toni Klemm, Ph.D.
Postdoctoral Research Associate
Department of Ecosystem Science and Management
College of Agriculture and Life Sciences
Texas A&M University, College Station, TX
Contributor at the Early Career Climate Forum <http://www.eccforum.org/>
www.toni-klemm.de <http://www.toni-klemm.de/> | @toniklemm <http://twitter.com/toniklemm>







> On Aug 8, 2018, at 12:41 PM, Mauro Morichetti <s1073755 at pm.univpm.it> wrote:
> 
> Hi all,
> 
> based on the object, I am comparing a different wrfout with some AEREONET AOD data.
> 
> The AOD data are in csv format:
> 
> yyyy mm dd hh1 AOD
> yyyy mm dd hh1 AOD
> yyyy mm dd hh1 AOD
> …..
> yyyy mm dd hh2 AOD
> yyyy mm dd hh2 AOD
> yyyy mm dd hh2 AOD
> ….
> 
> 
> The values per hour have not the same interval (one hour can have 3 values and the other one can have 5 values etc.). 
> 
> I should need the hourly mean in order to compare with the wrfoutput, but I can not figure out how to solve the problem, any suggestions?
> 
> Thanks,
> 
> Mauro
> ---------
> 
> Eng. Mauro Morichetti
> Department of Industrial Engineering and Mathematical Science
> Marche Polytechnic University, Ancona, Italy
> m.morichetti at pm.univpm.it <mailto:m.morichetti at pm.univpm.it>
> _______________________________________________
> 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/20180808/db11a4da/attachment.html>