<div dir="ltr"><div>Hi Steven,<br><br></div>Yes a quick test indicates that the missing_values is the culprit (for more, check out:<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 7, 2015 at 8:13 AM, BLIUJUS, STEVEN D CTR USAF AFMC AFLCMC/HBAW-OL <span dir="ltr"><<a href="mailto:steven.bliujus.3.ctr@us.af.mil" target="_blank">steven.bliujus.3.ctr@us.af.mil</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I am trying to take two files and add the variables together. One variable<br>
has data over water while the other has missing values, when I try adding<br>
the values, nothing prints out for the areas over the ocean. I tried setting<br>
the variables to 0 instead of missing values but that did not work either. I<br>
have attached three images: var1, var2, and var. The variable var is the<br>
addition, and the addition gets done over land but not over water. Not sure<br>
why this is occurring, my only guess is the struggle with missing data.<br>
However, when I look at var, it doesn't look as though it is adding the<br>
variables properly either since the max for var is 2.49807, but the max for<br>
var2 is 2.544, so there is an obvious discrepancy. Here is my info:<br>
<br>
1. I am using Version 6.1.2<br>
2. My Script:<br>
<br>
load "/home/bliujuss/ncl/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
load "/home/bliujuss/ncl/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br>
load "/home/bliujuss/ncl/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
begin<br>
<br>
a=addfile("/home/bliujuss/Plots/<a href="http://MYD08_D3.051.timeAverage.1.nc" target="_blank">MYD08_D3.051.timeAverage.1.nc</a>","r")<br>
b=addfile("/home/bliujuss/Plots/<a href="http://MOD08_D3.051.timeAverage.0.nc" target="_blank">MOD08_D3.051.timeAverage.0.nc</a>","r")<br>
wks = gsn_open_wks("x11","plt_aod")<br>
gsn_define_colormap(wks,"3gauss")<br>
<br>
var1 = a->Deep_Blue_Aerosol_Optical_Depth_550_Land_QA_Mean(:,:)<br>
var2 = b->Optical_Depth_Land_And_Ocean_Mean(:,:)<br>
<br>
printVarSummary(var1)<br>
printMinMax(var1, True)<br>
printVarSummary(var2)<br>
printMinMax(var2, True)<br>
<br>
<br>
var = var1<br>
var = 0.0<br>
var = var2+var1<br>
<br>
printVarSummary(var)<br>
printMinMax(var, True)<br>
<br>
pltres = True<br>
pltres@tiMainString = "AOD550"<br>
pltres@cnFillOn = True<br>
pltres@gsnSpreadColors = False<br>
pltres@gsnAddCyclic = False<br>
<br>
pltres@cnLineLabelsOn = False<br>
pltres@cnLineLabelBackgroundColor = "white"<br>
;plot = gsn_csm_contour_map_ce(wks,var({-46.2:46.5},{-4.5:184.5}),pltres)<br>
;plot = gsn_csm_contour_map_ce(wks,var1({-46.2:46.5},{-4.5:184.5}),pltres)<br>
;plot = gsn_csm_contour_map_ce(wks,var2({-46.2:46.5},{-4.5:184.5}),pltres)<br>
End<br>
<br>
3. Files are attached<br>
<br>
4. Images attached<br>
<br>
5. No error message<br>
<br>
6. gcc (GCC) 4.1.2 20080704 (Red Hat 4.1.2-44)<br>
<br>
<br>
<br>
Steven Bliujus, Contractor<br>
SEMS/16WS WXE<br>
557th Weather Wing<br>
DNS: 232-7151<br>
<br>
<br>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>