Re: vector/array writing

Ioana,

You can use ispan(1970,1999,1).

Best,
Ekin

On 04/28/2010 04:42 PM, Ioana Colfescu wrote:
> Hi,
>
> I have a very simple question: I want to write an array with years (1970 1971 .... 1999) and var0 = (/{1970:1:1999}/) doesn't work though {| : } works in selecting special parts of a dimension in a variable.Could someone tell me how should I write this correctly? Sorry the question is very simple but I didn't find in the documentation yet anything.
>
>
> ----- Original Message -----
> From: ncl-talk-request@ucar.edu
> To: ncl-talk@ucar.edu
> Sent: Tuesday, April 27, 2010 10:34:26 PM
> Subject: ncl-talk Digest, Vol 77, Issue 77
>
> Send ncl-talk mailing list submissions to
> ncl-talk@ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk or, via email, send a
> message with subject or body 'help' to
> ncl-talk-request@ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner@ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
>
> Today's Topics:
>
> 1. usage of get1Dindex function ... (u.utku.turuncoglu@be.itu.edu.tr)
> 2. Re: usage of get1Dindex function ... (Dennis Shea)
> 3. Conservative remapping - grid cell corners (Ashley Watson)
> 4. sorting out the files (p s)
> 5. Re: sorting out the files (David Brown)
> 6. Re: sorting out the files (Jonathan Vigh)
> 7. Re: missing_value problem about GODAS monthly data (moon)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 27 Apr 2010 21:22:33 +0300 (EEST)
> From: u.utku.turuncoglu@be.itu.edu.tr
> Subject: usage of get1Dindex function ...
> To: ncl-talk@ucar.edu
> Message-ID:
> <50178.160.75.177.22.1272392553.squirrel@webmail.be.itu.edu.tr>
> Content-Type: text/plain;charset=iso-8859-1
>
> Hi,
>
> I wrote the following procedure to sort x,y coordinate pair based on
> selected one. The opt controls the behavior of the sort procedure. If
> the value is equal to 1 then, it sorts the x values and match the
> correct y
> with it.
>
> undef("sort_xy")
> procedure sort_xy(x:numeric, y:numeric, opt:numeric)
> local x_old, y_old, x_ind, y_ind
> begin ;--- sort x ---
> if (opt .eq. 1) then
> x_old = x
> qsort(x) x_ind = get1Dindex(x_old, x)
> y_old = y
> y = y(x_ind)
> ;--- sort y ---
> else y_old = y
> qsort(y) y_ind = get1Dindex(y_old, y)
> x_old = x
> x = x(y_ind)
> end if
> end
>
> when i run this procedure,
>
> xxx = (/ 3, 5, 2, 10, 6, 4, 3 /)
> yyy = (/ 1, 8, 4, 3 , 7, 6, 2 /)
> sort_xy(xxx, yyy, 1)
>
> it gives the "fatal:Dimension sizes on right hand side of assignment do
> not match dimension sizes of left hand side" error in the "get1Dindex"
> call line. The dimension of the arrays y_old and y or x_old and x is
> same and the second argument is sorted. I could not find the source of
> problem. Is there any special requirement to use "get1Dindex" function?
> How can i
> solve it?
>
> Thanks,
>
> --ufuk
>
>
> -- This message has been scanned for viruses and
> dangerous content by MailScanner, and is
> believed to be clean.
>
>
>
> ------------------------------
>
> Message: 2
> Date: Tue, 27 Apr 2010 13:02:07 -0600
> From: Dennis Shea<shea@ucar.edu>
> Subject: Re: usage of get1Dindex function ...
> To: u.utku.turuncoglu@be.itu.edu.tr
> Cc: ncl-talk@ucar.edu
> Message-ID:<4BD734AF.9010403@ucar.edu>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi Ufuk
>
> Actually the "get1Dindex" fuction notes the following:
>
> "The input arguments should contain only unique values. Duplicate
> entries will lead to a fatal error."
>
>
> I ahve attached a script that contains an example.
> It also (I think) illustrates an error in your
> original script. The fact that subscripts are repeated means that
> the reurn size might be different ... thay are!!!.
>
> Good luck
>
>
> On 04/27/2010 12:22 PM, u.utku.turuncoglu@be.itu.edu.tr wrote:
>
>> Hi,
>>
>> I wrote the following procedure to sort x,y coordinate pair based on
>> selected one. The opt controls the behavior of the sort procedure. If
>> the value is equal to 1 then, it sorts the x values and match the
>> correct y
>> with it.
>>
>> undef("sort_xy")
>> procedure sort_xy(x:numeric, y:numeric, opt:numeric)
>> local x_old, y_old, x_ind, y_ind
>> begin
>> ;--- sort x ---
>> if (opt .eq. 1) then
>> x_old = x
>> qsort(x) x_ind = get1Dindex(x_old, x)
>> y_old = y
>> y = y(x_ind)
>> ;--- sort y ---
>> else
>> y_old = y
>> qsort(y) y_ind = get1Dindex(y_old, y)
>> x_old = x
>> x = x(y_ind)
>> end if
>> end
>>
>> when i run this procedure,
>>
>> xxx = (/ 3, 5, 2, 10, 6, 4, 3 /)
>> yyy = (/ 1, 8, 4, 3 , 7, 6, 2 /)
>> sort_xy(xxx, yyy, 1)
>>
>> it gives the "fatal:Dimension sizes on right hand side of assignment
>> do not match dimension sizes of left hand side" error in the
>> "get1Dindex" call line. The dimension of the arrays y_old and y or
>> x_old and x is same
>> and the second argument is sorted. I could not find the source of
>> problem. Is there any special requirement to use "get1Dindex"
>> function? How can i
>> solve it?
>>
>> Thanks,
>>
>> --ufuk
>>
>>
>>
>
> -- ======================================================
> Dennis J. Shea tel: 303-497-1361 |
> P.O. Box 3000 fax: 303-497-1333 |
> Climate Analysis Section |
> Climate& Global Dynamics Div. |
> National Center for Atmospheric Research |
> Boulder, CO 80307 |
> USA email: shea 'at' ucar.edu |
> ======================================================
> -------------- next part --------------
> An embedded and charset-unspecified text was scrubbed...
> Name: tst_get1Dindex_multiple.ncl
> Url:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100427/3488c4ef/attachment.pl
>
> ------------------------------
>
> Message: 3
> Date: Tue, 27 Apr 2010 22:11:11 +0200
> From: Ashley Watson<ashleywats@gmail.com>
> Subject: Conservative remapping - grid cell corners
> To: ncl-talk@ucar.edu
> Message-ID:
> <n2gf83b83e41004271311hf7e21d7cx58ce312f2864dad7@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear All,
>
> I want interpolate model output to regular grid using conservative
> remapping. As NCL function "area_conserve_remap" works only for regular
> grids, I have to use either SCRIP or CDO. I tried cdo but I found it
> requires corner points of grid cells. In my input file, I have only
> lat2d and lon2d ( grid cell centers). Does someone have ncl script to
> calculate grid cell corners and invoke cdo *remapcon* to do conservative
> remapping. Or
> any hint or function to calculate get corners from grid cell centers?
>
> Thanks in advance
> Regards Ashley
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100427/87ebd2ca/attachment.html
>
> ------------------------------
>
> Message: 4
> Date: Tue, 27 Apr 2010 17:26:30 -0400
> From: p s<prabshr@gmail.com>
> Subject: sorting out the files
> To: ncl-talk@ucar.edu
> Message-ID:
> <s2x8c65c1ab1004271426j7521e99fi92fdcd08295cba90@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi,
> I am trying to read multipled 3B42 hdf 3hourly files and write it out as
> a single netcdf file. The problem I am coming across is due to the way
> the HDF
> filenames are:
>
>
> 3B42.090609.0.6A.HDF
> 3B42.090609.12.6A.HDF 3B42.090609.15.6A.HDF
> 3B42.090609.18.6A.HDF 3B42.090609.21.6A.HDF
> 3B42.090609.3.6A.HDF 3B42.090609.6.6A.HDF
>
>
>
> When I use systemfunc to list the files and read in, they do not appear
> in sequential time order, like 0, 3, 6, 15, 18, 21 because of the size
> of hr
> array in filename
> so the data is stored in netcdf in different time order (0, 12, 15, 18,
> 21, 3, 6), which I want to correct.
>
> Could you please kindly help me to sort out this issue.
>
> Regards,
> Prabhakar
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100427/62392ddb/attachment.html
>
> ------------------------------
>
> Message: 5
> Date: Tue, 27 Apr 2010 16:06:54 -0600
> From: David Brown<dbrown@ucar.edu>
> Subject: Re: sorting out the files
> To: p s<prabshr@gmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID:<6ECFA61B-7A47-4758-AAE2-42864F25B23D@ucar.edu>
> Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes
>
> Hi Prabhakar,
>
> One solution would be to create a file with just the names of the HDF
> files arranged in the order you want:
>
> 3B42.090609.0.6A.HDF
> 3B42.090609.3.6A.HDF 3B42.090609.6.6A.HDF
> 3B42.090609.12.6A.HDF 3B42.090609.15.6A.HDF
> 3B42.090609.18.6A.HDF 3B42.090609.21.6A.HDF
>
> You could use
> ls -1 *.HDF> input_hdf_files
>
> to get the initial list, and then rearrange it manually into the
> correct order. Then instead of doing "ls" as the systemfunc command,
> you could say:
>
> infiles = systemfunc("cat input_hdf_files")
>
> Hope this helps,
> -dave
>
>
> On Apr 27, 2010, at 3:26 PM, p s wrote:
>
>
>> Hi,
>> I am trying to read multipled 3B42 hdf 3hourly files and write it
>> out as a single netcdf file. The problem I am coming across is due
>> to the way the HDF filenames are:
>>
>>
>> 3B42.090609.0.6A.HDF
>> 3B42.090609.12.6A.HDF 3B42.090609.15.6A.HDF
>> 3B42.090609.18.6A.HDF 3B42.090609.21.6A.HDF
>> 3B42.090609.3.6A.HDF 3B42.090609.6.6A.HDF
>>
>>
>>
>> When I use systemfunc to list the files and read in, they do not
>> appear in sequential time order, like 0, 3, 6, 15, 18, 21 because of
>> the size of hr array in filename
>> so the data is stored in netcdf in different time order (0, 12, 15,
>> 18, 21, 3, 6), which I want to correct.
>>
>> Could you please kindly help me to sort out this issue.
>>
>> Regards,
>> Prabhakar
>> _______________________________________________ ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> ------------------------------
>
> Message: 6
> Date: Tue, 27 Apr 2010 16:36:00 -0600
> From: Jonathan Vigh<jvigh@ucar.edu>
> Subject: Re: sorting out the files
> To: p s<prabshr@gmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID:<4BD766D0.70704@ucar.edu>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi Prabhakar,
> If you have a small number of files, then David Brown's solution is
> the easiest But if you have many files, or you are not able to
> physically rename them (for example, if they are on a remote server or
> you don't have write permission), then you have to be clever.
>
> What you can do is to read it the filenames, parse them into their
> separate parts, rewrite the file names in a way that can be sorted by
> NCL's string sort function (sqsort). If you first attach the original
> filenames as an attribute array to the list of new filenames, then you
> sort on the rehashed filenames and the attached attribute array of old
> names will be sorted into the correct order. You can then use that list
> of the old names sorted in the correct order.
>
> This might sounds hard, but it's really not too bad. I've attached a
> script that will do this for you.
>
> Best regards,
> Jonathan
>
>
>
>
>
>
>
>
> p s wrote:
>
>> Hi, I am trying to read multipled 3B42 hdf 3hourly files and write it
>> out as a single netcdf file. The problem I am coming across is due to
>> the
>> way the HDF filenames are:
>>
>>
>> 3B42.090609.0.6A.HDF
>> 3B42.090609.12.6A.HDF 3B42.090609.15.6A.HDF
>> 3B42.090609.18.6A.HDF 3B42.090609.21.6A.HDF
>> 3B42.090609.3.6A.HDF 3B42.090609.6.6A.HDF
>>
>>
>>
>> When I use systemfunc to list the files and read in, they do not
>> appear in sequential time order, like 0, 3, 6, 15, 18, 21 because of
>> the size of hr array in filename
>> so the data is stored in netcdf in different time order (0, 12, 15,
>> 18, 21, 3, 6), which I want to correct.
>>
>> Could you please kindly help me to sort out this issue.
>>
>> Regards,
>> Prabhakar
>> ------------------------------------------------------------------------
>>
>> _______________________________________________ ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
> -------------- next part --------------
> An embedded and charset-unspecified text was scrubbed...
> Name: rename.ncl
> Url:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100427/870f3335/attachment.pl
>
> ------------------------------
>
> Message: 7
> Date: Wed, 28 Apr 2010 10:34:19 +0800
> From: moon<nightmoon35@gmail.com>
> Subject: Re: missing_value problem about GODAS monthly data
> To: ncl-talk@ucar.edu
> Message-ID:<op.vbupfhb8gbhnh9@guanyue>
> Content-Type: text/plain; charset="gbk"
>
> Dear Dennis,
>
> Thanks, _fillvalue worked. But there is anonther problem about panel
> common bar.
> Still using this ncl script, I want to draw a panel figure using a
> common bar. So I add the command
>
> res@lbLabelBarOn = False ; turn off individual cb
>
> and creat the panel,but the common label bar didn't appear. Is there
> any problem?
> __________________________________________________________________________________________
> do nt=0,3
> plot(nt%2,nt/2)=gsn_csm_vector_scalar_map_ce(wks,u(nt,{minLat:maxLat},{minLon:maxLon})
> \ ,v(nt,{minLat:maxLat},{minLon:maxLon})
> \ ,dzdt(nt,{minLat:maxLat},{minLon:maxLon})
> , res)
> end do
> ;************************************************ ; create panel
> ;************************************************ resP = True ; modify
> the panel plot
> ; resP@gsnFrame = False ; don't advance panel
> plot resP@gsnPanelLabelBar = True ; add common colorbar
> resP@lbLabelStride = 2 ; automatically choose
> best stride
> resP@gsnSpreadColors = True ; use full range of color
> map resP@gsnMaximize = True ; fill up the page
> plotone=ndtooned(plot) gsn_panel(wks,plotone,(/2,2/),resP) ; now draw as
> one plot
> ___________________________________________________________________________________
>
> Thanks
>
> moon
>
>
>
>
>
> ? Wed, 28 Apr 2010 00:27:38 +0800?Dennis Shea<shea@ucar.edu> ??:
>
>
>> I see the issue. This is a good example of the "golden rule of
>> data processing" .... Look at your data before processing.
>>
>> Consider:
>>
>> b = addfile(diri + f2, "r")
>> us = b->ucur(:,{5},:,:) ; us is type short
>> printVarSummary(us)
>>
>> Variable: us
>> Type: short<=============================
>> Total Size: 3611520 bytes
>> 1805760 values
>> Number of Dimensions: 3
>> Dimensions and sizes: [time | 12] x [lat | 418] x [lon | 360]
>> Coordinates:
>> time: [35429..35763]
>> lat: [-74..64.999]
>> lon: [ 1..360]
>> Number Of Attributes: 20
>> level : 5
>> missing_value : 32767<=============================
>> [snip]
>> _FillValue : 32766<=============================
>> [snip]
>> add_offset : 0
>> scale_factor : 0.000122074
>> valid_range : ( -2, 2 )
>> actual_range : ( -1.51604, 1.53069 )
>>
>> --- Note the the missing_value and _FillValue attribute are different.
>> If you do a
>>
>> %> ncdump -v ucur ucur.1998.nc | less
>>
>> you will see many 32767 [land] and no 32766.
>>
>> NCL only recognizes (_FillValue) as recommended by Unidata
>> and the Climate and Forecast Convention.
>>
>> In this case, it is the user responsibility to 'manually'
>> make the change. For this dataset
>>
>> us@_FillValue = us@missing_value ; Force over write
>> u = short2flt(us) ; convert to float
>> delete(us) ; no longer needed
>>
>> The attached script shows the complete "fix"
>>
>> Good luck
>>
>> On 04/27/2010 01:10 AM, moon wrote:
>>
>>> Dear NCL users,
>>>
>>> I download GODAS monthly data from the website
>>> "http://www.esrl.noaa.gov/psd/data/gridded/data.godas.html". Draw the
>>> sea surface current over the dzdt, but the value over the land
>>> can't be masked.
>>> I think it's the problem of the missing value, delete the
>>> "missing_value", but it still doesn't work.
>>> Do you have any solution?
>>>
>>> Thank you very much.
>>>
>>> Regards
>>> moon
>>>
>> [SNIP]
>>
>
>
> -- ?? Opera ????????????: http://www.opera.com/mail/
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: godas_2_panel.ncl
> Type: application/octet-stream
> Size: 4965 bytes
> Desc: not available
> Url :
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100428/45ca6642/attachment.obj
>
> ------------------------------
>
> _______________________________________________ ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 77, Issue 77
> ****************************************
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 28 07:56:42 2010