<div dir="ltr"><div><div><div><div>Hello I'm having trouble doing a timeaverage of pm2.5 since it has 4 dimensions (time,e_we, s_sn, bot_top)<br><br></div>this is my code: <br>--------------------------------------------------------------------------------------------<br><br><br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br>;load "./WRFUserARW.ncl"<br><br><br> diri = "./"<br> jake = (/ "d03_JUN_P_11_YSU.nc", "d03_JUL_P_15_YSU.nc", "d03_AUG_P_15_YSU.nc", "d03_AUG_P_16_YSU.nc"/)<br> pthi = diri+jake<br> j = addfile("d03_DEC_MYJ_F.nc", "r")<br>fa = addfiles(pthi, "r")<br>pm25 = fa[:]->PM2_5_DRY<br>pm25_avg = dim_avg_n_Wrap(t2a, 0)<br>printVarSummary(pm25_avg) ; south_north | ??] x [west_east | ??]<br><br><br><br><br>; We generate plots, but what kind do we prefer?<br> type = "pdf"<br>; type = "pdf"<br>; type = "ps"<br>; type = "ncgm"<br><br>wks = gsn_open_wks(type,"avepm_d03")<br> <br> <br> <br> gsn_define_colormap(wks, "BlWhRe")<br> <br><br><br><br> pltres = True<br> mpres = True<br> res = True<br> mpres@mpDataBaseVersion = "HighRes"<br> res@MainTitle = "2 meter temperature Difference"<br> mpres@mpGridAndLimbOn = False<br> res@gsnMaximize = True<br> res@cnFillOn = True<br> res@gsnSpreadColors = True<br> res@cnLevelSelectionMode = "ManualLevels" <br> res@cnMinLevelValF = -1.8 ; set the minimum contour level<br> res@cnMaxLevelValF = 2. ; set the maximum contour level <br> res@cnLevelSpacingF = 0.1<br> res@cnLinesOn = False<br> res@cnLineLabelsOn = False<br> mpres@mpGeophysicalLineColor = "Black"<br> <br> <br> <br><br>plot = gsn_csm_contour_map(wks,pm25_avg,res)<br><br> contour_ave = wrf_contour(j,wks,pm25_avg,res) <br><br> plot = wrf_map_overlays(j,wks,(/contour_ave/),pltres,mpres)<br> <br> <br><br><br><br><br>--------------------------------------------------------------------------------------------------------<br><br>pm25 = fa[:]->PM2_5_DRY<br>pm25_avg = dim_avg_n_Wrap(t2a, 0)<br><br></div>maybe in these two lines i need to change something to just let it get thhe lowest level? I cannot find out how. Hoping for any help at all.<br><br></div>Cheers,<br></div>Jacob<br></div>
<br>
<a href="http://www.dlsu.edu.ph" target="_blank"><img src="http://www.dlsu.edu.ph/offices/mco/images/email-disclaimer/dlsu-logo.jpg"></a><a href="https://www.facebook.com/DLSU.Manila.100" target="_blank"><img src="http://www.dlsu.edu.ph/offices/mco/images/email-disclaimer/dlsu-fb.jpg"></a><a href="http://instagram.com/dlsu" target="_blank"><img src="http://www.dlsu.edu.ph/offices/mco/images/email-disclaimer/dlsu-instagram.jpg"></a><a href="https://twitter.com/dlsumanila" target="_blank"><img src="http://www.dlsu.edu.ph/offices/mco/images/email-disclaimer/dlsu-twitter.jpg"></a><a href="http://dlsumanila.tumblr.com/" target="_blank"><img src="http://www.dlsu.edu.ph/offices/mco/images/email-disclaimer/dlsu-tumbler.jpg"></a><a href="http://iblog.lasalle.ph/" target="_blank"><img src="http://www.dlsu.edu.ph/offices/mco/images/email-disclaimer/dlsu-blog.jpg"></a><font color="green" face="Arial" size="2"><br><br>DISCLAIMER AND CONFIDENTIALITY NOTICE </font><div><font color="green" face="Arial" size="2">The information contained in this e-mail, including those in its attachments, is confidential and intended only for the person(s) or entity(ies) to which it is addressed. If you are not an intended recipient, you must not read, copy, store, disclose, distribute this message, or act in reliance upon the information contained in it. If you received this e-mail in error, please contact the sender and delete the material from any computer or system. Any views expressed in this message are those of the individual sender and may not necessarily reflect the views of De La Salle University. </font></div>