PROGRAM GRADS_READ

IMPLICIT NONE

INTEGER, PARAMETER ::  NXP=137,NYP=134,NZP=50,NT=864,NG=12,NRECSIZE=NXP*NYP*4
INTEGER :: TT,ICOUNT,K,i,j
REAL, DIMENSION(NZP,NXP,NYP) :: PRESS,U_AVG_TN,V_AVG_TN,THETA,U_AVG,V_AVG,W_AVG,TEMPK,SGS_TKE
CHARACTER*80  filea, fileo


FILEA='Jezero_Mars2020_Ls0-atm-S-g07.dat'
FILEO='Jezero_press_Ls0_review.txt'

OPEN(unit=21,file=filea,access='direct',recl=nrecsize)
OPEN(unit=22,file=fileo)

ICOUNT=0
DO TT=1,NT
      print *,'Reading Time ',tt
      ! read temperature
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) press(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) u_avg_tN(k,:,:)
      ENDDO
      ! read  theta
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) v_avg_tN(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) theta(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) u_avg(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) v_avg(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) w_avg(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) tempk(k,:,:)
      ENDDO
      DO K=1,NZP
         ICOUNT=ICOUNT+1  !increase record for each z slab read
         read(21,rec=ICOUNT) sgs_tke(k,:,:)
      ENDDO
      IF (TT > 289 .AND. TT < 578) THEN
        !print *,tempk (2,73,38)
        !write (22,*) tempk(3,69,26)
        write (22,*) press(2,18,104)
        !write (22,*) v(2,69,26)      
      END IF
ENDDO
CLOSE(21)

END