sppad.f

Go to the documentation of this file.

C> @file
C> @brief Pad or truncate a spectral field.
C> @author Iredell @date 92-10-31
 
C> Pad or truncate a spectral field.
C>
C> @param I1 input spectral domain shape
C> (0 for triangular, 1 for rhomboidal)
C> @param M1 input spectral truncation
C> @param Q1 ((M+1)*((I+1)*M+2)) input field
C> @param I2 output spectral domain shape
C> (0 for triangular, 1 for rhomboidal)
C> @param M2 output spectral truncation
C> @param Q2 ((M+1)*((I+1)*M+2)) output field
C>
C> @author Iredell @date 92-10-31
      SUBROUTINE sppad(I1,M1,Q1,I2,M2,Q2)
 
      REAL Q1((M1+1)*((I1+1)*M1+2))
      REAL Q2((M2+1)*((I2+1)*M2+2))
 
      DO l=0,m2
        DO n=l,i2*l+m2
          ks2=l*(2*m2+(i2-1)*(l-1))+2*n
          IF(l.LE.m1.AND.n.LE.i1*l+m1) THEN
            ks1=l*(2*m1+(i1-1)*(l-1))+2*n
            q2(ks2+1)=q1(ks1+1)
            q2(ks2+2)=q1(ks1+2)
          ELSE
            q2(ks2+1)=0
            q2(ks2+2)=0
          ENDIF
        ENDDO
      ENDDO
      RETURN
      END