<div dir="ltr"><div>Hi,</div><div><br></div><div>I don't understand why you need to re-order the array dimensions, but I believe the re-ordering is indeed taking place. It is a requirement of gsn_csm_contour_map() that the leftmost dimension corresponds to latitude, while the rightmost in longitude., hence the map you are getting.</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, May 19, 2019 at 9:46 PM Atsuyoshi MANDA <<a href="mailto:am@bio.mie-u.ac.jp">am@bio.mie-u.ac.jp</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
I am trying to reorder the array dimensions in the monthly mean OISST<br>
version 2 dataset, which was downloaded from<br>
<a href="ftp://ftp.cdc.noaa.gov/Datasets/noaa.oisst.v2/sst.mnmean.nc" rel="noreferrer" target="_blank">ftp://ftp.cdc.noaa.gov/Datasets/noaa.oisst.v2/sst.mnmean.nc</a>.<br>
<br>
However, reordering using the named subscripting in my script does not work.<br>
<br>
Please see the attached files:<br>
OIV2LR.sst.1982JUL.png (original data)<br>
OIV2LR.sstre.1982JUL.png (reordered data)<br>
<br>
Any suggestions would be greatly appreciated.<br>
<br>
Here is my script:<br>
;=========================<br>
indir="./"<br>
infle="<a href="http://sst.mnmean.nc" rel="noreferrer" target="_blank">sst.mnmean.nc</a>"<br>
in=indir+"/"+infle<br>
<br>
script_name  = get_script_name()<br>
outdir="BIN_"+systemfunc("basename "+script_name+" .ncl")<br>
<br>
dset="OIV2LR"<br>
MM=7<br>
<br>
month_abbr = (/"","JAN","FEB","MAR","APR","MAY","JUN","JUL","AUG","SEP", \<br>
                    "OCT","NOV","DEC"/)<br>
MMM=month_abbr(MM)<br>
<br>
system("mkdir -vp "+outdir)<br>
<br>
<br>
f=addfile(in, "r")<br>
<br>
;printVarSummary(f)<br>
;print(f)<br>
<br>
lon=f->lon<br>
lat=f->lat<br>
sst=f->sst<br>
time=f->time<br>
<br>
utc_date = cd_calendar(time, 0)<br>
<br>
year   = tointeger(utc_date(:,0))    ; Convert to integer for<br>
month  = tointeger(utc_date(:,1))    ; use sprinti<br>
day    = tointeger(utc_date(:,2))<br>
hour   = tointeger(utc_date(:,3))<br>
minute = tointeger(utc_date(:,4))<br>
second = utc_date(:,5)<br>
date_str = sprinti("%0.2iZ ", hour) + sprinti("%0.2i ", day) + \<br>
           month_abbr(month) + " "  + sprinti("%0.4i", year)<br>
<br>
yyyys=1982<br>
yyyye=2017<br>
<br>
sst_scaled=sst*sst@scale_factor + sst@add_offset<br>
copy_VarAtts(sst,sst_scaled)<br>
copy_VarCoords(sst,sst_scaled)<br>
printVarSummary(sst_scaled)<br>
<br>
sstre=sst_scaled(time|:,lon|:,lat|:)<br>
copy_VarAtts(sst_scaled,sst)<br>
copy_VarCoords(sst_scaled,sstre)<br>
printVarSummary(sstre)<br>
<br>
<br>
<br>
dim=dimsizes(year)<br>
nt=dim(0)<br>
<br>
do n=7,7 ;0,nt-1<br>
<br>
if(year(n) .ge. yyyys .and. year(n) .le. yyyye .and. month(n) .eq. MM)then<br>
print(n+" "+date_str(n))<br>
<br>
<br>
figdir="FIG_"+systemfunc("basename "+script_name+" .ncl")+"_"+dset<br>
system("mkdir -vp "+figdir)<br>
type="png"<br>
<br>
<br>
<br>
figfle=dset+"."+"sstre."+year(n)+MMM+""<br>
fig=figdir+"/"+figfle<br>
<br>
print("figfle="+figfle+"."+type)<br>
<br>
wks = gsn_open_wks(type, fig)<br>
<br>
res=True<br>
  res@cnFillOn              = True               ; turn on color fill<br>
  res@cnLinesOn             = False              ; turn off contour lines<br>
plot1=gsn_csm_contour_map(wks,sstre(n,:,:),res)<br>
<br>
<br>
<br>
figfle=dset+"."+"sst."+year(n)+MMM+""<br>
fig=figdir+"/"+figfle<br>
<br>
print("figfle="+figfle+"."+type)<br>
<br>
wks = gsn_open_wks(type, fig)<br>
<br>
res=True<br>
  res@cnFillOn              = True               ; turn on color fill<br>
  res@cnLinesOn             = False              ; turn off contour lines<br>
plot1=gsn_csm_contour_map(wks,sst(n,:,:),res)<br>
<br>
<br>
print("sst(n,90-30,130)="+sst(n,90-30,130)+" "+lon(130)+"E "+lat(90-30)+"N")<br>
print("sstre(n,130,90+30-1)="+sstre(n,130,90+30-1)+" "+lon(130)+"E<br>
"+sstre&lat(90+30-1)+"N")<br>
<br>
end if<br>
<br>
end do<br>
;=========================<br>
<br>
Thank you,<br>
Atsuyoshi<br>
<br>
-<br>
Atsuyoshi Manda, Ph.D<br>
Associate Professor,<br>
Mie University<br>
1577 Kurimamachiya-cho Tsu city, Mie 514-8507 JAPAN<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>