<div dir="ltr">Hi<br><br>I am trying to add a variable from netcdf files, and make an average plot. My scripts is mentioned below, but the process gets killed. Maybe because of memory. <br><br>Can you please suggest a "efficient" way of writing this? Please help.<br><br>Script ----------<br>begin<br>a = addfile("../wrfout_d01_2012-05-01_00:00:<a href="http://00.nc">00.nc</a>","r")<br><br> it        = 0     ; first time step<br>  hgt       = wrf_user_getvar(a,"HGT_M",it)    ; Terrain elevation<br>  hgt@lat2d = wrf_user_getvar(a,"XLAT",it)   ; latitude/longitude<br>  hgt@lon2d = wrf_user_getvar(a,"XLONG",it)  ; required for plotting<br><br>  wks = gsn_open_wks("png","try1")<br><br>    FILES = systemfunc ("ls ../wrfout_d01_2012-05*") ; file paths<div><br> numFILES = dimsizes(FILES)<br> aa    = addfiles (FILES+".nc", "r")     ;add multiple files<br>;;;;;;;;;;;;<br><br>no = aa[:]->no<br>no_avg=dim_avg_n(no,0)<br>no2 = aa[:]->no2<br>no2_avg=dim_avg_n(no2,0)<br>oli = aa[:]->oli<br>oli_avg=dim_avg_n(oli,0)<br>csl = aa[:]->csl<br>csl_avg=dim_avg_n(csl,0)<br>olt = aa[:]->olt<br>olt_avg=dim_avg_n(olt,0)<br>xyl = aa[:]->xyl<br>xyl_avg=dim_avg_n(xyl,0)<br>ald = aa[:]->ald<br>ald_avg=dim_avg_n(ald,0)<br>hc8 = aa[:]->hc8<br>hc8_avg=dim_avg_n(hc8,0)<br>hcho = aa[:]->hcho<br>hcho_avg=dim_avg_n(hcho,0)<br>ol2 = aa[:]->ol2<br>ol2_avg=dim_avg_n(ol2,0)<br>tol = aa[:]->tol<br>tol_avg=dim_avg_n(tol,0)<br>hc5= aa[:]->hc5<br>hc5_avg=dim_avg_n(hc5,0)<br>hc3= aa[:]->hc3<br>hc3_avg=dim_avg_n(hc3,0)<br>ket= aa[:]->ket<br>ket_avg=dim_avg_n(ket,0)<br><br>Mpan_avg=no_avg+no2_avg+oli_avg+csl_avg+olt_avg+xyl_avg+ald_avg+hc8_avg+hcho_avg+ol2_avg+tol_avg+hc5_avg+hc3_avg+ket_avg<br>Mpan_avg=Mpan_avg*1000<br>Mpan_avg2D=Mpan_avg(0,:,:)<br>Mpan_avg2D@lat2d = wrf_user_getvar(a,"XLAT",it)   ; latitude/longitude<br>Mpan_avg2D@lon2d = wrf_user_getvar(a,"XLONG",it)  ; required for plotting<br>--some plotting options here--<br>contour = gsn_csm_contour_map(wks,Mpan_avg2D,res)<br>draw(contour)<br>  frame(wks)<br> <div>Thanks<div>Kim</div></div></div></div>