[ncl-talk] sub: mask-area specifier error...

Adam Phillips asphilli at ucar.edu
Fri Dec 9 14:53:03 MST 2016


Hi Adv,
Add the following to your script:
res at mpDataBaseVersion     = "Ncarg4_1"
That should allow you to specify the states in mpMaskAreaSpecifiers.

Also, the valid settings of mpAreaMaskingOn  are True or False. In your
case I believe you want to set it to True.

Hope that helps!
Adam


On Fri, Dec 9, 2016 at 12:57 PM, Adv <advita6 at gmail.com> wrote:

> Hi,
>
> could someone help me to fix this error?
>
> Thank you,
> Adv
>
>
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wyoming"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "idaho"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "utah"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wisconsin"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wyoming"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "idaho"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "utah"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wisconsin"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wyoming"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "idaho"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "utah"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wisconsin"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wyoming"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "idaho"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "utah"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "wisconsin"
> warning:MapV40DHUpdateDrawList: invalid boundary specification string:
> "iowa"
>
>
>
> load "./contributed.ncl_beta_640"
>
> begin
> ;***************************************************
> strll  = asciiread("stationmac.txt", -1, "string") ;
>   print(strll)   ; STATION","STATION_NAME", "LAT", "LON"
>   print("=====")
>
> sloc = str_get_field(strll , 2, ",")
>   lat  = tofloat( str_get_field(strll , 3, ",")  )
>   lon  = tofloat( str_get_field(strll , 4, ",")  )
>   print(sloc +" : "+lat+" "+lon)
> mnlat=(min(lat))
> mxlat=(max(lat))
> mnlon=(min(lon))
> mxlon=(max(lon))
> ;print(lon)
> print(mnlon)
> print(mxlon)
> print(mnlat)
> print(mxlat)
> ;iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
> mlon = 360
> dlon = -90.35 / mlon
> lonn  = fspan (-116, (mlon - 1) * dlon, mlon)
> print(lonn)
> mlat = 90
> dlat = 49.0 / mlat
> latt  = fspan (40.01, (mlat - 1) * dlat, mlat)
> print(latt)
> diri="./"
> z1 = asciiread(diri+"Temp_Mean1965-2005",-1,"float")
> print(z1)
> printVarSummary(z1)
> opt    = True
> rscan  = (/8,5,1 /)
> grid   = obj_anal_ic_deprecated(lon, lat, z1, lonn, latt, rscan, True)
> printVarSummary(grid)
> ;print(grid)
> print(min(grid))
> print(max(grid))
> ;************************************************
> ; Specify geographical region and time span (year-month start and end
> ;************************************************
> ;latS=-90
> ;latN     =  90
> ;lonL     =  0;360-116.5
> ;lonR     =  358;360-90
> ;  latS     = 40
> ;  latS     = 0
>  ; latN     =  49
>  ; lonL     =  244;360-116.5
>  ; lonR     =  270;360-90
>
>   ymStrt   = 196501
>   ymLast   = 200512
>
>   pltType  = "x11"
>   pltName  = "regress"
>   pltTitle = "Globe: "+(ymStrt/100)+"-"+(ymLast/100)
>
> ;************************************************
> ; Read from netCDF file: variable is type short...unpack
> ;************************************************
>    version= "v2"         ; "v2", "v2c"
>
>    diri   = "./"
> ;   fili   = "air.sig995.mon.mean."+version+".nc"
>    fili1   = "Netcdfmean.nc"
>    fili   = "air.mon.mean.nc"
>    ff      = addfile(diri+fili1,"r")
>    t      = ff->mean(:,:)
>    printVarSummary(t)                            ; [time| 720]x[lat|
> 91]x[lon| 180]
> printMinMax(t,0)
>    f      = addfile(diri+fili,"r")
> f1    = addfile("tas_Amon_CCSM4_historical_r3i1p1_185001-200512.nc",
> "r")   ; note the "s" of addfile
> f2    = addfile("Ecearth_tas_320x160x1961-2005.nc", "r")   ; note the "s"
> of addfile
> f3    = addfile("CanCM4_Temp_320x160x1961-2005.nc", "r")   ; note the "s"
> of addfile
>
>    YYYYMM = cd_calendar( f->time, -1)
>    YYYYMMa = cd_calendar( f1->time, -1)
>    YYYYMMb = cd_calendar( f2->time, -1)
>    YYYYMMc = cd_calendar( f3->time, -1)
>
>    iStrt  = ind(YYYYMM.eq.ymStrt)
>    iStrta  = ind(YYYYMMa.eq.ymStrt)
>    iStrtb  = ind(YYYYMMb.eq.ymStrt)
>    iStrtc  = ind(YYYYMMc.eq.ymStrt)
>    iLast  = ind(YYYYMM.eq.ymLast)
>    iLasta  = ind(YYYYMMa.eq.ymLast)
>    iLastb  = ind(YYYYMMb.eq.ymLast)
>    iLastc  = ind(YYYYMMc.eq.ymLast)
>           ;                                      ; monthly data
>    xt      = f->air(:,::-1,:)
>    printVarSummary(xt)                            ; [time| 720]x[lat|
> 91]x[lon| 180]
>   xt1      = f1->tas(:,:,:)
>    printVarSummary(xt1)
>   xt2      = f2->t2(:,:,:)
>    printVarSummary(xt2)
>   xt3      = f3->t2(:,:,:)
>    printVarSummary(xt3)
> ;   x      = xt(iStrt:iLast,{latN:latS},{lonL:lonR})
>    x      = xt(iStrt:iLast,:,:)
>    x1      = xt1(iStrta:iLasta,:,:)
>    x2      = xt2(iStrtb:iLastb,:,:)
>    x3      = xt3(iStrtc:iLastc,:,:)
>    delete(xt)
>    delete(xt1)
>    delete(xt2)
>    delete(xt3)
>    x1      = x1-273.15
>    x2      = x2-273.15
>    x3      = x3-273.15
>    x1 at units= "degC"
>    x2 at units= "degC"
>    x3 at units= "degC"
>    printVarSummary(x)                            ; [time| 720]x[lat|
> 91]x[lon| 180]
>    printVarSummary(x1)                            ; [time| 720]x[lat|
> 91]x[lon| 180]
>    printVarSummary(x2)                            ; [time| 720]x[lat|
> 91]x[lon| 180]
>    printVarSummary(x3)                            ; [time| 720]x[lat|
> 91]x[lon| 180]
> ;************************************************
> ; Calculate the Annual Mean (slopes)
> ;************************************************
> rc=  dim_avg_n(x,0)
> rc1=  dim_avg_n(x1,0)
> rc2=  dim_avg_n(x2,0)
> rc3=  dim_avg_n(x3,0)
>    printVarSummary(rc)
>    printVarSummary(rc)
>    printVarSummary(rc1)
>    printVarSummary(rc2)
>    printVarSummary(rc3)
> copy_VarCoords(x(0,:,:),rc)
> copy_VarCoords(x1(0,:,:),rc1)
> copy_VarCoords(x2(0,:,:),rc2)
> copy_VarCoords(x3(0,:,:),rc3)
>    printMinMax(rc,0)
>    printMinMax(rc1,0)
>    printMinMax(rc2,0)
>    printMinMax(rc3,0)
> ;==============================O===========;
> ; interpolate to new grid
> ;=========================================;
>   newgrid  = f2fsh_Wrap(grid(:,:),(/73,144/))
>   newgrid1  = f2fsh_Wrap(grid(:,:),(/192,288/))
>   newgrid2  = f2fsh_Wrap(grid(:,:),(/160,320/))
>   newgrid3  = f2fsh_Wrap(grid(:,:),(/64,128/))
> y= newgrid-rc
> y1= newgrid1-rc1
> y2= newgrid2-rc2
> y3= newgrid3-rc3
> copy_VarCoords(rc,y)
> copy_VarCoords(rc1,y1)
> copy_VarCoords(rc2,y2)
> copy_VarCoords(rc3,y3)
>  mask_specs = (/"Iowa","Wyoming","Idaho","Utah","Wisconsin","Iowa"/)
>
> ;************************************************
> ; plotting parameters
> ;************************************************
>    wks  = gsn_open_wks("x11","Tempdifference_Mean_"+ymStrt+"-"+ymLast)
> ; specifies a ps plot
> gsn_define_colormap(wks,"matlab_jet")
>     plot = new(4,graphic)
>    res                       = True
>    res at gsnMaximize           = True             ; make large
>    res at cnFillOn              = True             ; turn on color
>    res at cnLinesOn             = False            ; turn off contour lines
>    res at cnLineLabelsOn        = False            ; turn off contour line
> labels
>    res at cnInfoLabelOn       = False
> res at lbLabelBarOn        = False           ; turn off individual cb's
>
> res at gsnTickMarksOn = False
> res at mpGridAndLimbOn       = False
>
>    res at cnLevelSelectionMode  = "ManualLevels"   ; set manual contour
> levels
>    res at cnMinLevelValF        =  -5            ; set min contour level
>    res at cnMaxLevelValF        =   5.0            ; set max contour level
>    res at cnLevelSpacingF       =   0.5            ; set contour interval
>    res at gsnAddCyclic          =  False
>    res at mpFillOn              = False            ; turn off default
> background gray
>   res at mpCenterLonF          = 180
>   res at mpMinLatF             =   40.
>   res at mpMaxLatF             =   49.
>   res at mpMinLonF             = -116
>   res at mpMaxLonF             = -90
> res at gsnStringFont   = "helvetica-bold"             ; instead of using
> txFontHeightF or gsnStringFontHeightF
>
> res at mpGeophysicalLineThicknessF = 3.0
> res at mpGeophysicalLineColor = "Black"; (/22/)
> ;res at mpNationalLineThicknessF  = 3.0
> res at mpOutlineBoundarySets = "GeophysicalAndUSStates" ; add state
> boundaries
> res at mpNationalLineColor  = res at mpGeophysicalLineColor
> res at mpUSStateLineThicknessF = 3.0
> res at mpUSStateLineColor  = res at mpGeophysicalLineColor
> ; res at gsnCenterString       = year(0)+"-"+year(nyrs-1)
>   res at gsnLeftString   = ""               ; add the gsn titles
>   res at gsnRightString  = ""
> ;   res at tiMainString          = "Trend, Reanalysis-2m
> Temperature(~S~o~N~C)"    ; fili
> ;res at tiMainFontHeightF  = 0.02
> ;************************************************
> ; Turn on masking, and set resource indicating
> ; areas to mask.
> ;************************************************
>   res at mpAreaMaskingOn       = 1
>   res at mpMaskAreaSpecifiers  = mask_specs
>
> ;:::::::::::::::::::::::::::;;;;;:::::::::::
> res at gsnCenterStringFontHeightF     = 0.035
> res at gsnCenterString     = "Obs-Reanal"
>    plot(0) = gsn_csm_contour_map_ce(wks,y,res)
> res at gsnCenterString     = "Obs-CCSM4"
>    plot(1) = gsn_csm_contour_map_ce(wks,y1,res)
> res at gsnCenterString     = "Obs-EC-Earth"
>    plot(2) = gsn_csm_contour_map_ce(wks,y2,res)
> res at gsnCenterString     = "Obs-CanCM4"
>    plot(3) = gsn_csm_contour_map_ce(wks,y3,res)
> ;************************************************
> ; create panel
> ;************************************************
>   resP                  = True                   ; modify the panel plot
>   resP at gsnFrame         = False                  ; don't advance the
> frame, so we can use gsn_text_ndc
>   resP at gsnPanelLabelBar = True                   ; add common colorbar
>   resP at lbLabelFontHeightF = 0.015                ; set font height of
> Label Bar labels
> resP at gsnPanelBottom   = 0.2                    ; shrink panel plot by
> setting bottom edge of plot
>   resP at gsnPanelTop      = 0.9                    ; shrink panel plot by
> setting top edge of plot
>  resP at pmLabelBarWidthF  = 0.3                 ; make thinner
>   resP at pmLabelBarHeightF = 0.05
>   resP at gsnPanelYWhiteSpacePercent = 5.           ; increase spacing along
> Y-axis between panel plots
>  resP at cnLinesOn       = False                   ; turn off contour lines
> ;   resP at cnFillOn        = True               ; color plot desired
> ;  resP at lbTitleOn        = True                  ; turn on title
> resP at lbTitleFont    = "helvetica-bold"
>  resP at lbLabelFontQuality      = "high"     ; label bar font
>   resP at lbLabelFont      = "helvetica-bold"     ; label bar font
>
>   resP at lbTitleFontHeightF= .013                 ;
>   resP at lbTitleString    = "~S~o~N~C"
>    resP at cnLineLabelsOn  = False
>
>   resP at txString   = "Mean_Temp_"+ymStrt+"-"+ymLast+""
>   gsn_panel(wks,plot,(/4,1/),resP)               ; now draw as one plot
>  frame(wks)
>    end
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161209/488515e9/attachment.html 


More information about the ncl-talk mailing list