Subversion Repositories lagranto.wrf

Rev

Rev 2 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 2 Rev 11
Line 319... Line 319...
319 
c     Read ascii format (mode=1,2)
 319 
c     Read ascii format (mode=1,2)
320 
      if ( (mode.eq.1).or.(mode.eq.2) ) then
 320 
      if ( (mode.eq.1).or.(mode.eq.2) ) then
321 
 
 321 
 
322 
c        Read the time specification (old and new format)
 322 
c        Read the time specification (old and new format)
323 
         read(fid,'(a80)') linestr
 323 
         read(fid,'(a80)') linestr
324
324 
 
325 
         if ( linestr(1:14).eq.'Reference date' ) then
 325 
         if ( linestr(1:14).eq.'Reference date' ) then
326 
            read(linestr,'(a15,i4,i2,i2,a1,i2,i2,a13,i8,a4)')
326 
            read(linestr,'(a15,i4,i2,i2,a1,i2,i2,a13,i8,a4)')
327 
     >           str1,
 327 
     >           str1,
328 
     >           time(1),time(2),time(3),str2,time(4),time(5),
 328 
     >           time(1),time(2),time(3),str2,time(4),time(5),
329 
     >           str3,time(6),str4
 329 
     >           str3,time(6),str4
Line 609... Line 609...
609 
            elseif ( (isstr.eq.1).and.(ch.eq.' ') ) then
 609 
            elseif ( (isstr.eq.1).and.(ch.eq.' ') ) then
610 
               isstr=0
 610 
               isstr=0
611 
               ncol=ncol+1
 611 
               ncol=ncol+1
612 
            endif
 612 
            endif
613 
         enddo
 613 
         enddo
614 
 
 - 
 
- 
 
 614 
c           backspace(fid)
615 
c        Get the first data block
 615 
c        Get the first data block
616 
         nline0  = 5
 616 
         nline0  = 5
617 
         nline1  = 5
 617 
         nline1  = 5
618 
         read(fid,*)
 618 
         read(fid,*)
619 
 100     read(fid,'(a500)',end=110) str
 619 
 100     read(fid,'(a500)',end=110) str
620 
         if (str.ne.'') then
 620 
         if (str.ne.'') then
621 
            nline1 = nline1 + 1
 621 
            nline1 = nline1 + 1
622 
            goto 100
 622 
            goto 100
623 
         endif
 623 
         endif
624 
 110     continue
 624 
 110     continue
- 
 
 625 
           backspace(fid)
625
626
626 
c        Get the total numbers of lines in the data block
 627 
c        Get the total numbers of lines in the data block
627 
         nline2 = nline1
 628 
         nline2 = nline1
628 
 120     read(fid,*,end=130)
 629 
 120     read(fid,*,end=130)
629 
         nline2 = nline2 + 1
 630 
         nline2 = nline2 + 1