program testread_john
  use wrf_data_pnc
  implicit none
  include 'wrf_status_codes.h'
  include 'netcdf.inc'
  character (80) FileName
  integer Comm
  character (80) SysDepInfo
  integer     :: DataHandle
  integer Status
  integer NCID
  real data(200)
  integer idata(200)
  real*8 ddata(200)
  logical ldata(200)
  character (80) cdata
  integer OutCount
  integer i,j,k

  integer, parameter ::    pad = 3
  integer, parameter ::    jds=1       , jde=6      , &
                           ids=1       , ide=9      , &
                           kds=1       , kde=5
  integer, parameter ::    jms=jds-pad , jme=jde+pad , &
                           ims=ids-pad , ime=ide+pad , &
                           kms=kds     , kme=kde
  integer, parameter ::    jps=jds     , jpe=jde    , &
                           ips=ids     , ipe=ide    , &
                           kps=kds     , kpe=kde

  real u( ims:ime , kms:kme , jms:jme )
  real v( ims:ime , kms:kme , jms:jme )
  real rho( ims:ime , kms:kme , jms:jme )
  real u2( ims:ime , jms:jme )
  real u1( ims:ime )

  integer int( ims:ime , kms:kme , jms:jme )
  real*8  r8 ( ims:ime , kms:kme , jms:jme )

  integer Dom
  character*3 MemOrd
  integer , Dimension(3) :: DomS,DomE,MemS,MemE,PatS,PatE
  integer , Dimension(2) :: Dom2S,Dom2E,Mem2S,Mem2E,Pat2S,Pat2E
  integer                   Dom1S,Dom1E,Mem1S,Mem1E,Pat1S,Pat1E
  character (19) Time, DateStr
  character (31) VarName
  character (19) Date

  print *, 'Testing wrf read'
  Date = '2000-09-18_16:42:01'
  call ext_init(Status)
  print *,'After call ext_init, Status =',Status
  FileName = 'foo.nc'
  Comm = 1
  SysDepInfo = 'sys info'
  call ext_open_for_read( FileName, Comm, SysDepInfo, DataHandle, Status)
  print *, 'Status = ',Status,DataHandle

  MemOrd = "XZY"

  DomS(1) = ids
  DomE(1) = ide
  DomS(2) = kds
  DomE(2) = kde
  DomS(3) = jds
  DomE(3) = jde

  PatS(1) = ips
  PatE(1) = ipe
  PatS(2) = kps
  PatE(2) = kpe
  PatS(3) = jps
  PatE(3) = jpe

  MemS(1) = ims
  MemE(1) = ime
  MemS(2) = kms
  MemE(2) = kme
  MemS(3) = jms
  MemE(3) = jme

  Dom2S(1) = ids
  Dom2S(2) = jds
  Dom2E(1) = ide
  Dom2E(2) = jde
  Mem2S(1) = ims
  Mem2S(2) = jms
  Mem2E(1) = ime
  Mem2E(2) = jme
  Pat2S(1) = ips
  Pat2S(2) = jps
  Pat2E(1) = ipe
  Pat2E(2) = jpe

  Dom1S = ids
  Dom1E = ide
  Mem1S = ims
  Mem1E = ime
  Pat1S = ips
  Pat1E = ipe

  call ext_get_next_time(DataHandle, Time, Status)
  print *, Time, Status

  call ext_read_field(DataHandle,Time,'u',u,WRF_REAL,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'ext_read_field Status = ',Status, ' u(2,3,4) ', u(2,3,4) 
  call ext_read_field(DataHandle,Time,'v',v,WRF_REAL,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'ext_read_field Status = ',Status, ' v(4,3,2) ', v(4,3,2)
  call ext_read_field(DataHandle,Time,'rho',rho,WRF_REAL,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'ext_read_field Status = ',Status, ' rho(3,4,5) ' , rho(3,4,5)
  call ext_read_field(DataHandle,Date,'u2',u2,WRF_REAL,Comm,Dom,'XY',Dom2S,Dom2E,Mem2S,Mem2E,Pat2S,Pat2E,Status)
  print *,'ext_read_field Status = ',Status, ' u2(6,5) ', u2(6,5) 
  call ext_read_field(DataHandle,Date,'ud2',u,WRF_REAL,Comm,Dom,"XzY",DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'ext_read_field Status = ',Status, ' u(2,3,4) ', u(2,3,4) 
  call ext_read_field(DataHandle,Date,'u1',u1,WRF_REAL,Comm,Dom,'Z',Dom1S,Dom1E,Mem1S,Mem1E,Pat1S,Pat1E,Status)
  print *,'ext_read_field Status = ',Status, ' u1(9) ', u1(9)

  call ext_read_field(DataHandle,Time,'int',int,WRF_INTEGER,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'ext_read_field Status = ',Status, ' int(8,5,6) ', int(8,5,6) 
  call ext_read_field(DataHandle,Time,'double',r8,WRF_DOUBLE,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'ext_read_field Status = ',Status, ' r8(7,4,5) ', r8(7,4,5) 

  call ext_get_next_time(DataHandle, Time, Status)
  print *, Time, Status

  call ext_read_field(DataHandle,Time,'u',u,WRF_REAL,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'> ext_read_field Status = ',Status, ' u(3,3,3) ' ,u(3,3,3)
  call ext_read_field(DataHandle,Time,'v',v,WRF_REAL,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'> ext_read_field Status = ',Status, ' v(4,4,4) ' ,v(4,4,4)
  call ext_read_field(DataHandle,Time,'rho',rho,WRF_REAL,Comm,Dom,MemOrd,DomS,DomE,MemS,MemE,PatS,PatE,Status)
  print *,'> ext_read_field Status = ',Status, ' rho(3,4,5) ' ,rho(3,4,5)

  call ext_close( DataHandle, Status)
  print *, 'After ext_close, Status = ',Status
  call ext_exit(Status)
  print *,'End of test program',Status
  stop
  end program testread_john
