Subversion Repositories lagranto.ecmwf

Rev

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

Rev 3 Rev 11
Line 401... Line 401...
401 
          klo=k
 401 
          klo=k
402 
        endif
 402 
        endif
403 
      goto 1
 403 
      goto 1
404 
      endif
 404 
      endif
405 
      h=xa(khi)-xa(klo)
 405 
      h=xa(khi)-xa(klo)
- 
 
 406 
      if (h.eq.0.) then
406 
      if (h.eq.0.) pause 'bad xa input in splint'
 407 
         print*, 'bad xa input in splint'
- 
 
 408 
         stop
- 
 
 409 
      endif
407 
      a=(xa(khi)-x)/h
 410 
      a=(xa(khi)-x)/h
408 
      b=(x-xa(klo))/h
 411 
      b=(x-xa(klo))/h
409 
      y=a*ya(klo)+b*ya(khi)+((a**3-a)*y2a(klo)+(b**3-b)*y2a(khi))*(h**
 412 
      y=a*ya(klo)+b*ya(khi)+((a**3-a)*y2a(klo)+(b**3-b)*y2a(khi))*(h**
410 
     *2)/6.
 413 
     *2)/6.
411 
      return
 414 
      return