WebSVN – lagranto.ecmwf – Diff – /trunk/select/select.f

             var = 'INCIRCLE'
          elseif ( ivar.eq. -6 ) then
             var = 'INREGION'
          elseif ( ivar.eq. -7 ) then
             var = 'TRIGGER'
+         elseif ( ivar.eq. -8 ) then
+            var = 'VERT0'
          else
             var=vars(ivar)
          endif
          write(*,'(i5,f15.4,10x,a10)') i,out(i),trim(var)
          if (imode.eq.4) mode='MIN'
          if (imode.eq.5) mode='VAR'
          if (imode.eq.6) mode='SUM'
          if (imode.eq.7) mode='CHANGE'
          if (imode.eq.8) mode='DIFF'
+         if (imode.eq.9) mode='RANGE'
          write(*,'(i5,f15.4,10x,a10)') i,out(i),trim(mode)
 c        Get values
          i=i+1
          nval=nint(out(i))
             else
                write(*,'(i5,f15.4,a2)') i,out(i),char(nint(out(i)))
             endif
          enddo
-c        Get times
+c        Get the number of times
          i=i+1
          ntim=nint(out(i))
+c        the number of times is variable - depending on TRIGGER
+         if ( ntim .eq. -993 ) then
+            write(*,'(i5,f15.4,7x,7x,a15)') i,out(i),'TIMES @ TRIGGER'
+c        Get the detailed list of times
+         else
-         write(*,'(i5,f15.4,7x,7x,a6)') i,out(i),'#TIMES'
+            write(*,'(i5,f15.4,7x,7x,a6)') i,out(i),'#TIMES'
-         do j=1,ntim
+            do j=1,ntim
-            i=i+1
+               i=i+1
-            write(*,'(i5,f15.4,f7.0)') i,out(i),times(nint(out(i)))
+               write(*,'(i5,f15.4,f7.0)') i,out(i),times(nint(out(i)))
+            enddo
-         enddo
+         endif
 c        Get time mode
          i=i+1
          itime=nint(out(i))
          if (itime.eq.1) strtim='ALL'
 c        OUT      4     MIN      4       TRIGGER -i (i the trigger index)
 c        EQ       5     VAR      5
 c        BEGIN   -1     SUM      6
 c        END     -2     CHANGE   7
 c        AND     -3     DIFF     8
-c        OR      -4
+c        OR      -4     RANGE    9
 c        TRUE     6
 c        FALSE    7
 c        SPECIAL  0
 c        ALL      8 (TRIGGER)
 c        ANY      9 (TRIGGER)
 c        INPOLYGON -3 : Specified polygon region
 c        INBOX     -4 : Longitude/latitude rectangle
 c        INCIRCLE  -5 : Within a specified radius
 c        INREGION  -6 : Within a specified rehion on the region file
 c        TRIGGER   -7 : Trigger field
+c        VERT0     -8 : Vertical distance from starting position
+c
 c     For the special commands BEGIN, END, AND and OR, only one field
 c     in <out> is used.
       implicit none
       integer      i,j
       integer      j1,j2
       integer      flag(ntimes)
       integer      count
       integer      ok
-      integer      itrigger
+      integer      itrigger,ttrigger
 c     Common block for initialisation of polygon check
       real    tlonv(1000),vlat_c(1000),vlon_c(1000)
       real    xlat_c,xlon_c
       integer ibndry,nv_c
         else
            read(fid,*) nval
            read(fid,*) (val(i),i=1,nval)
         endif
-c       Read times
+c       Read times (on request, change to special trigger times)
         read(fid,*) ntim
+        if ( ntim.eq.-993 ) then
+           ttrigger = 1
+           ntim     = 1
+           tim(1)   = -999.
+        else
+           ttrigger = 0
-        read(fid,*) (tim(i),i=1,ntim)
+           read(fid,*) (tim(i),i=1,ntim)
+        endif
         read(fid,*) tmode
 c       Bring CAPITAL "TIME,LAT,LON,P" into "time,lat,lon,p"
         if (var.eq.'TIME') var='time'
         if (var.eq.'LAT' ) var='lat'
               ok = -5
            elseif (var.eq.'INREGION') then
               ok = -6
            elseif (var.eq.'TRIGGER') then
               ok = -7
