Subroutine Write_2D_SenSlp (Fname)
! Real, Dimension (Nst,Nmon) :: Sen_Mon,  Pro_Mon   ! Sen's slopes and p-values  

Integer :: ist, imon
Real :: Xlon(Nmon), Xlat(Nst), X(Nst,Nmon)
Character (Len=*) :: Fname
Character (Len=255) :: Var_ln

Do imon=1, Nmon; Xlat(imon)=imon; Enddo   ! Mang Lon
Do ist=1, Nst; Xlon(ist)=ist; Enddo ! Mang Lat

!------------------------------------------------------
Var_ln="Sens Slopes of PDSI for each month"
Do imon=1, Nmon; X(1:Nst,imon)=Sen_Mon(1:Nst,imon); Enddo
Call wr_ncdf_2d (trim(Fname), Xlon, Xlat, Nst, Nmon, X, "SenSlp", trim(Var_ln), "-/year", "Stn", "mon", Miss)

Do imon=1, Nmon; X(1:Nst,imon)=Pro_Mon(1:Nst,imon); Enddo
Var_ln="P-values of PDSI trend"
Call wr_ncdf_2d_apd (trim(Fname), X, "P_value", Nst, Nmon, trim(Var_ln), "-", "Stn", "mon", Miss)

End Subroutine Write_2D_SenSlp