units

I am using an open source FEM platform, which requires you to convert your equation system to non-dimensional form. So, there are no units specified for the parameters in the problem. If you use compatible input units, you should expect the output in matching unit - just like one would do in a normal python or C++ code for solving physical problems. I have an elasto-plastic dynamic problem to solve, and these are the parameters for the Drucker-Prager Soil Model for soft clay.

E = 5000000;                 // Young's Modulus in N/m^2 or Pa

I'm using the ABAQUS DLOAD subroutines, and I'm curious what the dimensions of the 

volume distributed body load should be. 

I understand ABAQUS uses [LT-2] for GRAV body loads. However, I'm curious if the 

units for BZNU acting over a volume should be [density]*[acceleration] or just

[acceleration].

Thank you,

Anirban