Subversion Repositories lagranto.ecmwf

Rev

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

Rev 5 Rev 9
Line 248... Line 248...
248 
      endif
 248 
      endif
249 
 
 249 
 
250 
c     Read the constant grid parameters (nx,ny,nz,xmin,xmax,ymin,ymax,pollon,pollat)
 250 
c     Read the constant grid parameters (nx,ny,nz,xmin,xmax,ymin,ymax,pollon,pollat)
251 
c     The negative <-fid> of the file identifier is used as a flag for parameter retrieval
251 
c     The negative <-fid> of the file identifier is used as a flag for parameter retrieval
252 
      filename = charp//dat(1)
 252 
      filename = charp//dat(1)
253 
      varname  = 'U'
 253 
      varname  = tvar(1)
254 
      call input_open (fid,filename)
 254 
      call input_open (fid,filename)
255 
      call input_grid (-fid,varname,xmin,xmax,ymin,ymax,dx,dy,nx,ny,
 255 
      call input_grid (-fid,varname,xmin,xmax,ymin,ymax,dx,dy,nx,ny,
256 
     >                 tstart,pollon,pollat,rd,rd,nz,rd,rd,rd,timecheck)
 256 
     >                 tstart,pollon,pollat,rd,rd,nz,rd,rd,rd,timecheck)
257 
      call input_close(fid)
 257 
      call input_close(fid)
258 
 
 258 
 
Line 638... Line 638...
638 
                         y0=89.99
 638 
                         y0=89.99
639 
                      endif
639 
                      endif
640 
 
 640 
 
641 
c                 If requested, dump the lidar coordinates
 641 
c                 If requested, dump the lidar coordinates
642 
                  if ( (dumpcoord.eq.'yes').and.(i.eq.1) ) then
 642 
                  if ( (dumpcoord.eq.'yes').and.(i.eq.1) ) then
643 
                     write(10,'3f10.2') x0,y0,trainp(k,j,1)
 643 
                     write(10,'(3f10.2)') x0,y0,trainp(k,j,1)
644 
                     write(10,'3f10.2') x0_tra,y0_tra,5.
 644 
                     write(10,'(3f10.2)') x0_tra,y0_tra,5.
645 
                  endif
 645 
                  endif
646 
 
 646 
 
647 
C                 Get the index where to interpolate (x0,y0,p0)
 647 
C                 Get the index where to interpolate (x0,y0,p0)
648 
                  if ( (abs(x0-mdv).gt.eps).and.
 648 
                  if ( (abs(x0-mdv).gt.eps).and.
649 
     >                 (abs(y0-mdv).gt.eps) )
 649 
     >                 (abs(y0-mdv).gt.eps) )