Subversion Repositories lagranto.um

Rev

Rev 7 | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 7 Rev 16
Line 117... Line 117...
117 
      real              tmp1,tmp2
 117 
      real              tmp1,tmp2
118 
      real              tfrac,frac
 118 
      real              tfrac,frac
119 
      real              radius,dist
 119 
      real              radius,dist
120 
      character*80      string
 120 
      character*80      string
121 
      character*80      selectstr
121 
      character*80      selectstr
- 
 
 122 
      integer           noutside
122 
      real,allocatable, dimension (:,:,:) :: fld0
 123 
      real,allocatable, dimension (:,:,:) :: fld0
123 
      real,allocatable, dimension (:,:,:) :: fld1
 124 
      real,allocatable, dimension (:,:,:) :: fld1
124 
      real,allocatable, dimension (:,:  ) :: sfc0
 125 
      real,allocatable, dimension (:,:  ) :: sfc0
125 
      real,allocatable, dimension (:,:)   :: sfc1
 126 
      real,allocatable, dimension (:,:)   :: sfc1
126 
 
 127 
 
Line 378... Line 379...
378 
 
 379 
 
379 
c     ------------------------------------------------------------------
 380 
c     ------------------------------------------------------------------
380 
c     Read grid parameters from inital files
 381 
c     Read grid parameters from inital files
381 
c     ------------------------------------------------------------------
 382 
c     ------------------------------------------------------------------
382 
 
 383 
 
- 
 
 384 
c     Init missing data value
- 
 
 385 
      mdv = -999.0
- 
 
 386 
 
383 
c     Get the time of the first and second data file
 387 
c     Get the time of the first and second data file
384 
      tstart = -timeshift                              ! Format hh.mm
 388 
      tstart = -timeshift                              ! Format hh.mm
385 
      call hhmm2frac(tstart,frac)
 389 
      call hhmm2frac(tstart,frac)
386 
      frac   = frac + timeinc
390 
      frac   = frac + timeinc
387 
      call frac2hhmm(frac,tend)                        ! Format hh.mm
391 
      call frac2hhmm(frac,tend)                        ! Format hh.mm
Line 1099... Line 1103...
1099 
 
 1103 
 
1100 
      enddo
 1104 
      enddo
1101
1105
1102 
c     Remove all starting points outside the domain
 1106 
c     Remove all starting points outside the domain
1103 
      i        = 1
 1107 
      i        = 1
- 
 
 1108 
      noutside = 0
1104 
      do while ( i.le.start_n )
 1109 
      do while ( i.le.start_n )
1105
1110
1106 
         if ( abs(start_pre(i)-mdv).lt.eps ) then
 1111 
         if ( abs(start_pre(i)-mdv).lt.eps ) then
1107
1112
- 
 
 1113 
            noutside = noutside +1
- 
 
 1114 
            if ( noutside.lt.10 ) then
- 
 
 1115 
                print*,'  Outside ',
1108 
            print*,'  Outside ', start_lon(i),start_lat(i),start_lev(i)
 1116 
     >           start_lon(i),start_lat(i),start_lev(i)
- 
 
 1117 
            else if ( noutside.eq.10 ) then
- 
 
 1118 
                print*,'  Outside > 10 starting points'
- 
 
 1119 
            endif
1109 
 
 1120 
 
1110 
            do j=i,start_n
 1121 
            do j=i,start_n
1111 
               start_lon(j) = start_lon(j+1)
 1122 
               start_lon(j) = start_lon(j+1)
1112 
               start_lat(j) = start_lat(j+1)
 1123 
               start_lat(j) = start_lat(j+1)
1113 
               start_pre(j) = start_pre(j+1)
 1124 
               start_pre(j) = start_pre(j+1)