<div dir="ltr"><div>Hi all, <br><br></div>I try to plot t-test figure, with significant change &gt;95% between 01-January and 02-Janaury. But i found an error as shown below<br><div><div><br>&quot;fatal:Number of dimensions in parameter (2) of (wrf_contour) is (1), (2) dimensions were expected                                                                       <br>fatal:[&quot;Execute.c&quot;:8578]:Execute: Error occurred at or near line 56 in file t-test2.ncl  &quot;<br><br>Can anyone help me on this? Your assistance is highly appreciated.<br>Thank you.<br><br><br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;<br><br>begin<br>   a = addfile(&quot;./JAN2010/01-05/wrfout_d02_2010-01-01_00:00:00&quot;,&quot;r&quot;)<br>   b = addfile(&quot;./JAN2010/01-05/wrfout_d02_2010-01-02_00:00:00&quot;,&quot;r&quot;)<br><br>  wks = gsn_open_wks(&quot;x11&quot;,&quot;plt_geo2_6&quot;)<br>        gsn_define_colormap(wks,&quot;BlAqGrYeOrReVi200&quot;)            ; Create a plot workstation<br><br>  opts = True                                      ; Set some Basic Plot options<br>  opts@MainTitle = &quot;GEOGRID FIELDS&quot;<br>  opts@InitTime = False                            ; Do not plot time or footers<br>  opts@Footer = False<br><br>;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 01 January;;;;;;;;;;;;;;<br><br>temp1=a-&gt;T2(0,:,:)<br>temp1=0<br>  do i = 0,23<br>    temp1 = temp1 + a-&gt;T2(i,:,:)<br>  end do<br><br>  mean1 = temp1/24<br><br>;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 02 January;;;;;;;;;;;;;;<br><br>temp2=b-&gt;T2(0,:,:)<br>temp2=0<br>  do i = 0,23<br>    temp2 = temp2 + b-&gt;T2(i,:,:)<br>  end do<br><br>  mean2 = temp2/24<br><br>;------------------T-Test: Area with 95% significant difference----------------------------------------------<br><br>aveX = avg(mean1)<br>varX = variance(mean1)<br><br>aveY = avg(mean2)<br>varY = variance(mean2)<br><br>alphat = 100.*(1. - ttest(aveX,varX,24, aveY,varY,24, True, False))<br>aveX = where(alphat.lt.95.,aveX@_FillValue, aveX)<br><br>  res = opts                                       ; Use basic options for this field<br>  res@cnFillOn = True                              ; Create a color fill plot<br>;  res@ContourParameters = (/ -50., 50., 5. /)    ; Set the levels<br><br>  contour = wrf_contour(a,wks,aveX,res)<br><br>  pltres = True                                    ; Set plot options<br>  mpres = True                                     ; Set map options<br>  mpres@mpGeophysicalLineColor      = &quot;Black&quot;      ; Overwrite basic map settings<br>  mpres@mpGridLineColor             = &quot;Black&quot;<br>  mpres@mpLimbLineColor             = &quot;Black&quot;<br>  mpres@mpNationalLineColor         = &quot;Black&quot;<br>  mpres@mpPerimLineColor            = &quot;Black&quot;<br>  mpres@mpUSStateLineColor          = &quot;Black&quot;<br>  plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background<br><br>end<br><br></div></div></div>