Subversion Repositories lagranto.wrf

Rev

Rev 17 | Rev 21 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 17 Rev 19
Line 657... Line 657...
657 
      enddo
 657 
      enddo
658 
 
 658 
 
659 
c     Transform trajectory position from WRF grid do grid index
 659 
c     Transform trajectory position from WRF grid do grid index
660 
      do i=1,ntra
 660 
      do i=1,ntra
661 
        do j=1,ntim
 661 
        do j=1,ntim
- 
 
 662 
           if ( abs(traout(i,j,2)-mdv).gt.eps ) then
662 
           traout(i,j,2) = ( traout(i,j,2) - xmin ) / dx + 1.
 663 
              traout(i,j,2) = ( traout(i,j,2) - xmin ) / dx + 1.
- 
 
 664 
           else
- 
 
 665 
              traout(i,j,2) = mdv
- 
 
 666 
           endif
- 
 
 667 
           if ( abs(traout(i,j,3)-mdv).gt.eps ) then
663 
           traout(i,j,3) = ( traout(i,j,3) - ymin ) / dy + 1.
 668 
              traout(i,j,3) = ( traout(i,j,3) - ymin ) / dy + 1.
- 
 
 669 
           else
- 
 
 670 
              traout(i,j,3) = mdv
- 
 
 671 
           endif
664 
        enddo
 672 
        enddo
665 
      enddo
 673 
      enddo
666 
 
 674 
 
667 
c     Write trajectory file
 675 
c     Write trajectory file
668 
      vars(1)  ='time'
 676 
      vars(1)  ='time'