00001 subroutine iteration(iter_count)
00002 use phys_constant, only : long
00003 use grid_parameter
00004 use coordinate_grav_r
00005 use coordinate_grav_phi
00006 use coordinate_grav_theta
00007 use weight_midpoint_grav
00008 use make_array_2d
00009 use make_array_3d
00010 use make_array_4d
00011 use def_metric
00012 use def_matter
00013 use interface_interpo_fl2gr
00014 use interface_poisson_solver
00015 use interface_sourceterm_HaC
00016 use interface_sourceterm_trG
00017 use interface_sourceterm_MoC
00018 use interface_compute_shift
00019 use interface_compute_alps2alph
00020 use interface_hydrostatic_eq
00021 use interface_calc_surface
00022 use interface_update_grfield
00023 use interface_update_matter
00024 use interface_update_parameter
00025 use interface_update_surface
00026 use interface_error_matter
00027 implicit none
00028 real(long), pointer :: sou(:,:,:), pot(:,:,:), sou2(:,:,:)
00029 real(long), pointer :: potf(:,:,:), emd_bak(:,:,:)
00030 real(long), pointer :: potrs(:,:)
00031 real(long), pointer :: potx(:,:,:), poty(:,:,:), potz(:,:,:)
00032 real(long), pointer :: souvec(:,:,:,:)
00033 real(long) :: error_emd, count
00034 integer :: iter_count, flag = 0, hydro_iter = 4
00035 integer :: irf, itf, ipf, irg, itg, ipg, ihy
00036
00037 call alloc_array3d(sou,0,nrg,0,ntg,0,npg)
00038 call alloc_array3d(sou2,0,nrg,0,ntg,0,npg)
00039 call alloc_array3d(potx,0,nrg,0,ntg,0,npg)
00040 call alloc_array3d(poty,0,nrg,0,ntg,0,npg)
00041 call alloc_array3d(potz,0,nrg,0,ntg,0,npg)
00042 call alloc_array3d(pot,0,nrg,0,ntg,0,npg)
00043 call alloc_array3d(potf,0,nrf,0,ntf,0,npf)
00044 call alloc_array3d(emd_bak,0,nrf,0,ntf,0,npf)
00045 call alloc_array2d(potrs,0,ntf,0,npf)
00046 call alloc_array4d(souvec,0,nrg,0,ntg,0,npg,1,3)
00047
00048 iter_count = 0
00049 do
00050 iter_count = iter_count + 1
00051 count = dble(iter_count)
00052 conv_gra = dmin1(conv0_gra,conv_ini*count)
00053 conv_den = dmin1(conv0_den,conv_ini*count)
00054
00055 emdg = 0.0d0
00056 call interpo_fl2gr(emd, emdg)
00057 call calc_vector_x_grav(1)
00058 call calc_vector_x_matter(1)
00059 call calc_vector_phi_grav(1)
00060 call calc_vector_phi_matter(1)
00061 call excurve
00062
00063 call sourceterm_HaC(sou)
00064 call poisson_solver(sou,pot)
00065 pot = pot + 1.0d0
00066 call update_grfield(pot,psi,conv_gra)
00067 call update_parameter(conv_den)
00068
00069
00070 call sourceterm_trG(sou)
00071 call poisson_solver(sou,pot)
00072 pot = pot + 1.0d0
00073 call compute_alps2alph(pot,psi)
00074 call update_grfield(pot,alph,conv_gra)
00075 call update_parameter(conv_den)
00076
00077 call sourceterm_MoC(souvec,sou)
00078 sou2(1:nrg, 1:ntg, 1:npg) = souvec(1:nrg, 1:ntg, 1:npg, 1)
00079 call poisson_solver(sou2,potx)
00080 sou2(1:nrg, 1:ntg, 1:npg) = souvec(1:nrg, 1:ntg, 1:npg, 2)
00081 call poisson_solver(sou2,poty)
00082 sou2(1:nrg, 1:ntg, 1:npg) = souvec(1:nrg, 1:ntg, 1:npg, 3)
00083 call poisson_solver(sou2,potz)
00084 call poisson_solver(sou,pot)
00085 souvec(0:nrg, 0:ntg, 0:npg, 1) = potx(0:nrg, 0:ntg, 0:npg)
00086 souvec(0:nrg, 0:ntg, 0:npg, 2) = poty(0:nrg, 0:ntg, 0:npg)
00087 souvec(0:nrg, 0:ntg, 0:npg, 3) = potz(0:nrg, 0:ntg, 0:npg)
00088 call compute_shift(potx, poty, potz, souvec, pot)
00089 call update_grfield(potx,bvxd,conv_gra)
00090 call update_grfield(poty,bvyd,conv_gra)
00091 call update_grfield(potz,bvzd,conv_gra)
00092 call update_parameter(conv_den)
00093
00094 do ihy = 1, hydro_iter
00095 call hydrostatic_eq(potf)
00096 call calc_surface(potrs, potf)
00097 emd_bak = emd
00098 call update_matter(potf,emd,conv_den)
00099 call update_surface(potrs,rs,conv_den)
00100 call reset_fluid
00101 call update_parameter(conv_den)
00102 call calc_vector_x_matter(1)
00103 call calc_vector_phi_matter(1)
00104
00105 if (ihy.eq.1) then
00106 call error_matter(emd,emd_bak,error_emd,flag)
00107 call printout_error(iter_count,error_emd)
00108 end if
00109 end do
00110
00111 if (flag == 0) exit
00112 if (iter_count >= iter_max) exit
00113 flag = 0
00114
00115 end do
00116
00117 deallocate(sou)
00118 deallocate(sou2)
00119 deallocate(potx)
00120 deallocate(poty)
00121 deallocate(potz)
00122 deallocate(pot)
00123 deallocate(potf)
00124 deallocate(emd_bak)
00125 deallocate(potrs)
00126 deallocate(souvec)
00127 end subroutine iteration
