anchor.diff:              
               iSub = mod(I-1,6)+1                         ! Finds the 6x6 sub-matrix indexing for the SttcDim multiplier matrix
               jSub = mod(J-1,6)+1  
  -!TODO need to transform this when PtfmRefztRot is nonzero per plan
  +
               p%HdroSttc (I,J) = TmpData1*SttcDim(iSub,jSub)    ! Redimensionalize the data and place it at the appropriate location within the array
   
            ELSE                                           ! We must have reached the end of the file, so stop reading in data
  
positive.code:         
     
  CONTAINS
  !----------------------------------------------------------------------------------------------------------------------------------
  !> This routine is called at the start of the simulation to perform initialization steps. 
  !! The parameters are set here and not changed during the simulation.
  
negative.code: 
        CLOSE ( UnWh ) ! Close file.
  
           ! need to transform p%HdroSttc when PtfmRefztRot is nonzero per plan
        call TransformWAMITMatrices( p%NBody, InitInp%PtfmRefztRot, p%HdroSttc )
  
           ! Linear, frequency-dependent hydrodynamic added mass and damping from the
           !   radiation problem:
  