A student's question about Slepian in GRACE

[aNannan6]

      Excuse me ,I am a graduate student of geoscience in China. I had the honor to listen to F. J. Simons's lecture —— "Study of galaxy lithosphere using Slepian function" at the University of Chinese Academy of Sciences. I am very interested in the application of Slepian function method in earth science. I also want to use this method to calculate gravity changes and restore the time-varying gravity field model in local areas. Therefore, I am currently learning the software package (slepian_delta) for GRACE products for free, and I have read your article (doi:10.1073/pnas.1206785109), but there are some questions in this learning process. I would like to ask you. I hope you can take time out of your busy schedule to read the email and reply to the email.
    Problem content is mainly, the regional gravity change spatial distribution map, such as your article Greenland glacier total change figure (as shown in figure 2), I would like to know if I input the GRACE product, calculated the base function (G), and the corresponding slepian expansion coefficient (slepcoffs), the gravity value f is how to calculate, also is the specification mentioned the fifth point (as shown)
    Here's my idea to see if it's correct. If I want to calculate the gravity value, then the gravity value f=slepcoffs' * G, G is the matrix of addmout form, then f is also the matrix of addmout form, convert f of addmout form using coef2lmcosi() to F of lmcosi form, and then draw with plm2xyz() and plotplm().

Then I want to ask the unit of calculation is how much, how to pay attention to the problem of the unit?
If you could take time out of your busy schedule to reply to my email and answer my stupid question, I would really appreciate it！