other/tign_history_api/tign_test.F90

   1 ! This program tests the lfn_history api.  It assumes a valid wrfinput_d01
   2 ! file exists in the current path.
   3
   4 program tign_test
   5 use tign_history
   6 implicit none
   7
   8 integer,parameter::idom=1
   9 integer::nx,ny,i,mtime,sr_x,sr_y
  10 real,allocatable,dimension(:,:)::tign
  11 real::dx,dy,dt
  12 real,parameter::testvalue=11.5
  13
  14 call get_grid_info(idom,nx=nx,ny=ny,sr_x=sr_x,sr_y=sr_y,dx=dx,dy=dy,dt=dt)
  15 allocate(tign(nx,ny))
  16
  17 tign(:,:)=testvalue
  18
  19 call write_tign_history(idom,nx,ny,tign)
  20
  21 tign(:,:)=-1
  22
  23 call read_tign_history(idom,nx,ny,tign)
  24
  25 if(any(tign(:,:).ne.testvalue))then
  26   print*,'ERROR tign set incorrectly'
  27   call abort()
  28 endif
  29 print*,'TEST SUCCESSFUL'
  30 end program tign_test