[cig-commits] r17016 - in seismo/3D/SPECFEM3D/trunk/meshfem3D: . DATA
pieyre at geodynamics.org
pieyre at geodynamics.org
Tue Jun 22 06:32:42 PDT 2010
Author: pieyre
Date: 2010-06-22 06:32:42 -0700 (Tue, 22 Jun 2010)
New Revision: 17016
Modified:
seismo/3D/SPECFEM3D/trunk/meshfem3D/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/meshfem3D/check_mesh_quality.f90
seismo/3D/SPECFEM3D/trunk/meshfem3D/create_regions_mesh.f90
seismo/3D/SPECFEM3D/trunk/meshfem3D/meshfem3D.f90
seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90
Log:
suppressed Delta_t parameter in meshfem3D/DATA/Par_file used to check stability during mesh generation. Instead the program suggests a good time step in output.
Modified: seismo/3D/SPECFEM3D/trunk/meshfem3D/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/meshfem3D/DATA/Par_file 2010-06-22 09:08:08 UTC (rev 17015)
+++ seismo/3D/SPECFEM3D/trunk/meshfem3D/DATA/Par_file 2010-06-22 13:32:42 UTC (rev 17016)
@@ -33,9 +33,6 @@
CREATE_ABAQUS_FILES = .false.
CREATE_DX_FILES = .false.
-# anticipated time step for simulation (in order to check the stability condition)
-DT = 0.01
-
# path to store the databases files
LOCAL_PATH = OUTPUT_FILES
Modified: seismo/3D/SPECFEM3D/trunk/meshfem3D/check_mesh_quality.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/meshfem3D/check_mesh_quality.f90 2010-06-22 09:08:08 UTC (rev 17015)
+++ seismo/3D/SPECFEM3D/trunk/meshfem3D/check_mesh_quality.f90 2010-06-22 13:32:42 UTC (rev 17016)
@@ -36,7 +36,7 @@
-subroutine check_mesh_quality(myrank,VP_MAX,delta_t,NPOIN,NSPEC,x,y,z,ibool)
+subroutine check_mesh_quality(myrank,VP_MAX,NPOIN,NSPEC,x,y,z,ibool)
implicit none
@@ -44,7 +44,6 @@
integer :: NPOIN ! number of nodes
integer :: NSPEC
- double precision :: delta_t ! arbitrary, initial guess
double precision :: VP_MAX ! maximum vp in volume block id 3
!------------------------------------------------------------------------------------------------
@@ -71,6 +70,7 @@
double precision :: distmin_MPI,distmax_MPI
! for stability
+ double precision :: dt_suggested,dt_suggested_max,dt_suggested_max_MPI
double precision :: stability,stability_min,stability_max,max_CFL_stability_limit
double precision :: stability_MPI,stability_min_MPI,stability_max_MPI,max_CFL_stability_limit_MPI
@@ -113,6 +113,7 @@
diagonal_aspect_ratio_max = - HUGEVAL
stability_max = - HUGEVAL
distance_max = - HUGEVAL
+ dt_suggested_max = HUGEVAL
ispec_min_edge_length = -1
ispec_max_edge_length = -1
@@ -120,7 +121,7 @@
! loop on all the elements
do ispec = 1,NSPEC
- call create_mesh_quality_data_3D(x,y,z,ibool,ispec,NSPEC,NPOIN,VP_MAX,delta_t, &
+ call create_mesh_quality_data_3D(x,y,z,ibool,ispec,NSPEC,NPOIN,VP_MAX,dt_suggested, &
equiangle_skewness,edge_aspect_ratio,diagonal_aspect_ratio,stability,distmin,distmax)
! store element number in which the edge of minimum or maximum length is located
@@ -139,6 +140,7 @@
edge_aspect_ratio_max = max(edge_aspect_ratio_max,edge_aspect_ratio)
diagonal_aspect_ratio_max = max(diagonal_aspect_ratio_max,diagonal_aspect_ratio)
stability_max = max(stability_max,stability)
+ dt_suggested_max = min(dt_suggested_max,dt_suggested)
distance_max = max(distance_max,distmax)
enddo
@@ -155,7 +157,9 @@
skewness_max_rank = int(buf_maxloc_recv(2))
- call max_all_dp(stability_max,stability_max_MPI)
+ call max_all_dp(stability,stability_max_MPI)
+ call min_all_dp(dt_suggested_max,dt_suggested_max_MPI)
+
call max_all_dp(edge_aspect_ratio_max,edge_aspect_ratio_max_MPI)
call max_all_dp(diagonal_aspect_ratio_max,diagonal_aspect_ratio_max_MPI)
@@ -192,8 +196,10 @@
write(IMAIN,*)
write(IMAIN,*) 'maximum length of an edge in the whole mesh (m) = ',distance_max_MPI!,' in element ',ispec_max_edge_length
write(IMAIN,*)
- write(IMAIN,*) 'max equiangle skewness = ',equiangle_skewness_max_MPI,' in element ',ispec_max_skewness_MPI, &
+ write(IMAIN,*) '***'
+ write(IMAIN,*) '*** max equiangle skewness = ',equiangle_skewness_max_MPI,' in element ',ispec_max_skewness_MPI, &
' of slice ',skewness_max_rank
+ write(IMAIN,*) '***'
! write(IMAIN,*) 'min equiangle skewness = ',equiangle_skewness_min
write(IMAIN,*)
write(IMAIN,*) 'max deviation angle from a right angle (90 degrees) is therefore = ',90.*equiangle_skewness_max_MPI
@@ -207,8 +213,13 @@
write(IMAIN,*) 'max diagonal aspect ratio = ',diagonal_aspect_ratio_max_MPI
! write(IMAIN,*) 'min diagonal aspect ratio = ',diagonal_aspect_ratio_min
write(IMAIN,*)
- write(IMAIN,*) 'max stability = ',stability_max_MPI
- write(IMAIN,*) 'computed using VP_MAX = ',VP_MAX
+ !write(IMAIN,*) 'max stability = ',stability_max_MPI
+ write(IMAIN,*) '***'
+ write(IMAIN,'(a50,f13.8)') ' *** Maximum suggested time step for simulation = ',dt_suggested_max_MPI
+ write(IMAIN,*) '***'
+ write(IMAIN,*) '*** max stability = ',stability_max_MPI
+ write(IMAIN,*) '*** computed using VP_MAX = ',VP_MAX
+ write(IMAIN,*) '***'
! write(IMAIN,*) 'min stability = ',stability_min
! max stability CFL value is different in 2D and in 3D
@@ -247,7 +258,7 @@
! loop on all the elements
do ispec = 1,NSPEC
- call create_mesh_quality_data_3D(x,y,z,ibool,ispec,NSPEC,NPOIN,VP_MAX,delta_t, &
+ call create_mesh_quality_data_3D(x,y,z,ibool,ispec,NSPEC,NPOIN,VP_MAX,dt_suggested, &
equiangle_skewness,edge_aspect_ratio,diagonal_aspect_ratio,stability,distmin,distmax)
! store skewness in histogram
@@ -324,7 +335,7 @@
! create mesh quality data for a given 3D spectral element
-subroutine create_mesh_quality_data_3D(x,y,z,ibool,ispec,NSPEC,NPOIN,VP_MAX,delta_t, &
+subroutine create_mesh_quality_data_3D(x,y,z,ibool,ispec,NSPEC,NPOIN,VP_MAX,dt_suggested, &
equiangle_skewness,edge_aspect_ratio,diagonal_aspect_ratio,stability,distmin,distmax)
implicit none
@@ -348,7 +359,7 @@
double precision equiangle_skewness,edge_aspect_ratio,diagonal_aspect_ratio
! for stability
- double precision :: stability,VP_MAX,delta_t
+ double precision :: stability,VP_MAX,dt_suggested
! maximum polynomial degree for which we can compute the stability condition
integer, parameter :: NGLL_MAX_STABILITY = 15
@@ -470,8 +481,11 @@
! compute edge aspect ratio
edge_aspect_ratio = distmax / distmin
- stability = delta_t * VP_MAX / (distmin * percent_GLL(true_NGLLX))
+ !stability = delta_t * VP_MAX / (distmin * percent_GLL(true_NGLLX))
+ dt_suggested = ((1.d0 - 0.02d0)*0.48d0) * (distmin * percent_GLL(true_NGLLX)) / VP_MAX
+ stability = dt_suggested * VP_MAX / (distmin * percent_GLL(true_NGLLX))
+
! compute diagonal aspect ratio
dist1 = sqrt((xelm(5) - xelm(4))**2 + (yelm(5) - yelm(4))**2 + (zelm(5) - zelm(4))**2)
dist2 = sqrt((xelm(1) - xelm(8))**2 + (yelm(1) - yelm(8))**2 + (zelm(1) - zelm(8))**2)
Modified: seismo/3D/SPECFEM3D/trunk/meshfem3D/create_regions_mesh.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/meshfem3D/create_regions_mesh.f90 2010-06-22 09:08:08 UTC (rev 17015)
+++ seismo/3D/SPECFEM3D/trunk/meshfem3D/create_regions_mesh.f90 2010-06-22 13:32:42 UTC (rev 17016)
@@ -39,7 +39,7 @@
nsubregions,subregions,nblayers,ner_layer,NMATERIALS,material_properties, &
myrank,LOCAL_PATH,UTM_X_MIN,UTM_X_MAX,UTM_Y_MIN,UTM_Y_MAX,Z_DEPTH_BLOCK, &
CREATE_ABAQUS_FILES,CREATE_DX_FILES,&
- USE_REGULAR_MESH,NDOUBLINGS,ner_doublings,delta_t)
+ USE_REGULAR_MESH,NDOUBLINGS,ner_doublings)
! create the different regions of the mesh
@@ -66,9 +66,6 @@
integer NDOUBLINGS
integer, dimension(2) :: ner_doublings
-
- double precision delta_t
-
double precision UTM_X_MIN,UTM_X_MAX,UTM_Y_MIN,UTM_Y_MAX,Z_DEPTH_BLOCK
double precision horiz_size,vert_size
@@ -197,147 +194,6 @@
allocate(yp(npointot))
allocate(zp(npointot))
-
-
- !--- apply heuristic rule to modify doubling regions to balance angles
-
-
- ! if(APPLY_HEURISTIC_RULE .and. .not. USE_REGULAR_MESH) then
-
- ! stop 'pas encore implemente'
-
- ! ! define number of subregions affected by heuristic rule in doubling regions
- ! nsubregions = 8
-
- ! do isubregion = 1,nsubregions
-
- ! ! define shape of elements for heuristic
- ! call define_subregions_heuristic(myrank,isubregion,iaddx,iaddy,iaddz, &
- ! nblayers,ner_layer, &
- ! ix1,ix2,dix,iy1,iy2,diy,ir1,ir2,dir,iax,iay,iar, &
- ! itype_element,npx,npy)
- ! ! NER_BOTTOM_MOHO,NER_MOHO_16,NER_16_BASEMENT,NER_BASEMENT_SEDIM)
-
- ! ! loop on all the mesh points in current subregion
- ! do ir = ir1,ir2,dir
- ! do iy = iy1,iy2,diy
- ! do ix = ix1,ix2,dix
-
- ! ! this heuristic rule is only valid for 8-node elements
- ! ! it would not work in the case of 27 nodes
-
- ! !----
- ! if(itype_element == ITYPE_UNUSUAL_1) then
-
- ! ! side 1
- ! horiz_size = xgrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2)) &
- ! - xgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! xgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) = &
- ! xgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + horiz_size * MAGIC_RATIO
-
- ! vert_size = zgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) &
- ! - zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! zgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) = &
- ! zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + vert_size * MAGIC_RATIO / 0.50
-
- ! ! side 2
- ! horiz_size = xgrid(ir+iar*iaddz(3),ix+iax*iaddx(3),iy+iay*iaddy(3)) &
- ! - xgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4))
- ! xgrid(ir+iar*iaddz(8),ix+iax*iaddx(8),iy+iay*iaddy(8)) = &
- ! xgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4)) + horiz_size * MAGIC_RATIO
-
- ! vert_size = zgrid(ir+iar*iaddz(8),ix+iax*iaddx(8),iy+iay*iaddy(8)) &
- ! - zgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4))
- ! zgrid(ir+iar*iaddz(8),ix+iax*iaddx(8),iy+iay*iaddy(8)) = &
- ! zgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4)) + vert_size * MAGIC_RATIO / 0.50
-
- ! !----
- ! else if(itype_element == ITYPE_UNUSUAL_1p) then
-
- ! ! side 1
- ! horiz_size = xgrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2)) &
- ! - xgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! xgrid(ir+iar*iaddz(6),ix+iax*iaddx(6),iy+iay*iaddy(6)) = &
- ! xgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + horiz_size * (1. - MAGIC_RATIO)
-
- ! vert_size = zgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) &
- ! - zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! zgrid(ir+iar*iaddz(6),ix+iax*iaddx(6),iy+iay*iaddy(6)) = &
- ! zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + vert_size * MAGIC_RATIO / 0.50
-
- ! ! side 2
- ! horiz_size = xgrid(ir+iar*iaddz(3),ix+iax*iaddx(3),iy+iay*iaddy(3)) &
- ! - xgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4))
- ! xgrid(ir+iar*iaddz(7),ix+iax*iaddx(7),iy+iay*iaddy(7)) = &
- ! xgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4)) + horiz_size * (1. - MAGIC_RATIO)
-
- ! vert_size = zgrid(ir+iar*iaddz(8),ix+iax*iaddx(8),iy+iay*iaddy(8)) &
- ! - zgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4))
- ! zgrid(ir+iar*iaddz(7),ix+iax*iaddx(7),iy+iay*iaddy(7)) = &
- ! zgrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4)) + vert_size * MAGIC_RATIO / 0.50
-
- ! !----
- ! else if(itype_element == ITYPE_UNUSUAL_4) then
-
- ! ! side 1
- ! horiz_size = ygrid(ir+iar*iaddz(3),ix+iax*iaddx(3),iy+iay*iaddy(3)) &
- ! - ygrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2))
- ! ygrid(ir+iar*iaddz(7),ix+iax*iaddx(7),iy+iay*iaddy(7)) = &
- ! ygrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2)) + horiz_size * (1. - MAGIC_RATIO)
-
- ! vert_size = zgrid(ir+iar*iaddz(6),ix+iax*iaddx(6),iy+iay*iaddy(6)) &
- ! - zgrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2))
- ! zgrid(ir+iar*iaddz(7),ix+iax*iaddx(7),iy+iay*iaddy(7)) = &
- ! zgrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2)) + vert_size * MAGIC_RATIO / 0.50
-
- ! ! side 2
- ! horiz_size = ygrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4)) &
- ! - ygrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! ygrid(ir+iar*iaddz(8),ix+iax*iaddx(8),iy+iay*iaddy(8)) = &
- ! ygrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + horiz_size * (1. - MAGIC_RATIO)
-
- ! vert_size = zgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) &
- ! - zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! zgrid(ir+iar*iaddz(8),ix+iax*iaddx(8),iy+iay*iaddy(8)) = &
- ! zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + vert_size * MAGIC_RATIO / 0.50
-
- ! !----
- ! else if(itype_element == ITYPE_UNUSUAL_4p) then
-
- ! ! side 1
- ! horiz_size = ygrid(ir+iar*iaddz(3),ix+iax*iaddx(3),iy+iay*iaddy(3)) &
- ! - ygrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2))
- ! ygrid(ir+iar*iaddz(6),ix+iax*iaddx(6),iy+iay*iaddy(6)) = &
- ! ygrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2)) + horiz_size * MAGIC_RATIO
-
- ! vert_size = zgrid(ir+iar*iaddz(6),ix+iax*iaddx(6),iy+iay*iaddy(6)) &
- ! - zgrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2))
- ! zgrid(ir+iar*iaddz(6),ix+iax*iaddx(6),iy+iay*iaddy(6)) = &
- ! zgrid(ir+iar*iaddz(2),ix+iax*iaddx(2),iy+iay*iaddy(2)) + vert_size * MAGIC_RATIO / 0.50
-
- ! ! side 2
- ! horiz_size = ygrid(ir+iar*iaddz(4),ix+iax*iaddx(4),iy+iay*iaddy(4)) &
- ! - ygrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! ygrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) = &
- ! ygrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + horiz_size * MAGIC_RATIO
-
- ! vert_size = zgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) &
- ! - zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1))
- ! zgrid(ir+iar*iaddz(5),ix+iax*iaddx(5),iy+iay*iaddy(5)) = &
- ! zgrid(ir+iar*iaddz(1),ix+iax*iaddx(1),iy+iay*iaddy(1)) + vert_size * MAGIC_RATIO / 0.50
-
- ! endif
-
- ! enddo
- ! enddo
- ! enddo
-
- ! enddo
-
- ! endif
-
- !---
-
! generate the elements in all the regions of the mesh
ispec = 0
@@ -528,7 +384,7 @@
NSPEC2DMAX_XMIN_XMAX,NSPEC2DMAX_YMIN_YMAX)
VP_MAX = maxval(material_properties(:,2))
- call check_mesh_quality(myrank,VP_MAX,delta_t,nglob,nspec,nodes_coords(:,1),nodes_coords(:,2),nodes_coords(:,3),ibool)
+ call check_mesh_quality(myrank,VP_MAX,nglob,nspec,nodes_coords(:,1),nodes_coords(:,2),nodes_coords(:,3),ibool)
call save_databases(prname,nspec,nglob,iproc_xi,iproc_eta,NPROC_XI,NPROC_ETA,addressing,iMPIcut_xi,iMPIcut_eta,&
ibool,nodes_coords,true_material_num,nspec2D_xmin,nspec2D_xmax,nspec2D_ymin,nspec2D_ymax,NSPEC2D_BOTTOM,NSPEC2D_TOP,&
Modified: seismo/3D/SPECFEM3D/trunk/meshfem3D/meshfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/meshfem3D/meshfem3D.f90 2010-06-22 09:08:08 UTC (rev 17015)
+++ seismo/3D/SPECFEM3D/trunk/meshfem3D/meshfem3D.f90 2010-06-22 13:32:42 UTC (rev 17016)
@@ -196,8 +196,6 @@
integer, dimension(:,:,:,:), allocatable :: ibool
- double precision DT
-
! arrays with the mesh in double precision
double precision, dimension(:,:,:,:), allocatable :: xstore,ystore,zstore
@@ -340,7 +338,7 @@
LOCAL_PATH,SUPPRESS_UTM_PROJECTION,&
INTERFACES_FILE,NSUBREGIONS,subregions,NMATERIALS,material_properties, &
CREATE_ABAQUS_FILES,CREATE_DX_FILES,&
- USE_REGULAR_MESH,NDOUBLINGS,ner_doublings,DT)
+ USE_REGULAR_MESH,NDOUBLINGS,ner_doublings)
if (sizeprocs == 1 .and. (NPROC_XI /= 1 .or. NPROC_ETA /= 1)) then
stop 'must have NPROC_XI = NPROC_ETA = 1 for a serial run'
@@ -731,7 +729,7 @@
NSUBREGIONS,subregions,number_of_layers,ner_layer,NMATERIALS,material_properties, &
myrank,LOCAL_PATH,UTM_X_MIN,UTM_X_MAX,UTM_Y_MIN,UTM_Y_MAX,Z_DEPTH_BLOCK,&
CREATE_ABAQUS_FILES,CREATE_DX_FILES,&
- USE_REGULAR_MESH,NDOUBLINGS,ner_doublings,DT)
+ USE_REGULAR_MESH,NDOUBLINGS,ner_doublings)
if(myrank == 0) then
Modified: seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90 2010-06-22 09:08:08 UTC (rev 17015)
+++ seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90 2010-06-22 13:32:42 UTC (rev 17016)
@@ -32,7 +32,7 @@
LOCAL_PATH,SUPPRESS_UTM_PROJECTION,&
INTERFACES_FILE,NSUBREGIONS,subregions,NMATERIALS,material_properties,&
CREATE_ABAQUS_FILES,CREATE_DX_FILES,&
- USE_REGULAR_MESH,NDOUBLINGS,ner_doublings,DT)
+ USE_REGULAR_MESH,NDOUBLINGS,ner_doublings)
implicit none
@@ -81,9 +81,6 @@
integer i,ireg,imat
-! anticipated time step
- double precision DT
-
! open parameter file
open(unit=IIN,file='DATA/Par_file',status='old',action='read')
@@ -144,10 +141,6 @@
call read_value_logical(IIN,IGNORE_JUNK,CREATE_DX_FILES, 'mesher.CREATE_DX_FILES')
if(err_occurred() /= 0) return
-
- call read_value_double_precision(IIN,IGNORE_JUNK,DT, 'mesher.NPROC_DT')
- if(err_occurred() /= 0) return
-
! file in which we store the databases
call read_value_string(IIN,IGNORE_JUNK,LOCAL_PATH, 'LOCAL_PATH')
if(err_occurred() /= 0) return
More information about the CIG-COMMITS
mailing list