subroutine read_nu(Nb,nu_low,nu_hi) implicit none include 'max.inc' include 'formats.inc' integer strlen,ios,dim,Nb,i,band character*(Nchar_mx) file_nu character*(Nchar_mx) strtmp1,strtmp2 double precision nu_low(1:Nbmx) double precision nu_hi(1:Nbmx) file_nu='./venus_nu.dat' open(10,file=file_nu(1:strlen(file_nu)), & status='old',iostat=ios) if (ios.ne.0) then ! file not found write(*,*) 'Error from routine read_nu:' write(*,*) 'File could not ne found:' write(*,*) file_nu(1:strlen(file_nu)) stop endif do i=1,4 read(10,*) enddo read(10,33) strtmp1,dim,strtmp2,Nb do i=1,3 read(10,*) enddo do band=1,Nb read(10,*) nu_low(band),nu_hi(band) enddo close(10) return end