Subversion Repositories lagranto.ecmwf

Rev

Rev 27 | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 27 Rev 46
Line 557... Line 557...
557 
             p1 = pp0(i)
 557 
             p1 = pp0(i)
558 
             call get_index4 (xind,yind,pind,x1,y1,p1,0.,
 558 
             call get_index4 (xind,yind,pind,x1,y1,p1,0.,
559 
     >                 p3t1,p3t1,spt1,spt1,3,
 559 
     >                 p3t1,p3t1,spt1,spt1,3,
560 
     >                 nx,ny,nz,xmin,ymin,dx,dy,mdv)
 560 
     >                 nx,ny,nz,xmin,ymin,dx,dy,mdv)
561 
             zindex(i) = pind
 561 
             zindex(i) = pind
- 
 
 562 
             if ( i.lt.20 ) then
- 
 
 563 
                  print*,pp0(i),' hPa -> IND ',zindex(i)
- 
 
 564 
             endif
562 
          enddo
 565 
          enddo
563 
 
 566 
 
564 
          do i=1,nz
 - 
 
565 
            print*,i,p3t1(189+(141-1)*nx+(i-1)*nx*ny)
 - 
 
566 
          enddo
 - 
 
567 
             print*,x1,y1,p1
 - 
 
568 
             print*,xind,yind,pind
 - 
 
569 
 
 - 
 
570 
c         Write info about zindex range of starting positions
 567 
c         Write info about zindex range of starting positions
571 
          zindexmin = zindex(1)
 568 
          zindexmin = zindex(1)
572 
          zindexmax = zindex(1)
 569 
          zindexmax = zindex(1)
573 
          do i=2,ntra
 570 
          do i=2,ntra
574 
            if ( zindex(i).gt.zindexmax ) zindexmax = zindex(i)
 571 
            if ( zindex(i).gt.zindexmax ) zindexmax = zindex(i)
Line 888... Line 885...
888 
           endif
 885 
           endif
889 
 
 886 
 
890 
          enddo
 887 
          enddo
891 
 
 888 
 
892 
C         Save positions only every deltout minutes
 889 
C         Save positions only every deltout minutes
- 
 
 890 
          delta =
893 
          delta = aint(iloop*60*ts/deltout)-iloop*60*ts/deltout
 891 
     > aint(real(iloop)*60.*ts/real(deltout))-
- 
 
 892 
     > real(iloop)*60.*ts/real(deltout)
894 
          if (abs(delta).lt.eps) then
 893 
       if (abs(delta).lt.eps) then
895 
c          wstep = wstep + abs(ts)
 894 
c          wstep = wstep + abs(ts)
896 
c          if ( mod(wstep,deltout).eq.0 ) then
 895 
c          if ( mod(wstep,deltout).eq.0 ) then
897 
            time = time0+reltpos1*timeinc*fbflag
 896 
            time = time0+reltpos1*timeinc*fbflag
898 
            itim = itim + 1
 897 
            itim = itim + 1