<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=utf-8">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page Section1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EL link=blue vlink=purple style='word-wrap: break-word;-webkit-nbsp-mode: space;
-webkit-line-break: after-white-space'>
<div class=Section1>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Dear Jack,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>thanks again for your attention….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I have already made a script which the dust_tot for the zero
level (level=0) was estimated and wrote down the lat, lon and the dust_tot for
each grid point (see attached files at the previous email readwrite_dust_all.ncl,
with the output file All_Dust.txt (<a
href="https://www.dropbox.com/s/g8yscu2qvj3muuk/All_Dust.txt?dl=0">https://www.dropbox.com/s/g8yscu2qvj3muuk/All_Dust.txt?dl=0</a>
))….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I have to clarify that I want the average value of the produced
columns of each file that is calculated from your latest version and the script
readwrite_dust_all.ncl of mine….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>So in order to be clear, the revised script calculates the dust_tot
value for each grid point and for the first file (wrfoutput) which is written
into the 4<sup>th</sup> column.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>The next dust_tot value for the next file is written into the 5<sup>th</sup>
column, and etc<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I would like to calculate the average value of the 4<sup>th</sup>
and 5<sup>th</sup> columns (for 2 wrfoutput files) for each grid point….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>The problem is that I cannot use the dim_avg_n because the
dust_tot variable is 2d, one value for lat, lon of the dimensions….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>May be a new do loop or the modification of the existence loop
is the solution….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I hope to clarify the aim of the desired script…<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I apologise for the inconvenience and the possible waste of your
time….<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Best regards <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Ioannis <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<div>
<div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'>
<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span
style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>
ncl-talk-bounces@ucar.edu [mailto:ncl-talk-bounces@ucar.edu] <b>On Behalf Of </b>Jack
Ogaja<br>
<b>Sent:</b> Monday, November 23, 2015 8:07 PM<br>
<b>To:</b> Ioannis Koletsis<br>
<b>Cc:</b> ncl-talk [ncl-talk]<br>
<b>Subject:</b> Re: [ncl-talk] Calculate Average value for a number of files<o:p></o:p></span></p>
</div>
</div>
<p class=MsoNormal><o:p> </o:p></p>
<div>
<p class=MsoNormal>Ioannis,<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<p class=MsoNormal>Sorry, I oversaw one point in your script, the "density
weighted sum”. However, <o:p></o:p></p>
<div>
<blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'>
<div>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>I would like to calculate the average value for big number of
files wrfoutput which it seems to be infeasible with the revised version of the
script (readwrite_dust_all_2_Jack….</span><o:p></o:p></p>
</div>
</blockquote>
<div>
<p class=MsoNormal>for which dimension do you want to calculate the average?
time? level? It was not clear in your first email and unfortunately it is still
not clear in your recent email what you mean by “average”. If the rest of your
files have been written at each time step as it seems to be with the two files
you uploaded, then the only possible average will be in “levels” if you want
the output as indicated in the last email.<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>Sure, you wouldn’t use “readwrite_dust_all_2_Jack…” script
to read multiple files. In this case you need the “addfiles” function. Please,
find attached a modification of your latest script. Unlike the earlier script,
I have stayed as close as possible to the original coding and I have only added
“addfiles” and “dim_sum_n” functions to your script. You can read more about
the functions in the documentation.<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>Assuming you mean to calculate average in the vertical
direction, just switch between<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>;==<o:p></o:p></p>
</div>
<div>
<div>
<p class=MsoNormal> dust_tot = dust_J(level,:,:)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>;==<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>and<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>;==<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> ;+dust_tot = dim_avg_n(dust_J,0)<o:p></o:p></p>
</div>
</div>
<div>
<p class=MsoNormal>;==<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>in lines 54 and 56 respectively, otherwise for the two data
files you uploaded, the modified script should give you the exact numbers as
your original script(you should confirm it). If you add more files, the
modified script should print to the output file multiple columns representing
values from each file.<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>If you need further clarification/help, don’t hesitate to
write back to the ncl-talk. <o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>Jack<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</div>
</div>
<br /><br />
<hr style='border:none; color:#909090; background-color:#B0B0B0; height: 1px; width: 99%;' />
<table style='border-collapse:collapse;border:none;'>
<tr>
<td style='border:none;padding:0px 15px 0px 8px'>
<a href="https://www.avast.com/antivirus">
<img border=0 src="http://static.avast.com/emails/avast-mail-stamp.png" alt="Avast logo" />
</a>
</td>
<td>
<p style='color:#3d4d5a; font-family:"Calibri","Verdana","Arial","Helvetica"; font-size:12pt;'>
Αυτό το e-mail ελέγχθηκε για ιούς από το πρόγραμμα Avast antivirus.
<br><a href="https://www.avast.com/antivirus">www.avast.com</a>
</p>
</td>
</tr>
</table>
<br />
</body>
</html>