Use netcdf
Integer :: ncid, id_lat_dim, id_lon_dim, id_t_dim, varid, Nlat, Nlon, Ntime, ifile
Integer :: i,j,it, i0, i1, j0, j1, ic, Nx, Ny
Real, Allocatable :: Lat(:), Lon(:), Time(:), X(:,:,:)
Real, Allocatable :: Lat2(:), Lon2(:), Time2(:), Y(:,:,:)
Integer :: Nlat2, Nlon2, Ntime2
Character (Len=255) :: FileName, VarName, LatName, LonName, TName
Real :: Lat0, Lat1, Lon0, Lon1
Real, Allocatable :: XLat(:), XLon(:), XDat(:,:)
!------------------------------------------------------------------------------------

FileName = "../Data/VnGP_0.25_1980_2010_Mon.nc"

Call rd_ncdf_3d (trim(FileName),lon,lat,Time, Nlon, Nlat, Ntime, X, "rain")
print*,">>>>>>>>>>>>>>>>"
Where (X==-99.9) X=-1
Do it=1, Ntime
   Do j=Nlat, 1, -1
      print "(40F4.1)", (X(i,j,it),i=1, Nlon)
   enddo
   pause
Enddo

!------------------------------------------------------------------------------------
CONTAINS

!------------------------------------------------------------------------------------
Include "sub_rd_ncdf.F90"

END

!Include "sub_rd_ncdf.F90"