00001 subroutine allocate_emfield
00002   use phys_constant, only : long
00003   use grid_parameter
00004   use def_emfield
00005   use def_emfield_derivatives
00006   use def_faraday_tensor
00007   use make_array_3d
00008   use make_array_4d
00009   use make_array_5d
00010   implicit none
00011 
00012   call alloc_array3d(va , 0, nrg, 0, ntg, 0, npg)
00013   call alloc_array3d(alva,0, nrg, 0, ntg, 0, npg)
00014   call alloc_array3d(vaxd,0, nrg, 0, ntg, 0, npg)
00015   call alloc_array3d(vayd,0, nrg, 0, ntg, 0, npg)
00016   call alloc_array3d(vazd,0, nrg, 0, ntg, 0, npg)
00017   call alloc_array3d(vaxu,0, nrg, 0, ntg, 0, npg)
00018   call alloc_array3d(vayu,0, nrg, 0, ntg, 0, npg)
00019   call alloc_array3d(vazu,0, nrg, 0, ntg, 0, npg)
00020   call alloc_array3d(jtd, 0, nrg, 0, ntg, 0, npg)
00021   call alloc_array3d(jxd, 0, nrg, 0, ntg, 0, npg)
00022   call alloc_array3d(jyd, 0, nrg, 0, ntg, 0, npg)
00023   call alloc_array3d(jzd, 0, nrg, 0, ntg, 0, npg)
00024   call alloc_array3d(jtu, 0, nrg, 0, ntg, 0, npg)
00025   call alloc_array3d(jxu, 0, nrg, 0, ntg, 0, npg)
00026   call alloc_array3d(jyu, 0, nrg, 0, ntg, 0, npg)
00027   call alloc_array3d(jzu, 0, nrg, 0, ntg, 0, npg)
00028   call alloc_array3d(jtuf, 0, nrf, 0, ntf, 0, npf)
00029   call alloc_array3d(jxuf, 0, nrf, 0, ntf, 0, npf)
00030   call alloc_array3d(jyuf, 0, nrf, 0, ntf, 0, npf)
00031   call alloc_array3d(jzuf, 0, nrf, 0, ntf, 0, npf)
00032 
00033   call alloc_array3d(fxd, 1, nrg, 1, ntg, 1, npg)
00034   call alloc_array3d(fyd, 1, nrg, 1, ntg, 1, npg)
00035   call alloc_array3d(fzd, 1, nrg, 1, ntg, 1, npg)
00036   call alloc_array3d(fxu, 1, nrg, 1, ntg, 1, npg)
00037   call alloc_array3d(fyu, 1, nrg, 1, ntg, 1, npg)
00038   call alloc_array3d(fzu, 1, nrg, 1, ntg, 1, npg)
00039   call alloc_array4d(fijd,1, nrg, 1, ntg, 1, npg, 1, 3) 
00040   call alloc_array4d(fiju,1, nrg, 1, ntg, 1, npg, 1, 3) 
00041   call alloc_array5d(fijdu, 1, nrg, 1, ntg, 1, npg, 1, 3, 1, 3) 
00042   call alloc_array3d(fidfiu,1, nrg, 1, ntg, 1, npg)
00043   call alloc_array3d(fijfij,1, nrg, 1, ntg, 1, npg)
00044   call alloc_array3d(fxd_grid, 0, nrg, 0, ntg, 0, npg)
00045   call alloc_array3d(fyd_grid, 0, nrg, 0, ntg, 0, npg)
00046   call alloc_array3d(fzd_grid, 0, nrg, 0, ntg, 0, npg)
00047   call alloc_array3d(fxu_grid, 0, nrg, 0, ntg, 0, npg)
00048   call alloc_array3d(fyu_grid, 0, nrg, 0, ntg, 0, npg)
00049   call alloc_array3d(fzu_grid, 0, nrg, 0, ntg, 0, npg)
00050   call alloc_array4d(fijd_grid,0, nrg, 0, ntg, 0, npg, 1, 3) 
00051   call alloc_array4d(fiju_grid,0, nrg, 0, ntg, 0, npg, 1, 3) 
00052   call alloc_array5d(fijdu_grid, 0, nrg, 0, ntg, 0, npg, 1, 3, 1, 3) 
00053   call alloc_array3d(fidfiu_grid,0, nrg, 0, ntg, 0, npg)
00054   call alloc_array3d(fijfij_grid,0, nrg, 0, ntg, 0, npg)
00055   call alloc_array3d(Lie_bFxd, 1, nrg, 1, ntg, 1, npg)
00056   call alloc_array3d(Lie_bFyd, 1, nrg, 1, ntg, 1, npg)
00057   call alloc_array3d(Lie_bFzd, 1, nrg, 1, ntg, 1, npg)
00058 
00059   call alloc_array4d(pdvaxd, 1, nrg, 1, ntg, 1, npg, 1, 3) 
00060   call alloc_array4d(pdvayd, 1, nrg, 1, ntg, 1, npg, 1, 3) 
00061   call alloc_array4d(pdvazd, 1, nrg, 1, ntg, 1, npg, 1, 3) 
00062   call alloc_array4d(cdvaxd, 1, nrg, 1, ntg, 1, npg, 1, 3) 
00063   call alloc_array4d(cdvayd, 1, nrg, 1, ntg, 1, npg, 1, 3) 
00064   call alloc_array4d(cdvazd, 1, nrg, 1, ntg, 1, npg, 1, 3) 
00065   call alloc_array3d(Lie_bAxd, 1, nrg, 1, ntg, 1, npg)
00066   call alloc_array3d(Lie_bAyd, 1, nrg, 1, ntg, 1, npg)
00067   call alloc_array3d(Lie_bAzd, 1, nrg, 1, ntg, 1, npg)
00068   call alloc_array3d(Lie_bAxd_grid, 0, nrg, 0, ntg, 0, npg)
00069   call alloc_array3d(Lie_bAyd_grid, 0, nrg, 0, ntg, 0, npg)
00070   call alloc_array3d(Lie_bAzd_grid, 0, nrg, 0, ntg, 0, npg)
00071   call alloc_array3d(gLie_bAxu_grid, 0, nrg, 0, ntg, 0, npg)
00072   call alloc_array3d(gLie_bAyu_grid, 0, nrg, 0, ntg, 0, npg)
00073   call alloc_array3d(gLie_bAzu_grid, 0, nrg, 0, ntg, 0, npg)
00074 end subroutine allocate_emfield