[ncl-talk] Array reordering in OISST

Dennis Shea shea at ucar.edu
Mon May 20 11:16:53 MDT 2019


 I *think* I understand what you want to do. You use:

[A]  This following uses *standard subscripting (indexing):*

print("sst(n,*90-30,130*)="+sst(n,*90-30,130)*+" "+lon(130)+"E
"+lat(90-30)+"N")
print("sstre(n,130,90+30-1)="+sstre(n,130,90+30-1)+" "+lon(130)+"E
"+sstre&lat(90+30-1)+"N")


[B] I think you intended to use *coordinate subscripting
<https://www.ncl.ucar.edu/Document/Language/subscript_coord.shtml> *which
uses the following syntax:
*{...}*
A 'best practices' programming style is to use variables rather than
hard-wiring numbers. Generally, this makes ithe intended use more clear.

   LATS = 30
   LATN = 90
   LON = 130

    SST = sst(n,*{*LATS:LATN*}*,*{*LON*}*)
    printVarSummary(SST)

or, say
   LATS = 30
   LATN = 90
   LONW = 110
   LONE  = 150

    SST = sst(n,*{*LATS:LATN*}*,*{*LON*W:LONE}*)
    printVarSummary(SST)

=================================
[1] Your file variable is:
*     short *sst(time, lat, lon)

     The coordinates are:

       Coordinates:
            time: [66443..80078]
            *lat: [89.5..-89.5]               ; <=== North-to-South
ordering*
            lon: [0.5..359.5]

[2] NCL has the *short2flt *
<https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml>function
that unpacks an maintains meta data.

    sst =* short2flt*( f->sst )
*; import and unpack the variable maintaining meta data*
*    printVarSummary
<https://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml>(sst)*
    *printMinMax*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml>
(sst,0)

[3]  It is not necessary to explicitly read the following coordinate
variables from the file. NCL imports the variable as a *data structure*
['object'] that contains the array values and meta data including
coordinate information associated with the variable.

time = sst*&*time
lat    = sst*&*lat
lon   = sst*&*lon

[4] There are different methods for accessing subsets of a variable:
classic indexing and coordinate variable indexing. See above.

[5] Generally, *NCL graphics* *are aware of coordinate order.* Hence, no
reordering is necessary.

---
[6]




On Sun, May 19, 2019 at 9:46 PM Atsuyoshi MANDA <am at bio.mie-u.ac.jp> wrote:

> Hi,
>
> I am trying to reorder the array dimensions in the monthly mean OISST
> version 2 dataset, which was downloaded from
> ftp://ftp.cdc.noaa.gov/Datasets/noaa.oisst.v2/sst.mnmean.nc.
>
> However, reordering using the named subscripting in my script does not
> work.
>
> Please see the attached files:
> OIV2LR.sst.1982JUL.png (original data)
> OIV2LR.sstre.1982JUL.png (reordered data)
>
> Any suggestions would be greatly appreciated.
>
> Here is my script:
> ;=========================
> indir="./"
> infle="sst.mnmean.nc"
> in=indir+"/"+infle
>
> script_name  = get_script_name()
> outdir="BIN_"+systemfunc("basename "+script_name+" .ncl")
>
> dset="OIV2LR"
> MM=7
>
> month_abbr = (/"","JAN","FEB","MAR","APR","MAY","JUN","JUL","AUG","SEP", \
>                     "OCT","NOV","DEC"/)
> MMM=month_abbr(MM)
>
> system("mkdir -vp "+outdir)
>
>
> f=addfile(in, "r")
>
> ;printVarSummary(f)
> ;print(f)
>
> lon=f->lon
> lat=f->lat
> sst=f->sst
> time=f->time
>
> utc_date = cd_calendar(time, 0)
>
> year   = tointeger(utc_date(:,0))    ; Convert to integer for
> month  = tointeger(utc_date(:,1))    ; use sprinti
> day    = tointeger(utc_date(:,2))
> hour   = tointeger(utc_date(:,3))
> minute = tointeger(utc_date(:,4))
> second = utc_date(:,5)
> date_str = sprinti("%0.2iZ ", hour) + sprinti("%0.2i ", day) + \
>            month_abbr(month) + " "  + sprinti("%0.4i", year)
>
> yyyys=1982
> yyyye=2017
>
> sst_scaled=sst*sst at scale_factor + sst at add_offset
> copy_VarAtts(sst,sst_scaled)
> copy_VarCoords(sst,sst_scaled)
> printVarSummary(sst_scaled)
>
> sstre=sst_scaled(time|:,lon|:,lat|:)
> copy_VarAtts(sst_scaled,sst)
> copy_VarCoords(sst_scaled,sstre)
> printVarSummary(sstre)
>
>
>
> dim=dimsizes(year)
> nt=dim(0)
>
> do n=7,7 ;0,nt-1
>
> if(year(n) .ge. yyyys .and. year(n) .le. yyyye .and. month(n) .eq. MM)then
> print(n+" "+date_str(n))
>
>
> figdir="FIG_"+systemfunc("basename "+script_name+" .ncl")+"_"+dset
> system("mkdir -vp "+figdir)
> type="png"
>
>
>
> figfle=dset+"."+"sstre."+year(n)+MMM+""
> fig=figdir+"/"+figfle
>
> print("figfle="+figfle+"."+type)
>
> wks = gsn_open_wks(type, fig)
>
> res=True
>   res at cnFillOn              = True               ; turn on color fill
>   res at cnLinesOn             = False              ; turn off contour lines
> plot1=gsn_csm_contour_map(wks,sstre(n,:,:),res)
>
>
>
> figfle=dset+"."+"sst."+year(n)+MMM+""
> fig=figdir+"/"+figfle
>
> print("figfle="+figfle+"."+type)
>
> wks = gsn_open_wks(type, fig)
>
> res=True
>   res at cnFillOn              = True               ; turn on color fill
>   res at cnLinesOn             = False              ; turn off contour lines
> plot1=gsn_csm_contour_map(wks,sst(n,:,:),res)
>
>
> print("sst(n,90-30,130)="+sst(n,90-30,130)+" "+lon(130)+"E
> "+lat(90-30)+"N")
> print("sstre(n,130,90+30-1)="+sstre(n,130,90+30-1)+" "+lon(130)+"E
> "+sstre&lat(90+30-1)+"N")
>
> end if
>
> end do
> ;=========================
>
> Thank you,
> Atsuyoshi
>
> -
> Atsuyoshi Manda, Ph.D
> Associate Professor,
> Mie University
> 1577 Kurimamachiya-cho Tsu city, Mie 514-8507 JAPAN
> _______________________________________________
> 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/20190520/bdedeab1/attachment.html>


More information about the ncl-talk mailing list