Subversion Repositories lagranto.ecmwf

Rev

Rev 13 | 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