Subversion Repositories lagranto.ecmwf

Rev

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

Rev 13 Rev 15
Line 439... Line 439...
439
439
440 
         ppo0 = int_index3(vert,nx,ny,nz,rid,rjd,real(1),0.)
 440 
         ppo0 = int_index3(vert,nx,ny,nz,rid,rjd,real(1),0.)
441 
         rkd=0
 441 
         rkd=0
442 
         do i=1,nz-1
 442 
         do i=1,nz-1
443 
            ppo1 = int_index3(vert,nx,ny,nz,rid,rjd,real(i+1),0.)
 443 
            ppo1 = int_index3(vert,nx,ny,nz,rid,rjd,real(i+1),0.)
444 
            if ( (ppo0.lt.ppo).and.(ppo1.ge.ppo) ) then
 444 
            if ( (ppo.gt.ppo0).and.(ppo.le.ppo1) ) then
445 
               rkd=real(i)+(ppo0-ppo)/(ppo0-ppo1)
 445 
               rkd=real(i)+(ppo-ppo0)/(ppo1-ppo0)
446 
               goto 100
 446 
               goto 100
447 
            endif
 447 
            endif
448 
            ppo0 = ppo1
 448 
            ppo0 = ppo1
449 
         enddo
 449 
         enddo
450 
 
 450