[ncl-talk] efficient way of reading numbered directories?

[Tabish Ansari]

Hi

I am currently reading in files in a very inefficient way from 60
directories:

DATADir1 = "/mnt/rasher/60sensruns/run1/"
DATADir2 = "/mnt/rasher/60sensruns/run2/"
DATADir3 = "/mnt/rasher/60sensruns/run3/"
DATADir4 = "/mnt/rasher/60sensruns/run4/"
DATADir5 = "/mnt/rasher/60sensruns/run5/"
DATADir6 = "/mnt/rasher/60sensruns/run6/"
DATADir7 = "/mnt/rasher/60sensruns/run7/"
DATADir8 = "/mnt/rasher/60sensruns/run8/"
DATADir9 = "/mnt/rasher/60sensruns/run9/"
DATADir10 = "/mnt/rasher/60sensruns/run10/"
DATADir11 = "/mnt/rasher/60sensruns/run11/"
DATADir12 = "/mnt/rasher/60sensruns/run12/"
DATADir13 = "/mnt/rasher/60sensruns/run13/"
DATADir14 = "/mnt/rasher/60sensruns/run14/"
DATADir15 = "/mnt/rasher/60sensruns/run15/"
DATADir16 = "/mnt/rasher/60sensruns/run16/"
DATADir17 = "/mnt/rasher/60sensruns/run17/"
DATADir18 = "/mnt/rasher/60sensruns/run18/"
DATADir19 = "/mnt/rasher/60sensruns/run19/"
DATADir20 = "/mnt/rasher/60sensruns/run20/"
DATADir21 = "/mnt/rasher/60sensruns/run21/"
DATADir22 = "/mnt/rasher/60sensruns/run22/"
DATADir23 = "/mnt/rasher/60sensruns/run23/"
DATADir24 = "/mnt/rasher/60sensruns/run24/"
DATADir25 = "/mnt/rasher/60sensruns/run25/"
DATADir26 = "/mnt/rasher/60sensruns/run26/"
DATADir27 = "/mnt/rasher/60sensruns/run27/"
DATADir28 = "/mnt/rasher/60sensruns/run28/"
DATADir29 = "/mnt/rasher/60sensruns/run29/"
DATADir30 = "/mnt/rasher/60sensruns/run30/"
DATADir31 = "/mnt/rasher/60sensruns/run31/"
DATADir32 = "/mnt/rasher/60sensruns/run32/"
DATADir33 = "/mnt/rasher/60sensruns/run33/"
DATADir34 = "/mnt/rasher/60sensruns/run34/"
DATADir35 = "/mnt/rasher/60sensruns/run35/"
DATADir36 = "/mnt/rasher/60sensruns/run36/"
DATADir37 = "/mnt/rasher/60sensruns/run37/"
DATADir38 = "/mnt/rasher/60sensruns/run38/"
DATADir39 = "/mnt/rasher/60sensruns/run39/"
DATADir40 = "/mnt/rasher/60sensruns/run40/"
DATADir41 = "/mnt/rasher/60sensruns/run41/"
DATADir42 = "/mnt/rasher/60sensruns/run42/"
DATADir43 = "/mnt/rasher/60sensruns/run43/"
DATADir44 = "/mnt/rasher/60sensruns/run44/"
DATADir45 = "/mnt/rasher/60sensruns/run45/"
DATADir46 = "/mnt/rasher/60sensruns/run46/"
DATADir47 = "/mnt/rasher/60sensruns/run47/"
DATADir48 = "/mnt/rasher/60sensruns/run48/"
DATADir49 = "/mnt/rasher/60sensruns/run49/"
DATADir50 = "/mnt/rasher/60sensruns/run50/"
DATADir51 = "/mnt/rasher/60sensruns/run51/"
DATADir52 = "/mnt/rasher/60sensruns/run52/"
DATADir53 = "/mnt/rasher/60sensruns/run53/"
DATADir54 = "/mnt/rasher/60sensruns/run54/"
DATADir55 = "/mnt/rasher/60sensruns/run55/"
DATADir56 = "/mnt/rasher/60sensruns/run56/"
DATADir57 = "/mnt/rasher/60sensruns/run57/"
DATADir58 = "/mnt/rasher/60sensruns/run58/"
DATADir59 = "/mnt/rasher/60sensruns/run59/"
DATADir60 = "/mnt/rasher/60sensruns/run60/"

For each directory I will further have to do this:
FILES1 = systemfunc (" ls -1 " + DATADir1 + "subset_d03_2014-* ")
a1 = addfiles(FILES1,"r")
pm25_sim_a1    = a1[:]->PM2_5_DRY(:,0,145,79)

(I will have to type this another 59 times for the other directories)

Is there a way to do this more efficiently perhaps in a loop?

Thanks

Tabish

Tabish U Ansari
PhD student, Lancaster Environment Center
Lancaster Univeristy
Bailrigg, Lancaster,
LA1 4YW, United Kingdom
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190221/19090a5b/attachment.html>