00001 subroutine allocate_metric_and_matter_WL
00002 use phys_constant, only : long
00003 use grid_parameter
00004 use def_metric
00005 use def_matter
00006 use def_gamma_crist
00007 use def_gamma_crist_grid
00008 use def_metric_hij
00009 use def_metric_rotshift
00010 use def_ricci_tensor
00011 use def_metric_excurve_grid
00012 use def_shift_derivatives
00013 use def_shift_derivatives_grid
00014 use def_Lie_derivatives
00015 use def_Lie_derivatives_grid
00016 use def_matter_parameter
00017 use def_cristoffel
00018 use def_cristoffel_grid
00019 use def_matter
00020 use def_vector_x
00021 use def_vector_phi
00022 use make_array_2d
00023 use make_array_3d
00024 use make_array_4d
00025 use make_array_5d
00026 implicit none
00027
00028 call alloc_array2d(rs, 0, ntf, 0, npf)
00029 call alloc_array3d(emdg, 0, nrg, 0, ntg, 0, npg)
00030 call alloc_array3d(emd, 0, nrf, 0, ntf, 0, npf)
00031 call alloc_array3d(utg, 0, nrg, 0, ntg, 0, npg)
00032 call alloc_array3d(utf, 0, nrf, 0, ntf, 0, npf)
00033 call alloc_array3d(lambda, 0, nrf, 0, ntf, 0, npf)
00034 call alloc_array3d(psi, 0, nrg, 0, ntg, 0, npg)
00035 call alloc_array3d(alph, 0, nrg, 0, ntg, 0, npg)
00036 call alloc_array3d(alps, 0, nrg, 0, ntg, 0, npg)
00037 call alloc_array3d(tfkijkij, 1, nrg, 1, ntg, 1, npg)
00038 call alloc_array3d(bvxd, 0, nrg, 0, ntg, 0, npg)
00039 call alloc_array3d(bvyd, 0, nrg, 0, ntg, 0, npg)
00040 call alloc_array3d(bvzd, 0, nrg, 0, ntg, 0, npg)
00041 call alloc_array5d(tfkij, 1, nrg, 1, ntg, 1, npg, 1, 3, 1, 3)
00042
00043 call alloc_array3d(alps2, 0, nrg, 0, ntg, 0, npg)
00044 call alloc_array3d(bvxu, 0, nrg, 0, ntg, 0, npg)
00045 call alloc_array3d(bvyu, 0, nrg, 0, ntg, 0, npg)
00046 call alloc_array3d(bvzu, 0, nrg, 0, ntg, 0, npg)
00047 call alloc_array3d(tfkijkij_grid, 0, nrg, 0, ntg, 0, npg)
00048 call alloc_array5d(tfkij_grid, 0, nrg, 0, ntg, 0, npg, 1, 3, 1, 3)
00049 call alloc_array5d(cri, 1, nrg, 1, ntg, 1, npg, 1, 3, 1, 6)
00050 call alloc_array5d(cri_grid, 0, nrg, 0, ntg, 0, npg, 1, 3, 1, 6)
00051 call alloc_array3d(ovxu, 0, nrg, 0, ntg, 0, npg)
00052 call alloc_array3d(ovyu, 0, nrg, 0, ntg, 0, npg)
00053 call alloc_array3d(ovzu, 0, nrg, 0, ntg, 0, npg)
00054 call alloc_array3d(ovxd, 0, nrg, 0, ntg, 0, npg)
00055 call alloc_array3d(ovyd, 0, nrg, 0, ntg, 0, npg)
00056 call alloc_array3d(ovzd, 0, nrg, 0, ntg, 0, npg)
00057 call alloc_array3d(hxxd, 0, nrg, 0, ntg, 0, npg)
00058 call alloc_array3d(hxyd, 0, nrg, 0, ntg, 0, npg)
00059 call alloc_array3d(hxzd, 0, nrg, 0, ntg, 0, npg)
00060 call alloc_array3d(hyyd, 0, nrg, 0, ntg, 0, npg)
00061 call alloc_array3d(hyzd, 0, nrg, 0, ntg, 0, npg)
00062 call alloc_array3d(hzzd, 0, nrg, 0, ntg, 0, npg)
00063 call alloc_array3d(hxxu, 0, nrg, 0, ntg, 0, npg)
00064 call alloc_array3d(hxyu, 0, nrg, 0, ntg, 0, npg)
00065 call alloc_array3d(hxzu, 0, nrg, 0, ntg, 0, npg)
00066 call alloc_array3d(hyyu, 0, nrg, 0, ntg, 0, npg)
00067 call alloc_array3d(hyzu, 0, nrg, 0, ntg, 0, npg)
00068 call alloc_array3d(hzzu, 0, nrg, 0, ntg, 0, npg)
00069 call alloc_array3d(gmcrix, 0, nrg, 0, ntg, 0, npg)
00070 call alloc_array3d(gmcriy, 0, nrg, 0, ntg, 0, npg)
00071 call alloc_array3d(gmcriz, 0, nrg, 0, ntg, 0, npg)
00072 call alloc_array3d(gmcrix_grid, 0, nrg, 0, ntg, 0, npg)
00073 call alloc_array3d(gmcriy_grid, 0, nrg, 0, ntg, 0, npg)
00074 call alloc_array3d(gmcriz_grid, 0, nrg, 0, ntg, 0, npg)
00075 call alloc_array4d(rab, 0, nrg, 0, ntg, 0, npg, 1, 6)
00076 call alloc_array4d(rabnl, 0, nrg, 0, ntg, 0, npg, 1, 6)
00077 call alloc_array3d(elpxx, 0, nrg, 0, ntg, 0, npg)
00078 call alloc_array3d(elpxy, 0, nrg, 0, ntg, 0, npg)
00079 call alloc_array3d(elpxz, 0, nrg, 0, ntg, 0, npg)
00080 call alloc_array3d(elpyy, 0, nrg, 0, ntg, 0, npg)
00081 call alloc_array3d(elpyz, 0, nrg, 0, ntg, 0, npg)
00082 call alloc_array3d(elpzz, 0, nrg, 0, ntg, 0, npg)
00083 call alloc_array3d(rlpxx, 0, nrg, 0, ntg, 0, npg)
00084 call alloc_array3d(rlpxy, 0, nrg, 0, ntg, 0, npg)
00085 call alloc_array3d(rlpxz, 0, nrg, 0, ntg, 0, npg)
00086 call alloc_array3d(rlpyy, 0, nrg, 0, ntg, 0, npg)
00087 call alloc_array3d(rlpyz, 0, nrg, 0, ntg, 0, npg)
00088 call alloc_array3d(rlpzz, 0, nrg, 0, ntg, 0, npg)
00089 call alloc_array3d(rlbxx, 0, nrg, 0, ntg, 0, npg)
00090 call alloc_array3d(rlbxy, 0, nrg, 0, ntg, 0, npg)
00091 call alloc_array3d(rlbxz, 0, nrg, 0, ntg, 0, npg)
00092 call alloc_array3d(rlbyy, 0, nrg, 0, ntg, 0, npg)
00093 call alloc_array3d(rlbyz, 0, nrg, 0, ntg, 0, npg)
00094 call alloc_array3d(rlbzz, 0, nrg, 0, ntg, 0, npg)
00095 call alloc_array3d(elpxx_grid, 0, nrg, 0, ntg, 0, npg)
00096 call alloc_array3d(elpxy_grid, 0, nrg, 0, ntg, 0, npg)
00097 call alloc_array3d(elpxz_grid, 0, nrg, 0, ntg, 0, npg)
00098 call alloc_array3d(elpyy_grid, 0, nrg, 0, ntg, 0, npg)
00099 call alloc_array3d(elpyz_grid, 0, nrg, 0, ntg, 0, npg)
00100 call alloc_array3d(elpzz_grid, 0, nrg, 0, ntg, 0, npg)
00101 call alloc_array3d(rlpxx_grid, 0, nrg, 0, ntg, 0, npg)
00102 call alloc_array3d(rlpxy_grid, 0, nrg, 0, ntg, 0, npg)
00103 call alloc_array3d(rlpxz_grid, 0, nrg, 0, ntg, 0, npg)
00104 call alloc_array3d(rlpyy_grid, 0, nrg, 0, ntg, 0, npg)
00105 call alloc_array3d(rlpyz_grid, 0, nrg, 0, ntg, 0, npg)
00106 call alloc_array3d(rlpzz_grid, 0, nrg, 0, ntg, 0, npg)
00107 call alloc_array3d(rlbxx_grid, 0, nrg, 0, ntg, 0, npg)
00108 call alloc_array3d(rlbxy_grid, 0, nrg, 0, ntg, 0, npg)
00109 call alloc_array3d(rlbxz_grid, 0, nrg, 0, ntg, 0, npg)
00110 call alloc_array3d(rlbyy_grid, 0, nrg, 0, ntg, 0, npg)
00111 call alloc_array3d(rlbyz_grid, 0, nrg, 0, ntg, 0, npg)
00112 call alloc_array3d(rlbzz_grid, 0, nrg, 0, ntg, 0, npg)
00113 call alloc_array4d(cdbvxd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00114 call alloc_array4d(cdbvyd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00115 call alloc_array4d(cdbvzd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00116 call alloc_array4d(cdbvxd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00117 call alloc_array4d(cdbvyd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00118 call alloc_array4d(cdbvzd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00119 call alloc_array3d(cdivbv, 0, nrg, 0, ntg, 0, npg)
00120 call alloc_array3d(cdivbv_grid, 0, nrg, 0, ntg, 0, npg)
00121 call alloc_array4d(pdbvxd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00122 call alloc_array4d(pdbvyd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00123 call alloc_array4d(pdbvzd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00124 call alloc_array4d(pdbvxd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00125 call alloc_array4d(pdbvyd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00126 call alloc_array4d(pdbvzd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00127
00128 call allocate_vector_x
00129 call allocate_vector_phi
00130
00131 end subroutine allocate_metric_and_matter_WL