+           elseif (var.eq.'VERT0') then
+              ok = -8
            else
              goto 400
           endif
         endif
         n=n+1
         if (mode.eq.'MIN'    ) ok=4
         if (mode.eq.'VAR'    ) ok=5
         if (mode.eq.'SUM'    ) ok=6
         if (mode.eq.'CHANGE' ) ok=7
         if (mode.eq.'DIFF'   ) ok=8
+        if (mode.eq.'RANGE'  ) ok=9
         if (ok.eq.0) goto 400
         n=n+1
         out(n)=real(ok)
 c       Write the parameter values: INPOLYGON
               n=n+1
               out(n)=val(i)
            enddo
         endif
+c       Special time handling: only trigger times are cosidered
+        if ( ttrigger.eq.1 ) then
+           n = n+1
+           out(n)=-993.
 c       All times are selected
-        if ( abs(tim(1)+999.).lt.eps ) then
+        elseif ( abs(tim(1)+999.).lt.eps ) then
            n=n+1
            out(n)=real(ntimes)
            do i=1,ntimes
               n=n+1
               out(n)=real(i)
       integer   flag
       real      dist
       real      xcorner(4),ycorner(4)
       integer   iparam
       character ch
+      real      varmin,varmax
+      real      lev0,lev1
 c     Common block for initialisation of polygon check
       real    tlonv(1000),vlat_c(1000),vlon_c(1000)
       real    xlat_c,xlon_c
       integer ibndry,nv_c
          enddo
 c        Get times (<time(j)> gets the row indices of <tra>)
          i=i+1
          nsel=nint(cmd(i))
+         if ( nsel .eq. -993 ) then
+            nsel = 0
+            do k=1,ntim
+              if ( trigger(k).ne.0 ) then
+                nsel       = nsel + 1
+                time(nsel) = k
+              endif
+           enddo
+         else
-         do j=1,nsel
+           do j=1,nsel
-            i=i+1
+              i=i+1
-            time(j)=nint(cmd(i))
+              time(j)=nint(cmd(i))
-         enddo
+           enddo
+         endif
+c        If no times are selected, exit with non-select status
+         if ( nsel.eq.0 ) then
+             stack(1) = 0
+             goto 300
+         endif
 c        Get time mode
          i=i+1
          itime=nint(cmd(i))
 c        Implict variable: TRIGGER
          elseif ( ivar.eq. -7 ) then
             do j=1,nsel
                intvar(j) = trigger( time(j) )
             enddo
+c        Implicit variable: VERT0
+         elseif ( ivar.eq. -8 ) then
+            do j=1,nsel
+               lev0   = tra(1      ,4)
+               lev1   = tra(time(j),4)
+               var(j) = lev0 - lev1
+            enddo
 c        Explicit variable (column index <ivar>)
          else
             do j=1,nsel
                var(j) = tra(time(j),ivar)
                tmp=tmp+var(j)
             enddo
             var(1)=tmp
             nsel=1
-c        Take CHANGE of the variable (change between first and last time)
+c        Take CHANGE of the variable (absolute difference between first and last time)
          elseif (imode.eq.7) then
             var(1)=abs(var(1)-var(nsel))
             nsel=1
 c        Take DIFF of the variable (first minus last time)
          elseif (imode.eq.8) then
             var(1)=var(1)-var(nsel)
             nsel=1
+c        Take RANGE of the variable
+         elseif (imode.eq.9) then
+            varmax=var(1)
+            varmin=var(1)
+            do j=2,nsel
+               if (var(j).gt.varmax) varmax=var(j)
+               if (var(j).lt.varmin) varmin=var(j)
+            enddo
+            var(1) = varmax - varmin
+            nsel=1
          endif
 c        --- Apply the operators to the single values ---
          do j=1,nsel

Subversion Repositories lagranto.ecmwf

(root)/trunk/select/select.f @ 42 – Rev 3 → 9