<div dir="ltr">Hi NCL users, <div><br></div><div>I am using the ncl script below to plot temperature values from GFS data,</div><div>this script works great when i regrid my data to a lower resolution.</div><div>but when i want plotting the whole data (720 x 361 values) and after more than 3 hours, ncl still running without generate any image or error message.</div>
<div>this very strange since NCL is supposed to be a powerful tool. </div><div><br></div><div>Have you any suggest to this problem.</div><div><br></div><div>Thanks a lot.</div><div><br></div><div>Rabah</div><div><br></div>
<div><div>;************************************************</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot; </div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;  </div>
<div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;   </div><div>;************************************************</div><div> a = addfile(&quot;gfs.t00z.pgrb2f00.grib2&quot;,&quot;r&quot;)</div><div>
    lon = a-&gt;lon_0</div><div>    lat = a-&gt;lat_0</div><div>    tt      = short2flt(a-&gt;TMP_P0_L100_GLL0({100000},:,:))</div><div>    tt = (/tt-273/)</div><div><br></div><div>lat2d= new((/361,720/),float)</div><div>
;********* create 2D lat lon </div><div>do l=0,719</div><div> lat2d(:,l)=lat(:)</div><div>end do</div><div><br></div><div>lon2d= new((/361,720/),float)</div><div>do l=0,360</div><div> lon2d(l,:)=lon(:)</div><div>end do</div>
<div>;******************************************</div><div>wks_type=&quot;png&quot;</div><div><br></div><div>  wks_type@wkWidth = 4000</div><div>   wks_type@wkHeight = 4000</div><div><br></div><div>  wks = gsn_open_wks(wks_type,&quot;plot_val&quot;)         </div>
<div>  gsn_define_colormap(wks,&quot;gui_default&quot;)</div><div> </div><div>  res = True               ; plot mods desired</div><div>res@pmTickMarkDisplayMode = &quot;Always&quot;</div><div>  res@mpProjection      = &quot;Mercator&quot;       ; choose projection</div>
<div>  res@mpGridAndLimbOn   = False             ; turn on lat/lon lines</div><div>  res@mpPerimOn         = True             ; turn off box around plot</div><div>  res@mpGridLatSpacingF = 45.               ; spacing for lat lines</div>
<div>  res@mpGridLonSpacingF = 45.               ; spacing for lon lines</div><div>  res@mpFillOn          = False</div><div>   res@mpLimitMode  = &quot;LatLon&quot;</div><div>   res@mpMinLatF            =   85.0</div><div>
    res@mpMaxLatF            =85.0 </div><div>    res@mpMinLonF            = -180</div><div>    res@mpMaxLonF            = 180</div><div> </div><div>  res@cnFillOn          = False              ; color plot desired</div><div>
  res@cnLineLabelsOn    = False             ; turn off contour lines</div><div><br></div><div>  res@vpXF            = 0                 ; make plot bigger</div><div>  res@vpYF            = 1  </div><div>  res@vpWidthF        = 1</div>
<div>  res@vpHeightF       = 1</div><div>map = gsn_csm_map(wks,res)</div><div>  txres               = True</div><div>  txres@txFontQuality = True</div><div>  txres@txFontHeightF = 0.00125</div><div>  txres@txFont        = &quot;helvetica-bold&quot;</div>
<div>  txres@txFontColor =1</div><div>  txres@txFontThicknessF= 10</div><div>  txres@gsnDraw      =  False                   ; do not draw the plot</div><div>  txres@gsnFrame     =  False</div><div>text = gsn_add_text(wks,map,sprintf(&quot;%2.0f&quot;, tt) , lon2d, lat2d,txres)</div>
<div>draw(wks)      ; Now draw map with text strings and</div><div>frame(wks)     ; advance the frame</div><div>end</div>-- <br><div dir="ltr"><div>------------------------------</div>Cordialement,<br>Best regards,<br>Rabah Hachelaf <br>
</div>
</div></div>