00001 subroutine allocate_metric_WL
00002 use phys_constant, only : long
00003 use grid_parameter, only : nrg, ntg, npg
00004 use def_metric
00005 use def_gamma_crist
00006 use def_gamma_crist_grid
00007 use def_metric_hij
00008 use def_metric_hij_dirac
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_cristoffel
00017 use def_cristoffel_grid
00018 use def_transverse_part
00019 use make_array_2d
00020 use make_array_3d
00021 use make_array_4d
00022 use make_array_5d
00023 implicit none
00024
00025 call alloc_array3d(psi, 0, nrg, 0, ntg, 0, npg)
00026 call alloc_array3d(alph, 0, nrg, 0, ntg, 0, npg)
00027 call alloc_array3d(alps, 0, nrg, 0, ntg, 0, npg)
00028 call alloc_array3d(alps2, 0, nrg, 0, ntg, 0, npg)
00029 call alloc_array3d(bvxd, 0, nrg, 0, ntg, 0, npg)
00030 call alloc_array3d(bvyd, 0, nrg, 0, ntg, 0, npg)
00031 call alloc_array3d(bvzd, 0, nrg, 0, ntg, 0, npg)
00032 call alloc_array3d(bvxu, 0, nrg, 0, ntg, 0, npg)
00033 call alloc_array3d(bvyu, 0, nrg, 0, ntg, 0, npg)
00034 call alloc_array3d(bvzu, 0, nrg, 0, ntg, 0, npg)
00035 call alloc_array3d(ovxu, 0, nrg, 0, ntg, 0, npg)
00036 call alloc_array3d(ovyu, 0, nrg, 0, ntg, 0, npg)
00037 call alloc_array3d(ovzu, 0, nrg, 0, ntg, 0, npg)
00038 call alloc_array3d(ovxd, 0, nrg, 0, ntg, 0, npg)
00039 call alloc_array3d(ovyd, 0, nrg, 0, ntg, 0, npg)
00040 call alloc_array3d(ovzd, 0, nrg, 0, ntg, 0, npg)
00041 call alloc_array3d(hxxd, 0, nrg, 0, ntg, 0, npg)
00042 call alloc_array3d(hxyd, 0, nrg, 0, ntg, 0, npg)
00043 call alloc_array3d(hxzd, 0, nrg, 0, ntg, 0, npg)
00044 call alloc_array3d(hyyd, 0, nrg, 0, ntg, 0, npg)
00045 call alloc_array3d(hyzd, 0, nrg, 0, ntg, 0, npg)
00046 call alloc_array3d(hzzd, 0, nrg, 0, ntg, 0, npg)
00047 call alloc_array3d(hxxu, 0, nrg, 0, ntg, 0, npg)
00048 call alloc_array3d(hxyu, 0, nrg, 0, ntg, 0, npg)
00049 call alloc_array3d(hxzu, 0, nrg, 0, ntg, 0, npg)
00050 call alloc_array3d(hyyu, 0, nrg, 0, ntg, 0, npg)
00051 call alloc_array3d(hyzu, 0, nrg, 0, ntg, 0, npg)
00052 call alloc_array3d(hzzu, 0, nrg, 0, ntg, 0, npg)
00053
00054 call alloc_array3d(gaugex, 0, nrg, 0, ntg, 0, npg)
00055 call alloc_array3d(gaugey, 0, nrg, 0, ntg, 0, npg)
00056 call alloc_array3d(gaugez, 0, nrg, 0, ntg, 0, npg)
00057
00058 call alloc_array3d(tfkijkij, 1, nrg, 1, ntg, 1, npg)
00059 call alloc_array5d(tfkij, 1, nrg, 1, ntg, 1, npg, 1, 3, 1, 3)
00060 call alloc_array3d(tfkijkij_grid, 0, nrg, 0, ntg, 0, npg)
00061 call alloc_array5d(tfkij_grid, 0, nrg, 0, ntg, 0, npg, 1, 3, 1, 3)
00062 call alloc_array3d(trk, 1, nrg, 1, ntg, 1, npg)
00063 call alloc_array3d(trk_grid, 0, nrg, 0, ntg, 0, npg)
00064
00065 call alloc_array5d(cri, 1, nrg, 1, ntg, 1, npg, 1, 3, 1, 6)
00066 call alloc_array5d(cri_grid, 0, nrg, 0, ntg, 0, npg, 1, 3, 1, 6)
00067 call alloc_array5d(crid, 1, nrg, 1, ntg, 1, npg, 1, 3, 1, 6)
00068 call alloc_array5d(crid_grid, 0, nrg, 0, ntg, 0, npg, 1, 3, 1, 6)
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(dagmabu, 1, nrg, 1, ntg, 1, npg, 1, 3)
00076
00077 call alloc_array3d(Ftvx, 0, nrg, 0, ntg, 0, npg)
00078 call alloc_array3d(Ftvy, 0, nrg, 0, ntg, 0, npg)
00079 call alloc_array3d(Ftvz, 0, nrg, 0, ntg, 0, npg)
00080 call alloc_array3d(Ftvx_grid, 0, nrg, 0, ntg, 0, npg)
00081 call alloc_array3d(Ftvy_grid, 0, nrg, 0, ntg, 0, npg)
00082 call alloc_array3d(Ftvz_grid, 0, nrg, 0, ntg, 0, npg)
00083
00084 call alloc_array4d(rab, 0, nrg, 0, ntg, 0, npg, 1, 6)
00085 call alloc_array4d(rabnl, 0, nrg, 0, ntg, 0, npg, 1, 6)
00086 call alloc_array4d(rabDF, 0, nrg, 0, ntg, 0, npg, 1, 6)
00087 call alloc_array3d(elpxx, 0, nrg, 0, ntg, 0, npg)
00088 call alloc_array3d(elpxy, 0, nrg, 0, ntg, 0, npg)
00089 call alloc_array3d(elpxz, 0, nrg, 0, ntg, 0, npg)
00090 call alloc_array3d(elpyy, 0, nrg, 0, ntg, 0, npg)
00091 call alloc_array3d(elpyz, 0, nrg, 0, ntg, 0, npg)
00092 call alloc_array3d(elpzz, 0, nrg, 0, ntg, 0, npg)
00093 call alloc_array3d(rlpxx, 0, nrg, 0, ntg, 0, npg)
00094 call alloc_array3d(rlpxy, 0, nrg, 0, ntg, 0, npg)
00095 call alloc_array3d(rlpxz, 0, nrg, 0, ntg, 0, npg)
00096 call alloc_array3d(rlpyy, 0, nrg, 0, ntg, 0, npg)
00097 call alloc_array3d(rlpyz, 0, nrg, 0, ntg, 0, npg)
00098 call alloc_array3d(rlpzz, 0, nrg, 0, ntg, 0, npg)
00099 call alloc_array3d(rlbxx, 0, nrg, 0, ntg, 0, npg)
00100 call alloc_array3d(rlbxy, 0, nrg, 0, ntg, 0, npg)
00101 call alloc_array3d(rlbxz, 0, nrg, 0, ntg, 0, npg)
00102 call alloc_array3d(rlbyy, 0, nrg, 0, ntg, 0, npg)
00103 call alloc_array3d(rlbyz, 0, nrg, 0, ntg, 0, npg)
00104 call alloc_array3d(rlbzz, 0, nrg, 0, ntg, 0, npg)
00105 call alloc_array3d(elpxx_grid, 0, nrg, 0, ntg, 0, npg)
00106 call alloc_array3d(elpxy_grid, 0, nrg, 0, ntg, 0, npg)
00107 call alloc_array3d(elpxz_grid, 0, nrg, 0, ntg, 0, npg)
00108 call alloc_array3d(elpyy_grid, 0, nrg, 0, ntg, 0, npg)
00109 call alloc_array3d(elpyz_grid, 0, nrg, 0, ntg, 0, npg)
00110 call alloc_array3d(elpzz_grid, 0, nrg, 0, ntg, 0, npg)
00111 call alloc_array3d(rlpxx_grid, 0, nrg, 0, ntg, 0, npg)
00112 call alloc_array3d(rlpxy_grid, 0, nrg, 0, ntg, 0, npg)
00113 call alloc_array3d(rlpxz_grid, 0, nrg, 0, ntg, 0, npg)
00114 call alloc_array3d(rlpyy_grid, 0, nrg, 0, ntg, 0, npg)
00115 call alloc_array3d(rlpyz_grid, 0, nrg, 0, ntg, 0, npg)
00116 call alloc_array3d(rlpzz_grid, 0, nrg, 0, ntg, 0, npg)
00117 call alloc_array3d(rlbxx_grid, 0, nrg, 0, ntg, 0, npg)
00118 call alloc_array3d(rlbxy_grid, 0, nrg, 0, ntg, 0, npg)
00119 call alloc_array3d(rlbxz_grid, 0, nrg, 0, ntg, 0, npg)
00120 call alloc_array3d(rlbyy_grid, 0, nrg, 0, ntg, 0, npg)
00121 call alloc_array3d(rlbyz_grid, 0, nrg, 0, ntg, 0, npg)
00122 call alloc_array3d(rlbzz_grid, 0, nrg, 0, ntg, 0, npg)
00123 call alloc_array4d(cdbvxd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00124 call alloc_array4d(cdbvyd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00125 call alloc_array4d(cdbvzd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00126 call alloc_array4d(cdbvxd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00127 call alloc_array4d(cdbvyd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00128 call alloc_array4d(cdbvzd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00129 call alloc_array3d(cdivbv, 0, nrg, 0, ntg, 0, npg)
00130 call alloc_array3d(cdivbv_grid, 0, nrg, 0, ntg, 0, npg)
00131 call alloc_array4d(pdbvxd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00132 call alloc_array4d(pdbvyd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00133 call alloc_array4d(pdbvzd_grid, 0, nrg, 0, ntg, 0, npg, 1, 3)
00134 call alloc_array4d(pdbvxd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00135 call alloc_array4d(pdbvyd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00136 call alloc_array4d(pdbvzd, 1, nrg, 1, ntg, 1, npg, 1, 3)
00137
00138 end subroutine allocate_metric_WL
