Subversion Repositories lagranto.wrf

Rev

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

Rev 2 Rev 11
Line 720... Line 720...
720 
     >		           spt0,spt1,p3d0,p3d1,uut0,uut1,vvt0,vvt1,wwt0,wwt1,
 720 
     >		           spt0,spt1,p3d0,p3d1,uut0,uut1,vvt0,vvt1,wwt0,wwt1,
721 
     >                 xmin,ymin,dx,dy,per,hem,nx,ny,nz,mpx,mpy)
 721 
     >                 xmin,ymin,dx,dy,per,hem,nx,ny,nz,mpx,mpy)
722 
 
 722 
 
723 
      implicit none
 723 
      implicit none
724 
 
 724 
 
- 
 
 725 
c     Flag for test mode
- 
 
 726 
      integer      test
- 
 
 727 
      parameter    (test=0)
- 
 
 728 
 
725 
c     Declaration of subroutine parameters
 729 
c     Declaration of subroutine parameters
726 
      integer      nx,ny,nz
 730 
      integer      nx,ny,nz
727 
      real         x1,y1,p1
 731 
      real         x1,y1,p1
728 
      integer      left
 732 
      integer      left
729 
      real	       x0,y0,p0
 733 
      real	       x0,y0,p0
Line 807... Line 811...
807 
C        Calculate new positions (adapted for cartesian grid)
 811 
C        Calculate new positions (adapted for cartesian grid)
808 
         x1 = x0 + fbflag * deltat * u * dx/mpsc_x
 812 
         x1 = x0 + fbflag * deltat * u * dx/mpsc_x
809 
         y1 = y0 + fbflag * deltat * v * dy/mpsc_y
 813 
         y1 = y0 + fbflag * deltat * v * dy/mpsc_y
810 
         p1 = p0 + fbflag * deltat * w * wfactor
 814 
         p1 = p0 + fbflag * deltat * w * wfactor
811 
 
 815 
 
- 
 
 816 
c        Write the update positions for tests
- 
 
 817 
         if ( (icount.eq.3).and.(test.eq.1) ) then
- 
 
 818 
            write(*,'(10f10.2)') x0,u,x1-x0,p0,w,p1-p0
- 
 
 819 
         endif
- 
 
 820 
 
812 
C       Interpolate surface geop.height to actual position
 821 
C       Interpolate surface geop.height to actual position
813 
        call get_index4 (xind,yind,pind,x1,y1,0.,reltpos1,
 822 
        call get_index4 (xind,yind,pind,x1,y1,0.,reltpos1,
814 
     >                   p3d0,p3d1,spt0,spt1,3,
 823 
     >                   p3d0,p3d1,spt0,spt1,3,
815 
     >                   nx,ny,nz,xmin,ymin,dx,dy,mdv)
 824 
     >                   nx,ny,nz,xmin,ymin,dx,dy,mdv)
816 
        sp = int_index4 (spt0,spt1,nx,ny,1,xind,yind,1.,reltpos1,mdv)
 825 
        sp = int_index4 (spt0,spt1,nx,ny,1,xind,yind,1.,reltpos1,mdv)