Our main program is F4.gp, the other two GAP programs are used for the enumeration of conjugacy classes of some permutation groups.
Reading F4.gp in PARI/GP, then you can use the function RepDim to calculate the dimensions of invariant subspaces:
Here we also provide a larger table containing the list of [w, d1(w), d2(w), d(w)], where w=[a, b, c, d] is a dominant weight of F4 such that the d(w) is nonzero and 2*a+3*b+2*c+d<=40.
Admitting a series of conjectures of Arthur and Langlands, we have the program CountPara.gp to calculate the number F4(w) in the paper (need to load F4.gp at first):
In this program we use some tables in Discrete series multiplicities for classical groups over Z and level 1 algebraic cusp forms so if you want to calculate for larger weights you need to ask them for larger tables.