Re: NARCCAP Resolution Problem

Check the lat / lon values of the new grid you are interpolating to. The values generated by the example you sent are outside of the NARCCAP domain.
On Sep 20, 2010, at 1:41 PM, ncl-talk-request@ucar.edu wrote:

>
> Message: 1
> Date: Mon, 20 Sep 2010 19:38:43 +0000
> From: "eric.holt@huskers.unl.edu" <eric.holt@huskers.unl.edu>
> Subject: FW: NARCCAP Resolution Problem
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <37DA92CD6F883B4FAC3AAC4CC1AAE569DD40@CH1PRD0104MB018.prod.exchangelabs.com>
>
> Content-Type: text/plain; charset="iso-8859-1"
>
> ________________________________
> From: eric.holt@huskers.unl.edu
> Sent: Saturday, September 18, 2010 1:24 PM
> To: ncl-users@ucar.edu
> Subject: NARCCAP Resolution Problem
>
> Hello,
>
> I am trying to simply plot winds from RCM3-cgcm3 and RCM3-gdfl solutions under NARCCAP. I can run my code at 2 degree resolution, and my image is fine (2reswind.pdf). However, when I run my code at 0.5 degrees, the plot displays choppy data, or none at all (pt5reswind.pdf). This has been bugging me for a while now, can anyone provide a tip on why this is, and what I can do? It would be much appreciated. Thank You!
>
> Here is the error when the 0.5 degree doesn't plot:
>
> warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible :[errno=1101]
>
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>
>
> begin
>
> file1 = "vas_RCM3_cgcm3_2038010103.nc"
> file2 = "uas_RCM3_cgcm3_2038010103.nc"
> file3 = "vas_RCM3_cgcm3_2066010103.nc"
> file4 = "uas_RCM3_cgcm3_2066010103.nc"
>
>
> print("Starting file 1")
> f1 = addfile(file1,"r")
> print("Starting file 2")
> f2 = addfile(file2,"r")
> print("Starting file 3")
> f3 = addfile(file1,"r")
> print("Starting file 4")
> f4 = addfile(file2,"r")
>
> ;===============================
> ;2038-2041
> print("Starting U")
> U10_o = f2->uas(0:5,:,:)
>
>
> print("Starting V")
> V10_o = f1->vas(0:5,:,:)
>
>
> ;Total Wind Calculation
> tot10 = V10_o ; Transfering meta data
> tot10 = sqrt(U10_o^2+V10_o^2)
>
> delete(V10_o)
> delete(U10_o)
>
> ;;2D Arrays of Lat and Lon
> lat2d = f1->lat
> lon2d = f1->lon
>
>
>
> y1 = dim_avg_Wrap( tot10(yc|:,xc|:,time|:))
> printVarSummary(y1)
>
> delete(tot10)
>
> ;;2D Arrays of Lat and Lon
> lat2d = f1->lat
> lon2d = f1->lon
> print(" min="+min(lat2d)+" max="+max(lat2d))
> print(" min="+min(lon2d)+" max="+max(lon2d))
>
> dimlc = dimsizes(lat2d)
> nlat = dimlc(0)
> mlon = dimlc(1)
>
> ; NLAT = 180
> ; MLON = 360
> NLAT = 046
> MLON = 090
>
> dlat = 0.5 ; output lat, in this case dlat=2
> lat = ispan ( 0,NLAT-1,1 )*dlat ; 0-90N, by dlat=0.5, so 0, 0.5N 1N, 1.5N. 90N
> lat!0 = "lat"
> lat@units = "degrees_north"
> lat&lat = lat
>
> dlon = 0.5 ; output lon, in this case dlon=2
> lon = ispan ( 0,MLON-1,1 )*dlon + 180. ; -180W-180E by dlon=2, so -180, -178W... 180E
> lon!0 = "lon"
> lon@units = "degrees_east"
> lon&lon = lon
>
>
> ; ydiff = y2
> ; ydiff = y2 - y1
>
> X2D = triple2grid(ndtooned(lon2d),ndtooned(lat2d),ndtooned(y1(:,:)), lon,lat, False)
> ; X2D = triple2grid2d(ndtooned(lon2d),ndtooned(lat2d),ndtooned(y1(:,:)), lon,lat, False)
> X2D!0 = "lat" ; triple2grid regrids based on a nearest neighbor routine.
> X2D&lat = lat
> X2D!1 = "lon"
> X2D&lon = lon
> X2D@units = y1@units
> printVarSummary(X2D)
> print(" min="+min(X2D)+" max="+max(X2D))
>
> print(" min="+min(lat2d)+" max="+max(lat2d))
> print(" min="+min(lon2d)+" max="+max(lon2d))
>
> ; print (lat)
> ; print (lon)
>
> ;********************************************
> ; create plots
> ;********************************************
> wks = gsn_open_wks ("pdf", "pt5reswind") ; open workstation
> gsn_define_colormap (wks,"gui_default") ; choose color map
>
> sres = True ; plot mods desired for 2x2 array
> sres@gsnAddCyclic = False
> sres@cnFillOn = True ; color fill
> sres@cnLinesOn = False ; no contour lines
> sres@cnLineLabelsOn = False ; no contour labels
> sres@gsnSpreadColors = True ; use total colormap
> sres@gsnSpreadColorStart = 4
> sres@gsnSpreadColorEnd = -1
> sres@cnInfoLabelOn = False ; no contour info label
>
> sres@mpMinLatF = min(lat2d)+1
> sres@mpMaxLatF = max(lat2d)-1
> sres@mpMinLonF = min(lon2d)+1
> sres@mpMaxLonF = max(lon2d)-1
> ; sres@mpMinLonF = 220. ; min(lon2d) - 1
> ; sres@mpMaxLonF = 300. ; max(lon2d) - 1
> sres@mpCenterLonF = 270 ; Centers the plot at 180
>
> sres@gsnFrame = False
> sres@gsnDraw = False
> sres@lbLabelBarOn = False
>
> panres = True ; set up panel resources
> panres@gsnPaperOrientation = "portrait"
> panres@gsnMaximize = True
> panres@gsnPanelLabelBar = True
> ; panres@gsnPanelYWhiteSpacePercent= 3.0
> panres@gsnPanelYWhiteSpacePercent= 0.5
>
> plot = new(1,graphic)
>
> sres@gsnLeftString = "Regridded to 0.5x0.5"
> ; sres@gsnCenterString = "Gridded Precipitation "
> ; plot(1) = gsn_csm_contour_map_ce(wks,X2D,sres) ; Draw new 2x2 grid
> plot = gsn_csm_contour_map_ce(wks,X2D,sres) ; Draw new 2x2 grid
>
> gsn_panel(wks,plot,(/1,1/),panres)
>
> end
>
>
>
> Eric Holt
> Department Of Geosciences
> Graduate Student
> University of Nebraska
> 119 Bessey Hall
> eric.holt@huskers.unl.edu
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100920/1c6d30ca/attachment.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: 2reswind.pdf
> Type: application/pdf
> Size: 165480 bytes
> Desc: 2reswind.pdf
> Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100920/1c6d30ca/attachment.pdf
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: pt5reswind.pdf
> Type: application/pdf
> Size: 87407 bytes
> Desc: pt5reswind.pdf
> Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20100920/1c6d30ca/attachment-0001.pdf
>
> ------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 82, Issue 32
> ****************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Sep 20 14:03:53 2010