Subversion Repositories lagranto.ecmwf

Rev

Rev 29 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 29 Rev 46
Line 19... Line 19...
19 
      integer   nlevmax
 19 
      integer   nlevmax
20 
      parameter (nlevmax=100)
 20 
      parameter (nlevmax=100)
21 
 
 21 
 
22 
c     Maximum number of input files (dates, length of trajectories)
 22 
c     Maximum number of input files (dates, length of trajectories)
23 
      integer   ndatmax
 23 
      integer   ndatmax
24 
      parameter (ndatmax=5000)
 24 
      parameter (ndatmax=10000)
25 
 
 25 
 
26 
c     Numerical epsilon (for float comparison)
 26 
c     Numerical epsilon (for float comparison)
27 
      real      eps
 27 
      real      eps
28 
      parameter (eps=0.001)
 28 
      parameter (eps=0.001)
29 
 
 29 
 
Line 105... Line 105...
105 
      real                                   stagz           ! Vertical staggering
 105 
      real                                   stagz           ! Vertical staggering
106 
      real                                   mdv             ! Missing data value
 106 
      real                                   mdv             ! Missing data value
107 
 
 107 
 
108 
c	  Output grid  and fields
 108 
c	  Output grid  and fields
109 
      real                                   levels(1000)    ! Ouput levels
 109 
      real                                   levels(1000)    ! Ouput levels
110 
      real                                   times (5000)    ! Output times
 110 
      real                                   times (10000)    ! Output times
111 
	  real,allocatable, dimension (:,:)   :: out_pos         ! Position of trajectories
 111 
	  real,allocatable, dimension (:,:)   :: out_pos         ! Position of trajectories
112 
	  real,allocatable, dimension (:,:)   :: out_val         ! Output lidar field
 112 
	  real,allocatable, dimension (:,:)   :: out_val         ! Output lidar field
113 
	  real,allocatable, dimension (:,:)   :: out_cnt         ! # output lidar sum ups
 113 
	  real,allocatable, dimension (:,:)   :: out_cnt         ! # output lidar sum ups
114 
 
 114 
 
115 
 
 115 
 
Line 485... Line 485...
485 
                p3t1     = p3t0
 485 
                p3t1     = p3t0
486 
                spt1     = spt0
 486 
                spt1     = spt0
487 
                iloaded1 = itime1
 487 
                iloaded1 = itime1
488 
             endif
 488 
             endif
489 
 
 489 
 
490 
c            Load manager:  Load first time (tracing variable and grid)
 490 
c            Load manager:  Load first time (tracing variable and grid)
491 
             if ( itime0.ne.iloaded0 ) then
 491 
             if ( itime0.ne.iloaded0 ) then
492 
 
 492 
 
493 
                filename = tfil(i)//dat(itime0)
 493 
                filename = tfil(i)//dat(itime0)
494 
                call frac2hhmm(time0,tload)
 494 
                call frac2hhmm(time0,tload)
495 
                varname  = tvar(i)
495 
                varname  = tvar(i)
Line 498... Line 498...
498 
     >               trim(filename),'  ',trim(varname),tload
 498 
     >               trim(filename),'  ',trim(varname),tload
499 
                call input_open (fid,filename)
499 
                call input_open (fid,filename)
500 
                call input_wind
500 
                call input_wind
501 
     >               (fid,varname,f3t0,tload,stagz,mdv,
 501 
     >               (fid,varname,f3t0,tload,stagz,mdv,
502 
     >               xmin,xmax,ymin,ymax,dx,dy,nx,ny,nz,timecheck)
502 
     >               xmin,xmax,ymin,ymax,dx,dy,nx,ny,nz,timecheck)
503
- 
 
504 
                call input_grid
503 
                call input_grid
505 
     >               (fid,varname,xmin,xmax,ymin,ymax,dx,dy,nx,ny,
 504 
     >               (fid,varname,xmin,xmax,ymin,ymax,dx,dy,nx,ny,
506 
     >                tload,pollon,pollat,p3t0,spt0,nz,ak,bk,stagz,
 505 
     >                tload,pollon,pollat,p3t0,spt0,nz,ak,bk,stagz,
507 
     >                timecheck)
 506 
     >                timecheck)
508 
                call input_close(fid)
 507 
                call input_close(fid)