added writing out of general anisotropy parameters as 21 binary format files

Author: aakash10gupta

src/generate_databases/save_arrays_solver.F90

@@ -443,7 +443,7 @@ subroutine save_arrays_solver_files()
     myrank
 
   use shared_parameters, only: ACOUSTIC_SIMULATION, ELASTIC_SIMULATION, POROELASTIC_SIMULATION, &
-    NPROC
+    NPROC, ANISOTROPY
 
   ! global indices
   use generate_databases_par, only: nspec => NSPEC_AB, ibool
@@ -553,6 +553,115 @@ subroutine save_arrays_solver_files()
   write(IOUT) v_tmp
   close(IOUT)
 
+  if (ANISOTROPY) then
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c11.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c11.bin'
+    write(IOUT) c11store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c12.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c12.bin'
+    write(IOUT) c12store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c13.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c13.bin'
+    write(IOUT) c13store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c14.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c14.bin'
+    write(IOUT) c14store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c15.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c15.bin'
+    write(IOUT) c15store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c16.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c16.bin'
+    write(IOUT) c16store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c22.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c22.bin'
+    write(IOUT) c22store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c23.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c23.bin'
+    write(IOUT) c23store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c24.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c24.bin'
+    write(IOUT) c24store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c25.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c25.bin'
+    write(IOUT) c25store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c26.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c26.bin'
+    write(IOUT) c26store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c33.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c33.bin'
+    write(IOUT) c33store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c34.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c34.bin'
+    write(IOUT) c34store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c35.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c35.bin'
+    write(IOUT) c35store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c36.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c36.bin'
+    write(IOUT) c36store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c44.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c44.bin'
+    write(IOUT) c44store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c45.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c45.bin'
+    write(IOUT) c45store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c46.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c46.bin'
+    write(IOUT) c46store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c55.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c55.bin'
+    write(IOUT) c55store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c56.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c56.bin'
+    write(IOUT) c56store
+    close(IOUT)
+
+    open(unit=IOUT,file=prname(1:len_trim(prname))//'c66.bin',status='unknown',form='unformatted',iostat=ier)
+    if (ier /= 0) stop 'error opening file c66.bin'
+    write(IOUT) c66store
+    close(IOUT)
+
+  endif
+
   ! attenuation
   ! shear attenuation Qmu
   open(unit=IOUT,file=prname(1:len_trim(prname))//'qmu.bin',status='unknown',form='unformatted',iostat=ier)