<div dir="ltr"><div><div><div><div><div><div><div><div>Please see the addfiles documentation<br><br></div><div> diri = "./"<br></div> fili = (/ "FOO_1", "FOO_2.nc", "FOO_3.nc" /)<br><br></div> pthi = diri+fili<br></div> f = addfiles(pthi, "r")<br><br></div> T2 = f[:]->T2<br></div> printVarSummary(T2)<br><br></div> T2_avg = dim_avg_n_Wrap(T2, 0)<br></div> printVarSummary(T2_avg)<br><br>;======<br> <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml</a><br><br></div>Good luck<br><div><div><div><div><div><div><div> <br><div>fa = addfile("d03_FEB_YSU_F.nc","r"<wbr>)<br>t2a = fa->T2<br>t2a_avg = dim_avg_n_Wrap(t2a, 0)<br>printVarSummary(t2a_avg) ; south_north | ??] x [west_east | ??]<br><br>fb = addfile("d03_Feb_MYJ_f.nc","r"<wbr>)<br>t2b = fb->T2<br>t2b_avg = dim_avg_n_Wrap(t2b, 0)<br>printVarSummary(t2b_avg) ; south_north |??] x [west_east | ??]<br><br></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Aug 6, 2016 at 10:13 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>Hello,<br><br></div>I want to take the three files and find the time average of the 2 meter temperature. I have this script here that calculates the mean difference of T2 for the two files given below. What I want to do is instead of just one month, I'd like to Average three months for the YSU and three months for the MYJ simulations. What do I change in this script?<br></div><br><br><div>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>fa = addfile("d03_FEB_YSU_F.nc","r"<wbr>)<br>t2a = fa->T2<br>t2a_avg = dim_avg_n_Wrap(t2a, 0)<br>printVarSummary(t2a_avg) ; south_north | ??] x [west_east | ??]<br><br>fb = addfile("d03_Feb_MYJ_f.nc","r"<wbr>)<br>t2b = fb->T2<br>t2b_avg = dim_avg_n_Wrap(t2b, 0)<br>printVarSummary(t2b_avg) ; south_north |??] x [west_east | ??]<br><br>t2_diff = t2a_avg-t2b_avg<br>copy_VarMeta(t2a_avg, t2_diff)<br>t2_diff@long_name = "T2 Difference"<br>t2_diff@units = t2a@units<br>printVarSummary(t2_diff) ; south_north |??] x [west_east | ??]<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,"plot_DIFF_<wbr>d03_FEB")<br> gsn_define_colormap(wks,"<wbr>BlueDarkRed18")<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@cnLinesOn = False<br> res@cnLineLabelsOn = False<br> mpres@mpGeophysicalLineColor = "Black"<br> <br> <br><br><br> plot = gsn_csm_contour_map(wks,t2_<wbr>diff,res)<br><br> contour_DIFF = wrf_contour(fa,wks,t2_diff,<wbr>res) <br><br> plot = wrf_map_overlays(fa,wks,(/<wbr>contour_DIFF/),pltres,mpres)<br> <br>==============================<wbr>=============================<br><br></div><div>Would something like this work?<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>fa = addfile("d03_DEC_YSU_F.nc","r"<wbr>)<br>faa=addfile("d03_JAN_YSU_F.nc"<wbr>,r)<br>faaa=addfile('d03_FEB_YSU_F.<wbr>nc",r)<br>t2a = fa->T2<br>t2aa = faa->T2<br>t2aaa = faaa-> T2<br>t2a_avg = dim_avg_n_Wrap(t2a, t2aa, t2aaa, 0)<br>printVarSummary(t2a_avg) ; south_north | ??] x [west_east | ??]<br><br>fb = addfile("d03_Dec_MYJ_f.nc","r"<wbr>)<br>fbb = addfile("d03_Jan_MYJ_f.nc","r"<wbr>)<br>fbbb = addfile("d03_Feb_MYJ_f.nc","r"<wbr>)<br>t2b = fb->T2<br>t2bb = fbb->T2<br>t2bbb = fbbb ->T2<br>t2b_avg = dim_avg_n_Wrap(t2b, t2bb, t2bbb, 0)<br>printVarSummary(t2b_avg) ; south_north |??] x [west_east | ??]<br><br>t2_diff = t2a_avg-t2b_avg<br>copy_VarMeta(t2a_avg, t2_diff)<br>t2_diff@long_name = "T2 Difference"<br>t2_diff@units = t2a@units<br>printVarSummary(t2_diff) ; south_north |??] x [west_east | ??]<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,"plot_DIFF_<wbr>d03_FEB")<br> gsn_define_colormap(wks,"<wbr>BlueDarkRed18")<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@cnLinesOn = False<br> res@cnLineLabelsOn = False<br> mpres@mpGeophysicalLineColor = "Black"<br> <br> <br><br><br> plot = gsn_csm_contour_map(wks,t2_<wbr>diff,res)<br><br> contour_DIFF = wrf_contour(fa,wks,t2_diff,<wbr>res) <br><br> plot = wrf_map_overlays(fa,wks,(/<wbr>contour_DIFF/),pltres,mpres)<br> <br> <br><br><br><br><br><br></div><div> <br><br><br><br><br><br></div></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>