I want to compute only specific subroutines in a group of SQUID vector nodes with VH.

You can perform calculations on VH by using VH CALL for specific subroutines. Detailed instructions are posted on the Aurora Forum, so please take a look.
libvhcall-fortran

A simple program for you to verify the operation is provided below.
VH CALL Verification Program

*Compilation Method
module load BaseVEC/2024
sh comp.sh
 
*Submission Method (edit run.sh to change the group setting before submitting)
qsub run.sh