[cig-commits] r12543 - seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Wed Aug 6 06:42:54 PDT 2008
Author: dkomati1
Date: 2008-08-06 06:42:53 -0700 (Wed, 06 Aug 2008)
New Revision: 12543
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/attenuation_model.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/compute_arrays_source.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_list_files_chunks.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_regions_mesh.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/exit_mpi.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/jp3d1994_model.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_sources.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/meshfem3D.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/s362ani.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/save_header_file.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/sea99_s_model.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.f90
Log:
added action='read' or action='write' to open() statements
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/attenuation_model.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/attenuation_model.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/attenuation_model.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -293,7 +293,7 @@
call MPI_COMM_RANK(MPI_COMM_WORLD, myrank, ier)
if(myrank == 0 .AND. iregion_code == IREGION_CRUST_MANTLE .AND. first_time_called == 1) then
first_time_called = 0
- open(unit=27,file=prname(1:len_trim(prname))//'1D_Q.bin',status='unknown',form='unformatted')
+ open(unit=27,file=prname(1:len_trim(prname))//'1D_Q.bin',status='unknown',form='unformatted',action='write')
write(27) AM_V%QT_c_source
write(27) AM_V%Qtau_s
write(27) AM_V%Qn
@@ -660,7 +660,7 @@
if(myrank == 0 .AND. iregion_code == IREGION_CRUST_MANTLE .AND. first_time_called == 1) then
first_time_called = 0
- open(unit=27, file=prname(1:len_trim(prname))//'1D_Q.bin', status='unknown', form='unformatted')
+ open(unit=27, file=prname(1:len_trim(prname))//'1D_Q.bin', status='unknown', form='unformatted',action='read')
read(27) AM_V%QT_c_source
read(27) tau_s
read(27) AM_V%Qn
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/compute_arrays_source.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/compute_arrays_source.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/compute_arrays_source.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -215,7 +215,7 @@
do icomp = 1, NDIM
filename = 'SEM/'//trim(adj_source_file) // '.'// comp(icomp) // '.adj'
- open(unit = IIN, file = trim(filename), iostat = ios)
+ open(unit = IIN, file = trim(filename), iostat = ios, status='old', action='read')
if (ios /= 0) call exit_MPI(myrank, ' file '//trim(filename)//' does not exist')
do itime = 1, NSTEP
read(IIN,*) junk, adj_src(itime,icomp)
@@ -298,7 +298,7 @@
do icomp = 1, NDIM
filename = 'SEM/'//trim(adj_source_file) // '.'// comp(icomp) // '.adj'
- open(unit = IIN, file = trim(filename), iostat = ios)
+ open(unit = IIN, file = trim(filename), iostat = ios, status='old', action='read')
if (ios /= 0) call exit_MPI(myrank, ' file '//trim(filename)//'does not exist')
do itime =1,iadjsrc(it_sub_adj,1)-1
read(IIN,*) junk,junk
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_list_files_chunks.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_list_files_chunks.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_list_files_chunks.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -161,7 +161,7 @@
call get_value_string(OUTPUT_FILES, 'OUTPUT_FILES', 'OUTPUT_FILES')
! file to store the list of processors for each message for faces
- open(unit=IOUT,file=trim(OUTPUT_FILES)//'/list_messages_faces.txt',status='unknown')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//'/list_messages_faces.txt',status='unknown',action='write')
endif
@@ -510,7 +510,7 @@
endif
! file to store the list of processors for each message for corners
- if(myrank == 0) open(unit=IOUT,file=trim(OUTPUT_FILES)//'/list_messages_corners.txt',status='unknown')
+ if(myrank == 0) open(unit=IOUT,file=trim(OUTPUT_FILES)//'/list_messages_corners.txt',status='unknown',action='write')
! loop over all the messages to create the addressing
do imsg = 1,NCORNERSCHUNKS
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_regions_mesh.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_regions_mesh.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/create_regions_mesh.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -2354,7 +2354,7 @@
call exit_mpi(myrank,'Not the same number of 670 surface elements')
! writing surface topology databases
- open(unit=27,file=prname(1:len_trim(prname))//'boundary_disc.bin',status='unknown',form='unformatted')
+ open(unit=27,file=prname(1:len_trim(prname))//'boundary_disc.bin',status='unknown',form='unformatted',action='write')
write(27) NSPEC2D_MOHO, NSPEC2D_400, NSPEC2D_670
write(27) ibelm_moho_top
write(27) ibelm_moho_bot
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/exit_mpi.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/exit_mpi.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/exit_mpi.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -54,7 +54,7 @@
! write error message to file
call get_value_string(OUTPUT_FILES, 'OUTPUT_FILES', 'OUTPUT_FILES')
write(outputname,"('/error_message',i6.6,'.txt')") myrank
- open(unit=IERROR,file=trim(OUTPUT_FILES)//outputname,status='unknown')
+ open(unit=IERROR,file=trim(OUTPUT_FILES)//outputname,status='unknown',action='write')
write(IERROR,*) error_msg(1:len(error_msg))
write(IERROR,*) 'Error detected, aborting MPI... proc ',myrank
close(IERROR)
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/jp3d1994_model.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/jp3d1994_model.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/jp3d1994_model.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -125,12 +125,12 @@
type (jp3d_model_variables) JP3DM_V
! jp3d_model_variables
- OPEN(2,FILE="DATA/Zhao_JP_model/m3d1341")
- OPEN(3,FILE="DATA/Zhao_JP_model/datadis")
+ open(unit=2,FILE="DATA/Zhao_JP_model/m3d1341",status='old',action='read')
+ open(unit=3,FILE="DATA/Zhao_JP_model/datadis",status='old',action='read')
- CALL INPUTJP(JP3DM_V)
- CALL INPUT1(JP3DM_V)
- CALL INPUT2(JP3DM_V)
+ call INPUTJP(JP3DM_V)
+ call INPUT1(JP3DM_V)
+ call INPUT2(JP3DM_V)
end subroutine read_iso3d_dpzhao_model
!==========================================================================
@@ -220,9 +220,9 @@
HE = (ONE - radius)*R_EARTH_KM
! calculate depths of the Conrad, the Moho and
! the plate boundary beneath the location (PHI,RAM)
- CALL HLAY(PE,RE,H1,1,JP3DM_V)
- CALL HLAY(PE,RE,H2,2,JP3DM_V)
- CALL HLAY(PE,RE,H3,3,JP3DM_V)
+ call HLAY(PE,RE,H1,1,JP3DM_V)
+ call HLAY(PE,RE,H2,2,JP3DM_V)
+ call HLAY(PE,RE,H3,3,JP3DM_V)
! when LAY = 1, the focus is in the upper crust;
! when LAY = 2, the focus is in the lower crust;
! when LAY = 3, the focus is in the mantle wedge;
@@ -238,10 +238,10 @@
ELSE
LAY = 4
END IF
- CALL VEL1D(HE,vp,LAY,1,JP3DM_V)
- CALL VEL1D(HE,vs,LAY,2,JP3DM_V)
+ call VEL1D(HE,vp,LAY,1,JP3DM_V)
+ call VEL1D(HE,vs,LAY,2,JP3DM_V)
- CALL VEL3(PE,RE,HE,dvp,LAY,JP3DM_V)
+ call VEL3(PE,RE,HE,dvp,LAY,JP3DM_V)
dvp = 0.01d0*dvp
dvs = 1.5d0*dvp
vp = vp*(1.0d0+dvp)
@@ -341,10 +341,10 @@
100 FORMAT(3I3)
READ(2,100) JP3DM_V%NPA,JP3DM_V%NRA,JP3DM_V%NHA
- CALL PUT1(JP3DM_V%NPA,JP3DM_V%NRA,JP3DM_V%NHA,JP3DM_V%PNA,JP3DM_V%RNA,JP3DM_V%HNA,JP3DM_V%VELAP)
+ call PUT1(JP3DM_V%NPA,JP3DM_V%NRA,JP3DM_V%NHA,JP3DM_V%PNA,JP3DM_V%RNA,JP3DM_V%HNA,JP3DM_V%VELAP)
READ(2,100) JP3DM_V%NPB,JP3DM_V%NRB,JP3DM_V%NHB
- CALL PUT1(JP3DM_V%NPB,JP3DM_V%NRB,JP3DM_V%NHB,JP3DM_V%PNB,JP3DM_V%RNB,JP3DM_V%HNB,JP3DM_V%VELBP)
- CALL BLDMAP(JP3DM_V)
+ call PUT1(JP3DM_V%NPB,JP3DM_V%NRB,JP3DM_V%NHB,JP3DM_V%PNB,JP3DM_V%RNB,JP3DM_V%HNB,JP3DM_V%VELBP)
+ call BLDMAP(JP3DM_V)
RETURN
END SUBROUTINE INPUT1
@@ -530,9 +530,9 @@
type (jp3d_model_variables) JP3DM_V
! jp3d_model_variables
- CALL LOCX(JP3DM_V%PNA,JP3DM_V%RNA,JP3DM_V%HNA,JP3DM_V%NPA,JP3DM_V%NRA,JP3DM_V%NHA,MKA, &
+ call LOCX(JP3DM_V%PNA,JP3DM_V%RNA,JP3DM_V%HNA,JP3DM_V%NPA,JP3DM_V%NRA,JP3DM_V%NHA,MKA, &
JP3DM_V%PLA,JP3DM_V%RLA,JP3DM_V%HLA,JP3DM_V%IPLOCA,JP3DM_V%IRLOCA,JP3DM_V%IHLOCA)
- CALL LOCX(JP3DM_V%PNB,JP3DM_V%RNB,JP3DM_V%HNB,JP3DM_V%NPB,JP3DM_V%NRB,JP3DM_V%NHB,MKB, &
+ call LOCX(JP3DM_V%PNB,JP3DM_V%RNB,JP3DM_V%HNB,JP3DM_V%NPB,JP3DM_V%NRB,JP3DM_V%NHB,MKB, &
JP3DM_V%PLB,JP3DM_V%RLB,JP3DM_V%HLB,JP3DM_V%IPLOCB,JP3DM_V%IRLOCB,JP3DM_V%IHLOCB)
RETURN
END
@@ -653,10 +653,10 @@
JP3DM_V%R = RE/DEGREES_TO_RADIANS
JP3DM_V%H = HE
IF(LAY.LE.3) THEN
- CALL PRHF(JP3DM_V%IPLOCA,JP3DM_V%IRLOCA,JP3DM_V%IHLOCA,JP3DM_V%PLA,JP3DM_V%RLA,JP3DM_V%HLA, &
+ call PRHF(JP3DM_V%IPLOCA,JP3DM_V%IRLOCA,JP3DM_V%IHLOCA,JP3DM_V%PLA,JP3DM_V%RLA,JP3DM_V%HLA, &
JP3DM_V%PNA,JP3DM_V%RNA,JP3DM_V%HNA,MPA,MRA,MHA,MKA,JP3DM_V)
ELSE IF(LAY.EQ.4) THEN
- CALL PRHF(JP3DM_V%IPLOCB,JP3DM_V%IRLOCB,JP3DM_V%IHLOCB,JP3DM_V%PLB,JP3DM_V%RLB,JP3DM_V%HLB, &
+ call PRHF(JP3DM_V%IPLOCB,JP3DM_V%IRLOCB,JP3DM_V%IHLOCB,JP3DM_V%PLB,JP3DM_V%RLB,JP3DM_V%HLB, &
JP3DM_V%PNB,JP3DM_V%RNB,JP3DM_V%HNB,MPB,MRB,MHB,MKB,JP3DM_V)
ELSE
END IF
@@ -670,9 +670,9 @@
JP3DM_V%WV(8) = JP3DM_V%PF*JP3DM_V%RF*JP3DM_V%HF
! calculate velocity
IF(LAY.LE.3) THEN
- CALL VABPS(MPA,MRA,MHA,JP3DM_V%VELAP,V,JP3DM_V)
+ call VABPS(MPA,MRA,MHA,JP3DM_V%VELAP,V,JP3DM_V)
ELSE IF(LAY.EQ.4) THEN
- CALL VABPS(MPB,MRB,MHB,JP3DM_V%VELBP,V,JP3DM_V)
+ call VABPS(MPB,MRB,MHB,JP3DM_V%VELBP,V,JP3DM_V)
ELSE
END IF
@@ -848,12 +848,12 @@
integer :: IPLOCX(MKX),IRLOCX(MKX),IHLOCX(MKX)
double precision :: PNX(MPX),RNX(MRX),HNX(MHX)
double precision :: PLX,RLX,HLX
- CALL LIMIT(PNX(1),PNX(MPX),JP3DM_V%P)
- CALL LIMIT(RNX(1),RNX(MRX),JP3DM_V%R)
- CALL LIMIT(HNX(1),HNX(MHX),JP3DM_V%H)
- CALL INTMAP(JP3DM_V%P*100.0,IPLOCX,MKX,PLX,JP3DM_V%IP)
- CALL INTMAP(JP3DM_V%R*100.0,IRLOCX,MKX,RLX,JP3DM_V%JP)
- CALL INTMAP(JP3DM_V%H,IHLOCX,MKX,HLX,JP3DM_V%KP)
+ call LIMIT(PNX(1),PNX(MPX),JP3DM_V%P)
+ call LIMIT(RNX(1),RNX(MRX),JP3DM_V%R)
+ call LIMIT(HNX(1),HNX(MHX),JP3DM_V%H)
+ call INTMAP(JP3DM_V%P*100.0,IPLOCX,MKX,PLX,JP3DM_V%IP)
+ call INTMAP(JP3DM_V%R*100.0,IRLOCX,MKX,RLX,JP3DM_V%JP)
+ call INTMAP(JP3DM_V%H,IHLOCX,MKX,HLX,JP3DM_V%KP)
JP3DM_V%IP1 = JP3DM_V%IP+1
JP3DM_V%JP1 = JP3DM_V%JP+1
JP3DM_V%KP1 = JP3DM_V%KP+1
@@ -944,8 +944,8 @@
integer :: IJK,J,J1,I,I1
P = 90.0-PE/DEGREES_TO_RADIANS
R = RE/DEGREES_TO_RADIANS
- CALL LIMIT(JP3DM_V%PN(1),JP3DM_V%PN(51),P)
- CALL LIMIT(JP3DM_V%RRN(1),JP3DM_V%RRN(63),R)
+ call LIMIT(JP3DM_V%PN(1),JP3DM_V%PN(51),P)
+ call LIMIT(JP3DM_V%RRN(1),JP3DM_V%RRN(63),R)
DO 1 I = 1,50
I1 = I+1
IF(P.GE.JP3DM_V%PN(I).AND.P.LT.JP3DM_V%PN(I1)) GO TO 11
@@ -1069,10 +1069,10 @@
ELSE IF(LAY.GE.3) THEN
HM = 40.0
IF(HE.LT.HM) THEN
- CALL JPMODEL(IPS,HM,VM,JP3DM_V)
+ call JPMODEL(IPS,HM,VM,JP3DM_V)
V = VM-(HM-HE)*0.003
ELSE
- CALL JPMODEL(IPS,HE,V,JP3DM_V)
+ call JPMODEL(IPS,HE,V,JP3DM_V)
END IF
ELSE
END IF
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_receivers.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -349,7 +349,7 @@
! Harvard format does not support the network name
! therefore only the station name is included below
! compute total number of samples for normal modes with 1 sample per second
- open(unit=1,file=trim(OUTPUT_FILES)//'/RECORDHEADERS',status='unknown')
+ open(unit=1,file=trim(OUTPUT_FILES)//'/RECORDHEADERS',status='unknown',action='write')
nsamp = nint(dble(NSTEP-1)*DT)
do irec = 1,nrec
@@ -615,7 +615,7 @@
epidist(1:nrec) = epidist_found(1:nrec)
! write the list of stations and associated epicentral distance
- open(unit=27,file=trim(OUTPUT_FILES)//'/output_list_stations.txt',status='unknown')
+ open(unit=27,file=trim(OUTPUT_FILES)//'/output_list_stations.txt',status='unknown',action='write')
write(27,*)
write(27,*) 'total number of stations: ',nrec
write(27,*)
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_sources.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_sources.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/locate_sources.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -611,6 +611,12 @@
write(IMAIN,*)
write(IMAIN,*) 'printing the source-time function'
+ scalar_moment = 0.
+ do i = 1,6
+ scalar_moment = scalar_moment + moment_tensor(i,isource)**2
+ enddo
+ scalar_moment = dsqrt(scalar_moment/2.)
+
! print the source-time function
if(NSOURCES == 1) then
plot_file = '/plot_source_time_function.txt'
@@ -623,14 +629,8 @@
write(plot_file,"('/plot_source_time_function',i3,'.txt')") isource
endif
endif
- open(unit=27,file=trim(OUTPUT_FILES)//plot_file,status='unknown')
- scalar_moment = 0.
- do i = 1,6
- scalar_moment = scalar_moment + moment_tensor(i,isource)**2
- enddo
- scalar_moment = dsqrt(scalar_moment/2.)
-
+ open(unit=27,file=trim(OUTPUT_FILES)//plot_file,status='unknown',action='write')
do it=1,NSTEP
time_source = dble(it-1)*DT-t0-t_cmt(isource)
write(27,*) sngl(dble(it-1)*DT-t0),sngl(scalar_moment*comp_source_time_function(time_source,hdur_gaussian(isource)))
@@ -652,8 +652,8 @@
write(plot_file,"('/plot_source_spectrum',i3,'.txt')") isource
endif
endif
- open(unit=27,file=trim(OUTPUT_FILES)//plot_file,status='unknown')
+ open(unit=27,file=trim(OUTPUT_FILES)//plot_file,status='unknown',action='write')
do iom=1,NSAMP_PLOT_SOURCE
om=TWO_PI*(1.0d0/8.0d0)*(iom-1)/dble(NSAMP_PLOT_SOURCE-1)
write(27,*) sngl(om/TWO_PI),sngl(scalar_moment*om*comp_source_spectrum(om,hdur(isource)))
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/meshfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/meshfem3D.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/meshfem3D.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -654,7 +654,7 @@
! open main output file, only written to by process 0
if(myrank == 0 .and. IMAIN /= ISTANDARD_OUTPUT) &
- open(unit=IMAIN,file=trim(OUTPUT_FILES)//'/output_mesher.txt',status='unknown')
+ open(unit=IMAIN,file=trim(OUTPUT_FILES)//'/output_mesher.txt',status='unknown',action='write')
! get MPI starting time
time_start = MPI_WTIME()
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/s362ani.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/s362ani.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/s362ani.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -747,7 +747,7 @@
integer :: ncoef,i,ihor,ifst,ilst,ifst1,ios,lstr,nmodkern,idummy,nhorpar,lmax
- open(lu,file=filename,iostat=ios)
+ open(unit=lu,file=filename,iostat=ios,status='old',action='read')
if(ios /= 0) then
stop 'error opening 3-d model'
endif
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/save_header_file.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/save_header_file.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/save_header_file.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -102,7 +102,7 @@
! copy number of elements and points in an include file for the solver
call get_value_string(HEADER_FILE, 'solver.HEADER_FILE', 'OUTPUT_FILES/values_from_mesher.h')
- open(unit=IOUT,file=HEADER_FILE,status='unknown')
+ open(unit=IOUT,file=HEADER_FILE,status='unknown',action='write')
write(IOUT,*)
write(IOUT,*) '!'
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/sea99_s_model.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/sea99_s_model.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/sea99_s_model.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -54,9 +54,8 @@
!----------------------- choose input file: ------------------
! relative anomaly
+ open(unit=1,file='DATA/Lebedev_sea99/sea99_dvsvs',status='old',action='read')
- open(1,file='DATA/Lebedev_sea99/sea99_dvsvs')
-
!----------------------- read input file: ------------------
do i = 1, 6
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -904,7 +904,7 @@
! open main output file, only written to by process 0
if(myrank == 0 .and. IMAIN /= ISTANDARD_OUTPUT) &
- open(unit=IMAIN,file=trim(OUTPUT_FILES)//'/output_solver.txt',status='unknown')
+ open(unit=IMAIN,file=trim(OUTPUT_FILES)//'/output_solver.txt',status='unknown',action='write')
if(myrank == 0) then
@@ -2000,7 +2000,7 @@
! create an empty file to monitor the start of the simulation
if(myrank == 0) then
- open(unit=IOUT,file=trim(OUTPUT_FILES)//'/starttimeloop.txt',status='unknown')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//'/starttimeloop.txt',status='unknown',action='write')
write(IOUT,*) 'hello, starting time loop'
close(IOUT)
endif
@@ -2121,7 +2121,7 @@
! write time stamp file to give information about progression of simulation
write(outputname,"('/timestamp',i6.6)") it
- open(unit=IOUT,file=trim(OUTPUT_FILES)//outputname,status='unknown')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//outputname,status='unknown',action='write')
write(IOUT,*) 'Time step # ',it
write(IOUT,*) 'Time: ',sngl(((it-1)*DT-t0)/60.d0),' minutes'
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.f90 2008-08-06 13:08:18 UTC (rev 12542)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.f90 2008-08-06 13:42:53 UTC (rev 12543)
@@ -112,17 +112,17 @@
if(USE_BINARY_FOR_LARGE_FILE) then
if (seismo_offset==0) then
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin',status='unknown',form='unformatted')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin',status='unknown',form='unformatted',action='write')
else
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin',status='old',&
- form='unformatted',position='append')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin', &
+ status='old',form='unformatted',position='append',action='write')
endif
else
if (seismo_offset==0) then
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii',status='unknown',form='formatted')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii',status='unknown',form='formatted',action='write')
else
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii',status='old',&
- form='formatted',position='append')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii', &
+ status='old',form='formatted',position='append',action='write')
endif
endif
endif
@@ -176,17 +176,17 @@
if(USE_BINARY_FOR_LARGE_FILE) then
if (seismo_offset==0) then
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin',status='unknown',form='unformatted')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin',status='unknown',form='unformatted',action='write')
else
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin',status='old',&
- form='unformatted',position='append')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.bin', &
+ status='old',form='unformatted',position='append',action='write')
endif
else
if (seismo_offset==0) then
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii',status='unknown',form='formatted')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii',status='unknown',form='formatted',action='write')
else
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii',status='old',&
- form='formatted',position='append')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname)//'.ascii', &
+ status='old',form='formatted',position='append',action='write')
endif
endif
@@ -649,9 +649,9 @@
endif
else
if (seismo_offset==0) then
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname_2),status='unknown')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname_2),status='unknown',action='write')
else
- open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname_2),status='old',position='append')
+ open(unit=IOUT,file=trim(OUTPUT_FILES)//trim(sisname_2),status='old',position='append',action='write')
endif
endif
More information about the cig-commits
mailing list