[ncl-talk] Extracting all values of rh2 from WRF file and writing it to csv file

Dennis Shea shea at ucar.edu
Wed Jan 2 21:17:16 MST 2019 

This series of questions is a bit frustrating.

You have used NCL quite a bit and you have sent numerous questions.
As you know, ncl-talk always recommends  that users use '*printVarSummary*'
to look at variables. You should do that.... *ALWAYS*.   I don't have your
file but

rh2 = wrf_user_getvar(a,"rh2",-1)       ; *[Time | 673] x [south_north |
71] x [west_east | 82]*
lat=wrf_user_getvar(a,"XLAT",-1)       ;                                 "
lon=wrf_user_getvar(a,"XLONG",-1)  ;                                 "

Your initial question:

lat and lon are *variables* with a '*Time*' dimension.
You can not use multidimensional variables. The following is not correct
syntax.

allRH = rh2(:, *lon, lat)   ; ??????????????????*

Please read the User Manual on *subscripting *and '*do*' loops
  *http://www.ncl.ucar.edu/Document/Manuals/*
<http://www.ncl.ucar.edu/Document/Manuals/>

also

*  http://www.ncl.ucar.edu/Applications/write_ascii.shtml*
<http://www.ncl.ucar.edu/Applications/write_ascii.shtml>

*========*
  rh2 = wrf_user_getvar(a,"rh2",-1)
  *printVarSummary*(rh2)

  lat = wrf_user_getvar(a,"XLAT",-1)
  lon = wrf_user_getvar(a,"XLONG",-1)
  *printVarSummary*(lat)


  dimrh2 = dimsizes(rh2)
  ntim   = dimrh2(0)
  nlat   = dimrh2(1)
  mlon   = dimrh2(2)

  wrftxt = new(nlat, "string")
  wrftxt = ""

  do nl=0,nlat-1

    do nt=0,ntim-1
      do ml=0,mlon-1
         wrftxt(nl) = wrftxt(nl) +sprintf("%9.2f", rh2(nt,nl,ml))
      end do
    end do

  end do


  filo = "rh2.csv"
  system("/bin/rm -f "+filo)
  asciiwrite(filo,wrftxt)


On Wed, Jan 2, 2019 at 8:28 PM Muhammad Omer Mughal <
m.mughal1 at graduate.curtin.edu.au> wrote:

> Hi Rick
>
>
> Thank you for the reply.
>
>
> Yes , I want to write rh2 to the ascii format. The size of the rh2 in this
> particular case is [Time|673]*[south_north|81]*[west_east|72]. I want to
> save the file in the format so that it would have 72 rows and 54513 columns
> i.e. (673*81) for a particular code to read it.
>
>
> Regards
>
>
> Omer
>
>
>
> 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 <m.lynch at 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:* Rick Brownrigg <brownrig at ucar.edu>
> *Sent:* Thursday, 3 January 2019 10:59:30 AM
> *To:* Muhammad Omer Mughal
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] Extracting all values of rh2 from WRF file and
> writing it to csv file
>
> HI,
>
> That error message makes sense, since in a WRF file, the variables XLAT
> and XLONG are 2D (curvilinear).
>
> I guess I'm confused as to what you are trying to do.  If you just simply
> want the variable "rh2" written to an ascii file, instead of this:
>
>     allRH = rh2(:, lon, lat)
>     asciiwrite("rh2.csv",allRH)
>
> what about just:
>
>    asciiwrite("rh2.csv",rh2)
>
> A simple experiment revealed that this will write the 3D array, 1 value
> per line, in "row major" format. This means to read such a file later, you
> must have some way of knowing the shape and dimensionality.
>
> Hope that helps....
> Rick
>
>
> On Wed, Jan 2, 2019 at 7:42 PM Muhammad Omer Mughal <
> m.mughal1 at graduate.curtin.edu.au> wrote:
>
> Dear Rick
>
> Thanks for the reply. I changed the script to remove the SEGV error but
> now I am getting the following error messages
>
>
>
> Please see the script below
>
>
> a=addfile("wrfout_d03_2017-02-01_00:00:00.nc","r")
>
> rh2 = wrf_user_getvar(a,"rh2",-1)
>
> lat=wrf_user_getvar(a,"XLAT",-1)
> lon=wrf_user_getvar(a,"XLONG",-1)
>
>
>   allRH = rh2(:, lon, lat)
>
>   asciiwrite("rh2.csv",allRH)
>
> fatal:Illegal subscript. Subscripts must be scalar or one dimensional
> vectors
>
>
>
>
>
> 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 <m.lynch at 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:* Rick Brownrigg <brownrig at ucar.edu>
> *Sent:* Wednesday, 2 January 2019 10:18:33 PM
> *To:* Muhammad Omer Mughal
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] Extracting all values of rh2 from WRF file and
> writing it to csv file
>
> Hi,
>
> Can you share your script and dataset with me?  I would like to
> investigate why there's a SEGV.
>
> Rick
>
>
> On Wed, Jan 2, 2019 at 4:49 AM Muhammad Omer Mughal <
> m.mughal1 at graduate.curtin.edu.au> wrote:
>
> Dear NCL team
>
> Is there a way to extract all values from WRF output files from all grind
> points and write to an ascii file.
> In my case grid size is 71x82 and the number of times is 673. I tried to
> write using a table format but I am getting a segmentation fault.
>
> I will appreciate help. Kindly note that I have used Lambert conformal
> projection method.
>
> Regards
>
>
>
> 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
> Web | 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
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> _______________________________________________
> 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/20190102/19f340ae/attachment.html>

More information about the ncl-talk mailing list