Skip to content

Commit 74cf68c

Browse files
hiromatsuihiromatsui
committed
Define a_r_1d_rj_r and a_r_1d_rj_r with radius_1d_rj_r
1 parent f1f2606 commit 74cf68c

File tree

13 files changed

+126
-88
lines changed

13 files changed

+126
-88
lines changed

‎src/Fortran_libraries/MHD_src/sph_MHD/Makefile.depends‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -132,7 +132,7 @@ dup_fields_on_circle_to_IO.o: $(MHD_SPH_DIR)/dup_fields_on_circle_to_IO.f90 m_pr
132132
$(F90) -c $(F90OPTFLAGS) $<
133133
explicit_scalars_sph.o: $(MHD_SPH_DIR)/explicit_scalars_sph.f90 m_precision.o t_spheric_parameter.o t_spheric_rj_data.o t_boundary_data_sph_MHD.o t_physical_property.o t_base_field_labels.o t_base_force_labels.o t_diffusion_term_labels.o t_explicit_term_labels.o t_phys_data.o select_diff_adv_source.o cal_inner_core_rotation.o
134134
$(F90) -c $(F90OPTFLAGS) $<
135-
field_at_mid_equator.o: $(MHD_SPH_DIR)/field_at_mid_equator.f90 m_precision.o m_constants.o m_machine_parameter.o t_field_on_circle.o t_circle_transform.o t_spheric_parameter.o calypso_mpi.o calypso_mpi_real.o t_spheric_rj_data.o t_phys_data.o t_phys_address.o t_sph_circle_parameters.o transfer_to_long_integers.o
135+
field_at_mid_equator.o: $(MHD_SPH_DIR)/field_at_mid_equator.f90 m_precision.o m_constants.o m_machine_parameter.o t_field_on_circle.o t_circle_transform.o t_spheric_parameter.o t_spheric_rj_data.o calypso_mpi.o calypso_mpi_real.o t_phys_data.o t_phys_address.o t_sph_circle_parameters.o transfer_to_long_integers.o
136136
$(F90) -c $(F90OPTFLAGS) $<
137137
fill_scalar_field.o: $(MHD_SPH_DIR)/fill_scalar_field.f90 m_precision.o m_constants.o t_boundary_params_sph_MHD.o
138138
$(F90) -c $(F90OPTFLAGS) $<
@@ -146,7 +146,7 @@ init_external_magne_sph.o: $(MHD_SPH_DIR)/init_external_magne_sph.f90 m_precisio
146146
$(F90) -c $(F90OPTFLAGS) $<
147147
init_radial_infos_sph_mhd.o: $(MHD_SPH_DIR)/init_radial_infos_sph_mhd.f90 m_precision.o calypso_mpi.o m_constants.o m_spheric_constants.o m_machine_parameter.o t_control_parameter.o t_spheric_parameter.o t_spheric_group.o t_poloidal_rotation.o t_radial_reference_field.o t_fdm_coefs.o t_sph_boundary_input_data.o t_bc_data_list.o t_boundary_data_sph_MHD.o t_phys_address.o t_phys_data.o t_work_4_sph_trans.o t_physical_property.o second_fdm_node_coefs.o third_fdm_node_to_ele.o first_fdm_ele_to_node.o material_property.o init_sph_radius_variations.o set_bc_sph_mhd.o set_radius_func_noequi.o sph_mhd_rst_IO_control.o reference_sources_from_d_rj.o init_reference_scalar.o init_external_magne_sph.o radial_reference_field_IO.o m_base_field_labels.o
148148
$(F90) -c $(F90OPTFLAGS) $<
149-
init_sph_radius_variations.o: $(MHD_SPH_DIR)/init_sph_radius_variations.f90 m_precision.o m_constants.o t_control_parameter.o t_spheric_parameter.o t_phys_data.o t_sph_radial_interpolate.o t_field_data_IO.o t_fdm_coefs.o calypso_mpi_real.o transfer_to_long_integers.o set_sph_radial_variations.o m_base_field_labels.o m_diffusion_term_labels.o calypso_mpi.o field_file_IO.o copy_rj_phys_data_4_IO.o set_sph_extensions.o
149+
init_sph_radius_variations.o: $(MHD_SPH_DIR)/init_sph_radius_variations.f90 m_precision.o m_constants.o t_control_parameter.o t_spheric_rj_data.o t_phys_data.o t_sph_radial_interpolate.o t_field_data_IO.o t_fdm_coefs.o calypso_mpi_real.o transfer_to_long_integers.o set_sph_radial_variations.o m_base_field_labels.o m_diffusion_term_labels.o calypso_mpi.o field_file_IO.o copy_rj_phys_data_4_IO.o set_sph_extensions.o
150150
$(F90) -c $(F90OPTFLAGS) $<
151151
init_sphrical_transform_MHD.o: $(MHD_SPH_DIR)/init_sphrical_transform_MHD.f90 m_precision.o m_constants.o m_machine_parameter.o m_work_time.o calypso_mpi.o t_SPH_MHD_model_data.o t_SPH_mesh_field_data.o t_sph_trans_comm_tbl.o t_poloidal_rotation.o t_sph_trans_arrays_MHD.o t_schmidt_poly_on_rtm.o t_work_4_sph_trans.o t_legendre_trans_select.o t_sph_FFT_selector.o t_const_wz_coriolis_rtp.o t_coriolis_terms_rlm.o t_gaunt_coriolis_rlm.o t_boundary_data_sph_MHD.o t_solver_SR.o set_address_sph_trans_MHD.o set_address_sph_trans_snap.o check_sph_mhd_openmp_size.o init_sph_trans.o init_FFT_4_MHD.o pole_sph_transform.o skip_comment_f.o m_legendre_transform_list.o test_legendre_transforms.o
152152
$(F90) -c $(F90OPTFLAGS) $<

