ncl批量读取站点资料

tSSkill

1.我现在有如下的数据，请问如何写一个批量读取的循环，年月的循环都好写，但是日有28，29，30，31天的不知道怎么解决
do i=0,dimsizes(year)-1
        do j=1,12
          do k=1,30
          nrow = numAsciiRow("/gpfshome/fenggl1/tsk/data/19510101-20180930/temp_daily20-20_"+year(i)+sprinti("%0.2i",j)+sprinti("%0.2i",k)+".txt")
          data=asciiread("/gpfshome/fenggl1/tsk/data/19510101-20180930/temp_daily20-20_"+year(i)+sprinti("%0.2i",j)+".txt",(/nrow,13/),"integer")

2.数据里有几千多站点，如何截取特定区域（比如江淮地区）的站点呢
新人，请大神赐教

尽头的尽头

1、（1）借助systemfunc命令，读取这个文件夹下的数据就好了；（2）如果非要分月读可以加一个if判断，不过二月还得判断闰年；（3）还有一种办法就是找一个包含你数据整个时间段的nc逐日数据，比如ncep1的数据，用cd_calender将你需要的时间段截取出来去读，大概思路就这样，根据你的实际情况吧
2、知道江淮地区的经纬度范围，然后判断站点的经纬度，在这个范围内的就读取或者重新写出来