Jose:
      van los datos de presiones (en MPa) y saturaciones de CO2 para 1 y 3  anios.

sco2_1year.dat       CO2 saturation after 1 year injection 


pco2_1year.dat        CO2 pressure (MPa)  after 1 year injection


sco2_3years.dat       CO2 saturation after 3 years  injection 


pco2_3years.dat        CO2 pressure (MPa)  after 3 years injection


El dominio para el simulador es de  1200m en x  \times  10 Km en y  \times 400 m en z
   
 La malla que se uso para el simulador de flujo BOAST es de 

300   en   x  (hx= 4m)

 5    en   y  (hy=2 km)

400   en z    (hz=1 m)

y el punto de inyeccion esta en la 

celda 3   en la direccion  y

celda  150 en la direccion x 

celda  285 en la direccion z



Por lo tanto representa una inyeccion  en una linea 

de 2 km   en la direccion   y

 
El grafico correponde a un corte vertical en el nodo 3, donde esta el punto de inyeccion


Van tambien graficos .eps  de presion (en MPa)  y saturacion de CO2  para 1 y 3  anios:

 pco2_gnu.dat_1year.eps 

 pco2_gnu.dat_3years.eps 

 sco2_gnu.dat_1year.eps 

 sco2_gnu.dat_3years.eps


 y los files  

pco2_gnu_1year.dat
 
pco2_gnu_3years.dat 

sco2_gnu_1year.dat

sco2_gnu_3years.dat


usados  para generar los

.eps 

con gnuplot


tambien hay unas salidas del simulador sismico para onda plana:


gather_vz_line_600x200_1a_chimney.ps

uz_*.eps : snapshots de la comp. vertical del desplaz. (se ve el pushdown aqui)


vp*.eps, vs*.eps, qp*.eps qs*.eps  : vp vs qp qs a 50 Hz

Como se cambiaron los datos de porosidad y Clay content:

a mano,  se eliminaron las variaciones fractales de porosidad y  Clay content 
en todo el dominio y ademas en la zona de la chimenea 
se cambio el Clay content y la porosidad 

DEbajo va la parte del codigo

gen_petro_model_for_BOAST.f 

donde se hace, alli estan los valores......


     			do j=1,nz

				do i=1,nx
* lee el clay content

				read(21,*) k, l, C(k,l)
* Juan removes fractal from Clay content within the Utsira	

                             if(C(k,l).lt.10.d0) C(k,l)=3.d0 		
                              C(k,l)=C(k,l)*1.d-2
* reads porosity

				read(22,*) k, l, phi0(k,l) 
* Juan removes fractal from  porosity within the Utsira

		 if(phi0(k,l).gt.30.d0) phi0(k,l)=37.d0
				phi0(k,l)=phi0(k,l)*1.d-2



* Juan : chimney construction

		do j=135,165   ! ancho de chimenea en x

		   do k=54,280  ! ancho de chimenea en z 
*  54 es la celda tope de la chimenea, justo debajo del sello arriba en la Utsira
			
* 280  el la celda mas profunda, justo arriba del punto de inyeccion.

* cambio a Clay content = .65 y porosity =  .25 en la chimenea

                         
                      if(dabs(C(j,k)-0.85).lt.1.d-5) then
                                 C(j,k)=0.65d0

				 phi0(j,k)=0.25d0

* ademas, debajo hay cambios  a manos celda por celda
* artwork......

			if(j.eq.150)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			if(j.eq.151)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			if(j.eq.152)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			
			if(j.eq.153)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			if(j.eq.154)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			if(j.eq.155)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			
			if(j.eq.156)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
			
			
			if(j.eq.157)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			
			if(j.eq.158)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
						
			
			if(j.eq.160)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif
			
			if(j.eq.161)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif
			
			if(j.eq.162)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif
			
			if(j.eq.163)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif			
			
			if(j.eq.165)then
			C(j,k)=0.50d0
			phi0(j,k)=0.25d0
			endif

* I tried the ones below, but gave me a wider chimney 

					
c			if(j.eq.166)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif
			
c			if(j.eq.168)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif			
			
			
c			if(j.eq.169)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif			
						
			
			if(j.eq.149)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
			

			
			if(j.eq.148)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
						
			
			
			if(j.eq.147)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
			
			
			if(j.eq.146)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
			

			
			if(j.eq.145)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
						
			
			
			if(j.eq.144)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif
						
			
			if(j.eq.143)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
						
			
			
			if(j.eq.142)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif			
			
			if(j.eq.141)then
			C(j,k)=0.35d0
			phi0(j,k)=0.30d0
			endif	
			
			
			if(j.eq.139)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif
			
			if(j.eq.138)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif

			if(j.eq.137)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif
			
			if(j.eq.136)then
			C(j,k)=0.40d0
			phi0(j,k)=0.28d0
			endif
			
			
* I tried the ones below, but gave me a wider chimney 
			
c			if(j.eq.134)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif
			
c			if(j.eq.133)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif
			
c			if(j.eq.131)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif
			
c			if(j.eq.130)then
c			C(j,k)=0.50d0
c			phi0(j,k)=0.25d0
c			endif			
			
			
			
                      end if
                end do

                end do		



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The data and codes used to get this data come from the directories


[santos@banach]:~/santos2013/codes/co2_flow_and_seismics/test20

[santos@banach]:~/santos2013/codes/co2_flow_and_seismics/test20_3a