femwind/fortran/netcdf_test.f90

   1 program netcdf_test
   2
   3 use module_utils
   4 use module_netcdf
   5 use module_wrfout
   6
   7 !*** local
   8 character(len=128)::filename
   9 integer::ncid,frame
  10 real, pointer, dimension(:,:,:)::u0_fmw,v0_fmw,w0_fmw
  11
  12 !*** executable
  13 filename = "wrf.nc"
  14
  15 call ncopen(filename,nf90_nowrite,ncid)
  16
  17 do frame = 1,9999
  18
  19   if(read_initial_wind(ncid,u0_fmw,v0_fmw,w0_fmw,frame=frame).ne.0)then
  20     print *,'check sum does not agree'
  21     stop 1
  22   endif
  23
  24 enddo
  25
  26 call ncclose(ncid)
  27
  28 end program netcdf_test