<div dir="ltr">Hi Dave,<div><br></div><div style>I was unaware of this ":=" issue. I simply thought "fout" was a file handle. After your recommended changes I still got the error. I even passed fout to the function and returned 0 to a dummy variable but I still got the same error. In order to speed up the script I moved the fout->time(itime) = (/f->time(itime)/)/60.0 into the loop. The speeds things up tremendously when the time dimension is set to unlimited. But there is something going on that I cannot grasp about NCL.</div>
<div style><br></div><div style><div>print("Regridding the variables...") </div><div>do itime = 0, 3;endtime </div>
<div> print("Processing time step: "+itime) </div><div> fout->time(itime) = (/f->time(itime)/)/60.0 ; change the time to hours </div>
<div> do n = 0, nsize ; loop over the different variables </div><div> print("Processing variable:-> "+name(n)+" mapping to:->"+Sname(n)) </div>
<div> if(name(n).eq."var26" .or.\ </div><div> name(n).eq."var142" .or.\ </div>
<div> name(n).eq."var143" .or.\ </div><div> name(n).eq."var178" .or.\ </div>
<div> name(n).eq."var179" .or.\ </div><div> name(n).eq."var184" .or.\ </div>
<div> name(n).eq."var245" .or.\ </div><div> name(n).eq."var246" .or.\ </div>
<div> name(n).eq."var247" .or.\ </div><div> name(n).eq."var248" .or.\ </div>
<div> name(n).eq."var249" .or.\ </div><div> name(n).eq."var251" .or.\ </div>
<div> name(n).eq."var252" .or.\ </div><div> name(n).eq."var253") then </div>
<div> fout := Pro2D(f,fout,Sname(n),name(n),itime,LAT1x1,LON1x1,LAT2x2,LON2x2,n) </div><div> end if </div>
<div> </div><div> if(name(n).eq."var21" .or.\ </div>
<div> name(n).eq."var22" .or.\ </div><div> name(n).eq."var23" .or.\ </div>
<div> name(n).eq."var24" .or.\ </div><div> name(n).eq."var25" .or.\ </div>
<div> name(n).eq."var223") then </div><div> do level = 0, nlevs </div>
<div> fout := Pro3D(f,fout,Sname(n),name(n),itime,level,f->lev_2(level),LAT1x1,LON1x1,LAT2x2,LON2x2,n) </div><div> end do </div>
<div> end if </div><div> end do </div>
<div> ;print("Time in file: "+fout->time(itime)) </div><div style>end do</div><div style><br></div><div>@a7 ncl_codes]$ ncl echam_regrid.ncl</div><div> Copyright (C) 1995-2013 - All Rights Reserved</div><div>
University Corporation for Atmospheric Research</div><div> NCAR Command Language Version 6.1.2</div><div> The use of this software is governed by a License Agreement.</div><div> See <a href="http://www.ncl.ucar.edu/">http://www.ncl.ucar.edu/</a> for more details.</div>
<div>(0) Names in nc file: var26</div><div>(1) Names in nc file: var142</div><div>(2) Names in nc file: var143</div><div>(3) Names in nc file: var178</div><div>(4) Names in nc file: var179</div><div>(5) Names in nc file: var184</div>
<div>(6) Names in nc file: var245</div><div>(7) Names in nc file: var246</div><div>(8) Names in nc file: var247</div><div>(9) Names in nc file: var248</div><div>(10) Names in nc file: var249</div><div>(11) Names in nc file: var251</div>
<div>(12) Names in nc file: var252</div><div>(13) Names in nc file: var253</div><div>(14) Names in nc file: var21</div><div>(15) Names in nc file: var22</div><div>(16) Names in nc file: var23</div><div>(17) Names in nc file: var24</div>
<div>(18) Names in nc file: var25</div><div>(19) Names in nc file: var223</div><div>(20) Names in nc file: var243</div><div>(21) Names in nc file: var244</div><div>(0) Create a large file to hold the data</div>
<div>(0) Regridding the variables...</div><div>(0) Processing time step: 0</div><div>(0) Processing variable:-> var26 mapping to:->SP</div><div>(0) Processing variable:-> var142 mapping to:->LSP</div>
<div>(0) Processing variable:-> var143 mapping to:->CP</div><div>(0) Processing variable:-> var178 mapping to:->TSIR</div><div>(0) Processing variable:-> var179 mapping to:->TTR</div><div>(0) Processing variable:-> var184 mapping to:->SI</div>
<div>(0) Processing variable:-> var245 mapping to:->sunlit_ISCCP</div><div>(0) Processing variable:-> var246 mapping to:->TCC_ISCCP</div><div>(0) Processing variable:-> var247 mapping to:->CTP_ISCCP</div>
<div>(0) Processing variable:-> var248 mapping to:->CALB_ISCCP</div><div>(0) Processing variable:-> var249 mapping to:->TAU_ISCCP</div><div>(0) Processing variable:-> var251 mapping to:->DCFLAG</div>
<div>(0) Processing variable:-> var252 mapping to:->SCFLAG</div><div>(0) Processing variable:-> var253 mapping to:->MCFLAG</div><div>(0) Calculate the ALBEDO...</div><div>(0) Processing variable:-> var21 mapping to:->R</div>
<div>(0) Processing variable:-> var22 mapping to:->Q</div><div>(0) Processing variable:-> var23 mapping to:->CLWC</div><div>(0) Processing variable:-> var24 mapping to:->CIWC</div><div>(0) Processing variable:-> var25 mapping to:->T</div>
<div>(0) Processing variable:-> var223 mapping to:->CC</div><div>(0) Processing variable:-> var243 mapping to:->CC_ISCCP_3D</div><div>(0) Processing variable:-> var244 mapping to:->TAU_ISCCP_3D</div>
<div>(0) Processing time step: 1</div><div>(0) Processing variable:-> var26 mapping to:->SP</div><div>(0) Processing variable:-> var142 mapping to:->LSP</div><div>(0) Processing variable:-> var143 mapping to:->CP</div>
<div>(0) Processing variable:-> var178 mapping to:->TSIR</div><div>(0) Processing variable:-> var179 mapping to:->TTR</div><div>(0) Processing variable:-> var184 mapping to:->SI</div><div>(0) Processing variable:-> var245 mapping to:->sunlit_ISCCP</div>
<div>(0) Processing variable:-> var246 mapping to:->TCC_ISCCP</div><div>(0) Processing variable:-> var247 mapping to:->CTP_ISCCP</div><div>(0) Processing variable:-> var248 mapping to:->CALB_ISCCP</div>
<div>(0) Processing variable:-> var249 mapping to:->TAU_ISCCP</div><div>(0) Processing variable:-> var251 mapping to:->DCFLAG</div><div>(0) Processing variable:-> var252 mapping to:->SCFLAG</div>
<div>(0) Processing variable:-> var253 mapping to:->MCFLAG</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div>
<div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>(0) Calculate the ALBEDO...</div>
<div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div>
<div>fatal:Subscript out of range, error in subscript #0</div><div style> </div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Oct 22, 2013 at 8:50 PM, David Brown <span dir="ltr"><<a href="mailto:dbrown@ucar.edu" target="_blank">dbrown@ucar.edu</a>></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">Hi Marston,<div>You are using 'fout' as the variable for both the 2D and the 3D variables. Once the variable is initialized, any attempt to assign different dimensions or sizes will cause an error.</div>
<div>The fix, most likely, is either to call</div><div>delete(fout) </div><div>at the end of each loop iteration, or (recommended) to use the new 'reassignment' operator ':=', e.g.:</div><div><br></div><div>
fout : = Pro2D(f,fout,Sname(n),name(n),itime,LAT1x1,LON1x1,LAT2x2,LON2x2,n)</div><div><br></div><div>and later</div><div><br></div><div>fout := Pro3D(f,fout,Sname(n),name(n),itime,level,f->lev_2(level),LAT1x1,LON1x1,LAT2x2,LON2x2,n)</div>
<div><br></div><div> -dave</div><div><br></div><div><br><div><div><div class="h5"><div>On Oct 22, 2013, at 12:09 PM, Marston Johnston <<a href="mailto:shejo284@gmail.com" target="_blank">shejo284@gmail.com</a>> wrote:</div>
<br></div></div><blockquote type="cite"><div><div class="h5"><div dir="ltr">Hi all,<div><br></div><div>I know this error: fatal:Subscript out of range, error in subscript #0, seems like a simple enough error message to understand but I cannot seem to find the cause of the problem. I'm writing data to a file and after some variables are written, I read the data from the newly written file and create a new variable, then I write the new variable to the outfile. It is when I write this new variable that I get this error which frankly boggles my mind. The error occurs on the 2nd time step. I would love some assistance in understanding and solving this error.</div>
<div><br></div><div>I have attached the script and below is a copy of the output for two time steps:</div><div><br></div><div><div>%@a7 ncl_codes]$ ncl echam_regrid.ncl</div><div> Copyright (C) 1995-2013 - All Rights Reserved</div>
<div> University Corporation for Atmospheric Research</div><div> NCAR Command Language Version 6.1.2</div><div> The use of this software is governed by a License Agreement.</div><div> See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.</div>
<div>(0) Names in nc file: var26</div><div>(1) Names in nc file: var142</div><div>(2) Names in nc file: var143</div><div>(3) Names in nc file: var178</div><div>(4) Names in nc file: var179</div><div>(5) Names in nc file: var184</div>
<div>(6) Names in nc file: var245</div><div>(7) Names in nc file: var246</div><div>(8) Names in nc file: var247</div><div>(9) Names in nc file: var248</div><div>(10) Names in nc file: var249</div><div>
(11) Names in nc file: var251</div>
<div>(12) Names in nc file: var252</div><div>(13) Names in nc file: var253</div><div>(14) Names in nc file: var21</div><div>(15) Names in nc file: var22</div><div>(16) Names in nc file: var23</div><div>(17) Names in nc file: var24</div>
<div>(18) Names in nc file: var25</div><div>(19) Names in nc file: var223</div><div>(20) Names in nc file: var243</div><div>(21) Names in nc file: var244</div><div>(0) Create a large file to hold the data</div>
<div>(0) Regridding the variables...</div><div>(0) Processing time step: 0</div><div>(0) Processing variable:-> var26 mapping to:->SP</div><div>(0) Processing variable:-> var142 mapping to:->LSP</div>
<div>(0) Processing variable:-> var143 mapping to:->CP</div><div>(0) Processing variable:-> var178 mapping to:->TSIR</div><div>(0) Processing variable:-> var179 mapping to:->TTR</div><div>(0) Processing variable:-> var184 mapping to:->SI</div>
<div>(0) Processing variable:-> var245 mapping to:->sunlit_ISCCP</div><div>(0) Processing variable:-> var246 mapping to:->TCC_ISCCP</div><div>(0) Processing variable:-> var247 mapping to:->CTP_ISCCP</div>
<div>(0) Processing variable:-> var248 mapping to:->CALB_ISCCP</div><div>(0) Processing variable:-> var249 mapping to:->TAU_ISCCP</div></div><div><div>(0) Processing variable:-> var251 mapping to:->DCFLAG</div>
<div>(0) Processing variable:-> var252 mapping to:->SCFLAG</div><div>(0) Processing variable:-> var253 mapping to:->MCFLAG</div><div>(0) Calculate the ALBEDO...</div><div>(0) Processing variable:-> var21 mapping to:->R</div>
<div>(0) Processing variable:-> var22 mapping to:->Q</div><div>(0) Processing variable:-> var23 mapping to:->CLWC</div><div>(0) Processing variable:-> var24 mapping to:->CIWC</div><div>(0) Processing variable:-> var25 mapping to:->T</div>
<div><br></div><div>(0) Processing variable:-> var223 mapping to:->CC</div><div>(0) Processing variable:-> var243 mapping to:->CC_ISCCP_3D</div><div>(0) Processing variable:-> var244 mapping to:->TAU_ISCCP_3D</div>
<div>(0) Processing time step: 1</div><div>(0) Processing variable:-> var26 mapping to:->SP</div><div>(0) Processing variable:-> var142 mapping to:->LSP</div><div>(0) Processing variable:-> var143 mapping to:->CP</div>
<div>(0) Processing variable:-> var178 mapping to:->TSIR</div><div>(0) Processing variable:-> var179 mapping to:->TTR</div><div>(0) Processing variable:-> var184 mapping to:->SI</div><div>(0) Processing variable:-> var245 mapping to:->sunlit_ISCCP</div>
<div>(0) Processing variable:-> var246 mapping to:->TCC_ISCCP</div><div>(0) Processing variable:-> var247 mapping to:->CTP_ISCCP</div><div>(0) Processing variable:-> var248 mapping to:->CALB_ISCCP</div>
<div>(0) Processing variable:-> var249 mapping to:->TAU_ISCCP</div><div>(0) Processing variable:-> var251 mapping to:->DCFLAG</div><div>(0) Processing variable:-> var252 mapping to:->SCFLAG</div>
<div>(0) Processing variable:-> var253 mapping to:->MCFLAG</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div>
<div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:Subscript out of range, error in subscript #0</div><div><br></div><div>The input file is very large but if you have room on the NCL server I could upload it.</div>
<div><br></div><div>Thanks in advance,</div></div><div>/Marston</div></div>
</div></div><span><echam_regrid.ncl></span>_______________________________________________<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>
</blockquote></div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Only the fruitful thing is true!
</div>