Dear Mary,
My problem for extracting grid values is solved.
The Second one ( ascciwrite( "griddata.txt",precip) i cleared just after i
send you the mail itself.
But the array declaration was the problem. I was facing problem with array
manipulation in NCL.
Now its also solved. I got the desired output.
Thanks a ton Madame.
Am attaching the corrected code.
-- A.R. @ +919960462099 Sent from Maharashtra, India On Mon, Oct 12, 2009 at 12:18, Mary Haley <haley_at_ucar.edu> wrote: > I think the problem is that you've set "precip" to an > integer, and then you tried to assign floats to it. > > You need to do this: > > precip =0.0 > > However you are then trying to treate "precip" as an array, > but as far as I can tell, you've made it a scalar. > I think you will need to do this: > > precip = new(9,float) > precip = 0.0 > > Finally, you are doing this: > > asciiwrite("griddata.txt",precip(i)) > > What's going to happen is that it's going to keep writing one value to > the same file, overwriting the original value. > > You want to move "asciiwrite" outside the do loop and call it with: > > asciiwrite("griddata.txt",precip) > > --Mary > > > On Mon, 12 Oct 2009, A.R Ragi wrote: > > *Dear Mary,** >> * >> *Thanks.* >> >> *Still some probs are there.* >> * >> * >> *Now i changed my code like this:* >> >> >> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" >> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" >> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" >> begin >> >> wrf_out = addfile("./wrfout_d01_2000-01-24_12:00:00.nc","r") >> >> rain_con = wrf_user_getvar(wrf_out,"RAINC",-1) >> rain_exp = wrf_user_getvar(wrf_out,"RAINNC",-1) >> rainfall_tot = rain_con + rain_exp >> rainfall_tot_at_description = "Total Precipitation" >> >> printVarSummary(rainfall_tot) >> >> grid = asciiread("data.txt",(/9,3/),"integer") >> print(grid) >> printVarSummary(grid) >> >> ntime = grid(:,0) >> nlat = grid(:,1) >> mlon = grid(:,2) >> >> print(ntime) >> print(nlat) >> print(mlon) >> >> precip =0 >> >> do i = 0, 8 >> >> precip(i) = rainfall_tot(ntime(i),nlat(i),mlon(i)) >> asciiwrite("griddata.txt",precip(i)) >> end do >> >> end >> >> >> >> *The output am getting with the error showing is:* >> >> >> Variable: rainfall_tot >> Type: float >> Total Size: 87600 bytes >> 21900 values >> Number of Dimensions: 3 >> Dimensions and sizes: [5] x [60] x [73] >> Coordinates: >> Number Of Attributes: 1 >> description : Total Precipitation >> >> >> Variable: grid >> Type: integer >> Total Size: 108 bytes >> 27 values >> Number of Dimensions: 2 >> Dimensions and sizes: [9] x [3] >> Coordinates: >> Number Of Attributes: 1 >> _FillValue : -999 >> (0,0) 1 >> (0,1) 22 >> (0,2) 33 >> (1,0) 1 >> (1,1) 23 >> (1,2) 34 >> (2,0) 2 >> (2,1) 12 >> (2,2) 45 >> (3,0) 2 >> (3,1) 15 >> (3,2) 56 >> (4,0) 2 >> (4,1) 15 >> (4,2) 67 >> (5,0) 2 >> (5,1) 15 >> (5,2) 56 >> (6,0) 3 >> (6,1) 23 >> (6,2) 48 >> (7,0) 3 >> (7,1) 22 >> (7,2) 56 >> (8,0) 4 >> (8,1) 13 >> (8,2) 21 >> >> >> Variable: grid >> Type: integer >> Total Size: 108 bytes >> 27 values >> Number of Dimensions: 2 >> Dimensions and sizes: [9] x [3] >> Coordinates: >> Number Of Attributes: 1 >> _FillValue : -999 >> >> >> Variable: ntime >> Type: integer >> Total Size: 36 bytes >> 9 values >> Number of Dimensions: 1 >> Dimensions and sizes: [9] >> Coordinates: >> Number Of Attributes: 1 >> _FillValue : -999 >> (0) 1 >> (1) 1 >> (2) 2 >> (3) 2 >> (4) 2 >> (5) 2 >> (6) 3 >> (7) 3 >> (8) 4 >> >> >> Variable: nlat >> Type: integer >> Total Size: 36 bytes >> 9 values >> Number of Dimensions: 1 >> Dimensions and sizes: [9] >> Coordinates: >> Number Of Attributes: 1 >> _FillValue : -999 >> (0) 22 >> (1) 23 >> (2) 12 >> (3) 15 >> (4) 15 >> (5) 15 >> (6) 23 >> (7) 22 >> (8) 13 >> >> >> Variable: mlon >> Type: integer >> Total Size: 36 bytes >> 9 values >> Number of Dimensions: 1 >> Dimensions and sizes: [9] >> Coordinates: >> Number Of Attributes: 1 >> _FillValue : -999 >> (0) 33 >> (1) 34 >> (2) 45 >> (3) 56 >> (4) 67 >> (5) 56 >> (6) 48 >> (7) 56 >> (8) 21 >> fatal:Assignment type mismatch, right hand side can't be coerced to type >> of >> left hand side >> fatal:Execute: Error occurred at or near line 32 in file newextract.ncl >> >> >> >> Thanks again >> -- >> A.R. @ +919960462099 >> >> >> Sent from Maharashtra, India >> >> On Fri, Oct 9, 2009 at 18:55, Mary Haley <haley_at_ucar.edu> wrote: >> >> >>> On Oct 9, 2009, at 4:15 AM, A.R Ragi wrote: >>> >>> *Dear Mary,* >>> * >>> * >>> *I have one more problem to be solved.* >>> * >>> * >>> *That is for extraxting station data from the given output.* >>> * >>> * >>> *I wrote the script for one grid.* >>> * >>> * >>> >>> precip = rainfall_tot(1,27,33) >>> >>> >>> *At that time I was geeting the output.* >>> *Now i would like to found data for more than 10 station grids (or for >>> 100 >>> grids).* >>> * >>> * >>> *I saved these values of time, lat, lon in a file named "data.txt" >>> 10*3.* >>> * >>> * >>> *Now i tried to read the data using ascii. * >>> *Using:* >>> >>> grid = asciiread("data.txt",(/10,3/),"integer") >>> >>> >>> *Am getting thw whole data in the form:* >>> >>> (0,0) 1 >>> (0,1) 27 >>> (0,2) 33 >>> . >>> . >>> . >>> . >>> (9,2) 56 >>> >>> *Depends on the value i saved in the "data.txt" file.* >>> * >>> * >>> *How can i take (0,0),(0,1),(0,2) values for first time and so on??* >>> * >>> * >>> *Will you please help me to do so???* >>> >>> >>> Do something like this: >>> >>> >>> time = grid(0,:) >>> >>> The other two dimensions are similar, if you need that data too: >>> >>> x1 = grid(1,:) >>> x2 = grid(2,:) >>> >>> Of course, you'll probably want to use something more appropriate than >>> "x1" >>> and "x2". >>> >>> --Mary >>> >>> >>> * >>> * >>> *Thanks in advance.* >>> * >>> * >>> * >>> * >>> *The code i attempted is:* >>> >>> >>> >>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" >>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" >>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" >>> begin >>> >>> wrf_out = addfile("./wrfout_d01_2000-01-24_12:00:00.nc","r") >>> >>> rain_con = wrf_user_getvar(wrf_out,"RAINC",-1) >>> rain_exp = wrf_user_getvar(wrf_out,"RAINNC",-1) >>> rainfall_tot = rain_con + rain_exp >>> rainfall_tot_at_description = "Total Precipitation" >>> >>> printVarSummary(rainfall_tot) >>> >>> >>> ;gridvalues = (/(/1,22,34/),(/1,24,33/),(/3,12,56/)/) >>> ;print(gridvalues) >>> >>> >>> grid = asciiread("data.txt",(/9,3/),"integer") >>> print(grid) >>> >>> precip = 0 >>> >>> ;do i = 1, 3 >>> >>> precip = rainfall_tot(grid(0,0)) >>> >>> >>> asciiwrite("griddata.txt",precip) >>> >>> ;end do >>> end >>> -- >>> A.R. @ +919960462099 >>> >>> >>> Sent from Maharashtra, India >>> >>> >>> >>> >>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.2.0 : Mon Oct 12 2009 - 03:45:17 MDT