<div dir="ltr"><div><div><div>Hi<br><br></div>I am trying to read SEVIRI data with a lot of variables and the dimension of each variable is 3712*3712. I know the data are pretty large. But the computer should read them smoothly. Since the memory that I specified is about 50GB. Unfortunately, the code is becoming slower and slower while it do the time loop. Here is a sample of my code.<br><br>;  SET UP THE START TIME AND END TIME<br>   start_year = 2008<br>   end_year   = 2008<br>   start_month= 6<br>   end_month  = 6<br>   start_day  = 1<br>   start_hour = 0<br>   end_hour   = 23<br>   start_min  = 0<br>   end_min    = 45<br>   min_stride = 15<br>   start_ind_lat = 1400<br>   end_ind_lat   = 3000<br>   start_ind_lon = 1100<br>   end_ind_lon   = 2600<br><br>;  DO YEAR LOOP<br>   do iyear = start_year, end_year<br><br>;  DO MONTH LOOP<br>      do imonth = start_month, end_month<br><br>;  CALCULATE THE NUMBER OF DAYS IN THIS MONTH<br>         nday_month = days_in_month(iyear, imonth)<br>;  DO DAY LOOP<br>         do iday = start_day, 10;nday_month<br>;  DO HOUR LOOP<br>            do ihour = start_hour, end_hour<br>;  DO MINUTE LOOP<br>               do imin = start_min, end_min, min_stride<br>;  READ VARIABLES FROM HDF FILE<br>                     a     = addfile(dir + siyear + &quot;/&quot; + symd1 + &quot;/&quot; + filename, &quot;r&quot;)<br>                     lat   = (/a-&gt;MSG_Latitude(start_ind_lat:end_ind_lat, start_ind_lon:end_ind_lon)/)<br>                     lon   = (/a-&gt;MSG_Longitude(start_ind_lat:end_ind_lat, start_ind_lon:end_ind_lon)/)<br>                     Cloud_Optical_Thickness_16 = a-&gt;Cloud_Optical_Thickness_16(start_ind_lat:end_ind_lat, start_ind_lon:end_ind_lon)<br><br>               end do ;imin<br>            end do ;ihour<br>         end do ;iday<br>      end do ;imonth<br>   end do ;iyear<br><br><br></div>Thank you<br></div>Zhifeng<br></div>