The extraction and conversion of FIA data to css and pss files for ED2 model works.

This section is incorrect, you should not be filling in the soil C and N variables with zero's. The previous code, which used matrix(soil,...) was correct and you should revert to the old code. The only problem with the previous code was aesthetic, since the soil matrix didn't have named columns, but the model doesn't read those column headers. This problem could be solved more elegantly by just pre-computing the soil matrix and then giving it column names before the cbind.