<div dir="ltr"><div><div><div><div><div>Did you run out of memory?<b><br><br>print(12419l*22*145*512*4)  ; *4 for bytes<br>                              <wbr>                 ; the appended  &#39;l&#39; [ small L ] make a long integer<br>(0)    81134817280 </b><br><br></div>81+GB  for the variable &#39;x&#39;<br></div><br></div>Another 81GB for the variable &#39;xAnom&#39;<br><br></div>Do you need the original &#39;x&#39; ? If not, maybe the following will work:<br><br>Replace<br><br><span style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px" class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span"><span style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px" class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span"><div style="word-wrap:break-word"> 
xAnom = new((/ntim,nlev,nlat,mlon/),fl<wbr>oat) ; This line causes the segmentation fault. </div>
<div style="word-wrap:break-word">
<div style="word-wrap:break-word">
 do n = 0,nlev-1</div>
<div style="word-wrap:break-word">
    xAnom(:,n,:,:)  = calcDayAnomTLL (x(:,n,:,:), yyyyddd, xClmDay(:,n,:,:))</div>
<div style="word-wrap:break-word">
  end d0</div></div></span></span><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px"><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px"><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px"><br><br></span></span></span></div><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px"><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px"><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px">With<br>   </span></span></span><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px"><span class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px">
<div style="word-wrap:break-word">
<div style="word-wrap:break-word">
 do n = 0,nlev-1</div>
<div style="word-wrap:break-word">
    x(:,n,:,:)  = calcDayAnomTLL (x(:,n,:,:), yyyyddd, xClmDay(:,n,:,:))</div>
<div style="word-wrap:break-word">
  end do</div></div></span></span><br><span style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px" class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span"><span style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px" class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span"><span style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-indent:0px;border-spacing:0px" class="m_3972741238158577484gmail-m_5216042728470290617Apple-style-span"> <br></span></span></span></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 25, 2017 at 1:43 PM, Cristiana Stan <span dir="ltr">&lt;<a href="mailto:cstan@gmu.edu" target="_blank">cstan@gmu.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div style="word-wrap:break-word">
<br>
<div>
<div style="text-align:start;text-indent:0px;word-wrap:break-word">
<span class="m_5471710558388831061Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;border-spacing:0px">
<div style="word-wrap:break-word">
<span class="m_5471710558388831061Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;border-spacing:0px">
<div style="word-wrap:break-word">
<span class="m_5471710558388831061Apple-style-span" style="border-collapse:separate;font-variant-ligatures:normal;font-variant-numeric:normal;font-variant-alternates:normal;font-variant-east-asian:normal;line-height:normal;text-align:-webkit-auto;text-indent:0px;border-spacing:0px">
<div style="word-wrap:break-word">
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px">
<div style="word-wrap:break-word">
Hi,</div>
</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
I would like to compute daily anomalies of a large 4D variable that exceeds the NCL size limit and causes segmentation fault. I’m looking for a suggestion to work around this issue. </div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
My code looks like:</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<br>
</div>
<div style="word-wrap:break-word">
 x       =  f-&gt;$var$(iStrt:iLast,:,:,:)</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
ntim = 12419</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
nlev = 22</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
nlat = 145</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
mlon = 512</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<span style="text-align:-webkit-auto"><br>
</span></div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<span style="text-align:-webkit-auto"><br>
</span></div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<span style="text-align:-webkit-auto">xClmDay = clmDayTLLL((/x/), yyyyddd)    </span></div>
<div style="word-wrap:break-word">
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px">
<br>
</div>
</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<br>
</div>
<div style="word-wrap:break-word">
xAnom = new((/ntim,nlev,nlat,mlon/),<wbr>float) ; This line causes the segmentation fault. </div>
<div style="word-wrap:break-word">
<div style="word-wrap:break-word">
 do n = 0,nlev-1</div>
<div style="word-wrap:break-word">
    xAnom(:,n,:,:)  = calcDayAnomTLL (x(:,n,:,:), yyyyddd, xClmDay(:,n,:,:))</div>
<div style="word-wrap:break-word">
  end do</div>
<div>
<div style="word-wrap:break-word">
<br>
</div>
</div>
<div style="word-wrap:break-word">
<div style="word-wrap:break-word">
;*****************************<wbr>******************************</div>
<div style="word-wrap:break-word">
; Create netCDF using &#39;simple&#39; method</div>
<div style="word-wrap:break-word">
;*****************************<wbr>******************************</div>
</div>
<div style="word-wrap:break-word">
<br>
</div>
</div>
</div>
</span></div>
</span></div>
</span><br class="m_5471710558388831061Apple-interchange-newline">
</div>
<div style="text-align:start;text-indent:0px;word-wrap:break-word">
Thanks,</div>
<div style="text-align:start;text-indent:0px;word-wrap:break-word">
<br>
</div>
<div style="text-align:start;text-indent:0px;word-wrap:break-word">
Cristiana</div>
<br class="m_5471710558388831061Apple-interchange-newline">
<br class="m_5471710558388831061Apple-interchange-newline">
</div>
<br>
</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>