program readg5
! data over the domain 180W-180E and 50-90 E
! there are 4 time segments in each day: 0000z, 0600z, 1200z, 1800z
        integer, parameter :: ix = 540, iy = 81, iz = 72, it = 4
        real*8 lat(iy)
        real*8 lon(ix)
        real*8 lev(iz)
        real*8 time
        real*8 delp(ix,iy,iz,it)
        real*8 ad(ix,iy,iz,it)
        real*8 co(ix,iy,iz,it)
        real*8 pp(ix,iy,iz)
        character*50 dir
        character*50 fl_in

        dir ='/misc/mcc21/bian/chu/'
        fl_in = 'GEOS5_ARCTAS_co_20080701.dat'
        open(1,file=TRIM(dir)//TRIM(fl_in))
        print *, 'file=',TRIM(dir)//TRIM(fl_in)
        do n = 1,4 
        if (n .eq. 1) then
            read(1,'(A)') 
            read(1,'(10f8.2)') lat
        !    print *,'lat=',lat
            read(1,'(A)') 
            read(1,'(10f8.2)') lon
        !    print *,'lon=',lon
            read(1,'(A)') 
            read(1,'(10f8.2)') lev
        !    print *,'lev=',lev
        endif
            read(1,'(A)') 
          read(1,'(8e10.3)') delp(:,:,:,n)
            read(1,'(A)') 
          read(1,'(8e10.3)') ad(:,:,:,n)
            read(1,'(A)') 
          read(1,'(8e10.3)') co(:,:,:,n)
        enddo
        close(1)
        print *,'delp=',delp(10,10,:,3)
        do i=1,ix
        do j=1,iy
        do k=1,iz
           pp(i,j,k)=sum(delp(i,j,1:k,3))*0.01  ! from Pa to hPa
        enddo
        enddo
        enddo
        !print *,'co=',co(1,1,1,:)
        print *,'pp=',pp(10,10,:)
        stop
        end