PVSH_RAD_ALL_DP Subroutine

private subroutine PVSH_RAD_ALL_DP(OUT, LMAX, THETA, PHI)

All radial polar VSH members PVSH_RAD for , from one VSH_CORE pass.

Warning

Inherits ASSOC_LEGENDRE_NORM_ALL's -dependent stability limit (safe for unconditionally) -- see STABILITY_FINDINGS.md at the repo root.

1==, 2= (always 0), 3= (always 0); mode axis indexed by YLM_INDEX(l,m).

Arguments

Type IntentOptional Attributes Name
complex(kind=dp), intent(out) :: OUT(3,(LMAX+1)**2)
integer(kind=i4), intent(in) :: LMAX

Maximum degree, .
real(kind=dp), intent(in) :: THETA

Colatitude in radians, .
real(kind=dp), intent(in) :: PHI

Longitude in radians, .

Calls

proc~~pvsh_rad_all_dp~~CallsGraph proc~pvsh_rad_all_dp PVSH_RAD_ALL_DP interface~vsh_core VSH_CORE proc~pvsh_rad_all_dp->interface~vsh_core proc~vsh_core_dp VSH_CORE_DP interface~vsh_core->proc~vsh_core_dp dcos dcos proc~vsh_core_dp->dcos dsin dsin proc~vsh_core_dp->dsin interface~assoc_legendre_norm_all ASSOC_LEGENDRE_NORM_ALL proc~vsh_core_dp->interface~assoc_legendre_norm_all interface~ddx_assoc_legendre_norm_all DDX_ASSOC_LEGENDRE_NORM_ALL proc~vsh_core_dp->interface~ddx_assoc_legendre_norm_all proc~plm_index PLM_INDEX proc~vsh_core_dp->proc~plm_index proc~ylm_index YLM_INDEX proc~vsh_core_dp->proc~ylm_index proc~assoc_legendre_norm_all_dp ASSOC_LEGENDRE_NORM_ALL_DP interface~assoc_legendre_norm_all->proc~assoc_legendre_norm_all_dp proc~ddx_assoc_legendre_norm_all_dp DDX_ASSOC_LEGENDRE_NORM_ALL_DP interface~ddx_assoc_legendre_norm_all->proc~ddx_assoc_legendre_norm_all_dp proc~assoc_legendre_norm_all_dp->proc~plm_index proc~ddx_assoc_legendre_norm_all_dp->proc~plm_index

Called by

proc~~pvsh_rad_all_dp~~CalledByGraph proc~pvsh_rad_all_dp PVSH_RAD_ALL_DP interface~pvsh_rad_all PVSH_RAD_ALL interface~pvsh_rad_all->proc~pvsh_rad_all_dp proc~batch_pvsh_rad_cons BATCH_PVSH_RAD_CONS proc~batch_pvsh_rad_cons->interface~pvsh_rad_all proc~decompose DECOMPOSE proc~decompose->interface~pvsh_rad_all proc~decompose~2 DECOMPOSE proc~decompose~2->interface~pvsh_rad_all proc~reconstruct RECONSTRUCT proc~reconstruct->interface~pvsh_rad_all proc~reconstruct~2 RECONSTRUCT proc~reconstruct~2->interface~pvsh_rad_all proc~vsh_pol_inversion VSH_POL_INVERSION proc~vsh_pol_inversion->interface~pvsh_rad_all proc~run_part_a RUN_PART_A proc~run_part_a->proc~decompose proc~run_part_a~2 RUN_PART_A proc~run_part_a~2->proc~decompose~2 proc~run_part_b RUN_PART_B proc~run_part_b->proc~decompose proc~run_part_b->proc~reconstruct proc~run_part_b~2 RUN_PART_B proc~run_part_b~2->proc~decompose~2 proc~run_part_b~2->proc~reconstruct~2 program~main MAIN program~main->proc~batch_pvsh_rad_cons program~main->proc~vsh_pol_inversion program~vsh_decomposition VSH_DECOMPOSITION program~vsh_decomposition->proc~run_part_a program~vsh_decomposition->proc~run_part_b program~vsh_decomposition_tor VSH_DECOMPOSITION_TOR program~vsh_decomposition_tor->proc~run_part_a~2 program~vsh_decomposition_tor->proc~run_part_b~2

Source Code

  SUBROUTINE PVSH_RAD_ALL_DP(OUT, LMAX, THETA, PHI)
  IMPLICIT NONE
  INTEGER(KIND=i4), INTENT(IN)  :: LMAX
  REAL(KIND=dp),    INTENT(IN)  :: THETA
  REAL(KIND=dp),    INTENT(IN)  :: PHI
  COMPLEX(KIND=dp), INTENT(OUT) :: OUT(3,(LMAX+1)**2)
  INTEGER(KIND=i4) :: NYLM
  COMPLEX(KIND=dp), ALLOCATABLE :: YLM(:), GTH(:), GPH(:)
  NYLM = (LMAX+1)**2
  ALLOCATE(YLM(NYLM), GTH(NYLM), GPH(NYLM))
  CALL VSH_CORE(YLM, GTH, GPH, LMAX, THETA, PHI)
  OUT(1,:) = YLM
  OUT(2,:) = DCMPLX(0.d0, 0.d0)
  OUT(3,:) = DCMPLX(0.d0, 0.d0)
  DEALLOCATE(YLM, GTH, GPH)
  END SUBROUTINE PVSH_RAD_ALL_DP