00001 subroutine allocate_matter_emdrsrho
00002   use phys_constant, only : long
00003   use grid_parameter
00004   use def_matter, only : rs, emd, emdg, rhof, rhog
00005   use make_array_2d
00006   use make_array_3d
00007   implicit none
00008 
00009   call alloc_array2d(rs, 0, ntf, 0, npf)
00010   call alloc_array3d(emd, 0, nrf, 0, ntf, 0, npf)
00011   call alloc_array3d(rhof, 0, nrf, 0, ntf, 0, npf)
00012   call alloc_array3d(emdg, 0, nrg, 0, ntg, 0, npg)
00013   call alloc_array3d(rhog, 0, nrg, 0, ntg, 0, npg)
00014 
00015 end subroutine allocate_matter_emdrsrho