Test 2's target quantity computed directly from DOT products of PVSH_POL/PVSH_TOR/PVSH_RAD values -- must agree with TEST2_ANALYTIC to machine precision.
Returns: Complex array of size NTH, matching TEST2_ANALYTIC's
grid.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=i4), | intent(in) | :: | NTH |
Number of colatitude samples over . |
FUNCTION TEST2_VSH(NTH) IMPLICIT NONE INTEGER(KIND=i4), INTENT(IN) :: NTH INTEGER(KIND=i4) :: I REAL(KIND=dp) :: TH,PH COMPLEX(KIND=dp), DIMENSION(NTH) :: TEST2_VSH PH = 1.006d0 !- arbitrary choice DO I=1,NTH TH = (I-1)*pi/(NTH-1) TEST2_VSH(I) = DOT(PVSH_POL(2,0,TH,PH),PVSH_TOR(3,2,TH,PH))* & DOT(PVSH_RAD(1,0,TH,PH),CONJG(PVSH_RAD(2,-2,TH,PH))) ENDDO RETURN END FUNCTION TEST2_VSH