‎src/Fortran_libraries/MHD_src/sph_MHD/cal_write_sph_monitor_data.f90‎

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -150,8 +150,9 @@ subroutine init_rms_sph_mhd_control(MHD_prop, sph_MHD_bc, &
150150
if(monitor%bench%iflag_dynamobench .gt. 0) then
151151
call init_circle_field_name_dbench(SPH_MHD%ipol, &
152152
& monitor%circ_mid_eq%d_circle, monitor%bench)
153-
call init_mid_equator_point_global(SPH_MHD%sph, &
154-
& monitor%circ_mid_eq)
153+
call init_mid_equator_point_global &
154+
& (SPH_MHD%sph%sph_params, SPH_MHD%sph%sph_rj, &
155+
& monitor%circ_mid_eq)
155156
call init_circle_point_global &
156157
& (SPH_MHD%sph, SPH_MHD%comms, trans_p, &
157158
& monitor%circ_mid_eq, SR_sig, SR_r)

‎src/Fortran_libraries/MHD_src/sph_MHD/field_at_mid_equator.f90‎

Lines changed: 13 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -7,8 +7,10 @@
77
!>@brief data at mid-depth of the shell at equator for dynamo benchmark
88
!!
99
!!@verbatim
10-
!! subroutine init_mid_equator_point_global(sph, cdat)
11-
!! type(sph_grids), intent(in) :: sph
10+
!! subroutine init_mid_equator_point_global(sph_params, sph_rj, &
11+
!! & cdat)
12+
!! type(sph_shell_parameters), intent(in) :: sph_params
13+
!! type(sph_rj_grid), intent(in) :: sph_rj
1214
!! type(phys_address), intent(in) :: ipol
1315
!! type(circle_fld_maker), intent(inout) :: cdat
1416
!!
@@ -46,21 +48,24 @@ module field_at_mid_equator
4648
!
4749
! ----------------------------------------------------------------------
4850
!
49-
subroutine init_mid_equator_point_global(sph, cdat)
51+
subroutine init_mid_equator_point_global(sph_params, sph_rj, &
52+
& cdat)
5053
!
5154
use t_spheric_parameter
55+
use t_spheric_rj_data
5256
!
53-
type(sph_grids), intent(in) :: sph
57+
type(sph_shell_parameters), intent(in) :: sph_params
58+
type(sph_rj_grid), intent(in) :: sph_rj
5459
type(circle_fld_maker), intent(inout) :: cdat
5560
!
5661
integer(kind = kint) :: kr_ICB, kr_CMB
5762
real(kind = kreal) :: r_MID
5863
!
5964
!
60-
kr_ICB = sph%sph_params%nlayer_ICB
61-
kr_CMB = sph%sph_params%nlayer_CMB
62-
r_MID = half * (sph%sph_rj%radius_1d_rj_r(kr_ICB) &
63-
& + sph%sph_rj%radius_1d_rj_r(kr_CMB))
65+
kr_ICB = sph_params%nlayer_ICB
66+
kr_CMB = sph_params%nlayer_CMB
67+
r_MID = half * (sph_rj%radius_1d_rj_r(kr_ICB) &
68+
& + sph_rj%radius_1d_rj_r(kr_CMB))
6469
!
6570
cdat%circle%s_circle = r_MID
6671
cdat%circle%z_circle = zero

‎src/Fortran_libraries/MHD_src/sph_MHD/init_radial_infos_sph_mhd.f90‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -122,7 +122,7 @@ subroutine init_r_infos_sph_mhd_evo(bc_IO, sph_grps, MHD_BC, &
122122
call const_third_fdm_node_to_ele(sph%sph_rj, r_n2e_3rd)
123123
!
124124
!* ---------- Radial variations of density and diffusivities -------
125-
call init_radius_variations_sph_mhd(sph, r_2nd, MHD_prop, &
125+
call init_radius_variations_sph_mhd(sph%sph_rj, r_2nd, MHD_prop, &
126126
& radial_variation)
127127
!
128128
!* ---------- boundary conditions ---------------

‎src/Fortran_libraries/MHD_src/sph_MHD/init_sph_radius_variations.f90‎

Lines changed: 15 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -8,9 +8,9 @@
88
!>@brief Set radial variations for densityuy and diffusivities
99
!!
1010
!!@verbatim
11-
!! subroutine init_radius_variations_sph_mhd(sph, r_2nd, MHD_prop, &
12-
!! & radial_variation)
13-
!! type(sph_grids), intent(in) :: sph
11+
!! subroutine init_radius_variations_sph_mhd &
12+
!! & (sph_rj, r_2nd, MHD_prop, radial_variation)
13+
!! type(sph_rj_grid), intent(in) :: sph_rj
1414
!! type(fdm_matrices), intent(in) :: r_2nd
1515
!! type(MHD_evolution_param), intent(inout) :: MHD_prop
1616
!! type(phys_data), intent(inout) :: radial_variation
@@ -22,7 +22,7 @@ module init_sph_radius_variations
2222
use m_precision
2323
use m_constants
2424
use t_control_parameter
25-
use t_spheric_parameter
25+
use t_spheric_rj_data
2626
use t_phys_data
2727
use t_sph_radial_interpolate
2828
use t_field_data_IO
@@ -45,16 +45,16 @@ module init_sph_radius_variations
4545
!
4646
! -------------------------------------------------------------------
4747
!
48-
subroutine init_radius_variations_sph_mhd(sph, r_2nd, MHD_prop, &
49-
& radial_variation)
48+
subroutine init_radius_variations_sph_mhd &
49+
& (sph_rj, r_2nd, MHD_prop, radial_variation)
5050
!
5151
use calypso_mpi_real
5252
use transfer_to_long_integers
5353
use set_sph_radial_variations
5454
use m_base_field_labels
5555
use m_diffusion_term_labels
5656
!
57-
type(sph_grids), intent(in) :: sph
57+
type(sph_rj_grid), intent(in) :: sph_rj
5858
type(fdm_matrices), intent(in) :: r_2nd
5959

6060
type(MHD_evolution_param), intent(inout) :: MHD_prop
@@ -71,7 +71,7 @@ subroutine init_radius_variations_sph_mhd(sph, r_2nd, MHD_prop, &
7171
radial_variation%num_phys_viz = radial_variation%num_phys
7272
radial_variation%ntot_phys = radial_variation%num_phys
7373
call alloc_phys_name(radial_variation)
74-
call alloc_phys_data((sph%sph_rj%nidx_rj(1)+1), radial_variation)
74+
call alloc_phys_data((sph_rj%nidx_rj(1)+1), radial_variation)
7575
!
7676
call set_r_variation_data_address(MHD_prop, radial_variation)
7777

@@ -81,30 +81,30 @@ subroutine init_radius_variations_sph_mhd(sph, r_2nd, MHD_prop, &
8181
!
8282
if(my_rank .eq. 0) then
8383
radial_variation%d_fld(1,1) = 0.0d0
84-
do k = 1, sph%sph_rj%nidx_rj(1)
85-
radial_variation%d_fld(k+1,1) = sph%sph_rj%radius_1d_rj_r(k)
84+
do k = 1, sph_rj%nidx_rj(1)
85+
radial_variation%d_fld(k+1,1) = sph_rj%radius_1d_rj_r(k)
8686
end do
8787
!
8888
call set_sph_radial_density(my_rank, radius_name, &
8989
& density%name, MHD_prop%fl_prop%ir_rho, &
90-
& sph%sph_rj, r_2nd, MHD_prop%polytrope_param, &
90+
& sph_rj, r_2nd, MHD_prop%polytrope_param, &
9191
& radial_variation, r_itp, fld_IO)
9292
call set_sph_radial_diffusivity(my_rank, radius_name, &
9393
& kinetic_viscosity%name, MHD_prop%fl_prop%ir_nu, &
94-
& sph%sph_rj, r_2nd, MHD_prop%val_viscous_param, &
94+
& sph_rj, r_2nd, MHD_prop%val_viscous_param, &
9595
& radial_variation, r_itp, fld_IO)
9696
!
9797
call set_sph_radial_diffusivity(my_rank, radius_name, &
9898
& magnetic_diffusivity%name, MHD_prop%cd_prop%ir_eta, &
99-
& sph%sph_rj, r_2nd, MHD_prop%val_mag_diffuse_param, &
99+
& sph_rj, r_2nd, MHD_prop%val_mag_diffuse_param, &
100100
& radial_variation, r_itp, fld_IO)
101101
call set_sph_radial_diffusivity(my_rank, radius_name, &
102102
& thermal_diffusivity%name, MHD_prop%ht_prop%ir_kappa, &
103-
& sph%sph_rj, r_2nd, MHD_prop%val_thermal_diffuse_param, &
103+
& sph_rj, r_2nd, MHD_prop%val_thermal_diffuse_param, &
104104
& radial_variation, r_itp, fld_IO)
105105
call set_sph_radial_diffusivity(my_rank, radius_name, &
106106
& chemical_diffusivity%name, MHD_prop%cp_prop%ir_kappa, &
107-
& sph%sph_rj, r_2nd, MHD_prop%val_comp_diffuse_param, &
107+
& sph_rj, r_2nd, MHD_prop%val_comp_diffuse_param, &
108108
& radial_variation, r_itp, fld_IO)
109109
end if
110110
!

‎src/Fortran_libraries/PARALLEL_src/COMM_src/t_element_double_number.f90‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -116,14 +116,14 @@ subroutine find_belonged_pe_4_ele &
116116
iele_dbl%k_ref(1:numele) = 0
117117
!$omp end parallel workshare
118118
!
119-
!%omp parallel do private(iele,ie_one)
119+
!$omp parallel do private(iele,ie_one)
120120
do iele = 1, numele
121121
ie_one = ie(iele,1)
122122
call find_belonged_pe_each_ele &
123123
& (inod_dbl%num_dbl, inod_dbl%irank, ie_one, &
124124
& iele_dbl%irank(iele), iele_dbl%k_ref(iele))
125125
end do
126-
!%omp end parallel do
126+
!$omp end parallel do
127127
!
128128
end subroutine find_belonged_pe_4_ele
129129
!
@@ -154,7 +154,7 @@ subroutine find_belonged_pe_4_surf(my_rank, inod_dbl, &
154154
isurf_dbl%k_ref(1:numsurf) = 0
155155
!$omp end parallel workshare
156156
!
157-
!%omp parallel do private(isurf,ie_surf_one,nnod_same)
157+
!$omp parallel do private(isurf,ie_surf_one,nnod_same)
158158
do isurf = 1, numsurf
159159
ie_surf_one(1:num_linear_sf) = ie_surf(isurf,1:num_linear_sf)
160160
call find_belonged_pe_each_surf &
@@ -164,14 +164,14 @@ subroutine find_belonged_pe_4_surf(my_rank, inod_dbl, &
164164
interior_surf(isurf) &
165165
& = set_each_interior_flag(my_rank, isurf_dbl%irank(isurf))
166166
end do
167-
!%omp end parallel do
167+
!$omp end parallel do
168168
!
169169
icou = 0
170-
!%omp parallel do private(isurf) reduction(+:icou)
170+
!$omp parallel do private(isurf) reduction(+:icou)
171171
do isurf = 1, numsurf
172172
if(isurf_dbl%irank(isurf) .eq. my_rank) icou = icou + 1
173173
end do
174-
!%omp end parallel do
174+
!$omp end parallel do
175175
internal_surf = icou
176176
!
177177
end subroutine find_belonged_pe_4_surf
@@ -202,7 +202,7 @@ subroutine find_belonged_pe_4_edge(my_rank, inod_dbl, &
202202
iedge_dbl%k_ref(1:numedge) = 0
203203
!$omp end parallel workshare
204204
!
205-
!%omp parallel do private(iedge,ie_edge_one,nnod_same)
205+
!$omp parallel do private(iedge,ie_edge_one,nnod_same)
206206
do iedge = 1, numedge
207207
ie_edge_one(1) = ie_edge(iedge,1)
208208
ie_edge_one(2) = ie_edge(iedge,2)
@@ -213,14 +213,14 @@ subroutine find_belonged_pe_4_edge(my_rank, inod_dbl, &
213213
interior_edge(iedge) &
214214
& = set_each_interior_flag(my_rank, iedge_dbl%irank(iedge))
215215
end do
216-
!%omp end parallel do
216+
!$omp end parallel do
217217
!
218218
icou = 0
219-
!%omp parallel do private(iedge) reduction(+:icou)
219+
!$omp parallel do private(iedge) reduction(+:icou)
220220
do iedge = 1, numedge
221221
if(iedge_dbl%irank(iedge) .eq. my_rank) icou = icou + 1
222222
end do
223-
!%omp end parallel do
223+
!$omp end parallel do
224224
internal_edge = icou
225225
!
226226
end subroutine find_belonged_pe_4_edge

‎src/Fortran_libraries/PARALLEL_src/CONST_SPH_GRID/copy_sph_1d_global_index.f90‎

Lines changed: 16 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -187,18 +187,29 @@ subroutine copy_sph_1d_gl_idx_rj &
187187
! & write(*,*) 'Baka!', my_rank, i, sph_rj%idx_gl_1d_rj_r(i)
188188
! end do
189189
! end if
190-
191-
do i = 1, sph_rj%nidx_rj(1)
192-
j = sph_rj%idx_gl_1d_rj_r(i)
193-
sph_rj%radius_1d_rj_r(i) = s3d_radius%radius_1d_gl(j)
194-
end do
195190
!
196191
do i = 1, sph_rj%nidx_rj(2)
197192
j = i - 1 + sph_rj%ist_rj(2)
198193
sph_rj%idx_gl_1d_rj_j(i,1) = sph_gl1d%idx_global_rj_j(j,1)
199194
sph_rj%idx_gl_1d_rj_j(i,2) = sph_gl1d%idx_global_rj_j(j,2)
200195
sph_rj%idx_gl_1d_rj_j(i,3) = sph_gl1d%idx_global_rj_j(j,3)
201196
end do
197+
!
198+
do i = 1, sph_rj%nidx_rj(1)
199+
j = sph_rj%idx_gl_1d_rj_r(i)
200+
sph_rj%radius_1d_rj_r(i) = s3d_radius%radius_1d_gl(j)
201+
end do
202+
!
203+
!$omp parallel do private(i)
204+
do i = 1, sph_rj%nidx_rj(1)
205+
sph_rj%a_r_1d_rj_r(i) = one / sph_rj%radius_1d_rj_r(i)
206+
sph_rj%ar_1d_rj(i,1) = sph_rj%a_r_1d_rj_r(i)
207+
sph_rj%ar_1d_rj(i,2) = sph_rj%ar_1d_rj(i,1) &
208+
& * sph_rj%a_r_1d_rj_r(i)
209+
sph_rj%ar_1d_rj(i,3) = sph_rj%ar_1d_rj(i,2) &
210+
& * sph_rj%a_r_1d_rj_r(i)
211+
end do
212+
!$omp end parallel do
202213
!
203214
end subroutine copy_sph_1d_gl_idx_rj
204215
!

‎src/Fortran_libraries/UTILS_src/MERGE/Makefile.depends‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -10,15 +10,15 @@ new_SPH_restart.o: $(MERGE_DIR)/new_SPH_restart.f90 m_precision.o calypso_mpi.o
1010
$(F90) -c $(F90OPTFLAGS) $<
1111
parallel_assemble_sph.o: $(MERGE_DIR)/parallel_assemble_sph.f90 m_precision.o m_machine_parameter.o m_constants.o t_spheric_parameter.o t_phys_data.o t_spheric_rj_data.o t_sph_radial_interpolate.o
1212
$(F90) -c $(F90OPTFLAGS) $<
13-
r_interpolate_marged_sph.o: $(MERGE_DIR)/r_interpolate_marged_sph.f90 m_precision.o m_constants.o calypso_mpi.o t_spheric_parameter.o t_spheric_group.o t_sph_radial_interpolate.o calypso_mpi_int.o skip_comment_f.o radial_interpolation.o
13+
r_interpolate_marged_sph.o: $(MERGE_DIR)/r_interpolate_marged_sph.f90 m_precision.o m_constants.o calypso_mpi.o t_spheric_parameter.o t_spheric_group.o t_sph_radial_interpolate.o calypso_mpi_int.o skip_comment_f.o t_spheric_rj_data.o radial_interpolation.o
1414
$(F90) -c $(F90OPTFLAGS) $<
1515
set_control_newsph.o: $(MERGE_DIR)/set_control_newsph.f90 m_precision.o m_constants.o calypso_mpi.o t_file_IO_parameter.o t_control_param_assemble.o t_spectr_data_4_assemble.o calypso_mpi_real.o calypso_mpi_int.o bcast_file_IO_parameter.o m_error_IDs.o t_control_data_4_merge.o t_SPH_mesh_field_data.o t_ctl_params_gen_sph_shell.o m_file_format_switch.o set_control_platform_item.o set_control_platform_data.o set_ctl_4_shell_grids.o new_SPH_restart.o skip_comment_f.o
1616
$(F90) -c $(F90OPTFLAGS) $<
1717
set_merged_restart_data.o: $(MERGE_DIR)/set_merged_restart_data.f90 m_precision.o t_phys_data.o m_base_field_labels.o m_explicit_term_labels.o
1818
$(F90) -c $(F90OPTFLAGS) $<
1919
share_field_data.o: $(MERGE_DIR)/share_field_data.f90 m_precision.o m_constants.o calypso_mpi.o calypso_mpi_real.o calypso_mpi_int.o t_time_data.o t_phys_data.o calypso_mpi_logical.o transfer_to_long_integers.o cal_minmax_and_stacks.o t_field_data_IO.o calypso_mpi_char.o
2020
$(F90) -c $(F90OPTFLAGS) $<
21-
share_spectr_index_data.o: $(MERGE_DIR)/share_spectr_index_data.f90 m_precision.o m_constants.o calypso_mpi.o t_spheric_parameter.o calypso_mpi_real.o calypso_mpi_int.o transfer_to_long_integers.o
21+
share_spectr_index_data.o: $(MERGE_DIR)/share_spectr_index_data.f90 m_precision.o m_constants.o calypso_mpi.o t_spheric_rj_data.o calypso_mpi_real.o calypso_mpi_int.o transfer_to_long_integers.o
2222
$(F90) -c $(F90OPTFLAGS) $<
2323
t_control_data_4_merge.o: $(MERGE_DIR)/t_control_data_4_merge.f90 m_precision.o m_machine_parameter.o t_read_control_elements.o t_ctl_data_4_platforms.o t_ctl_data_4_fields.o t_ctl_data_4_time_steps.o t_control_array_real.o t_ctl_data_gen_sph_shell.o skip_comment_f.o ctl_data_platforms_IO.o ctl_file_gen_sph_shell_IO.o ctl_data_4_time_steps_IO.o
2424
$(F90) -c $(F90OPTFLAGS) $<

‎src/Fortran_libraries/UTILS_src/MERGE/assemble_sph_fields.f90‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -61,7 +61,7 @@ subroutine share_org_sph_rj_data(org_sph_array)
6161
!
6262
!
6363
do ip = 1, org_sph_array%num_pe
64-
call share_sph_rj_data(ip, org_sph_array%sph(ip))
64+
call share_sph_rj_data(ip, org_sph_array%sph(ip)%sph_rj)
6565
end do
6666
!
6767
end subroutine share_org_sph_rj_data

0 commit comments

Comments
 (0)