! $Id: dop2phs.f90 2022-10-19 09:51:28Z ychen $

!****s* fsi/dop2phs *
!
! NAME
!   dop2phs  - Reconstruct phase with smoothed Doppler
!
! SYNOPSIS
!    call dop2phs(phase, dop_sm, phase_sm)
!
! DESCRIPTION
!    This subroutine reconstruct phase with smoothed Doppler
!
! INPUTS
!    REAL(wp), DIMENSION(:,:),  :: phase     Excess phase (m) [ch, t]
!    REAL(wp), DIMENSION(:,:),  :: dop_sm    smoothed doppler shift (m) [ch, t]
!
! OUTPUT
!    REAL(wp), DIMENSION(:,:),  :: phase_sm  Reconstruct phase after smooth (m) [ch, t]
!
! AUTHOR
!   Yong Chen, NOAA/NESDIS/STAR, yong.chen@noaa.gov
!   Loknath Adhikari, UMD/CISESS, loknath.adhikari@noaa.gov
!
! COPYRIGHT
!   Copyright (c) 2022-2023 Yong Chen <yong.chen@noaa.gov>
!   For further details please refer to the file COPYRIGHT
!   which you should have received as part of this distribution.
!
!****

SUBROUTINE dop2phs(phase, dop_sm, phase_sm)

      USE typesizes,         ONLY: wp => EightByteReal

      IMPLICIT NONE
      REAL(wp), DIMENSION(:,:), INTENT(in)    :: phase    ! Excess phase (m) [ch, t]
      REAL(wp), DIMENSION(:,:), INTENT(in)    :: dop_sm   ! smoothed doppler shift (m) [ch, t]
      REAL(wp), DIMENSION(:,:), INTENT(inout) :: phase_sm ! Reconstruct phase after smooth (m) [ch, t]

      INTEGER  :: n, nc
      INTEGER  :: ic, i

      n  = SIZE(dop_sm, 2)   ! time dimension
      nc = SIZE(dop_sm, 1)   ! 2 frequency channels


      DO ic=1,nc
        phase_sm(ic, 1) = phase(ic, 1)
        DO i = 2, n
         phase_sm(ic, i) = phase_sm(ic, i-1) + dop_sm(ic, i)
        ENDDO
      ENDDO

END SUBROUTINE dop2phs