list the content of a directory specified by loop indices in fortran
list the content of a directory specified by loop indices in fortran
(OP)
Hello everybody
0
down vote
favorite
I need some help, well let's get started:
I'm working with numerical weather forecast named as follow:
sub_gfsanl_4_2011MMDD-IIII-FFFF.grb2
-MM stands for month from 01 to 12 -DD stands for days from 01 to 31 -IIII stands for initialisation time, the first and second digits are for hours the third and last are for minutes -FFFF stands for forecast hour , the first and second digits are for hours the third and last are for minutes
in my directory i have several files for a given days of a given months. A day has 4 data, one at every six hours IIII=0000 ,0600,1200,1800.
What i'm trying to do is to list all files of a given days, here the f90 code i wrote:
I would greatly appreciate any help of advises, for some reasons the execute_command_line doesn't seem to like the $ for variable...
Thanks!
0
down vote
favorite
I need some help, well let's get started:
I'm working with numerical weather forecast named as follow:
sub_gfsanl_4_2011MMDD-IIII-FFFF.grb2
-MM stands for month from 01 to 12 -DD stands for days from 01 to 31 -IIII stands for initialisation time, the first and second digits are for hours the third and last are for minutes -FFFF stands for forecast hour , the first and second digits are for hours the third and last are for minutes
in my directory i have several files for a given days of a given months. A day has 4 data, one at every six hours IIII=0000 ,0600,1200,1800.
What i'm trying to do is to list all files of a given days, here the f90 code i wrote:
program test_ec
implicit none
!==variable declaration==
integer :: mi,di,dil,mil
character*3 :: temp
!==Program instructions==
mil=1
write(temp,'(i2.2)') mil
read(temp,'(i2.2)') mi
!convert the month into a two digit value mi=01
! change to directory where the data are stored
CALL chdir('/media/Hello/ncfiles/GFS' )
do dil=1,31
!loop over days
write(temp,'(i2.2)') dil
read(temp,'(i2.2)') di
! converting day number into a two digit number, store this value into di. ex dil=9 then di=09
CALL execute_command_line( 'ls sub_gfsanl_4_2011${mi}${di}*.nc > yes.txt' )
!list all files with the correct month and days and store it to yes.txt
end do
end program test_ec
I would greatly appreciate any help of advises, for some reasons the execute_command_line doesn't seem to like the $ for variable...
Thanks!





RE: list the content of a directory specified by loop indices in fortran
program test_ec
implicit none
integer :: mi,di,iiii,fff
character(len=1024) :: filename
character(len=1024) :: format_string
logical exist
mi=1
CALL chdir('/media/Hello/ncfiles/GFS' )
!do di=1,10
do iiii=0,18,6
do fff=0,18,6
format_string = "(A17,i2.2,i2.2,A1,i4.4,A1,i3.3,A8)"
write (filename,format_string) "sub_gfsanl_4_2011",mi,di,"_",iiii,"_",fff,".grb2.nc"
inquire(file=trim(filename), exist=exist)
if (exist) then
write(*,*) 'file exists i can do do whatever i want with this file'
else
write(*,*) 'I did not find that file.'
end if
enddo
enddo
enddo
end program test_ec
RE: list the content of a directory specified by loop indices in fortran
You don't really need a format string - the code can be written as
CODE
CODE