! Maximum Covariance Analysis of Drought Inde:wqx

      Use NetCDF
      Real, Parameter :: Miss=-99.
      Integer, Parameter :: Nvar=6

      Real, Allocatable :: xlon(:), xlat(:), Lon(:), Lat(:), xtime(:)
      Integer :: Nlat, Nlon, ilat, ilon, Ntime,i,j

      Logical :: Status
      !---------------------------------------------------------------- 

      Character (Len=225), parameter :: Fname(Nvar)=(/"Diff_ENS_DD_ORIG.nc","Diff_ENS_DD_QM.nc","Diff_ENS_DS_ORIG.nc","Diff_ENS_DS_QM.nc", &
                                                      "AV_ENS_DD.nc","AV_ENS_DS.nc"/) 
      Character (Len=150) :: Fnm, t_unit, var_ln, var_unit, lon_nm, lat_nm, t_nm
      ! For MCA:
      Real, Allocatable :: X2d(:,:), X3d(:,:,:)        ! Uk in 3D array
      Integer :: Nt, Nx, Ny                 
      Real :: Sumx,Corr,prob
      !--------------------------------------------------------------------

      Do i=1,Nvar
         Call rd_ncdf_2d ("NC/"//trim(Fname(i)),xLon,xLat,Nlon,Nlat, X2d, "av")
         if (i==1) then
            Allocate (X3d(Nlon,Nlat,Nvar),Lon(Nlon),Lat(Nlat),xtime(Nvar))
            Lon=xLon; Lat=xLat
         endif
         Where (IsNaN(X2d)); X2d=Miss; Endwhere
         X3d(:,:,i)=X2d(:,:)
         Deallocate (xlon,xlat,X2d)
      Enddo
      Do i=1,Nvar; xtime(i)=i; Enddo
      t_unit="years since 1977-01-01 "
      var_ln="1: Diff_DD_ORIG; 2: Diff_DD_QM; " & 
           //"3: Diff_DS_ORIG; 4: Diff_DS_QM; 5: AV_DD; 6: AV_DS"
      
      Fnm="DD_DS_AV_ORG_QM.nc"
      Call wr_ncdf_3d (Fnm,lon,lat,xtime,nlon,nlat,Nvar,X3d,"ch", &
             t_unit, var_ln, "-", "lon", "lat", "cha", Miss)
  
      CONTAINS
      !--------------------------------------------------------------------
      Include "/home/tanpv/Mylibs/sub_rd_ncdf_Tan.F90"
      Include "/home/tanpv/Mylibs/sub_wr_ncdf_Tan.F90"
      !--------------------------------------------------------------------
      END