Subversion Repositories lagranto.um

Rev

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

Rev 7 Rev 16
Line 116... Line 116...
116 
      integer           stat,flag
 116 
      integer           stat,flag
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)