Use Netcdf

   Real, Allocatable :: slp(:,:,:), t2m(:,:,:)
   Real, Allocatable :: lon(:), lat(:), xtime(:)
   Integer :: Nlon, Nlat, Nt
   Character (Len=255) :: Fname

   Real, Allocatable :: pmsl(:,:), temp(:,:)
   !-------------------------------------------------

   Fname="wrfout/netcdf/regcm_25km_2019101600_6hourly.nc"
   Call rd_ncdf_3d (Fname,lon,lat,xtime,nlon,nlat,nt,slp,"slp")
   Deallocate (lon,lat,xtime)
   Call rd_ncdf_3d (Fname,lon,lat,xtime,nlon,nlat,nt,t2m,"t2")



   CONTAINS

   Include "/home/tanpv/Mylibs/sub_rd_ncdf_Tan.F90"
   Include "/home/tanpv/Mylibs/sub_wr_ncdf_Tan.F90"

   END