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