Dear @Kazuya.Nishiyama,
Regarding the tower-top mass, center of mass, and moments of inertia (BModes input ixx_tip, etc.), you can calculate these through a FAST / OpenFAST linearization analysis by only enabling the platform DOFs and only having mass in your ElastoDyn module pertaining to the rotor and nacelle, e.g., as discussed in the following forum topic: OpenFast 2nd order Linearization.
Regarding the setting of flp_iner, etc. in BModes when computing blade mode shapes for ElastoDyn, see my post dated March 18, 2015 in the following forum topic: NREL 5MW Blade Mode Shape. A similar process should be applied when using BMods to generate tower mode shapes for ElastoDyn, e.g., tor_stff should be set very high relative to the bending stiffness (100 to 1000 times more).
I’m not sure what you mean by a “new manual” for BModes. Which manual are you using?
Best regards,