When extracting tracers vectors i_run and j_run are used rather than min_i:max_i
I can't remember the reason why this is so, but it seems to work fine for the benchmark data. However, using @vleguenn setup the reading of data is painfully slow (even though the number of bdy points << benchmark)
A quick test using np.min(i_run):np.max(i_run)+1 appears to remove the issue.
Need to check a little deeper before implementing this - I have a nagging suspicion this was done for a reason (maybe to get around the east-west wrap or north fold?)
When extracting tracers vectors
i_run
andj_run
are used rather thanmin_i:max_i
I can't remember the reason why this is so, but it seems to work fine for the benchmark data. However, using @vleguenn setup the reading of data is painfully slow (even though the number of bdy points << benchmark)A quick test using
np.min(i_run):np.max(i_run)+1
appears to remove the issue.Need to check a little deeper before implementing this - I have a nagging suspicion this was done for a reason (maybe to get around the east-west wrap or north fold?)