C C This programme is to read Chinese soil moisture data. C Created by Mingquan Mu in April 20,2001. C character*12 name(43) real sst(43,11,3,12,1981:1991) real lat(43),lon(43) integer id(43) do iy=1981,1991 do im=1,12 do ida=1,3 do li=1,11 do i=1,43 sst(i,li,ida,im,iy)=-9.9 enddo enddo enddo enddo enddo C++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ open(11,file='soilm_china.dat',status='old') do i=1,9 read(11,*) enddo do 101 i=1,43 read(11,21) name(i),id(i),ibeg,iend, & lat(i),lon(i) write(*,21) name(i),id(i),ibeg,iend, & lat(i),lon(i) 101 continue do 201 i=1,43 print*,'i is',i do ii=1,5 read(11,*) enddo do 201 iy=ibeg,iend write(*,*) iy read(11,*) read(11,*) do 201 im=1,12 do 201 ida=1,3 read(11,31) (sst(i,li,ida,im,iy),li=1,11) write(*,31) (sst(i,li,ida,im,iy),li=1,11) 201 continue close(11) 21 format(4X,A12,4X,I5,3X,I4,1X,I4,2X,F6.2,3X,F6.2) 31 format(6X,11(1X,f6.2)) stop end