<div dir="ltr"><div><div><div>Hi ,<br></div>Could someone help me to fix this?<br><br></div>Thanks<br></div>Dz<br>fatal:gc_inout: the lat/lon arrays must have the same dimension sizes<br>fatal:[&quot;Execute.c&quot;:8575]:Execute: Error occurred at or near line 57 in file <br><br> f1       = addfile(&quot;<a href="http://mask_cfsv2.nc">mask_cfsv2.nc</a>&quot;, &quot;r&quot;)<br>printVarSummary (f1)<br>  mrb_lon = f1-&gt;lon<br>  mrb_lat = f1-&gt;lat<br>  nmrb    = dimsizes(mrb_lon)<br><br>  min_mrb_lat = min(mrb_lat)<br>  max_mrb_lat = max(mrb_lat)<br>  min_mrb_lon = min(mrb_lon)<br>  max_mrb_lon = max(mrb_lon)<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>f    = addfile(&quot;Concat_precSum_0Z_198203-201003.nc&quot;, &quot;r&quot;)   ; note the &quot;s&quot; of addfile<br>;T    = f[:]-&gt;precSum                ; read T from all files<br>T    = f-&gt;var(:,{min_mrb_lat:max_mrb_lat},{min_mrb_lon:max_mrb_lon})                ; read T from all files<br>lon1d    = f-&gt;lon_0                ; read T from all files<br>lat1d    = f-&gt;lat_0                ; read T from all files<br>printVarSummary (T)<br>printVarSummary (lat1d)<br>;printVarSummary (lon1d)<br>;printVarSummary (mrb_lon)<br>;print(lat1d-mrb_lat)<br>;print(lon1d-mrb_lon)<br>;return<br> dimp= dimsizes(T)<br>  ntim    = dimp(0)<br>  nlat    = dimp(1)<br>  mlon    = dimp(2)<br>;;;;;;;;;;;;;;;;;;;;;<br><br>;---Create an array and initialize to _FillValue<br>  pmask   = new(dimsizes(T), typeof(T), T@_FillValue)<br>printVarSummary(pmask)<br>  copy_VarCoords(T,pmask)<br>;;;;;;;;;;;;;;;;;;;;;;;;;;<br>;---Keep only data within the polygon<br>;   Use NCL array syntax (:) to propagate to all times<br><br>  do nl=0,nlat-1<br>    do ml=0,mlon-1<br>      if(gc_inout(lat1d(nl),lon1d(ml),mrb_lat,mrb_lon)) then<br>return<br><br></div>