Thanks for your reply, Travis.<div><br></div><div>I asked that because I´m searching possible errors in my runs. My simulations using MIROC dataset  as initial and boundary conditions always finish with CFL violation at specific days. I put timestep with 15 seconds, but CFL violation continues. I´m using dx with 50 km. I´m looking for errors in MIROC dataset, but no success until now. </div>
<div><br></div><div>An important aspect considering MIROC dataset is that it has many missing values in the bottom of atmosphere (1000, 975, 950 mb) of South America (because of Andes mountains). This occurs because MIROC is a spectral model. Is there any special treatment for avoiding missing values in input dataset? I´ve made a basic interpolation based on 4 neighborhood points, but I´m really don´t know if it has some influence on CFL violation.</div>
<div><br></div><div>I made the following in mod_miroc.F90 (I created it) after reading netcdf file.</div><div><br></div><div><div>! INTERPOLATION TO ELIMINATE MISSING VALUES (-999.)</div><div><div>! klev = 23 | form 1 to 1000 mb</div>
<div>!  ilon = 320 | ilon = 1 (0.00 E) | ilon = 320 (358.875 E)</div><div>!  jlat = 160 | ilat = 1 (89.1415 N) | ilat = 160 (89.1415 S)</div><div><br></div></div><div>          do k = 1 , klev</div><div>            do j = 2 , jlat-1</div>
<div>              do i = 2 , ilon-1</div><div>                if (work(i,j,k)==-999.0000)then</div><div>                   contador=0</div><div>                   if (work(i-1,j,k)==-999.0000) then</div><div>                      temp1=0.</div>
<div>                      contador=0</div><div>                   else</div><div>                      temp1=work(i-1,j,k)</div><div>                      contador=1</div><div>                  endif</div><div>                   if (work(i+1,j,k)==-999.0000) then</div>
<div>                      temp2=0.</div><div>                   else</div><div>                      temp2=work(i+1,j,k)</div><div>                      contador=contador+1</div><div>                   endif</div><div>                   if (work(i,j-1,k)==-999.0000) then</div>
<div>                      temp3=0.</div><div>                   else</div><div>                      temp3=work(i,j-1,k)</div><div>                      contador=contador+1</div><div>                   endif</div><div>                   if (work(i,j+1,k)==-999.0000) then</div>
<div>                      temp4=0.</div><div>                   else</div><div>                      temp4=work(i,j+1,k)</div><div>                      contador=contador+1</div><div>                   endif</div><div>                work(i,j,k)=(temp1+temp2+temp3+temp4)/contador</div>
<div>                endif</div><div>              end do</div><div>            end do</div><div>          end do</div><div><br></div><div>Anyone can give an opinion about this method? Is there anything better to solve missing value points?</div>
<div><br></div><div>Thanks to everybody!</div><div>Felipe.</div><br><div class="gmail_quote">2012/8/15 Travis O&#39;Brien <span dir="ltr">&lt;<a href="mailto:travis.obrien@gmail.com" target="_blank">travis.obrien@gmail.com</a>&gt;</span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Felipe,<br><br>This should not be anything to worry about.  This is a typical occurrence in CLM output.  To be certain, I have asked a CLM-developer colleague of mine (Dr. Jinyun Tang), and he indicated that this message often appears because of rounding errors.<br>


<br>Kind Regards,<br>-Travis-<br><br><div class="gmail_quote"><div><div class="h5">On Mon, Aug 13, 2012 at 8:00 AM, Felipe Roque <span dir="ltr">&lt;<a href="mailto:felipenrs@gmail.com" target="_blank">felipenrs@gmail.com</a>&gt;</span> wrote:<br>


</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">Hi Friends,<div><br></div><div>I´m using RegCM-4.1. Anyone knows what could cause an error in subroutine BalanceCheck (~/Main/clm3.5/src/biogeophys/BalanceCheckMod.F90)? This warning doesn´t stop regcm, but is it a critical error? Should I worry about it?</div>



<div><div>This subroutine accumulates the numerical truncation errors of the water and energy balance calculation. It is helpful to see the performance of the process of integration.</div><div>The error for energy balance = abs(Net radiation - change of internal energy - Sensible heat - Latent heat)</div>



</div><div><br></div><div><br></div><div><div> soil balance e nstep =     49500 point =  4105 imbalance =    0.00 W/m2</div><div> soil balance e nstep =     49500 point =  5456 imbalance =    0.00 W/m2</div><div> soil balance e nstep =     49500 point =  2304 imbalance =    0.00 W/m2</div>



<div>     at day =   68.7497, ktau =     198000 :  1st, 2nd time deriv of ps =  0.10248E-04 0.39686E-07,  no. of points w/convection =    3309</div><div>            *** Day      69.7500 solar declination angle =  -4.23964630 degrees.</div>



<div>  dectim =    397440.000000000</div><div> ATM variables written at   1989031018  0.000000000000000E+000</div><div> SRF variables written at   1989031018  0.000000000000000E+000</div><div>SEARCH BC data for 1989031006</div>



<div>READY  BC from     1989031000 to 1989031006</div><div> soil balance e nstep =     49501 point =   614 imbalance =    0.00 W/m2</div><div> soil balance e nstep =     49501 point =  4743 imbalance =    0.00 W/m2</div>


<div>
 soil balance e nstep =     49501 point =  2304 imbalance =    0.00 W/m2</div><div> soil balance e nstep =     49501 point =  5456 imbalance =    0.00 W/m2</div><div> soil balance e nstep =     49501 point =  4599 imbalance =    0.00 W/m2</div>



<div> soil balance e nstep =     49501 point =  2716 imbalance =    0.00 W/m2</div></div><div>.......</div><div>.......</div><div>.......</div><div><br></div><div>Thanks, Felipe.</div>
<br></div></div>_______________________________________________<br>
RegCNET mailing list<br>
<a href="mailto:RegCNET@lists.ictp.it" target="_blank">RegCNET@lists.ictp.it</a><br>
<a href="https://lists.ictp.it/mailman/listinfo.cgi/regcnet" target="_blank">https://lists.ictp.it/mailman/listinfo.cgi/regcnet</a><br></blockquote></div><br>
</blockquote></div><br></div>