GRAD_SSH_ALL_DP Subroutine

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

All angular gradients (GRAD_SSH) 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= (always 0), 2=, 3=; 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~~grad_ssh_all_dp~~CallsGraph proc~grad_ssh_all_dp GRAD_SSH_ALL_DP interface~vsh_core VSH_CORE proc~grad_ssh_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~~grad_ssh_all_dp~~CalledByGraph proc~grad_ssh_all_dp GRAD_SSH_ALL_DP interface~grad_ssh_all GRAD_SSH_ALL interface~grad_ssh_all->proc~grad_ssh_all_dp proc~batch_grad_ssh_cons BATCH_GRAD_SSH_CONS proc~batch_grad_ssh_cons->interface~grad_ssh_all program~main MAIN program~main->proc~batch_grad_ssh_cons

Source Code

  SUBROUTINE GRAD_SSH_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,:) = DCMPLX(0.d0, 0.d0)
  OUT(2,:) = GTH
  OUT(3,:) = GPH
  DEALLOCATE(YLM, GTH, GPH)
  END SUBROUTINE GRAD_SSH_ALL_DP