<div dir="ltr">; t2a(time,e_we, s_sn, bot_top)<br><br><div> pm25_avg = dim_avg_n_Wrap(t2a, 0) <br> printVarSummary(pm25_avg) ; pm25_avg(e_we, s_sn, bot_top)<br><br></div><div> kl = 0 ; index of lowest level<br><br></div><div> pm25_kl = pm25_<wbr>avg(:,:,kl) ; for clarity, explicitly extract the lowest level<br></div><div> printVarSummary(pm25_kl) ; pm25_kl(e_we, s_sn)<br><br></div><div>; I am not sure why you are using the gsn_csm_contour_map function.<br></div><div>; You are not used it correctly. You must set more resources.<br></div><div>; plot = gsn_csm_contour_map(wks,pm25_<wbr>kl,res)<br><br> contour_ave = wrf_contour(j,wks,pm25_kl,<wbr>res) <br> plot = wrf_map_overlays(j,wks,(/<wbr>contour_ave/),pltres,mpres)<br><div><br></div>---<br></div><div>Of course, you could just use the subscripted variable directly.<br></div><div><br> ;plot = gsn_csm_contour_map(wks,pm25_<wbr>avg(:,:,kl),res)<br> contour_ave = wrf_contour(j,wks,pm25_avg(:,:,kl),<wbr>res) <br><br>--<br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 8, 2016 at 3:59 AM, Jacob Alberto Garcia <span dir="ltr"><<a href="mailto:jacob_garcia@dlsu.edu.ph" target="_blank">jacob_garcia@dlsu.edu.ph</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--<br><br><br><br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>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 <wbr> = "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_<wbr>avg,res)<br><br> contour_ave = wrf_contour(j,wks,pm25_avg,<wbr>res) <br><br> plot = wrf_map_overlays(j,wks,(/<wbr>contour_ave/),pltres,mpres)<br> <br> <br><br><br><br><br>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------<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 face="Arial" color="green" size="2"><br><br>DISCLAIMER AND CONFIDENTIALITY NOTICE </font><div><font face="Arial" color="green" 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><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>