C*********************************************************************
C
C     PROGRAM TQDVLOP
C
C     水平平均温度・混合比鉛直分布作成プログラム(時間発展)
C
C     1998/09/22  小高正嗣
C     1998/10/02  小高正嗣
C     1999/01/17  小高正嗣
C
C*********************************************************************
      program tqdvlop
C---------------------------------------------------------------------  
      IMPLICIT REAL*8 ( A-H, O-Z )                                      
C---------------------------------------------------------------------  
#include "grid_size_M.f"
      parameter ( nzprof = nz )                              
      parameter ( nzp = nz+1, ntp = 120 )
C---------------------------------------------------------------------  
C                                                                       
      CHARACTER*10 NMZPRF( NZPROF )                                     
      character*5  nmrun
      character*4  from
      character*1  chara
      character*1  pz
      character*2  end
      real*4 tend

      real*4  dtpad0,dqvad0
      real*4  zgrid(nzp)
      real*4  tptend(nzp),tpmean(nzp),dtpadv(nzp),
     &        dtpdif(nzp),dtpdnl(nzp),
     &        dtpdi0(nzp),dtprad(nzp),tperror(nzp)
      real*8  qvtend(nzp),qvmean(nzp),dqvadv(nzp),
     &        dqvdif(nzp),dqvdnl(nzp),
     &        dqvdi0(nzp),dqvnfl(nzp),qverror(nzp)
      real*4  qvmean_4(nzp)

      integer fnum
      integer ibegin,iend,inter
      real*4 para1
C
C=====================================================================  

      open (11,file='exparam.1')
      read (11,*) irun
      read (11,*) time0, nloop1, nloop2, dtime, ntmoni, ntkubu
      close(11)

      write(6,*) 'data number = ', nloop1
      write(0,*) 'ibegin(1-', nloop1, 
     &           ') iend(1-', nloop1,  
     &           ') inter(1-', nloop1, ')'
      read (5,*) ibegin, iend, inter

      call sgpwsn
      read(*,*) iws

      call gropn(iws)
      call sglset( 'LCORNER', .FALSE. )
c      call sldiv( 'Y', 2, 1 )

      call grfrm
      call grswnd( 225., 270., 0., 10. )
c      call grswnd( 230., 276., 0., 10. )
      call grsvpt( 0.2, 0.8, 0.2, 0.8 )
      call grstrn(1)
      call grstrf
      call ussttl ( 'THETA', '[K]', '', 'Z[km]' )
      call usdaxs

      para1 = 0.95

      do i = ibegin, iend+1, inter

         fnum = 49 + i
         open(fnum)

         read(fnum,'(a2,a15,g15.7,a6,g15.7,a4,9a15)') 
     &        chara, nmrun, t_begin, from, t_end,
     &        pz, (nmzprf(j),j=40,48)

         do iz = 1, nzp
            read(fnum,'(10g15.7)') zgrid(iz), 
     &           tptend(iz),
     &           tpmean(iz),
     &           dtpadv(iz),
     &           dtpdif(iz),
     &           dtpdnl(iz),
     &           dtpad0,
     &           dtpdi0(iz),
     &           dtprad(iz),
     &           tperror(iz)
         end do

         close(fnum)

         tend = t_end / 3600.
         tend = tend - 48 + 6
         write(end,500) INT(tend)
 500     format(I2)

         do iz = 1, nzp
            zgrid(iz) = zgrid(iz) / 1000.
         end do

         call sglset( 'LCHAR', .TRUE. )
         call sgsplc( end ) 
         call sglset ( 'LROT', .TRUE. )
         call sgiset ( 'IROT', 90 )
         call sgrset ( 'FFCT', para1 )
         call uulin(nzp, tpmean, zgrid)
         call sglset( 'LCHAR', .FALSE. )
         call sglset ( 'LROT', .FALSE. )
         para1 = para1 - 0.08

      end do
*
*      

      call grcls


      end