diff --git a/experiments/SConstruct b/experiments/SConstruct new file mode 100644 index 0000000..25c05bd --- /dev/null +++ b/experiments/SConstruct @@ -0,0 +1,93 @@ +#!/usr/bin/env python +# +# coding: utf-8 +# +# SConstruct (Madagascar Script) +# +# Purpose: Tracked build of the experiments in this directory and +# automatic report. +# +# Site: https://dirack.github.io +# +# Version 1.0 +# +# Programer: Rodolfo A C Neves (Dirack) 10/09/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +from rsf.tex import * +from rsf.proj import Flow +import os +import logging +import tempfile + +env = Environment(ENV=os.environ) +env.AlwaysBuild('paper.tex') + +logfile = tempfile.mkstemp(suffix='.log')[1] + +logging.basicConfig(filename=logfile,level=logging.INFO, + format='%(asctime)s - %(levelname)s - %(message)s') + +logging.info('PSTA ***[Experiment Started]***') + +#experiments = ['fullInterpolationAndStack', +# 'fullInterpolationAndStackCDS', +experiments = [ 'multiLayerModel/cds', + 'multiLayerModel/cre'] +sconscripts = filter(os.path.isfile,map(lambda e: str(e)+'/SConstruct',experiments)) + +logging.debug('experiments = %s' % experiments) +logging.debug('sconscripts = %s' % sconscripts) + +succed=0 +fail=0 +total=len(experiments) + +for i in sconscripts: + try: + SConscript(i,exports='env') + logging.info("Experiment %s Succed!",i) + succed = succed + 1 + except: + logging.error("PERR Experiment %s Failed!",i) + fail = fail + 1 + currentExperimentDir = os.path.dirname(i) + logging.warning( + '''Experiment Failed - Show files list\n\t Files=> %s + ''', str(os.listdir(currentExperimentDir))) + + if os.path.isdir(currentExperimentDir+'/'+'Fig'): + logging.warning( + '''Experiment Failed - Show Fig list\n\t Fig=> %s + ''', str(os.listdir(currentExperimentDir+'/'+'Fig'))) + + # Try to run a BUGSCRIPT (This is a SConscript with another name + # to allow the experiment to do a set of commands if the build fails) + if os.path.isfile(currentExperimentDir+'/BUGSCRIPT'): + try: + SConscript(currentExperimentDir+'/BUGSCRIPT') + except: + logging.warning("BUGSCRIPT failed!") + +logging.info('PEND Total of '+str(fail+succed)+'/'+str(total)+' experiments executed, '+str(fail)+' failed and '+str(succed)+' succed') + +logging.disable(logging.CRITICAL) +Flow('paper.tex',None, + ''' + ./log2tex %s + ''' %(logfile)) + +End(options='reproduce', + include=r''' + \usepackage[utf8]{inputenc} + \usepackage{microtype} + \usepackage{listings} + \usepackage[left=2.5cm,top=2.5cm,right=2.5cm,bottom=2.5cm]{geometry} + \setlength{\parindent}{1.3cm} + \usepackage{color} + \definecolor{mygray}{rgb}{0.8,0.8,0.8} + \lstset{backgroundcolor=\color{mygray},breaklines=true} + ''') diff --git a/experiments/fullInterpolationAndStack/SConstruct b/experiments/fullInterpolationAndStack/SConstruct index cba23ea..4de1269 100644 --- a/experiments/fullInterpolationAndStack/SConstruct +++ b/experiments/fullInterpolationAndStack/SConstruct @@ -3,12 +3,12 @@ # # SConstruct (Madagascar Script) # -# Purpose: Build the interpolation process twice to improve CMP sampling -# and CRE stacking results. +# Purpose: Common Reflection Element (CRE) stack applied to a seismic data cube +# from Kirchhoff modeling of a gaussian reflector in a linear velocity model. # -# Site: http://www.dirackslounge.online +# Site: https://dirack.github.io # -# Version 1.0 +# Version 2.0 # # Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 # @@ -19,45 +19,23 @@ # Madagascar package from rsf.proj import * -# CRE recipe -from creStack import kirchoffModeling as kimod -from creStack import pefInterpolation as pefin - # Import glob python library import glob -# Generate Gaussian reflector model and data cube -kimod() - -# PEF interpolation of the data cube -pefin('dataCube', - 'interpolatedDataCube', - nm=401, - dm=0.025, - nt=1001, - dt=0.004, - nhi=161) - -# Do the PEF interpolation one more time -# to increase the CMP sampling -pefin('interpolatedDataCube', - 'interpolatedDataCube2', - nm=802, - dm=0.0125, - nt=1001, - dt=0.004, - nhi=161) - # CRE stacking # It uses Very Fast Simulated Aneeling and non hyperbolic CRS # to get zero offset CRS parameters (RN, RNIP and BETA) from data cube v0 = 1.5 ot0 = 1.0 dt0 = 0.004 -nt0 = 500 +nt0 = 500 om0 = 3 dm0 = 0.5 nm0 = 9 +inputDir = '../modelingAndPEFInterpolation/gaussianModel/' +dataCube = inputDir+'dataCube.rsf' +interpolatedDataCube = inputDir+'interpolatedDataCube2.rsf' +files = [] for i in range(nm0): @@ -82,22 +60,22 @@ for i in range(nm0): creGatherPlot = 'crePlot-m0-%g-t0-%g' % (i,j) # Very Fast Simulated Aneelling Global Optimization (VFSA) - Flow(crsParameters,'dataCube', + Flow(crsParameters,dataCube, ''' vfsacrsnh m0=%g v0=%g t0=%g verb=y repeat=3 ''' % (m0,v0,t0)) # Calculate CRE trajectory - Flow(creMhCoordinates,['interpolatedDataCube2',crsParameters], + Flow(creMhCoordinates,[interpolatedDataCube,crsParameters], ''' cretrajec verb=y m0=%g param=${SOURCES[1]} | put unit1="Offset" label1="Km" ''' % (m0)) #Get CRE Gather from interpolated Data Cube - Flow([creGather,creMcoordinate],['interpolatedDataCube2',creMhCoordinates], + Flow([creGather,creMcoordinate],[interpolatedDataCube,creMhCoordinates], ''' - getcregather verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} | + getcregather verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} aperture=20 | put label1="Time" unit1="s" label2="Offset" unit2="km" ''') @@ -155,9 +133,15 @@ for i in range(nm0): put label1=t0 unit1=s label2=m0 unit2=Km --out=stdout ''') + files.append(creStackedTrace) + +# Stacked section depends of the last CRE stacked trace +Depends('stackedSection.rsf',[dataCube, + 'creStackedTrace-m0-'+str(nm0-1)+'.rsf', + 'creGatherTrace-m0-'+str(nm0-1)+'.rsf']) + # Build the cre stacked section # throughout cre stacked traces sorting -files = glob.glob('creStackedTrace-*.rsf') length = len(files) sortedFiles = [] diff --git a/experiments/fullInterpolationAndStack/creStack.py b/experiments/fullInterpolationAndStack/creStack.py deleted file mode 100644 index f997200..0000000 --- a/experiments/fullInterpolationAndStack/creStack.py +++ /dev/null @@ -1,165 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# creStack (Python) -# -# Purpose: Recipe to the cre model, interpolation and Stacking. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Madagascar package -from rsf.proj import * - -def kirchoffModeling(filename='dataCube'): - ''' - Modeling function of a gaussian reflector - ''' - - # Modeling: Gaussian reflector in a velocity linear model - # velocity increases with depth with a 0.5 velocity gradient - Flow('gaussianReflector',None, - ''' - math d1=0.01 n1=2001 o1=-5 unit1=km label1=Offset - output="4-3*exp(-(x1-5)^2/9)" - ''') - - # Velocity Model - Flow('velocityModel','gaussianReflector', - ''' - window min1=0 max1=10 | - spray axis=1 n=451 d=0.01 o=0 label=Depth unit=km | - math output="1.5+0.5*x1+0.0*x2" - ''') - - Flow('reflectorDip','gaussianReflector','math output="2/3*(x1-5)*input" ') - - # Kirchoff Modeling - Flow(filename,'gaussianReflector reflectorDip', - ''' - kirmod cmp=y dip=${SOURCES[1]} - nh=161 dh=0.025 h0=0 - ns=401 ds=0.025 s0=0 - freq=10 dt=0.004 nt=1001 - vel=1.5 gradz=0.5 gradx=0.0 verb=y | - put d2=0.0125 label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s - ''') - -def pefInterpolation( - dataCube, - interpolated, - nm, - dm, - nt, - dt, - nhi=1 - ): - ''' - PEF interpolation of the data cube - :param dataCube: filename, Seismic data cube to interpolate - :param interpolated: filename, Interpolated seismic data cube - :param nm: integer, number of CMPs in the seismic data cube - :param dm: float, CMP sampling - :param nt: integer, number of time samples - :param dt: float, time sampling - :param nhi: integer, number of constant offsets gathers to interpolate - ''' - - # Divide CMP sampling - dm = dm/2 - - # Define mask file names using input filename - mask1 = dataCube+'-mask1' - mask = dataCube+'-mask' - aa = dataCube+'-aa' - bb = dataCube+'-bb' - a = dataCube+'-a' - b = dataCube+'-b' - zeroTraceGather = dataCube+'-zeroedGather' - mask0 = dataCube+'-mask0' - - - # Build a mask to interleave zero traces with original data traces - Flow(aa,None,'spike n1=%i d1=%g o1=0' %(nm,dm)) - Flow(bb,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask1,[bb, aa], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - - Flow(a,None,'spike n1=%i d1=%g o1=0' % (nm,dm)) - Flow(b,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask,[a, b], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - Flow(zeroTraceGather,b, - ''' - spray axis=2 n=%i d=%g | - transp | - put label2=Offset unit2=Km label1=Time unit1=s - ''' %(nt,dm)) - - # Data Mask with double of traces in CMP (half of CMP sampling) - # Keep the same Time and Offset original data sampling - Flow(mask0,mask, - ''' - spray axis=1 n=%i d=%g - ''' %(nt,dt)) - - totalPefIterations = 100 - totalInterpolationIterations = 20 - - offsetGathers = [] - for offsetGatherIndex in range(nhi): - - offsetGather = dataCube+"-offsetGather-%i" % offsetGatherIndex - resampledOffsetGather = dataCube+"-resampledGather-%i" % offsetGatherIndex - interpolatedOffsetGather = dataCube+"-interpolatedGather-%i" % offsetGatherIndex - pefCoeficients = dataCube+"-pefCoeficients-%i" % offsetGatherIndex - - Flow(offsetGather,dataCube, - ''' - window n2=1 f2=%i - ''' % (offsetGatherIndex)) - - Flow(resampledOffsetGather,[offsetGather,zeroTraceGather], - ''' - interleave axis=2 ${SOURCES[1]} - ''') - - # Calculate adaptive PEF coeficients - Flow(pefCoeficients,[resampledOffsetGather,mask0], - ''' - apef jump=2 a=10,2 rect1=50 rect2=2 niter=%g verb=y - maskin=${SOURCES[1]} - ''' % (totalPefIterations)) - - # Interpolation - Flow(interpolatedOffsetGather, [resampledOffsetGather,pefCoeficients,mask0,mask1], - ''' - miss4 exact=y filt=${SOURCES[1]} mask=${SOURCES[2]} niter=%g verb=y | - put d2=%g - ''' % (totalInterpolationIterations,dm)) - - offsetGathers.append(interpolatedOffsetGather) - - # Concatenate interpolated sections - Flow(interpolated,offsetGathers, - ''' - rcat axis=3 ${SOURCES[1:%d]} | - transp plane=23 - ''' % nhi) - - diff --git a/experiments/fullInterpolationAndStackCDS/SConstruct b/experiments/fullInterpolationAndStackCDS/SConstruct index 9975780..5179dac 100644 --- a/experiments/fullInterpolationAndStackCDS/SConstruct +++ b/experiments/fullInterpolationAndStackCDS/SConstruct @@ -3,17 +3,19 @@ # # SConstruct (Madagascar Script) # -# Purpose: Build the interpolation process twice to improve CMP sampling -# and CRE stacking results. +# Purpose: Common Reflection Element (CRE) stack applied to a seismic data cube +# from Kirchhoff modeling of a gaussian reflector in a linear velocity model. # -# IMPORTANT: This SConstruct uses the non-hyperbolic CRS traveltime approximation -# with the CDS condition applied (RN=RNIP) as a stacking curve. +# Important!: This experiment uses non-hyperbolic Common Reflection Surface (CRS) +# traveltime approximation as stacking curve, with CDS condition (RNIP=RN) applied. +# Instead of conventional CRE traveltime approximation. # -# Site: http://www.dirackslounge.online +# +# Site: https://dirack.github.io # -# Version 1.0 +# Version 2.0 # -# Programer: Rodolfo A. C. Neves (Dirack) 16/03/2020 +# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 # # Email: rodolfo_profissional@hotmail.com # @@ -22,35 +24,9 @@ # Madagascar package from rsf.proj import * -# CRE recipe -from creStack import kirchoffModeling as kimod -from creStack import pefInterpolation as pefin - # Import glob python library import glob -# Generate Gaussian reflector model and data cube -kimod() - -# PEF interpolation of the data cube -pefin('dataCube', - 'interpolatedDataCube', - nm=401, - dm=0.025, - nt=1001, - dt=0.004, - nhi=161) - -# Do the PEF interpolation one more time -# to increase the CMP sampling -pefin('interpolatedDataCube', - 'interpolatedDataCube2', - nm=802, - dm=0.0125, - nt=1001, - dt=0.004, - nhi=161) - # CRE stacking # It uses Very Fast Simulated Aneeling and non hyperbolic CRS # to get zero offset CRS parameters (RN, RNIP and BETA) from data cube @@ -61,6 +37,10 @@ nt0 = 500 om0 = 3 dm0 = 0.5 nm0 = 9 +inputDir = '../modelingAndPEFInterpolation/gaussianModel/' +dataCube = inputDir+'dataCube.rsf' +interpolatedDataCube = inputDir+'interpolatedDataCube2.rsf' +files = [] for i in range(nm0): @@ -85,22 +65,22 @@ for i in range(nm0): creGatherPlot = 'crePlot-m0-%g-t0-%g' % (i,j) # Very Fast Simulated Aneelling Global Optimization (VFSA) - Flow(crsParameters,'dataCube', + Flow(crsParameters,dataCube, ''' vfsacrsnh m0=%g v0=%g t0=%g verb=y repeat=3 ''' % (m0,v0,t0)) # Calculate CRE trajectory - Flow(creMhCoordinates,['interpolatedDataCube2',crsParameters], + Flow(creMhCoordinates,[interpolatedDataCube,crsParameters], ''' cretrajec verb=y m0=%g param=${SOURCES[1]} | put unit1="Offset" label1="Km" ''' % (m0)) #Get CRE Gather from interpolated Data Cube - Flow([creGather,creMcoordinate],['interpolatedDataCube2',creMhCoordinates], + Flow([creGather,creMcoordinate],[interpolatedDataCube,creMhCoordinates], ''' - getcregather verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} | + getcregather verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} aperture=20 | put label1="Time" unit1="s" label2="Offset" unit2="km" ''') @@ -158,9 +138,15 @@ for i in range(nm0): put label1=t0 unit1=s label2=m0 unit2=Km --out=stdout ''') + files.append(creStackedTrace) + +# Stacked section depends of the last CRE stacked trace +Depends('stackedSection.rsf',[dataCube, + 'creStackedTrace-m0-'+str(nm0-1)+'.rsf', + 'creGatherTrace-m0-'+str(nm0-1)+'.rsf']) + # Build the cre stacked section # throughout cre stacked traces sorting -files = glob.glob('creStackedTrace-*.rsf') length = len(files) sortedFiles = [] diff --git a/experiments/fullInterpolationAndStackCDS/creStack.py b/experiments/fullInterpolationAndStackCDS/creStack.py deleted file mode 100644 index f997200..0000000 --- a/experiments/fullInterpolationAndStackCDS/creStack.py +++ /dev/null @@ -1,165 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# creStack (Python) -# -# Purpose: Recipe to the cre model, interpolation and Stacking. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Madagascar package -from rsf.proj import * - -def kirchoffModeling(filename='dataCube'): - ''' - Modeling function of a gaussian reflector - ''' - - # Modeling: Gaussian reflector in a velocity linear model - # velocity increases with depth with a 0.5 velocity gradient - Flow('gaussianReflector',None, - ''' - math d1=0.01 n1=2001 o1=-5 unit1=km label1=Offset - output="4-3*exp(-(x1-5)^2/9)" - ''') - - # Velocity Model - Flow('velocityModel','gaussianReflector', - ''' - window min1=0 max1=10 | - spray axis=1 n=451 d=0.01 o=0 label=Depth unit=km | - math output="1.5+0.5*x1+0.0*x2" - ''') - - Flow('reflectorDip','gaussianReflector','math output="2/3*(x1-5)*input" ') - - # Kirchoff Modeling - Flow(filename,'gaussianReflector reflectorDip', - ''' - kirmod cmp=y dip=${SOURCES[1]} - nh=161 dh=0.025 h0=0 - ns=401 ds=0.025 s0=0 - freq=10 dt=0.004 nt=1001 - vel=1.5 gradz=0.5 gradx=0.0 verb=y | - put d2=0.0125 label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s - ''') - -def pefInterpolation( - dataCube, - interpolated, - nm, - dm, - nt, - dt, - nhi=1 - ): - ''' - PEF interpolation of the data cube - :param dataCube: filename, Seismic data cube to interpolate - :param interpolated: filename, Interpolated seismic data cube - :param nm: integer, number of CMPs in the seismic data cube - :param dm: float, CMP sampling - :param nt: integer, number of time samples - :param dt: float, time sampling - :param nhi: integer, number of constant offsets gathers to interpolate - ''' - - # Divide CMP sampling - dm = dm/2 - - # Define mask file names using input filename - mask1 = dataCube+'-mask1' - mask = dataCube+'-mask' - aa = dataCube+'-aa' - bb = dataCube+'-bb' - a = dataCube+'-a' - b = dataCube+'-b' - zeroTraceGather = dataCube+'-zeroedGather' - mask0 = dataCube+'-mask0' - - - # Build a mask to interleave zero traces with original data traces - Flow(aa,None,'spike n1=%i d1=%g o1=0' %(nm,dm)) - Flow(bb,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask1,[bb, aa], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - - Flow(a,None,'spike n1=%i d1=%g o1=0' % (nm,dm)) - Flow(b,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask,[a, b], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - Flow(zeroTraceGather,b, - ''' - spray axis=2 n=%i d=%g | - transp | - put label2=Offset unit2=Km label1=Time unit1=s - ''' %(nt,dm)) - - # Data Mask with double of traces in CMP (half of CMP sampling) - # Keep the same Time and Offset original data sampling - Flow(mask0,mask, - ''' - spray axis=1 n=%i d=%g - ''' %(nt,dt)) - - totalPefIterations = 100 - totalInterpolationIterations = 20 - - offsetGathers = [] - for offsetGatherIndex in range(nhi): - - offsetGather = dataCube+"-offsetGather-%i" % offsetGatherIndex - resampledOffsetGather = dataCube+"-resampledGather-%i" % offsetGatherIndex - interpolatedOffsetGather = dataCube+"-interpolatedGather-%i" % offsetGatherIndex - pefCoeficients = dataCube+"-pefCoeficients-%i" % offsetGatherIndex - - Flow(offsetGather,dataCube, - ''' - window n2=1 f2=%i - ''' % (offsetGatherIndex)) - - Flow(resampledOffsetGather,[offsetGather,zeroTraceGather], - ''' - interleave axis=2 ${SOURCES[1]} - ''') - - # Calculate adaptive PEF coeficients - Flow(pefCoeficients,[resampledOffsetGather,mask0], - ''' - apef jump=2 a=10,2 rect1=50 rect2=2 niter=%g verb=y - maskin=${SOURCES[1]} - ''' % (totalPefIterations)) - - # Interpolation - Flow(interpolatedOffsetGather, [resampledOffsetGather,pefCoeficients,mask0,mask1], - ''' - miss4 exact=y filt=${SOURCES[1]} mask=${SOURCES[2]} niter=%g verb=y | - put d2=%g - ''' % (totalInterpolationIterations,dm)) - - offsetGathers.append(interpolatedOffsetGather) - - # Concatenate interpolated sections - Flow(interpolated,offsetGathers, - ''' - rcat axis=3 ${SOURCES[1:%d]} | - transp plane=23 - ''' % nhi) - - diff --git a/experiments/modelingAndPEFInterpolation/gaussianModel/SConstruct b/experiments/modelingAndPEFInterpolation/gaussianModel/SConstruct new file mode 100644 index 0000000..9615f7d --- /dev/null +++ b/experiments/modelingAndPEFInterpolation/gaussianModel/SConstruct @@ -0,0 +1,73 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- +# +# SConstruct (Madagascar Script) +# +# Purpose: Kirchhoff modeling of a gaussian reflector in a linear velocity model +# and PEF interpolation of the seismic data cube (Time x Offset x CMP) to increase +# CMP sampling. +# +# Site: https://dirack.github.io +# +# Version 2.0 +# +# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +# Madagascar package +from rsf.proj import * + +# Recipes +from kimodel import kirchoffModeling as kimod +from rsf.recipes.pefInterpolation import pefInterpolation as pef + +# Establish building dependencies +Depends('interpolatedDataCube.rsf','dataCube.rsf') +Depends('interpolatedDataCube2.rsf','interpolatedDataCube.rsf') +Depends('crsParameters-m0-0-t0-0.rsf','interpolatedDataCube.rsf') + +# Generate Gaussian reflector model and data cube +kimod() + +# PEF coeficients +a1=10 +a2=2 + +# PEF smooth radius +rect1=50 +rect2=2 + +# PEF interpolation +pef('dataCube', + 'interpolatedDataCube', + nm=401, + dm=0.025, + nt=1001, + dt=0.004, + nhi=21, + a1=a1, + a2=a2, + rect1=rect1, + rect2=rect2 +) + +# PEF interpolation +pef('interpolatedDataCube', + 'interpolatedDataCube2', + nm=802, + dm=0.0125, + nt=1001, + dt=0.004, + nhi=21, + a1=a1, + a2=a2, + rect1=rect1, + rect2=rect2 +) + +AddPostAction('interpolatedDataCube2.rsf',"./clean.sh") + +End() diff --git a/experiments/modelingAndPEFInterpolation/gaussianModel/clean.sh b/experiments/modelingAndPEFInterpolation/gaussianModel/clean.sh new file mode 100755 index 0000000..a47151a --- /dev/null +++ b/experiments/modelingAndPEFInterpolation/gaussianModel/clean.sh @@ -0,0 +1,23 @@ +#!/bin/bash +# +# clean.sh (Shell Script) +# +# Purpose: Clean temporary files after build. +# +# Site: https://dirack.github.io +# +# Version 1.0 +# +# Programer: Rodolfo A C Neves (Dirack) 19/09/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +TMP='gaussianReflector.rsf reflectorDip.rsf dataCube-offsetGather-* dataCube-interpolatedGather-* dataCube-pefCoeficients-* dataCube-resampledGather-* interpolatedDataCube-offsetGather-* interpolatedDataCube-interpolatedGather-* interpolatedDataCube-pefCoeficients-* interpolatedDataCube-resampledGather-* dataCube-a* dataCube-b* *-mask* interpolatedDataCube-a* interpolatedDataCube-b* *-zeroedGather*' + +for i in $TMP +do + $RSFROOT/bin/sfrm "$i" + #echo "$i" +done diff --git a/experiments/modelingAndPEFInterpolation/gaussianModel/kimodel.py b/experiments/modelingAndPEFInterpolation/gaussianModel/kimodel.py new file mode 100644 index 0000000..2255f78 --- /dev/null +++ b/experiments/modelingAndPEFInterpolation/gaussianModel/kimodel.py @@ -0,0 +1,57 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- +# +# kimodel.py (Madagascar Recipe) +# +# Purpose: Recipe to Kirchhoff modeling of a gaussian reflector. +# +# Important!: It should be called from a SConstruct +# +# Site: https://dirack.github.io +# +# Version 2.0 +# +# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +# Madagascar package +from rsf.proj import * + +def kirchoffModeling(filename='dataCube'): + ''' + Modeling function of a gaussian reflector + ''' + + # Modeling: Gaussian reflector in a velocity linear model + # velocity increases with depth with a 0.5 velocity gradient + Flow('gaussianReflector',None, + ''' + math d1=0.01 n1=2001 o1=-5 unit1=km label1=Offset + output="4-3*exp(-(x1-5)^2/9)" + ''') + + # Velocity Model + Flow('velocityModel','gaussianReflector', + ''' + window min1=0 max1=10 | + spray axis=1 n=451 d=0.01 o=0 label=Depth unit=km | + math output="1.5+0.5*x1+0.0*x2" + ''') + + Flow('reflectorDip','gaussianReflector','math output="2/3*(x1-5)*input" ') + + # Kirchoff Modeling + Flow(filename,'gaussianReflector reflectorDip', + ''' + kirmod cmp=y dip=${SOURCES[1]} + nh=161 dh=0.025 h0=0 + ns=401 ds=0.025 s0=0 + freq=10 dt=0.004 nt=1001 + vel=1.5 gradz=0.5 gradx=0.0 verb=y | + put d2=0.0125 label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s + ''') + + diff --git a/experiments/multiLayerModel/cds/SConscript b/experiments/modelingAndPEFInterpolation/multiLayerModel/SConstruct similarity index 57% rename from experiments/multiLayerModel/cds/SConscript rename to experiments/modelingAndPEFInterpolation/multiLayerModel/SConstruct index 4d8215f..f992703 100644 --- a/experiments/multiLayerModel/cds/SConscript +++ b/experiments/modelingAndPEFInterpolation/multiLayerModel/SConstruct @@ -1,15 +1,13 @@ #!/usr/bin/env python # -*- coding: utf-8 -*- # -# SConscript (Madagascar Script) +# SConstruct (Madagascar Script) # -# Purpose: Generate model and data cubes to CRE stacking process. -# This script generates model, original and PEF interpolated data cubes. -# It depends of the 'creStack.py' script that defines functions. +# Purpose: Kirchhoff Newton modeling of a multi layer velocity model. # -# Site: http://www.dirackslounge.online +# Site: https://dirack.github.io # -# Version 1.0 +# Version 2.0 # # Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 # @@ -20,10 +18,15 @@ # Madagascar package from rsf.proj import * -# CRE recipe -from madagascarRecipes.pefInterpolation import pefInterpolation as pefin -from madagascarRecipes.kimodel import multiLayerModel as mlmod -from madagascarRecipes.kimodel import kirchoffNewtonModeling as kinewmod +# Recipes +from rsf.recipes.pefInterpolation import pefInterpolation as pef +from rsf.recipes.kimodel import multiLayerModel as mlmod +from rsf.recipes.kimodel import kirchhoffNewtonModeling as kinewmod + +# Establish building dependencies +Depends('interpolatedDataCube.rsf','multiLayerDataCube.rsf') +Depends('interpolatedDataCube2.rsf','interpolatedDataCube.rsf') +Depends('crsParameters-m0-0-t0-0.rsf','interpolatedDataCube2.rsf') xmax = 6.0 zmax = 2.0 @@ -44,10 +47,17 @@ mlmod(interfaces='interfaces', layers=layers, velocities=velocities) +# Generate a modeled data cube kinewmod(reflectors='interfaces', reflectorsDip='interfacesDip', filename='multiLayerDataCube', - velocities=velocities) + velocities=velocities, + nt=1001, + dt=0.004, + ns=201, + ds=0.025, + nh=161, + dh=0.025) # PEF interpolation of the data cube a1=10 @@ -55,13 +65,13 @@ a2=2 rect1=50 rect2=2 -pefin(dataCube='multiLayerDataCube', +pef(dataCube='multiLayerDataCube', interpolated='interpolatedDataCube', nm=201, dm=0.025, nt=1001, dt=0.004, - nhi=161, + nhi=21, a1=a1, a2=a2, rect1=rect1, @@ -69,16 +79,18 @@ pefin(dataCube='multiLayerDataCube', # Do the PEF interpolation one more time # to increase the CMP sampling -pefin('interpolatedDataCube', +pef('interpolatedDataCube', 'interpolatedDataCube2', nm=402, dm=0.0125, nt=1001, dt=0.004, - nhi=161, + nhi=21, a1=a1, a2=a2, rect1=rect1, rect2=rect2) +AddPostAction('interpolatedDataCube2.rsf',"./clean.sh") + End() diff --git a/experiments/modelingAndPEFInterpolation/multiLayerModel/clean.sh b/experiments/modelingAndPEFInterpolation/multiLayerModel/clean.sh new file mode 100755 index 0000000..5a5e86c --- /dev/null +++ b/experiments/modelingAndPEFInterpolation/multiLayerModel/clean.sh @@ -0,0 +1,25 @@ +#!/bin/bash +# +# clean.sh (Shell Script) +# +# Purpose: Clean temporary files after build. +# +# Site: https://dirack.github.io +# +# Version 1.0 +# +# Programer: Rodolfo A C Neves (Dirack) 19/09/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +TMP='interfaces.rsf interfacesDip.rsf mod1.rsf layers.rsf multiLayerDataCube-offsetGather-* multiLayerDataCube-interpolatedGather-* multiLayerDataCube-pefCoeficients-* multiLayerDataCube-resampledGather-* interpolatedDataCube-offsetGather-* interpolatedDataCube-interpolatedGather-* interpolatedDataCube-pefCoeficients-* interpolatedDataCube-resampledGather-* multiLayerDataCube-a* multiLayerDataCube-b* *-mask* interpolatedDataCube-a* interpolatedDataCube-b* *-zeroedGather*' + +for i in $TMP +do + $RSFROOT/bin/sfrm "$i" + #echo "$i" +done + +rm layers.asc diff --git a/experiments/modelingAndPEFInterpolation/multiLayerModel/kimodel.py b/experiments/modelingAndPEFInterpolation/multiLayerModel/kimodel.py new file mode 100644 index 0000000..2255f78 --- /dev/null +++ b/experiments/modelingAndPEFInterpolation/multiLayerModel/kimodel.py @@ -0,0 +1,57 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- +# +# kimodel.py (Madagascar Recipe) +# +# Purpose: Recipe to Kirchhoff modeling of a gaussian reflector. +# +# Important!: It should be called from a SConstruct +# +# Site: https://dirack.github.io +# +# Version 2.0 +# +# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +# Madagascar package +from rsf.proj import * + +def kirchoffModeling(filename='dataCube'): + ''' + Modeling function of a gaussian reflector + ''' + + # Modeling: Gaussian reflector in a velocity linear model + # velocity increases with depth with a 0.5 velocity gradient + Flow('gaussianReflector',None, + ''' + math d1=0.01 n1=2001 o1=-5 unit1=km label1=Offset + output="4-3*exp(-(x1-5)^2/9)" + ''') + + # Velocity Model + Flow('velocityModel','gaussianReflector', + ''' + window min1=0 max1=10 | + spray axis=1 n=451 d=0.01 o=0 label=Depth unit=km | + math output="1.5+0.5*x1+0.0*x2" + ''') + + Flow('reflectorDip','gaussianReflector','math output="2/3*(x1-5)*input" ') + + # Kirchoff Modeling + Flow(filename,'gaussianReflector reflectorDip', + ''' + kirmod cmp=y dip=${SOURCES[1]} + nh=161 dh=0.025 h0=0 + ns=401 ds=0.025 s0=0 + freq=10 dt=0.004 nt=1001 + vel=1.5 gradz=0.5 gradx=0.0 verb=y | + put d2=0.0125 label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s + ''') + + diff --git a/experiments/modelingAndPEFInterpolation/test_modelBuild.py b/experiments/modelingAndPEFInterpolation/test_modelBuild.py new file mode 100644 index 0000000..9ef0b4e --- /dev/null +++ b/experiments/modelingAndPEFInterpolation/test_modelBuild.py @@ -0,0 +1,58 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- +# +# test_modelBuild.py (Python) +# +# Purpose: Test model building with SCons. +# +# Rum with: python -m unittest test_modelBuild +# +# Site: https://dirack.github.io +# +# Version 2.0 +# +# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +import unittest +import subprocess + +class TestScons(unittest.TestCase): + ''' + Unit tests of the model building + ''' + + def test_gaussianReflectorModeling(self): + ''' + Test the building of the gaussian reflector in + a linear velocity model + ''' + result = subprocess.call('cd gaussianModel && scons dataCube.rsf',shell=True) + self.assertEqual(result,0,'FAILED: gaussian reflector modeling failed!') + + + def test_gaussianReflectorPefInterpolation(self): + ''' + Test of the PEF interpolation step of the seismic data cube + modeled with the gaussian reflector linear velocity model. + ''' + result = subprocess.call('cd gaussianModel && scons -j3',shell=True) + self.assertEqual(result,0,'FAILED: gaussian reflector PEF interpolation failed!') + + def test_multiLayerModeling(self): + ''' + Test of the building of the multi layer model + ''' + result = subprocess.call('cd multiLayerModel && scons multiLayerDataCube.rsf',shell=True) + self.assertEqual(result,0,'FAILED: multi Layer modeling failed!') + + def test_multiLayerPefInterpolation(self): + ''' + Test of the PEF interpolation step of the seimic data cube + modeled with the multi layer model + ''' + result = subprocess.call('cd multiLayerModel && scons -j3',shell=True) + self.assertEqual(result,0,'FAILED: multi layer PEF interpolation failed!') diff --git a/experiments/multiLayerModel/cds/SConscript2 b/experiments/multiLayerModel/cds/SConscript2 deleted file mode 100644 index 0c63c46..0000000 --- a/experiments/multiLayerModel/cds/SConscript2 +++ /dev/null @@ -1,61 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# SConscript (Madagascar Script) -# -# Purpose: Contatenate CRE stacked traces. -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Madagascar package -from rsf.proj import * - -# Error report function -import atexit -from madagascarRecipes.errorReport import print_build_failures - -import glob - -# Build the cre stacked section -# throughout cre stacked traces sorting -files = glob.glob('creStackedTrace-*.rsf') -length = len(files) - -sortedFiles = [] -for i in range(length): - string = 'creStackedTrace-m0-%i.rsf' % i - sortedFiles.append(string) - -Flow('stackedSection',sortedFiles, - ''' - rcat axis=2 ${SOURCES[1:%d]} --out=stdout - ''' % len(files)) - -Result('stackedSection','grey title="Stacked Section"') - -Flow('filtStackedSection','stackedSection', - ''' - bandpass fhi=20 --out=stdout - ''') - -Result('filtStackedSection','grey title="Filtered Stacked Section"') - -# Use aliases to split building -Alias('stack',['stackedSection.rsf','filtStackedSection.rsf']) - -# Show error message if fail - -message = ''' -SCosncript Failed build -''' - -atexit.register(print_build_failures,message) -End() diff --git a/experiments/multiLayerModel/cds/SConstruct b/experiments/multiLayerModel/cds/SConstruct index 582bf15..2ddea89 100644 --- a/experiments/multiLayerModel/cds/SConstruct +++ b/experiments/multiLayerModel/cds/SConstruct @@ -3,12 +3,16 @@ # # SConstruct (Madagascar Script) # -# Purpose: Build the interpolation process twice to improve CMP sampling -# and CRE stacking results (Use CDS condition). +# Purpose: Common Reflection Element (CRE) stack applied to a seismic data cube +# from Kirchhoff Newton modeling of a multi layer velocity model. # -# Site: http://www.dirackslounge.online +# Important!: This experiment uses non-hyperbolic Common Reflection Surface (CRS) +# traveltime approximation as stacking curve, with CDS condition (RNIP=RN) applied. +# Instead of conventional CRE traveltime approximation. +# +# Site: https://dirack.github.io # -# Version 1.0 +# Version 2.0 # # Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 # @@ -22,14 +26,6 @@ from rsf.proj import * # Import glob python library import glob -# Error report function -import atexit -from madagascarRecipes.errorReport import print_build_failures - -# Call SConscript to generate input: original and interpolated data cubes -SConscript('SConscript') - -# CRE stacking # It uses Very Fast Simulated Aneeling and non hyperbolic CRS # to get zero offset CRS parameters (RN, RNIP and BETA) from data cube v0 = 1.5 @@ -37,9 +33,12 @@ ot0 = 0.3 dt0 = 0.008 nt0 = 331 om0 = 2 -dm0 = 0.025 -nm0 = 108 -dataCube='multiLayerDataCube' +dm0 = 0.25 +nm0 = 8 +inputDir = '../../modelingAndPEFInterpolation/multiLayerModel/' +dataCube=inputDir+'multiLayerDataCube.rsf' +interpolatedDataCube=inputDir+'interpolatedDataCube2.rsf' +files = [] for i in range(nm0): @@ -64,29 +63,29 @@ for i in range(nm0): creGatherPlot = 'crePlot-m0-%g-t0-%g' % (i,j) # Very Fast Simulated Aneelling Global Optimization (VFSA) - Flow(crsParameters,'interpolatedDataCube2', + Flow(crsParameters,dataCube, ''' vfsacrsnh m0=%g v0=%g t0=%g verb=y repeat=3 ''' % (m0,v0,t0)) # Calculate CRE trajectory - Flow(creMhCoordinates,['interpolatedDataCube2',crsParameters], + Flow(creMhCoordinates,[interpolatedDataCube,crsParameters], ''' cretrajec verb=y m0=%g param=${SOURCES[1]} | put unit1="Offset" label1="Km" ''' % (m0)) #Get CRE Gather from interpolated Data Cube - Flow([creGather,creMcoordinate],['interpolatedDataCube2',creMhCoordinates], + Flow([creGather,creMcoordinate],[interpolatedDataCube,creMhCoordinates], ''' - getcregather aperture=50 verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} | + getcregather aperture=50 verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} aperture=20 | put label1="Time" unit1="s" label2="Offset" unit2="km" ''') # Calculate CRE traveltime curve t(m,h) Flow(creTimeCurve,[creMcoordinate, crsParameters], ''' - getcretimecurve cds=y param=${SOURCES[1]} t0=%g m0=%g v0=%g verb=y | + getcretimecurve param=${SOURCES[1]} cds=y t0=%g m0=%g v0=%g verb=y | put label1="Offset" unit1="Km" ''' % (t0,m0,v0)) @@ -133,13 +132,19 @@ for i in range(nm0): # CRE stacking Flow(creStackedTrace,[creGatherTrace,creTimeCurveTrace], ''' - crestack timeCurves=${SOURCES[1]} verb=y aperture=10 | + crestack timeCurves=${SOURCES[1]} verb=y | put label1=t0 unit1=s label2=m0 unit2=Km --out=stdout ''') + files.append(creStackedTrace) + +# Stacked section depends of the last CRE stacked trace +Depends('stackedSection.rsf',[dataCube, + 'creStackedTrace-m0-'+str(nm0-1)+'.rsf', + 'creGatherTrace-m0-'+str(nm0-1)+'.rsf']) + # Build the cre stacked section # throughout cre stacked traces sorting -files = glob.glob('creStackedTrace-*.rsf') length = len(files) sortedFiles = [] @@ -157,11 +162,4 @@ Flow('filtStackedSection','stackedSection', bandpass fhi=20 --out=stdout ''') -# Show error message if fail -message = ''' -SConstruct Failed build when m0={m0} and t0={t0} -'''.format(m0=(i*dm0+om0),t0=(j*dt0+ot0)) - -atexit.register(print_build_failures,message) - End() diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/LICENSE b/experiments/multiLayerModel/cds/madagascarRecipes/LICENSE deleted file mode 100644 index f288702..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/LICENSE +++ /dev/null @@ -1,674 +0,0 @@ - GNU GENERAL PUBLIC LICENSE - Version 3, 29 June 2007 - - Copyright (C) 2007 Free Software Foundation, Inc. - Everyone is permitted to copy and distribute verbatim copies - of this license document, but changing it is not allowed. - - Preamble - - The GNU General Public License is a free, copyleft license for -software and other kinds of works. - - The licenses for most software and other practical works are designed -to take away your freedom to share and change the works. By contrast, -the GNU General Public License is intended to guarantee your freedom to -share and change all versions of a program--to make sure it remains free -software for all its users. We, the Free Software Foundation, use the -GNU General Public License for most of our software; it applies also to -any other work released this way by its authors. You can apply it to -your programs, too. - - When we speak of free software, we are referring to freedom, not -price. Our General Public Licenses are designed to make sure that you -have the freedom to distribute copies of free software (and charge for -them if you wish), that you receive source code or can get it if you -want it, that you can change the software or use pieces of it in new -free programs, and that you know you can do these things. - - To protect your rights, we need to prevent others from denying you -these rights or asking you to surrender the rights. Therefore, you have -certain responsibilities if you distribute copies of the software, or if -you modify it: responsibilities to respect the freedom of others. - - For example, if you distribute copies of such a program, whether -gratis or for a fee, you must pass on to the recipients the same -freedoms that you received. You must make sure that they, too, receive -or can get the source code. And you must show them these terms so they -know their rights. - - Developers that use the GNU GPL protect your rights with two steps: -(1) assert copyright on the software, and (2) offer you this License -giving you legal permission to copy, distribute and/or modify it. - - For the developers' and authors' protection, the GPL clearly explains -that there is no warranty for this free software. For both users' and -authors' sake, the GPL requires that modified versions be marked as -changed, so that their problems will not be attributed erroneously to -authors of previous versions. - - Some devices are designed to deny users access to install or run -modified versions of the software inside them, although the manufacturer -can do so. This is fundamentally incompatible with the aim of -protecting users' freedom to change the software. The systematic -pattern of such abuse occurs in the area of products for individuals to -use, which is precisely where it is most unacceptable. Therefore, we -have designed this version of the GPL to prohibit the practice for those -products. If such problems arise substantially in other domains, we -stand ready to extend this provision to those domains in future versions -of the GPL, as needed to protect the freedom of users. - - Finally, every program is threatened constantly by software patents. -States should not allow patents to restrict development and use of -software on general-purpose computers, but in those that do, we wish to -avoid the special danger that patents applied to a free program could -make it effectively proprietary. To prevent this, the GPL assures that -patents cannot be used to render the program non-free. - - The precise terms and conditions for copying, distribution and -modification follow. - - TERMS AND CONDITIONS - - 0. Definitions. - - "This License" refers to version 3 of the GNU General Public License. - - "Copyright" also means copyright-like laws that apply to other kinds of -works, such as semiconductor masks. - - "The Program" refers to any copyrightable work licensed under this -License. Each licensee is addressed as "you". "Licensees" and -"recipients" may be individuals or organizations. - - To "modify" a work means to copy from or adapt all or part of the work -in a fashion requiring copyright permission, other than the making of an -exact copy. The resulting work is called a "modified version" of the -earlier work or a work "based on" the earlier work. - - A "covered work" means either the unmodified Program or a work based -on the Program. - - To "propagate" a work means to do anything with it that, without -permission, would make you directly or secondarily liable for -infringement under applicable copyright law, except executing it on a -computer or modifying a private copy. Propagation includes copying, -distribution (with or without modification), making available to the -public, and in some countries other activities as well. - - To "convey" a work means any kind of propagation that enables other -parties to make or receive copies. Mere interaction with a user through -a computer network, with no transfer of a copy, is not conveying. - - An interactive user interface displays "Appropriate Legal Notices" -to the extent that it includes a convenient and prominently visible -feature that (1) displays an appropriate copyright notice, and (2) -tells the user that there is no warranty for the work (except to the -extent that warranties are provided), that licensees may convey the -work under this License, and how to view a copy of this License. If -the interface presents a list of user commands or options, such as a -menu, a prominent item in the list meets this criterion. - - 1. Source Code. - - The "source code" for a work means the preferred form of the work -for making modifications to it. "Object code" means any non-source -form of a work. - - A "Standard Interface" means an interface that either is an official -standard defined by a recognized standards body, or, in the case of -interfaces specified for a particular programming language, one that -is widely used among developers working in that language. - - The "System Libraries" of an executable work include anything, other -than the work as a whole, that (a) is included in the normal form of -packaging a Major Component, but which is not part of that Major -Component, and (b) serves only to enable use of the work with that -Major Component, or to implement a Standard Interface for which an -implementation is available to the public in source code form. A -"Major Component", in this context, means a major essential component -(kernel, window system, and so on) of the specific operating system -(if any) on which the executable work runs, or a compiler used to -produce the work, or an object code interpreter used to run it. - - The "Corresponding Source" for a work in object code form means all -the source code needed to generate, install, and (for an executable -work) run the object code and to modify the work, including scripts to -control those activities. However, it does not include the work's -System Libraries, or general-purpose tools or generally available free -programs which are used unmodified in performing those activities but -which are not part of the work. For example, Corresponding Source -includes interface definition files associated with source files for -the work, and the source code for shared libraries and dynamically -linked subprograms that the work is specifically designed to require, -such as by intimate data communication or control flow between those -subprograms and other parts of the work. - - The Corresponding Source need not include anything that users -can regenerate automatically from other parts of the Corresponding -Source. - - The Corresponding Source for a work in source code form is that -same work. - - 2. Basic Permissions. - - All rights granted under this License are granted for the term of -copyright on the Program, and are irrevocable provided the stated -conditions are met. This License explicitly affirms your unlimited -permission to run the unmodified Program. The output from running a -covered work is covered by this License only if the output, given its -content, constitutes a covered work. This License acknowledges your -rights of fair use or other equivalent, as provided by copyright law. - - You may make, run and propagate covered works that you do not -convey, without conditions so long as your license otherwise remains -in force. You may convey covered works to others for the sole purpose -of having them make modifications exclusively for you, or provide you -with facilities for running those works, provided that you comply with -the terms of this License in conveying all material for which you do -not control copyright. Those thus making or running the covered works -for you must do so exclusively on your behalf, under your direction -and control, on terms that prohibit them from making any copies of -your copyrighted material outside their relationship with you. - - Conveying under any other circumstances is permitted solely under -the conditions stated below. Sublicensing is not allowed; section 10 -makes it unnecessary. - - 3. Protecting Users' Legal Rights From Anti-Circumvention Law. - - No covered work shall be deemed part of an effective technological -measure under any applicable law fulfilling obligations under article -11 of the WIPO copyright treaty adopted on 20 December 1996, or -similar laws prohibiting or restricting circumvention of such -measures. - - When you convey a covered work, you waive any legal power to forbid -circumvention of technological measures to the extent such circumvention -is effected by exercising rights under this License with respect to -the covered work, and you disclaim any intention to limit operation or -modification of the work as a means of enforcing, against the work's -users, your or third parties' legal rights to forbid circumvention of -technological measures. - - 4. Conveying Verbatim Copies. - - You may convey verbatim copies of the Program's source code as you -receive it, in any medium, provided that you conspicuously and -appropriately publish on each copy an appropriate copyright notice; -keep intact all notices stating that this License and any -non-permissive terms added in accord with section 7 apply to the code; -keep intact all notices of the absence of any warranty; and give all -recipients a copy of this License along with the Program. - - You may charge any price or no price for each copy that you convey, -and you may offer support or warranty protection for a fee. - - 5. Conveying Modified Source Versions. - - You may convey a work based on the Program, or the modifications to -produce it from the Program, in the form of source code under the -terms of section 4, provided that you also meet all of these conditions: - - a) The work must carry prominent notices stating that you modified - it, and giving a relevant date. - - b) The work must carry prominent notices stating that it is - released under this License and any conditions added under section - 7. This requirement modifies the requirement in section 4 to - "keep intact all notices". - - c) You must license the entire work, as a whole, under this - License to anyone who comes into possession of a copy. This - License will therefore apply, along with any applicable section 7 - additional terms, to the whole of the work, and all its parts, - regardless of how they are packaged. This License gives no - permission to license the work in any other way, but it does not - invalidate such permission if you have separately received it. - - d) If the work has interactive user interfaces, each must display - Appropriate Legal Notices; however, if the Program has interactive - interfaces that do not display Appropriate Legal Notices, your - work need not make them do so. - - A compilation of a covered work with other separate and independent -works, which are not by their nature extensions of the covered work, -and which are not combined with it such as to form a larger program, -in or on a volume of a storage or distribution medium, is called an -"aggregate" if the compilation and its resulting copyright are not -used to limit the access or legal rights of the compilation's users -beyond what the individual works permit. Inclusion of a covered work -in an aggregate does not cause this License to apply to the other -parts of the aggregate. - - 6. Conveying Non-Source Forms. - - You may convey a covered work in object code form under the terms -of sections 4 and 5, provided that you also convey the -machine-readable Corresponding Source under the terms of this License, -in one of these ways: - - a) Convey the object code in, or embodied in, a physical product - (including a physical distribution medium), accompanied by the - Corresponding Source fixed on a durable physical medium - customarily used for software interchange. - - b) Convey the object code in, or embodied in, a physical product - (including a physical distribution medium), accompanied by a - written offer, valid for at least three years and valid for as - long as you offer spare parts or customer support for that product - model, to give anyone who possesses the object code either (1) a - copy of the Corresponding Source for all the software in the - product that is covered by this License, on a durable physical - medium customarily used for software interchange, for a price no - more than your reasonable cost of physically performing this - conveying of source, or (2) access to copy the - Corresponding Source from a network server at no charge. - - c) Convey individual copies of the object code with a copy of the - written offer to provide the Corresponding Source. This - alternative is allowed only occasionally and noncommercially, and - only if you received the object code with such an offer, in accord - with subsection 6b. - - d) Convey the object code by offering access from a designated - place (gratis or for a charge), and offer equivalent access to the - Corresponding Source in the same way through the same place at no - further charge. You need not require recipients to copy the - Corresponding Source along with the object code. If the place to - copy the object code is a network server, the Corresponding Source - may be on a different server (operated by you or a third party) - that supports equivalent copying facilities, provided you maintain - clear directions next to the object code saying where to find the - Corresponding Source. Regardless of what server hosts the - Corresponding Source, you remain obligated to ensure that it is - available for as long as needed to satisfy these requirements. - - e) Convey the object code using peer-to-peer transmission, provided - you inform other peers where the object code and Corresponding - Source of the work are being offered to the general public at no - charge under subsection 6d. - - A separable portion of the object code, whose source code is excluded -from the Corresponding Source as a System Library, need not be -included in conveying the object code work. - - A "User Product" is either (1) a "consumer product", which means any -tangible personal property which is normally used for personal, family, -or household purposes, or (2) anything designed or sold for incorporation -into a dwelling. In determining whether a product is a consumer product, -doubtful cases shall be resolved in favor of coverage. For a particular -product received by a particular user, "normally used" refers to a -typical or common use of that class of product, regardless of the status -of the particular user or of the way in which the particular user -actually uses, or expects or is expected to use, the product. A product -is a consumer product regardless of whether the product has substantial -commercial, industrial or non-consumer uses, unless such uses represent -the only significant mode of use of the product. - - "Installation Information" for a User Product means any methods, -procedures, authorization keys, or other information required to install -and execute modified versions of a covered work in that User Product from -a modified version of its Corresponding Source. The information must -suffice to ensure that the continued functioning of the modified object -code is in no case prevented or interfered with solely because -modification has been made. - - If you convey an object code work under this section in, or with, or -specifically for use in, a User Product, and the conveying occurs as -part of a transaction in which the right of possession and use of the -User Product is transferred to the recipient in perpetuity or for a -fixed term (regardless of how the transaction is characterized), the -Corresponding Source conveyed under this section must be accompanied -by the Installation Information. But this requirement does not apply -if neither you nor any third party retains the ability to install -modified object code on the User Product (for example, the work has -been installed in ROM). - - The requirement to provide Installation Information does not include a -requirement to continue to provide support service, warranty, or updates -for a work that has been modified or installed by the recipient, or for -the User Product in which it has been modified or installed. Access to a -network may be denied when the modification itself materially and -adversely affects the operation of the network or violates the rules and -protocols for communication across the network. - - Corresponding Source conveyed, and Installation Information provided, -in accord with this section must be in a format that is publicly -documented (and with an implementation available to the public in -source code form), and must require no special password or key for -unpacking, reading or copying. - - 7. Additional Terms. - - "Additional permissions" are terms that supplement the terms of this -License by making exceptions from one or more of its conditions. -Additional permissions that are applicable to the entire Program shall -be treated as though they were included in this License, to the extent -that they are valid under applicable law. If additional permissions -apply only to part of the Program, that part may be used separately -under those permissions, but the entire Program remains governed by -this License without regard to the additional permissions. - - When you convey a copy of a covered work, you may at your option -remove any additional permissions from that copy, or from any part of -it. (Additional permissions may be written to require their own -removal in certain cases when you modify the work.) You may place -additional permissions on material, added by you to a covered work, -for which you have or can give appropriate copyright permission. - - Notwithstanding any other provision of this License, for material you -add to a covered work, you may (if authorized by the copyright holders of -that material) supplement the terms of this License with terms: - - a) Disclaiming warranty or limiting liability differently from the - terms of sections 15 and 16 of this License; or - - b) Requiring preservation of specified reasonable legal notices or - author attributions in that material or in the Appropriate Legal - Notices displayed by works containing it; or - - c) Prohibiting misrepresentation of the origin of that material, or - requiring that modified versions of such material be marked in - reasonable ways as different from the original version; or - - d) Limiting the use for publicity purposes of names of licensors or - authors of the material; or - - e) Declining to grant rights under trademark law for use of some - trade names, trademarks, or service marks; or - - f) Requiring indemnification of licensors and authors of that - material by anyone who conveys the material (or modified versions of - it) with contractual assumptions of liability to the recipient, for - any liability that these contractual assumptions directly impose on - those licensors and authors. - - All other non-permissive additional terms are considered "further -restrictions" within the meaning of section 10. If the Program as you -received it, or any part of it, contains a notice stating that it is -governed by this License along with a term that is a further -restriction, you may remove that term. If a license document contains -a further restriction but permits relicensing or conveying under this -License, you may add to a covered work material governed by the terms -of that license document, provided that the further restriction does -not survive such relicensing or conveying. - - If you add terms to a covered work in accord with this section, you -must place, in the relevant source files, a statement of the -additional terms that apply to those files, or a notice indicating -where to find the applicable terms. - - Additional terms, permissive or non-permissive, may be stated in the -form of a separately written license, or stated as exceptions; -the above requirements apply either way. - - 8. Termination. - - You may not propagate or modify a covered work except as expressly -provided under this License. Any attempt otherwise to propagate or -modify it is void, and will automatically terminate your rights under -this License (including any patent licenses granted under the third -paragraph of section 11). - - However, if you cease all violation of this License, then your -license from a particular copyright holder is reinstated (a) -provisionally, unless and until the copyright holder explicitly and -finally terminates your license, and (b) permanently, if the copyright -holder fails to notify you of the violation by some reasonable means -prior to 60 days after the cessation. - - Moreover, your license from a particular copyright holder is -reinstated permanently if the copyright holder notifies you of the -violation by some reasonable means, this is the first time you have -received notice of violation of this License (for any work) from that -copyright holder, and you cure the violation prior to 30 days after -your receipt of the notice. - - Termination of your rights under this section does not terminate the -licenses of parties who have received copies or rights from you under -this License. If your rights have been terminated and not permanently -reinstated, you do not qualify to receive new licenses for the same -material under section 10. - - 9. Acceptance Not Required for Having Copies. - - You are not required to accept this License in order to receive or -run a copy of the Program. Ancillary propagation of a covered work -occurring solely as a consequence of using peer-to-peer transmission -to receive a copy likewise does not require acceptance. However, -nothing other than this License grants you permission to propagate or -modify any covered work. These actions infringe copyright if you do -not accept this License. Therefore, by modifying or propagating a -covered work, you indicate your acceptance of this License to do so. - - 10. Automatic Licensing of Downstream Recipients. - - Each time you convey a covered work, the recipient automatically -receives a license from the original licensors, to run, modify and -propagate that work, subject to this License. You are not responsible -for enforcing compliance by third parties with this License. - - An "entity transaction" is a transaction transferring control of an -organization, or substantially all assets of one, or subdividing an -organization, or merging organizations. If propagation of a covered -work results from an entity transaction, each party to that -transaction who receives a copy of the work also receives whatever -licenses to the work the party's predecessor in interest had or could -give under the previous paragraph, plus a right to possession of the -Corresponding Source of the work from the predecessor in interest, if -the predecessor has it or can get it with reasonable efforts. - - You may not impose any further restrictions on the exercise of the -rights granted or affirmed under this License. For example, you may -not impose a license fee, royalty, or other charge for exercise of -rights granted under this License, and you may not initiate litigation -(including a cross-claim or counterclaim in a lawsuit) alleging that -any patent claim is infringed by making, using, selling, offering for -sale, or importing the Program or any portion of it. - - 11. Patents. - - A "contributor" is a copyright holder who authorizes use under this -License of the Program or a work on which the Program is based. The -work thus licensed is called the contributor's "contributor version". - - A contributor's "essential patent claims" are all patent claims -owned or controlled by the contributor, whether already acquired or -hereafter acquired, that would be infringed by some manner, permitted -by this License, of making, using, or selling its contributor version, -but do not include claims that would be infringed only as a -consequence of further modification of the contributor version. For -purposes of this definition, "control" includes the right to grant -patent sublicenses in a manner consistent with the requirements of -this License. - - Each contributor grants you a non-exclusive, worldwide, royalty-free -patent license under the contributor's essential patent claims, to -make, use, sell, offer for sale, import and otherwise run, modify and -propagate the contents of its contributor version. - - In the following three paragraphs, a "patent license" is any express -agreement or commitment, however denominated, not to enforce a patent -(such as an express permission to practice a patent or covenant not to -sue for patent infringement). To "grant" such a patent license to a -party means to make such an agreement or commitment not to enforce a -patent against the party. - - If you convey a covered work, knowingly relying on a patent license, -and the Corresponding Source of the work is not available for anyone -to copy, free of charge and under the terms of this License, through a -publicly available network server or other readily accessible means, -then you must either (1) cause the Corresponding Source to be so -available, or (2) arrange to deprive yourself of the benefit of the -patent license for this particular work, or (3) arrange, in a manner -consistent with the requirements of this License, to extend the patent -license to downstream recipients. "Knowingly relying" means you have -actual knowledge that, but for the patent license, your conveying the -covered work in a country, or your recipient's use of the covered work -in a country, would infringe one or more identifiable patents in that -country that you have reason to believe are valid. - - If, pursuant to or in connection with a single transaction or -arrangement, you convey, or propagate by procuring conveyance of, a -covered work, and grant a patent license to some of the parties -receiving the covered work authorizing them to use, propagate, modify -or convey a specific copy of the covered work, then the patent license -you grant is automatically extended to all recipients of the covered -work and works based on it. - - A patent license is "discriminatory" if it does not include within -the scope of its coverage, prohibits the exercise of, or is -conditioned on the non-exercise of one or more of the rights that are -specifically granted under this License. You may not convey a covered -work if you are a party to an arrangement with a third party that is -in the business of distributing software, under which you make payment -to the third party based on the extent of your activity of conveying -the work, and under which the third party grants, to any of the -parties who would receive the covered work from you, a discriminatory -patent license (a) in connection with copies of the covered work -conveyed by you (or copies made from those copies), or (b) primarily -for and in connection with specific products or compilations that -contain the covered work, unless you entered into that arrangement, -or that patent license was granted, prior to 28 March 2007. - - Nothing in this License shall be construed as excluding or limiting -any implied license or other defenses to infringement that may -otherwise be available to you under applicable patent law. - - 12. No Surrender of Others' Freedom. - - If conditions are imposed on you (whether by court order, agreement or -otherwise) that contradict the conditions of this License, they do not -excuse you from the conditions of this License. If you cannot convey a -covered work so as to satisfy simultaneously your obligations under this -License and any other pertinent obligations, then as a consequence you may -not convey it at all. For example, if you agree to terms that obligate you -to collect a royalty for further conveying from those to whom you convey -the Program, the only way you could satisfy both those terms and this -License would be to refrain entirely from conveying the Program. - - 13. Use with the GNU Affero General Public License. - - Notwithstanding any other provision of this License, you have -permission to link or combine any covered work with a work licensed -under version 3 of the GNU Affero General Public License into a single -combined work, and to convey the resulting work. The terms of this -License will continue to apply to the part which is the covered work, -but the special requirements of the GNU Affero General Public License, -section 13, concerning interaction through a network will apply to the -combination as such. - - 14. Revised Versions of this License. - - The Free Software Foundation may publish revised and/or new versions of -the GNU General Public License from time to time. Such new versions will -be similar in spirit to the present version, but may differ in detail to -address new problems or concerns. - - Each version is given a distinguishing version number. If the -Program specifies that a certain numbered version of the GNU General -Public License "or any later version" applies to it, you have the -option of following the terms and conditions either of that numbered -version or of any later version published by the Free Software -Foundation. If the Program does not specify a version number of the -GNU General Public License, you may choose any version ever published -by the Free Software Foundation. - - If the Program specifies that a proxy can decide which future -versions of the GNU General Public License can be used, that proxy's -public statement of acceptance of a version permanently authorizes you -to choose that version for the Program. - - Later license versions may give you additional or different -permissions. However, no additional obligations are imposed on any -author or copyright holder as a result of your choosing to follow a -later version. - - 15. Disclaimer of Warranty. - - THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY -APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT -HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY -OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, -THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR -PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM -IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF -ALL NECESSARY SERVICING, REPAIR OR CORRECTION. - - 16. Limitation of Liability. - - IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING -WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MODIFIES AND/OR CONVEYS -THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, INCLUDING ANY -GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE -USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED TO LOSS OF -DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY YOU OR THIRD -PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER PROGRAMS), -EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE POSSIBILITY OF -SUCH DAMAGES. - - 17. Interpretation of Sections 15 and 16. - - If the disclaimer of warranty and limitation of liability provided -above cannot be given local legal effect according to their terms, -reviewing courts shall apply local law that most closely approximates -an absolute waiver of all civil liability in connection with the -Program, unless a warranty or assumption of liability accompanies a -copy of the Program in return for a fee. - - END OF TERMS AND CONDITIONS - - How to Apply These Terms to Your New Programs - - If you develop a new program, and you want it to be of the greatest -possible use to the public, the best way to achieve this is to make it -free software which everyone can redistribute and change under these terms. - - To do so, attach the following notices to the program. It is safest -to attach them to the start of each source file to most effectively -state the exclusion of warranty; and each file should have at least -the "copyright" line and a pointer to where the full notice is found. - - - Copyright (C) - - This program is free software: you can redistribute it and/or modify - it under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - This program is distributed in the hope that it will be useful, - but WITHOUT ANY WARRANTY; without even the implied warranty of - MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the - GNU General Public License for more details. - - You should have received a copy of the GNU General Public License - along with this program. If not, see . - -Also add information on how to contact you by electronic and paper mail. - - If the program does terminal interaction, make it output a short -notice like this when it starts in an interactive mode: - - Copyright (C) - This program comes with ABSOLUTELY NO WARRANTY; for details type `show w'. - This is free software, and you are welcome to redistribute it - under certain conditions; type `show c' for details. - -The hypothetical commands `show w' and `show c' should show the appropriate -parts of the General Public License. Of course, your program's commands -might be different; for a GUI interface, you would use an "about box". - - You should also get your employer (if you work as a programmer) or school, -if any, to sign a "copyright disclaimer" for the program, if necessary. -For more information on this, and how to apply and follow the GNU GPL, see -. - - The GNU General Public License does not permit incorporating your program -into proprietary programs. If your program is a subroutine library, you -may consider it more useful to permit linking proprietary applications with -the library. If this is what you want to do, use the GNU Lesser General -Public License instead of this License. But first, please read -. diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/README.md b/experiments/multiLayerModel/cds/madagascarRecipes/README.md deleted file mode 100644 index 9448f2b..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/README.md +++ /dev/null @@ -1,2 +0,0 @@ -# madagascarRecipes -Personal library of Madagascar recipes diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/__init__.py b/experiments/multiLayerModel/cds/madagascarRecipes/__init__.py deleted file mode 100644 index e69de29..0000000 diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/errorReport.py b/experiments/multiLayerModel/cds/madagascarRecipes/errorReport.py deleted file mode 100644 index 8509563..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/errorReport.py +++ /dev/null @@ -1,54 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# errorReport.py (Madagascar Recipe) -# -# Purpose: SCons error report function. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 01/06/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Python function to report errors in SCons - -If build fail report specific message. Import this module and add the following lineis to your SConstruct: - -import atexit -from madagascarRecipes.errorReport import print_build_failures - -Call error message by the following: - -atexit.register(print_build_failures,message) - -It will call function print_build_failures to report message variable that should be a string with the error message to be reported -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -__author__="Rodolfo Dirack " -__version__="1.0" - -def print_build_failures(message): - ''' - Print error message when build fail - :param message: string, error message - ''' - from SCons.Script import GetBuildFailures - for bf in GetBuildFailures(): - print("BUILD ERROR REPORT".center(40,'*')) - print("FAILED BUILD: %s" % (bf.node)) - print("MESSAGE REPORT".center(40,'*')) - print("%s" % (message)) - diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/generateSelfdoc b/experiments/multiLayerModel/cds/madagascarRecipes/generateSelfdoc deleted file mode 100644 index 4f98839..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/generateSelfdoc +++ /dev/null @@ -1,27 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# generateSelfdoc (SConscript) -# -# Purpose: Generate selfdoc for Madagascar recipes using pydoc. -# -# Important!: It should be called from 'scons' with the command -# 'scons -f generateSelfdoc' -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -import pydoc - -RECIPES=['kimodel', - 'pefInterpolation'] - -for i in RECIPES: - pydoc.writedoc(i) diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/kimodel.py b/experiments/multiLayerModel/cds/madagascarRecipes/kimodel.py deleted file mode 100644 index 177b625..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/kimodel.py +++ /dev/null @@ -1,187 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# kimodel.py (Madagascar Recipe) -# -# Purpose: Recipe to Kirchhoff modeling. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Madagascar recipe to Kirchhoff modeling - -Define functions to generate subsurface models, apply Kirchhoff modeling -and Kirchhoff Newton modeling. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -# Madagascar package -from rsf.proj import * - -__author__="Rodolfo Dirack " -__version__="1.0" - -def arr2str(array,sep=' '): - ''' - Convert a tuple into a comma separeted string - ''' - return string.join(map(str,array),sep) - - -def multiLayerModel( - interfaces, - dipsfile, - modelfile, - xmax, - zmax, - layers, - velocities - ): - ''' - - Generate a multi layer model to use in the program sfkirmod_newton - - :out interfaces: RSF filename, interpolated interfaces file - :out dipsfile: RSF filename, dips of interfaces - :out modelfile: RSF filename, model for ploting - :param xmax: interger, max x axis model distance - :param zmax: interger, max z axis model depth - :param layers: tuple array, points describing interfaces - :param velocities: tuple, layers velocities in Km/s - ''' - vstr=arr2str(velocities,',') - - n1 = len(layers[0]) - n2 = len(layers) - - Flow('layers.asc',None, - ''' - echo %s - n1=%d n2=%d o1=0 d1=%g - data_format=ascii_float in=$TARGET - ''' % (string.join(map(arr2str,layers),' '), - n1,n2,xmax/(n1-1))) - Flow('layers','layers.asc','dd form=native') - - d = 0.0101 # non-round for reproducibility - - Flow(interfaces,'layers', - 'spline o1=0 d1=%g n1=%d' % (d,int(1.5+xmax/d))) - Flow(dipsfile,interfaces,'deriv scale=y') - - Flow(modelfile,interfaces, - ''' - unif2 d1=%g n1=%d v00=%s - ''' % (d,int(1.5+zmax/d),vstr)) - -def kirchoffNewtonModeling( - reflectors, - reflectorsDip, - filename, - velocities): - ''' - Kirchhoff modeling for a multi layer model - ''' - vstr=arr2str(velocities,',') - - # Kirchoff modeling for multi layer model - Flow(filename,[reflectors,reflectorsDip], - ''' - kirmod_newton nt=1001 dt=0.004 freq=10 - ns=201 ds=0.025 nh=161 dh=0.025 h0=0 s0=0 verb=y cmp=y - vstatus=0 velocity=%s debug=n fwdxini=y - xref=0 zref=0 dip=${SOURCES[1]} | - put d2=0.0125 label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s - ''' % (vstr)) - -def gaussianReflector(filename='gaussianReflector'): - ''' - - Generate a gaussian reflector to use in the program sfkirmod - - :param filename: RSF filename, gaussian reflector file - ''' - - # Modeling: Gaussian reflector in a velocity linear model - # velocity increases with depth with a 0.5 velocity gradient - Flow('gaussianReflector',None, - ''' - math d1=0.01 n1=2001 o1=-5 unit1=km label1=Offset - output="4-3*exp(-(x1-5)^2/9)" - ''') - - # Velocity Model - Flow('velocityModel','gaussianReflector', - ''' - window min1=0 max1=10 | - spray axis=1 n=451 d=0.01 o=0 label=Depth unit=km | - math output="1.5+0.5*x1+0.0*x2" - ''') - - Flow('reflectorDip','gaussianReflector','math output="2/3*(x1-5)*input" ') - - -def kirchoffModeling( - filename, - reflector, - reflectorDip, - nh, - dh, - h0, - ns, - ds, - s0, - freq, - dt, - nt, - v0, - gradz - ): - ''' - - Kirchhoff modeling function. Generate (time,cmp,offset) datacube for - a given v(z) model (Velocity increases with depth). - - :out filename: RSF filename, Modeled data cube - :param reflector: RSF filename, interfaces file - :param reflectorDip: RSF filename, interfaces dips file - :param nh: integer, number of offsets - :param dh: float, offset sampling - :param h0: float, first offset - :param ns: integer, number of shots/cmps - :param ds: float, shots/cmps sampling - :param s0: float, first shot/cmp position - :param freq: integer, ricker pulse frequency - :param dt: float, time sampling - :param nt: float, number of time samples - :param v0: float, near surface velocity - :param gradz: float, velocity gradient - ''' - - halfOffset = dh/2. - - # Kirchoff Modeling - Flow(filename,[reflector, reflectorDip], - ''' - kirmod cmp=y dip=${SOURCES[1]} - nh=%i dh=%g h0=%g - ns=%i ds=%g s0=%g - freq=%i dt=%g nt=%i - vel=%g gradz=%g gradx=0.0 verb=y | - put d2=%g label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s - ''' % (nh,dh,h0,ns,ds,s0,freq,dt,nt,v0,gradz,halfOffset)) - diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/pefInterpolation.py b/experiments/multiLayerModel/cds/madagascarRecipes/pefInterpolation.py deleted file mode 100644 index 495b635..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/pefInterpolation.py +++ /dev/null @@ -1,156 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# pefInterpolation.py (Madagascar Recipe) -# -# Purpose: Recipe to Preditive Adaptative Error Filters (PEF) interpolation -# to double CMP number of samples. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Madagascar recipe to PEF interpolation - -Interpolate and increase CMP sampling using Preditive Adaptative Error Filters (PEF) interpolation. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -__author__="Rodolfo Dirack " -__version__="1.0" - -# Madagascar package -from rsf.proj import * - -def pefInterpolation( - dataCube, - interpolated, - nm, - dm, - nt, - dt, - nhi, - a1, - a2, - rect1, - rect2 - ): - ''' - - PEF interpolation to double CMP sampling of the data cube - - :param dataCube: RSF filename, Seismic data cube to interpolate - :param interpolated: RSF filename, Interpolated seismic data cube - :param nm: integer, number of CMPs in the seismic data cube - :param dm: float, CMP sampling - :param nt: integer, number of time samples - :param dt: float, time sampling - :param nhi: integer, number of constant offsets gathers to interpolate - :param a1: integer, Number of PEF coeficients in time axis - :param a2: integer, Number of PEF coeficients in space axis - :param rect1: integer, Smooth radius in time - :param rect2: integer, Smooth radius in space - ''' - - # Divide CMP sampling - dm = dm/2 - - # Define mask file names using input filename - mask1 = dataCube+'-mask1' - mask = dataCube+'-mask' - aa = dataCube+'-aa' - bb = dataCube+'-bb' - a = dataCube+'-a' - b = dataCube+'-b' - zeroTraceGather = dataCube+'-zeroedGather' - mask0 = dataCube+'-mask0' - - - # Build a mask to interleave zero traces with original data traces - Flow(aa,None,'spike n1=%i d1=%g o1=0' %(nm,dm)) - Flow(bb,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask1,[bb, aa], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - - Flow(a,None,'spike n1=%i d1=%g o1=0' % (nm,dm)) - Flow(b,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask,[a, b], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - Flow(zeroTraceGather,b, - ''' - spray axis=2 n=%i d=%g | - transp | - put label2=Offset unit2=Km label1=Time unit1=s - ''' %(nt,dm)) - - # Data Mask with double of traces in CMP (half of CMP sampling) - # Keep the same Time and Offset original data sampling - Flow(mask0,mask, - ''' - spray axis=1 n=%i d=%g - ''' %(nt,dt)) - - totalPefIterations = 100 - totalInterpolationIterations = 50 - - offsetGathers = [] - for offsetGatherIndex in range(nhi): - - offsetGather = dataCube+"-offsetGather-%i" % offsetGatherIndex - resampledOffsetGather = dataCube+"-resampledGather-%i" % offsetGatherIndex - interpolatedOffsetGather = dataCube+"-interpolatedGather-%i" % offsetGatherIndex - pefCoeficients = dataCube+"-pefCoeficients-%i" % offsetGatherIndex - - Flow(offsetGather,dataCube, - ''' - window n2=1 f2=%i - ''' % (offsetGatherIndex)) - - Flow(resampledOffsetGather,[offsetGather,zeroTraceGather], - ''' - interleave axis=2 ${SOURCES[1]} - ''') - - # Calculate adaptive PEF coeficients - Flow(pefCoeficients,[resampledOffsetGather,mask0], - ''' - apef jump=2 a=%i,%i rect1=%i rect2=%i niter=%g verb=y - maskin=${SOURCES[1]} - ''' % (a1,a2,rect1,rect2,totalPefIterations)) - - # Interpolation - Flow(interpolatedOffsetGather, [resampledOffsetGather,pefCoeficients,mask0,mask1], - ''' - miss4 exact=y filt=${SOURCES[1]} mask=${SOURCES[2]} niter=%g verb=y | - put d2=%g - ''' % (totalInterpolationIterations,dm)) - - offsetGathers.append(interpolatedOffsetGather) - - # Concatenate interpolated sections - Flow(interpolated,offsetGathers, - ''' - rcat axis=3 ${SOURCES[1:%d]} | - transp plane=23 - ''' % nhi) - - diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/plot.py b/experiments/multiLayerModel/cds/madagascarRecipes/plot.py deleted file mode 100644 index 3715a31..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/plot.py +++ /dev/null @@ -1,52 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# plot.py (Madagascar Recipe) -# -# Purpose: Recipe to generate plot files (VPL). -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Madagascar recipe to plot - -Define functions to generate vpl files. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -# Madagascar package -from rsf.proj import * - -__author__="Rodolfo Dirack " -__version__="1.0" - -def velplot(title,label1='Depth',unit1='Km'): - ''' - Velocity ploting function - :param title: RSF filename, velocity file - :param label1: string, first axis label - :param unit1: string, first axis unit - ''' - - return ''' - grey color=j allpos=y title="%s" scalebar=y - barlabel=Velocity barunit=Km/s - label1="%s" unit1="%s" label2=Lateral unit2=Km - barreverse=y pclip=100 bias=2.8 - ''' % (title,label1,unit1) - - diff --git a/experiments/multiLayerModel/cds/madagascarRecipes/velocityAnalisys.py b/experiments/multiLayerModel/cds/madagascarRecipes/velocityAnalisys.py deleted file mode 100644 index 96a5826..0000000 --- a/experiments/multiLayerModel/cds/madagascarRecipes/velocityAnalisys.py +++ /dev/null @@ -1,100 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# velocityAnalisys.py () -# -# Purpose: Madagascar recipe to do velocity analisys. -# -# Site: https://dirack.github.io -# -# Version 1.0 -# -# Programer: Rodolfo A C Neves (Dirack) 01/04/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -from rsf.proj import * -import math - -# Selfdoc string -''' -Madagascar recipe to Velocity analisys with automatic velocity scan and picking - -Define functions to do a velocity analisys of a (time,offset,CMP) datacube. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -# Madagascar package -from rsf.proj import * - -__author__="Rodolfo Dirack " -__version__="1.0" - -def velocityAnalisys( - dataCube, - pick, - nmo, - vrms, - v0, - dv, - nv, - vel0, - rect1=15, - rect2=40, - rect3=3, - dt=0.004) - ''' - Velocity analisys and automatic picking and semblance - :param datacube: RSF filename, (time, offset, CMP) datacube - :out pick: RSF filename, picking velocity - :out stack: RSF filename, stacked section - :out vrms: RSF filename, vrms velocity - :param dt: float, time sampling - :param v0: float, inicial scan velocity - :param dv: float, scan velocity sampling - :param nv: int, number of scan velocities - :param vel0: float, initial velocity picking - :param rect1: int, smooth time radius to velocity scan - :param rect2: int, smooth space radius to velocity scan - :param rect3: int, smooth time radius post picking - ''' - - scan = 'scan'+dataCube - nmo = 'nmo'+dataCube - - # Velocity scan - Flow(scan,[dataCube], - ''' - vscan semblance=y - v0=%g dv=%g nv=%d - ''' % (v0,dv,nv)) - - # Velocity picking - Flow(pick,scan, - ''' - pick rect1=%i rect2=%i vel0=%g smooth=y | - smooth rect1=%i - ''' %(rect1,rect2,vel0,rect3)) - - # NMO and stack - Flow(nmo,[dataCube,pick], - 'nmo half=n velocity=${SOURCES[1]}') - - Flow(stack,nmo,'stack') - - # generate RMS velocity : - # sfmul multiply input data by itself to make v**2 - # sfcausint integrate v**2 (this does not include scaling by dt) - # sfmath scale by dt and divide by t - Flow(vrms,pick, - ''' - sfmul $SOURCE | - sfcausint | - sfmath output="sqrt(input*%g/(x1+%g))" | - sfput n3=1 d3=1 o3=0 - '''%(dt,dt)) diff --git a/experiments/multiLayerModel/cre/SConscript b/experiments/multiLayerModel/cre/SConscript deleted file mode 100644 index d006cfd..0000000 --- a/experiments/multiLayerModel/cre/SConscript +++ /dev/null @@ -1,100 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# SConscript (Madagascar Script) -# -# Purpose: Generate model and data cubes to CRE stacking process. -# This script generates model, original and PEF interpolated data cubes. -# It depends of the 'creStack.py' script that defines functions. -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Madagascar package -from rsf.proj import * - -# Error report function -import atexit -from madagascarRecipes.errorReport import print_build_failures - -# CRE recipe -from madagascarRecipes.pefInterpolation import pefInterpolation as pefin -from madagascarRecipes.kimodel import multiLayerModel as mlmod -from madagascarRecipes.kimodel import kirchoffNewtonModeling as kinewmod - -xmax = 6.0 -zmax = 2.0 - -layers = ((0.30,0.50,0.20,0.30), - (1.65,1.85,1.55,1.65)) - -velocities = (1.508, - 1.690, - 2.0) - -# Generate multi layer model and data cube -mlmod(interfaces='interfaces', - dipsfile='interfacesDip', - modelfile='mod1', - xmax=xmax, - zmax=zmax, - layers=layers, - velocities=velocities) - -kinewmod(reflectors='interfaces', - reflectorsDip='interfacesDip', - filename='multiLayerDataCube', - velocities=velocities) - -# PEF interpolation of the data cube -a1=10 -a2=2 -rect1=50 -rect2=2 - -pefin(dataCube='multiLayerDataCube', - interpolated='interpolatedDataCube', - nm=201, - dm=0.025, - nt=1001, - dt=0.004, - nhi=161, - a1=a1, - a2=a2, - rect1=rect1, - rect2=rect2) - -# Do the PEF interpolation one more time -# to increase the CMP sampling -pefin('interpolatedDataCube', - 'interpolatedDataCube2', - nm=402, - dm=0.0125, - nt=1001, - dt=0.004, - nhi=161, - a1=a1, - a2=a2, - rect1=rect1, - rect2=rect2) - -# Use aliases to split building -Alias('model','multiLayerDataCube.rsf') -Alias('pef1','interpolatedDataCube.rsf') -Alias('pef2','interpolatedDataCube2.rsf') - -# Show error message if fail - -message = ''' -SCosncript Failed build -''' - -atexit.register(print_build_failures,message) -End() diff --git a/experiments/multiLayerModel/cre/SConscript2 b/experiments/multiLayerModel/cre/SConscript2 deleted file mode 100644 index 0c63c46..0000000 --- a/experiments/multiLayerModel/cre/SConscript2 +++ /dev/null @@ -1,61 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# SConscript (Madagascar Script) -# -# Purpose: Contatenate CRE stacked traces. -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Madagascar package -from rsf.proj import * - -# Error report function -import atexit -from madagascarRecipes.errorReport import print_build_failures - -import glob - -# Build the cre stacked section -# throughout cre stacked traces sorting -files = glob.glob('creStackedTrace-*.rsf') -length = len(files) - -sortedFiles = [] -for i in range(length): - string = 'creStackedTrace-m0-%i.rsf' % i - sortedFiles.append(string) - -Flow('stackedSection',sortedFiles, - ''' - rcat axis=2 ${SOURCES[1:%d]} --out=stdout - ''' % len(files)) - -Result('stackedSection','grey title="Stacked Section"') - -Flow('filtStackedSection','stackedSection', - ''' - bandpass fhi=20 --out=stdout - ''') - -Result('filtStackedSection','grey title="Filtered Stacked Section"') - -# Use aliases to split building -Alias('stack',['stackedSection.rsf','filtStackedSection.rsf']) - -# Show error message if fail - -message = ''' -SCosncript Failed build -''' - -atexit.register(print_build_failures,message) -End() diff --git a/experiments/multiLayerModel/cre/SConstruct b/experiments/multiLayerModel/cre/SConstruct index 5e8513d..9395bdd 100644 --- a/experiments/multiLayerModel/cre/SConstruct +++ b/experiments/multiLayerModel/cre/SConstruct @@ -3,12 +3,12 @@ # # SConstruct (Madagascar Script) # -# Purpose: Build the interpolation process twice to improve CMP sampling -# and CRE stacking results. +# Purpose: Common Reflection Element (CRE) stack applied to a seismic data cube +# from Kirchhoff Newton modeling of a multi layer velocity model. # -# Site: http://www.dirackslounge.online +# Site: https://dirack.github.io # -# Version 1.0 +# Version 2.0 # # Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 # @@ -22,13 +22,6 @@ from rsf.proj import * # Import glob python library import glob -# Error report function -import atexit -from madagascarRecipes.errorReport import print_build_failures - -# Call SConscript to generate input: original and interpolated data cubes -# SConscript('SConscript') - # It uses Very Fast Simulated Aneeling and non hyperbolic CRS # to get zero offset CRS parameters (RN, RNIP and BETA) from data cube v0 = 1.5 @@ -36,9 +29,12 @@ ot0 = 0.3 dt0 = 0.008 nt0 = 331 om0 = 2 -dm0 = 0.025 -nm0 = 108 -dataCube='multiLayerDataCube' +dm0 = 0.25 +nm0 = 8 +inputDir= '../../modelingAndPEFInterpolation/multiLayerModel/' +dataCube=inputDir+'multiLayerDataCube.rsf' +interpolatedDataCube=inputDir+'interpolatedDataCube2.rsf' +files = [] for i in range(nm0): @@ -63,22 +59,22 @@ for i in range(nm0): creGatherPlot = 'crePlot-m0-%g-t0-%g' % (i,j) # Very Fast Simulated Aneelling Global Optimization (VFSA) - Flow(crsParameters,'interpolatedDataCube2', + Flow(crsParameters,dataCube, ''' vfsacrsnh m0=%g v0=%g t0=%g verb=y repeat=3 ''' % (m0,v0,t0)) # Calculate CRE trajectory - Flow(creMhCoordinates,['interpolatedDataCube2',crsParameters], + Flow(creMhCoordinates,[interpolatedDataCube,crsParameters], ''' cretrajec verb=y m0=%g param=${SOURCES[1]} | put unit1="Offset" label1="Km" ''' % (m0)) #Get CRE Gather from interpolated Data Cube - Flow([creGather,creMcoordinate],['interpolatedDataCube2',creMhCoordinates], + Flow([creGather,creMcoordinate],[interpolatedDataCube,creMhCoordinates], ''' - getcregather aperture=50 verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} | + getcregather aperture=50 verb=y cremh=${SOURCES[1]} m=${TARGETS[1]} aperture=20 | put label1="Time" unit1="s" label2="Offset" unit2="km" ''') @@ -136,9 +132,15 @@ for i in range(nm0): put label1=t0 unit1=s label2=m0 unit2=Km --out=stdout ''') + files.append(creStackedTrace) + +# Stacked section depends of the last CRE stacked trace +Depends('stackedSection.rsf',[dataCube, + 'creStackedTrace-m0-'+str(nm0-1)+'.rsf', + 'creGatherTrace-m0-'+str(nm0-1)+'.rsf']) + # Build the cre stacked section # throughout cre stacked traces sorting -files = glob.glob('creStackedTrace-*.rsf') length = len(files) sortedFiles = [] @@ -156,11 +158,4 @@ Flow('filtStackedSection','stackedSection', bandpass fhi=20 --out=stdout ''') -# Show error message if fail -message = ''' -SConstruct Failed build when m0={m0} and t0={t0} -'''.format(m0=(i*dm0+om0),t0=(j*dt0+ot0)) - -atexit.register(print_build_failures,message) - End() diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/LICENSE b/experiments/multiLayerModel/cre/madagascarRecipes/LICENSE deleted file mode 100644 index f288702..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/LICENSE +++ /dev/null @@ -1,674 +0,0 @@ - GNU GENERAL PUBLIC LICENSE - Version 3, 29 June 2007 - - Copyright (C) 2007 Free Software Foundation, Inc. - Everyone is permitted to copy and distribute verbatim copies - of this license document, but changing it is not allowed. - - Preamble - - The GNU General Public License is a free, copyleft license for -software and other kinds of works. - - The licenses for most software and other practical works are designed -to take away your freedom to share and change the works. By contrast, -the GNU General Public License is intended to guarantee your freedom to -share and change all versions of a program--to make sure it remains free -software for all its users. We, the Free Software Foundation, use the -GNU General Public License for most of our software; it applies also to -any other work released this way by its authors. You can apply it to -your programs, too. - - When we speak of free software, we are referring to freedom, not -price. Our General Public Licenses are designed to make sure that you -have the freedom to distribute copies of free software (and charge for -them if you wish), that you receive source code or can get it if you -want it, that you can change the software or use pieces of it in new -free programs, and that you know you can do these things. - - To protect your rights, we need to prevent others from denying you -these rights or asking you to surrender the rights. Therefore, you have -certain responsibilities if you distribute copies of the software, or if -you modify it: responsibilities to respect the freedom of others. - - For example, if you distribute copies of such a program, whether -gratis or for a fee, you must pass on to the recipients the same -freedoms that you received. You must make sure that they, too, receive -or can get the source code. And you must show them these terms so they -know their rights. - - Developers that use the GNU GPL protect your rights with two steps: -(1) assert copyright on the software, and (2) offer you this License -giving you legal permission to copy, distribute and/or modify it. - - For the developers' and authors' protection, the GPL clearly explains -that there is no warranty for this free software. For both users' and -authors' sake, the GPL requires that modified versions be marked as -changed, so that their problems will not be attributed erroneously to -authors of previous versions. - - Some devices are designed to deny users access to install or run -modified versions of the software inside them, although the manufacturer -can do so. This is fundamentally incompatible with the aim of -protecting users' freedom to change the software. The systematic -pattern of such abuse occurs in the area of products for individuals to -use, which is precisely where it is most unacceptable. Therefore, we -have designed this version of the GPL to prohibit the practice for those -products. If such problems arise substantially in other domains, we -stand ready to extend this provision to those domains in future versions -of the GPL, as needed to protect the freedom of users. - - Finally, every program is threatened constantly by software patents. -States should not allow patents to restrict development and use of -software on general-purpose computers, but in those that do, we wish to -avoid the special danger that patents applied to a free program could -make it effectively proprietary. To prevent this, the GPL assures that -patents cannot be used to render the program non-free. - - The precise terms and conditions for copying, distribution and -modification follow. - - TERMS AND CONDITIONS - - 0. Definitions. - - "This License" refers to version 3 of the GNU General Public License. - - "Copyright" also means copyright-like laws that apply to other kinds of -works, such as semiconductor masks. - - "The Program" refers to any copyrightable work licensed under this -License. Each licensee is addressed as "you". "Licensees" and -"recipients" may be individuals or organizations. - - To "modify" a work means to copy from or adapt all or part of the work -in a fashion requiring copyright permission, other than the making of an -exact copy. The resulting work is called a "modified version" of the -earlier work or a work "based on" the earlier work. - - A "covered work" means either the unmodified Program or a work based -on the Program. - - To "propagate" a work means to do anything with it that, without -permission, would make you directly or secondarily liable for -infringement under applicable copyright law, except executing it on a -computer or modifying a private copy. Propagation includes copying, -distribution (with or without modification), making available to the -public, and in some countries other activities as well. - - To "convey" a work means any kind of propagation that enables other -parties to make or receive copies. Mere interaction with a user through -a computer network, with no transfer of a copy, is not conveying. - - An interactive user interface displays "Appropriate Legal Notices" -to the extent that it includes a convenient and prominently visible -feature that (1) displays an appropriate copyright notice, and (2) -tells the user that there is no warranty for the work (except to the -extent that warranties are provided), that licensees may convey the -work under this License, and how to view a copy of this License. If -the interface presents a list of user commands or options, such as a -menu, a prominent item in the list meets this criterion. - - 1. Source Code. - - The "source code" for a work means the preferred form of the work -for making modifications to it. "Object code" means any non-source -form of a work. - - A "Standard Interface" means an interface that either is an official -standard defined by a recognized standards body, or, in the case of -interfaces specified for a particular programming language, one that -is widely used among developers working in that language. - - The "System Libraries" of an executable work include anything, other -than the work as a whole, that (a) is included in the normal form of -packaging a Major Component, but which is not part of that Major -Component, and (b) serves only to enable use of the work with that -Major Component, or to implement a Standard Interface for which an -implementation is available to the public in source code form. A -"Major Component", in this context, means a major essential component -(kernel, window system, and so on) of the specific operating system -(if any) on which the executable work runs, or a compiler used to -produce the work, or an object code interpreter used to run it. - - The "Corresponding Source" for a work in object code form means all -the source code needed to generate, install, and (for an executable -work) run the object code and to modify the work, including scripts to -control those activities. However, it does not include the work's -System Libraries, or general-purpose tools or generally available free -programs which are used unmodified in performing those activities but -which are not part of the work. For example, Corresponding Source -includes interface definition files associated with source files for -the work, and the source code for shared libraries and dynamically -linked subprograms that the work is specifically designed to require, -such as by intimate data communication or control flow between those -subprograms and other parts of the work. - - The Corresponding Source need not include anything that users -can regenerate automatically from other parts of the Corresponding -Source. - - The Corresponding Source for a work in source code form is that -same work. - - 2. Basic Permissions. - - All rights granted under this License are granted for the term of -copyright on the Program, and are irrevocable provided the stated -conditions are met. This License explicitly affirms your unlimited -permission to run the unmodified Program. The output from running a -covered work is covered by this License only if the output, given its -content, constitutes a covered work. This License acknowledges your -rights of fair use or other equivalent, as provided by copyright law. - - You may make, run and propagate covered works that you do not -convey, without conditions so long as your license otherwise remains -in force. You may convey covered works to others for the sole purpose -of having them make modifications exclusively for you, or provide you -with facilities for running those works, provided that you comply with -the terms of this License in conveying all material for which you do -not control copyright. Those thus making or running the covered works -for you must do so exclusively on your behalf, under your direction -and control, on terms that prohibit them from making any copies of -your copyrighted material outside their relationship with you. - - Conveying under any other circumstances is permitted solely under -the conditions stated below. Sublicensing is not allowed; section 10 -makes it unnecessary. - - 3. Protecting Users' Legal Rights From Anti-Circumvention Law. - - No covered work shall be deemed part of an effective technological -measure under any applicable law fulfilling obligations under article -11 of the WIPO copyright treaty adopted on 20 December 1996, or -similar laws prohibiting or restricting circumvention of such -measures. - - When you convey a covered work, you waive any legal power to forbid -circumvention of technological measures to the extent such circumvention -is effected by exercising rights under this License with respect to -the covered work, and you disclaim any intention to limit operation or -modification of the work as a means of enforcing, against the work's -users, your or third parties' legal rights to forbid circumvention of -technological measures. - - 4. Conveying Verbatim Copies. - - You may convey verbatim copies of the Program's source code as you -receive it, in any medium, provided that you conspicuously and -appropriately publish on each copy an appropriate copyright notice; -keep intact all notices stating that this License and any -non-permissive terms added in accord with section 7 apply to the code; -keep intact all notices of the absence of any warranty; and give all -recipients a copy of this License along with the Program. - - You may charge any price or no price for each copy that you convey, -and you may offer support or warranty protection for a fee. - - 5. Conveying Modified Source Versions. - - You may convey a work based on the Program, or the modifications to -produce it from the Program, in the form of source code under the -terms of section 4, provided that you also meet all of these conditions: - - a) The work must carry prominent notices stating that you modified - it, and giving a relevant date. - - b) The work must carry prominent notices stating that it is - released under this License and any conditions added under section - 7. This requirement modifies the requirement in section 4 to - "keep intact all notices". - - c) You must license the entire work, as a whole, under this - License to anyone who comes into possession of a copy. This - License will therefore apply, along with any applicable section 7 - additional terms, to the whole of the work, and all its parts, - regardless of how they are packaged. This License gives no - permission to license the work in any other way, but it does not - invalidate such permission if you have separately received it. - - d) If the work has interactive user interfaces, each must display - Appropriate Legal Notices; however, if the Program has interactive - interfaces that do not display Appropriate Legal Notices, your - work need not make them do so. - - A compilation of a covered work with other separate and independent -works, which are not by their nature extensions of the covered work, -and which are not combined with it such as to form a larger program, -in or on a volume of a storage or distribution medium, is called an -"aggregate" if the compilation and its resulting copyright are not -used to limit the access or legal rights of the compilation's users -beyond what the individual works permit. Inclusion of a covered work -in an aggregate does not cause this License to apply to the other -parts of the aggregate. - - 6. Conveying Non-Source Forms. - - You may convey a covered work in object code form under the terms -of sections 4 and 5, provided that you also convey the -machine-readable Corresponding Source under the terms of this License, -in one of these ways: - - a) Convey the object code in, or embodied in, a physical product - (including a physical distribution medium), accompanied by the - Corresponding Source fixed on a durable physical medium - customarily used for software interchange. - - b) Convey the object code in, or embodied in, a physical product - (including a physical distribution medium), accompanied by a - written offer, valid for at least three years and valid for as - long as you offer spare parts or customer support for that product - model, to give anyone who possesses the object code either (1) a - copy of the Corresponding Source for all the software in the - product that is covered by this License, on a durable physical - medium customarily used for software interchange, for a price no - more than your reasonable cost of physically performing this - conveying of source, or (2) access to copy the - Corresponding Source from a network server at no charge. - - c) Convey individual copies of the object code with a copy of the - written offer to provide the Corresponding Source. This - alternative is allowed only occasionally and noncommercially, and - only if you received the object code with such an offer, in accord - with subsection 6b. - - d) Convey the object code by offering access from a designated - place (gratis or for a charge), and offer equivalent access to the - Corresponding Source in the same way through the same place at no - further charge. You need not require recipients to copy the - Corresponding Source along with the object code. If the place to - copy the object code is a network server, the Corresponding Source - may be on a different server (operated by you or a third party) - that supports equivalent copying facilities, provided you maintain - clear directions next to the object code saying where to find the - Corresponding Source. Regardless of what server hosts the - Corresponding Source, you remain obligated to ensure that it is - available for as long as needed to satisfy these requirements. - - e) Convey the object code using peer-to-peer transmission, provided - you inform other peers where the object code and Corresponding - Source of the work are being offered to the general public at no - charge under subsection 6d. - - A separable portion of the object code, whose source code is excluded -from the Corresponding Source as a System Library, need not be -included in conveying the object code work. - - A "User Product" is either (1) a "consumer product", which means any -tangible personal property which is normally used for personal, family, -or household purposes, or (2) anything designed or sold for incorporation -into a dwelling. In determining whether a product is a consumer product, -doubtful cases shall be resolved in favor of coverage. For a particular -product received by a particular user, "normally used" refers to a -typical or common use of that class of product, regardless of the status -of the particular user or of the way in which the particular user -actually uses, or expects or is expected to use, the product. A product -is a consumer product regardless of whether the product has substantial -commercial, industrial or non-consumer uses, unless such uses represent -the only significant mode of use of the product. - - "Installation Information" for a User Product means any methods, -procedures, authorization keys, or other information required to install -and execute modified versions of a covered work in that User Product from -a modified version of its Corresponding Source. The information must -suffice to ensure that the continued functioning of the modified object -code is in no case prevented or interfered with solely because -modification has been made. - - If you convey an object code work under this section in, or with, or -specifically for use in, a User Product, and the conveying occurs as -part of a transaction in which the right of possession and use of the -User Product is transferred to the recipient in perpetuity or for a -fixed term (regardless of how the transaction is characterized), the -Corresponding Source conveyed under this section must be accompanied -by the Installation Information. But this requirement does not apply -if neither you nor any third party retains the ability to install -modified object code on the User Product (for example, the work has -been installed in ROM). - - The requirement to provide Installation Information does not include a -requirement to continue to provide support service, warranty, or updates -for a work that has been modified or installed by the recipient, or for -the User Product in which it has been modified or installed. Access to a -network may be denied when the modification itself materially and -adversely affects the operation of the network or violates the rules and -protocols for communication across the network. - - Corresponding Source conveyed, and Installation Information provided, -in accord with this section must be in a format that is publicly -documented (and with an implementation available to the public in -source code form), and must require no special password or key for -unpacking, reading or copying. - - 7. Additional Terms. - - "Additional permissions" are terms that supplement the terms of this -License by making exceptions from one or more of its conditions. -Additional permissions that are applicable to the entire Program shall -be treated as though they were included in this License, to the extent -that they are valid under applicable law. If additional permissions -apply only to part of the Program, that part may be used separately -under those permissions, but the entire Program remains governed by -this License without regard to the additional permissions. - - When you convey a copy of a covered work, you may at your option -remove any additional permissions from that copy, or from any part of -it. (Additional permissions may be written to require their own -removal in certain cases when you modify the work.) You may place -additional permissions on material, added by you to a covered work, -for which you have or can give appropriate copyright permission. - - Notwithstanding any other provision of this License, for material you -add to a covered work, you may (if authorized by the copyright holders of -that material) supplement the terms of this License with terms: - - a) Disclaiming warranty or limiting liability differently from the - terms of sections 15 and 16 of this License; or - - b) Requiring preservation of specified reasonable legal notices or - author attributions in that material or in the Appropriate Legal - Notices displayed by works containing it; or - - c) Prohibiting misrepresentation of the origin of that material, or - requiring that modified versions of such material be marked in - reasonable ways as different from the original version; or - - d) Limiting the use for publicity purposes of names of licensors or - authors of the material; or - - e) Declining to grant rights under trademark law for use of some - trade names, trademarks, or service marks; or - - f) Requiring indemnification of licensors and authors of that - material by anyone who conveys the material (or modified versions of - it) with contractual assumptions of liability to the recipient, for - any liability that these contractual assumptions directly impose on - those licensors and authors. - - All other non-permissive additional terms are considered "further -restrictions" within the meaning of section 10. If the Program as you -received it, or any part of it, contains a notice stating that it is -governed by this License along with a term that is a further -restriction, you may remove that term. If a license document contains -a further restriction but permits relicensing or conveying under this -License, you may add to a covered work material governed by the terms -of that license document, provided that the further restriction does -not survive such relicensing or conveying. - - If you add terms to a covered work in accord with this section, you -must place, in the relevant source files, a statement of the -additional terms that apply to those files, or a notice indicating -where to find the applicable terms. - - Additional terms, permissive or non-permissive, may be stated in the -form of a separately written license, or stated as exceptions; -the above requirements apply either way. - - 8. Termination. - - You may not propagate or modify a covered work except as expressly -provided under this License. Any attempt otherwise to propagate or -modify it is void, and will automatically terminate your rights under -this License (including any patent licenses granted under the third -paragraph of section 11). - - However, if you cease all violation of this License, then your -license from a particular copyright holder is reinstated (a) -provisionally, unless and until the copyright holder explicitly and -finally terminates your license, and (b) permanently, if the copyright -holder fails to notify you of the violation by some reasonable means -prior to 60 days after the cessation. - - Moreover, your license from a particular copyright holder is -reinstated permanently if the copyright holder notifies you of the -violation by some reasonable means, this is the first time you have -received notice of violation of this License (for any work) from that -copyright holder, and you cure the violation prior to 30 days after -your receipt of the notice. - - Termination of your rights under this section does not terminate the -licenses of parties who have received copies or rights from you under -this License. If your rights have been terminated and not permanently -reinstated, you do not qualify to receive new licenses for the same -material under section 10. - - 9. Acceptance Not Required for Having Copies. - - You are not required to accept this License in order to receive or -run a copy of the Program. Ancillary propagation of a covered work -occurring solely as a consequence of using peer-to-peer transmission -to receive a copy likewise does not require acceptance. However, -nothing other than this License grants you permission to propagate or -modify any covered work. These actions infringe copyright if you do -not accept this License. Therefore, by modifying or propagating a -covered work, you indicate your acceptance of this License to do so. - - 10. Automatic Licensing of Downstream Recipients. - - Each time you convey a covered work, the recipient automatically -receives a license from the original licensors, to run, modify and -propagate that work, subject to this License. You are not responsible -for enforcing compliance by third parties with this License. - - An "entity transaction" is a transaction transferring control of an -organization, or substantially all assets of one, or subdividing an -organization, or merging organizations. If propagation of a covered -work results from an entity transaction, each party to that -transaction who receives a copy of the work also receives whatever -licenses to the work the party's predecessor in interest had or could -give under the previous paragraph, plus a right to possession of the -Corresponding Source of the work from the predecessor in interest, if -the predecessor has it or can get it with reasonable efforts. - - You may not impose any further restrictions on the exercise of the -rights granted or affirmed under this License. For example, you may -not impose a license fee, royalty, or other charge for exercise of -rights granted under this License, and you may not initiate litigation -(including a cross-claim or counterclaim in a lawsuit) alleging that -any patent claim is infringed by making, using, selling, offering for -sale, or importing the Program or any portion of it. - - 11. Patents. - - A "contributor" is a copyright holder who authorizes use under this -License of the Program or a work on which the Program is based. The -work thus licensed is called the contributor's "contributor version". - - A contributor's "essential patent claims" are all patent claims -owned or controlled by the contributor, whether already acquired or -hereafter acquired, that would be infringed by some manner, permitted -by this License, of making, using, or selling its contributor version, -but do not include claims that would be infringed only as a -consequence of further modification of the contributor version. For -purposes of this definition, "control" includes the right to grant -patent sublicenses in a manner consistent with the requirements of -this License. - - Each contributor grants you a non-exclusive, worldwide, royalty-free -patent license under the contributor's essential patent claims, to -make, use, sell, offer for sale, import and otherwise run, modify and -propagate the contents of its contributor version. - - In the following three paragraphs, a "patent license" is any express -agreement or commitment, however denominated, not to enforce a patent -(such as an express permission to practice a patent or covenant not to -sue for patent infringement). To "grant" such a patent license to a -party means to make such an agreement or commitment not to enforce a -patent against the party. - - If you convey a covered work, knowingly relying on a patent license, -and the Corresponding Source of the work is not available for anyone -to copy, free of charge and under the terms of this License, through a -publicly available network server or other readily accessible means, -then you must either (1) cause the Corresponding Source to be so -available, or (2) arrange to deprive yourself of the benefit of the -patent license for this particular work, or (3) arrange, in a manner -consistent with the requirements of this License, to extend the patent -license to downstream recipients. "Knowingly relying" means you have -actual knowledge that, but for the patent license, your conveying the -covered work in a country, or your recipient's use of the covered work -in a country, would infringe one or more identifiable patents in that -country that you have reason to believe are valid. - - If, pursuant to or in connection with a single transaction or -arrangement, you convey, or propagate by procuring conveyance of, a -covered work, and grant a patent license to some of the parties -receiving the covered work authorizing them to use, propagate, modify -or convey a specific copy of the covered work, then the patent license -you grant is automatically extended to all recipients of the covered -work and works based on it. - - A patent license is "discriminatory" if it does not include within -the scope of its coverage, prohibits the exercise of, or is -conditioned on the non-exercise of one or more of the rights that are -specifically granted under this License. You may not convey a covered -work if you are a party to an arrangement with a third party that is -in the business of distributing software, under which you make payment -to the third party based on the extent of your activity of conveying -the work, and under which the third party grants, to any of the -parties who would receive the covered work from you, a discriminatory -patent license (a) in connection with copies of the covered work -conveyed by you (or copies made from those copies), or (b) primarily -for and in connection with specific products or compilations that -contain the covered work, unless you entered into that arrangement, -or that patent license was granted, prior to 28 March 2007. - - Nothing in this License shall be construed as excluding or limiting -any implied license or other defenses to infringement that may -otherwise be available to you under applicable patent law. - - 12. No Surrender of Others' Freedom. - - If conditions are imposed on you (whether by court order, agreement or -otherwise) that contradict the conditions of this License, they do not -excuse you from the conditions of this License. If you cannot convey a -covered work so as to satisfy simultaneously your obligations under this -License and any other pertinent obligations, then as a consequence you may -not convey it at all. For example, if you agree to terms that obligate you -to collect a royalty for further conveying from those to whom you convey -the Program, the only way you could satisfy both those terms and this -License would be to refrain entirely from conveying the Program. - - 13. Use with the GNU Affero General Public License. - - Notwithstanding any other provision of this License, you have -permission to link or combine any covered work with a work licensed -under version 3 of the GNU Affero General Public License into a single -combined work, and to convey the resulting work. The terms of this -License will continue to apply to the part which is the covered work, -but the special requirements of the GNU Affero General Public License, -section 13, concerning interaction through a network will apply to the -combination as such. - - 14. Revised Versions of this License. - - The Free Software Foundation may publish revised and/or new versions of -the GNU General Public License from time to time. Such new versions will -be similar in spirit to the present version, but may differ in detail to -address new problems or concerns. - - Each version is given a distinguishing version number. If the -Program specifies that a certain numbered version of the GNU General -Public License "or any later version" applies to it, you have the -option of following the terms and conditions either of that numbered -version or of any later version published by the Free Software -Foundation. If the Program does not specify a version number of the -GNU General Public License, you may choose any version ever published -by the Free Software Foundation. - - If the Program specifies that a proxy can decide which future -versions of the GNU General Public License can be used, that proxy's -public statement of acceptance of a version permanently authorizes you -to choose that version for the Program. - - Later license versions may give you additional or different -permissions. However, no additional obligations are imposed on any -author or copyright holder as a result of your choosing to follow a -later version. - - 15. Disclaimer of Warranty. - - THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY -APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT -HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY -OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, -THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR -PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM -IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF -ALL NECESSARY SERVICING, REPAIR OR CORRECTION. - - 16. Limitation of Liability. - - IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING -WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MODIFIES AND/OR CONVEYS -THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, INCLUDING ANY -GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE -USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED TO LOSS OF -DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY YOU OR THIRD -PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER PROGRAMS), -EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE POSSIBILITY OF -SUCH DAMAGES. - - 17. Interpretation of Sections 15 and 16. - - If the disclaimer of warranty and limitation of liability provided -above cannot be given local legal effect according to their terms, -reviewing courts shall apply local law that most closely approximates -an absolute waiver of all civil liability in connection with the -Program, unless a warranty or assumption of liability accompanies a -copy of the Program in return for a fee. - - END OF TERMS AND CONDITIONS - - How to Apply These Terms to Your New Programs - - If you develop a new program, and you want it to be of the greatest -possible use to the public, the best way to achieve this is to make it -free software which everyone can redistribute and change under these terms. - - To do so, attach the following notices to the program. It is safest -to attach them to the start of each source file to most effectively -state the exclusion of warranty; and each file should have at least -the "copyright" line and a pointer to where the full notice is found. - - - Copyright (C) - - This program is free software: you can redistribute it and/or modify - it under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - This program is distributed in the hope that it will be useful, - but WITHOUT ANY WARRANTY; without even the implied warranty of - MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the - GNU General Public License for more details. - - You should have received a copy of the GNU General Public License - along with this program. If not, see . - -Also add information on how to contact you by electronic and paper mail. - - If the program does terminal interaction, make it output a short -notice like this when it starts in an interactive mode: - - Copyright (C) - This program comes with ABSOLUTELY NO WARRANTY; for details type `show w'. - This is free software, and you are welcome to redistribute it - under certain conditions; type `show c' for details. - -The hypothetical commands `show w' and `show c' should show the appropriate -parts of the General Public License. Of course, your program's commands -might be different; for a GUI interface, you would use an "about box". - - You should also get your employer (if you work as a programmer) or school, -if any, to sign a "copyright disclaimer" for the program, if necessary. -For more information on this, and how to apply and follow the GNU GPL, see -. - - The GNU General Public License does not permit incorporating your program -into proprietary programs. If your program is a subroutine library, you -may consider it more useful to permit linking proprietary applications with -the library. If this is what you want to do, use the GNU Lesser General -Public License instead of this License. But first, please read -. diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/README.md b/experiments/multiLayerModel/cre/madagascarRecipes/README.md deleted file mode 100644 index 9448f2b..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/README.md +++ /dev/null @@ -1,2 +0,0 @@ -# madagascarRecipes -Personal library of Madagascar recipes diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/__init__.py b/experiments/multiLayerModel/cre/madagascarRecipes/__init__.py deleted file mode 100644 index e69de29..0000000 diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/errorReport.py b/experiments/multiLayerModel/cre/madagascarRecipes/errorReport.py deleted file mode 100644 index 8509563..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/errorReport.py +++ /dev/null @@ -1,54 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# errorReport.py (Madagascar Recipe) -# -# Purpose: SCons error report function. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 01/06/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Python function to report errors in SCons - -If build fail report specific message. Import this module and add the following lineis to your SConstruct: - -import atexit -from madagascarRecipes.errorReport import print_build_failures - -Call error message by the following: - -atexit.register(print_build_failures,message) - -It will call function print_build_failures to report message variable that should be a string with the error message to be reported -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -__author__="Rodolfo Dirack " -__version__="1.0" - -def print_build_failures(message): - ''' - Print error message when build fail - :param message: string, error message - ''' - from SCons.Script import GetBuildFailures - for bf in GetBuildFailures(): - print("BUILD ERROR REPORT".center(40,'*')) - print("FAILED BUILD: %s" % (bf.node)) - print("MESSAGE REPORT".center(40,'*')) - print("%s" % (message)) - diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/generateSelfdoc b/experiments/multiLayerModel/cre/madagascarRecipes/generateSelfdoc deleted file mode 100644 index 4f98839..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/generateSelfdoc +++ /dev/null @@ -1,27 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# generateSelfdoc (SConscript) -# -# Purpose: Generate selfdoc for Madagascar recipes using pydoc. -# -# Important!: It should be called from 'scons' with the command -# 'scons -f generateSelfdoc' -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -import pydoc - -RECIPES=['kimodel', - 'pefInterpolation'] - -for i in RECIPES: - pydoc.writedoc(i) diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/kimodel.py b/experiments/multiLayerModel/cre/madagascarRecipes/kimodel.py deleted file mode 100644 index 177b625..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/kimodel.py +++ /dev/null @@ -1,187 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# kimodel.py (Madagascar Recipe) -# -# Purpose: Recipe to Kirchhoff modeling. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Madagascar recipe to Kirchhoff modeling - -Define functions to generate subsurface models, apply Kirchhoff modeling -and Kirchhoff Newton modeling. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -# Madagascar package -from rsf.proj import * - -__author__="Rodolfo Dirack " -__version__="1.0" - -def arr2str(array,sep=' '): - ''' - Convert a tuple into a comma separeted string - ''' - return string.join(map(str,array),sep) - - -def multiLayerModel( - interfaces, - dipsfile, - modelfile, - xmax, - zmax, - layers, - velocities - ): - ''' - - Generate a multi layer model to use in the program sfkirmod_newton - - :out interfaces: RSF filename, interpolated interfaces file - :out dipsfile: RSF filename, dips of interfaces - :out modelfile: RSF filename, model for ploting - :param xmax: interger, max x axis model distance - :param zmax: interger, max z axis model depth - :param layers: tuple array, points describing interfaces - :param velocities: tuple, layers velocities in Km/s - ''' - vstr=arr2str(velocities,',') - - n1 = len(layers[0]) - n2 = len(layers) - - Flow('layers.asc',None, - ''' - echo %s - n1=%d n2=%d o1=0 d1=%g - data_format=ascii_float in=$TARGET - ''' % (string.join(map(arr2str,layers),' '), - n1,n2,xmax/(n1-1))) - Flow('layers','layers.asc','dd form=native') - - d = 0.0101 # non-round for reproducibility - - Flow(interfaces,'layers', - 'spline o1=0 d1=%g n1=%d' % (d,int(1.5+xmax/d))) - Flow(dipsfile,interfaces,'deriv scale=y') - - Flow(modelfile,interfaces, - ''' - unif2 d1=%g n1=%d v00=%s - ''' % (d,int(1.5+zmax/d),vstr)) - -def kirchoffNewtonModeling( - reflectors, - reflectorsDip, - filename, - velocities): - ''' - Kirchhoff modeling for a multi layer model - ''' - vstr=arr2str(velocities,',') - - # Kirchoff modeling for multi layer model - Flow(filename,[reflectors,reflectorsDip], - ''' - kirmod_newton nt=1001 dt=0.004 freq=10 - ns=201 ds=0.025 nh=161 dh=0.025 h0=0 s0=0 verb=y cmp=y - vstatus=0 velocity=%s debug=n fwdxini=y - xref=0 zref=0 dip=${SOURCES[1]} | - put d2=0.0125 label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s - ''' % (vstr)) - -def gaussianReflector(filename='gaussianReflector'): - ''' - - Generate a gaussian reflector to use in the program sfkirmod - - :param filename: RSF filename, gaussian reflector file - ''' - - # Modeling: Gaussian reflector in a velocity linear model - # velocity increases with depth with a 0.5 velocity gradient - Flow('gaussianReflector',None, - ''' - math d1=0.01 n1=2001 o1=-5 unit1=km label1=Offset - output="4-3*exp(-(x1-5)^2/9)" - ''') - - # Velocity Model - Flow('velocityModel','gaussianReflector', - ''' - window min1=0 max1=10 | - spray axis=1 n=451 d=0.01 o=0 label=Depth unit=km | - math output="1.5+0.5*x1+0.0*x2" - ''') - - Flow('reflectorDip','gaussianReflector','math output="2/3*(x1-5)*input" ') - - -def kirchoffModeling( - filename, - reflector, - reflectorDip, - nh, - dh, - h0, - ns, - ds, - s0, - freq, - dt, - nt, - v0, - gradz - ): - ''' - - Kirchhoff modeling function. Generate (time,cmp,offset) datacube for - a given v(z) model (Velocity increases with depth). - - :out filename: RSF filename, Modeled data cube - :param reflector: RSF filename, interfaces file - :param reflectorDip: RSF filename, interfaces dips file - :param nh: integer, number of offsets - :param dh: float, offset sampling - :param h0: float, first offset - :param ns: integer, number of shots/cmps - :param ds: float, shots/cmps sampling - :param s0: float, first shot/cmp position - :param freq: integer, ricker pulse frequency - :param dt: float, time sampling - :param nt: float, number of time samples - :param v0: float, near surface velocity - :param gradz: float, velocity gradient - ''' - - halfOffset = dh/2. - - # Kirchoff Modeling - Flow(filename,[reflector, reflectorDip], - ''' - kirmod cmp=y dip=${SOURCES[1]} - nh=%i dh=%g h0=%g - ns=%i ds=%g s0=%g - freq=%i dt=%g nt=%i - vel=%g gradz=%g gradx=0.0 verb=y | - put d2=%g label3="CMP" unit3="Km" label2="Offset" unit2="Km" label1=Time unit1=s - ''' % (nh,dh,h0,ns,ds,s0,freq,dt,nt,v0,gradz,halfOffset)) - diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/pefInterpolation.py b/experiments/multiLayerModel/cre/madagascarRecipes/pefInterpolation.py deleted file mode 100644 index 495b635..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/pefInterpolation.py +++ /dev/null @@ -1,156 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# pefInterpolation.py (Madagascar Recipe) -# -# Purpose: Recipe to Preditive Adaptative Error Filters (PEF) interpolation -# to double CMP number of samples. -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Madagascar recipe to PEF interpolation - -Interpolate and increase CMP sampling using Preditive Adaptative Error Filters (PEF) interpolation. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -__author__="Rodolfo Dirack " -__version__="1.0" - -# Madagascar package -from rsf.proj import * - -def pefInterpolation( - dataCube, - interpolated, - nm, - dm, - nt, - dt, - nhi, - a1, - a2, - rect1, - rect2 - ): - ''' - - PEF interpolation to double CMP sampling of the data cube - - :param dataCube: RSF filename, Seismic data cube to interpolate - :param interpolated: RSF filename, Interpolated seismic data cube - :param nm: integer, number of CMPs in the seismic data cube - :param dm: float, CMP sampling - :param nt: integer, number of time samples - :param dt: float, time sampling - :param nhi: integer, number of constant offsets gathers to interpolate - :param a1: integer, Number of PEF coeficients in time axis - :param a2: integer, Number of PEF coeficients in space axis - :param rect1: integer, Smooth radius in time - :param rect2: integer, Smooth radius in space - ''' - - # Divide CMP sampling - dm = dm/2 - - # Define mask file names using input filename - mask1 = dataCube+'-mask1' - mask = dataCube+'-mask' - aa = dataCube+'-aa' - bb = dataCube+'-bb' - a = dataCube+'-a' - b = dataCube+'-b' - zeroTraceGather = dataCube+'-zeroedGather' - mask0 = dataCube+'-mask0' - - - # Build a mask to interleave zero traces with original data traces - Flow(aa,None,'spike n1=%i d1=%g o1=0' %(nm,dm)) - Flow(bb,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask1,[bb, aa], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - - Flow(a,None,'spike n1=%i d1=%g o1=0' % (nm,dm)) - Flow(b,None,'spike n1=%i d1=%g o1=0 mag=0' % (nm,dm)) - Flow(mask,[a, b], - ''' - interleave axis=1 ${SOURCES[1]} | - dd type=int - ''') - Flow(zeroTraceGather,b, - ''' - spray axis=2 n=%i d=%g | - transp | - put label2=Offset unit2=Km label1=Time unit1=s - ''' %(nt,dm)) - - # Data Mask with double of traces in CMP (half of CMP sampling) - # Keep the same Time and Offset original data sampling - Flow(mask0,mask, - ''' - spray axis=1 n=%i d=%g - ''' %(nt,dt)) - - totalPefIterations = 100 - totalInterpolationIterations = 50 - - offsetGathers = [] - for offsetGatherIndex in range(nhi): - - offsetGather = dataCube+"-offsetGather-%i" % offsetGatherIndex - resampledOffsetGather = dataCube+"-resampledGather-%i" % offsetGatherIndex - interpolatedOffsetGather = dataCube+"-interpolatedGather-%i" % offsetGatherIndex - pefCoeficients = dataCube+"-pefCoeficients-%i" % offsetGatherIndex - - Flow(offsetGather,dataCube, - ''' - window n2=1 f2=%i - ''' % (offsetGatherIndex)) - - Flow(resampledOffsetGather,[offsetGather,zeroTraceGather], - ''' - interleave axis=2 ${SOURCES[1]} - ''') - - # Calculate adaptive PEF coeficients - Flow(pefCoeficients,[resampledOffsetGather,mask0], - ''' - apef jump=2 a=%i,%i rect1=%i rect2=%i niter=%g verb=y - maskin=${SOURCES[1]} - ''' % (a1,a2,rect1,rect2,totalPefIterations)) - - # Interpolation - Flow(interpolatedOffsetGather, [resampledOffsetGather,pefCoeficients,mask0,mask1], - ''' - miss4 exact=y filt=${SOURCES[1]} mask=${SOURCES[2]} niter=%g verb=y | - put d2=%g - ''' % (totalInterpolationIterations,dm)) - - offsetGathers.append(interpolatedOffsetGather) - - # Concatenate interpolated sections - Flow(interpolated,offsetGathers, - ''' - rcat axis=3 ${SOURCES[1:%d]} | - transp plane=23 - ''' % nhi) - - diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/plot.py b/experiments/multiLayerModel/cre/madagascarRecipes/plot.py deleted file mode 100644 index 3715a31..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/plot.py +++ /dev/null @@ -1,52 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# plot.py (Madagascar Recipe) -# -# Purpose: Recipe to generate plot files (VPL). -# -# Important!: It should be called from a SConstruct -# -# Site: http://www.dirackslounge.online -# -# Version 1.0 -# -# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -# Selfdoc string -''' -Madagascar recipe to plot - -Define functions to generate vpl files. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -# Madagascar package -from rsf.proj import * - -__author__="Rodolfo Dirack " -__version__="1.0" - -def velplot(title,label1='Depth',unit1='Km'): - ''' - Velocity ploting function - :param title: RSF filename, velocity file - :param label1: string, first axis label - :param unit1: string, first axis unit - ''' - - return ''' - grey color=j allpos=y title="%s" scalebar=y - barlabel=Velocity barunit=Km/s - label1="%s" unit1="%s" label2=Lateral unit2=Km - barreverse=y pclip=100 bias=2.8 - ''' % (title,label1,unit1) - - diff --git a/experiments/multiLayerModel/cre/madagascarRecipes/velocityAnalisys.py b/experiments/multiLayerModel/cre/madagascarRecipes/velocityAnalisys.py deleted file mode 100644 index 96a5826..0000000 --- a/experiments/multiLayerModel/cre/madagascarRecipes/velocityAnalisys.py +++ /dev/null @@ -1,100 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- -# -# velocityAnalisys.py () -# -# Purpose: Madagascar recipe to do velocity analisys. -# -# Site: https://dirack.github.io -# -# Version 1.0 -# -# Programer: Rodolfo A C Neves (Dirack) 01/04/2020 -# -# Email: rodolfo_profissional@hotmail.com -# -# License: GPL-3.0 . - -from rsf.proj import * -import math - -# Selfdoc string -''' -Madagascar recipe to Velocity analisys with automatic velocity scan and picking - -Define functions to do a velocity analisys of a (time,offset,CMP) datacube. -''' - -if __name__ == "__main__": - print(__doc__) - exit() - -# Madagascar package -from rsf.proj import * - -__author__="Rodolfo Dirack " -__version__="1.0" - -def velocityAnalisys( - dataCube, - pick, - nmo, - vrms, - v0, - dv, - nv, - vel0, - rect1=15, - rect2=40, - rect3=3, - dt=0.004) - ''' - Velocity analisys and automatic picking and semblance - :param datacube: RSF filename, (time, offset, CMP) datacube - :out pick: RSF filename, picking velocity - :out stack: RSF filename, stacked section - :out vrms: RSF filename, vrms velocity - :param dt: float, time sampling - :param v0: float, inicial scan velocity - :param dv: float, scan velocity sampling - :param nv: int, number of scan velocities - :param vel0: float, initial velocity picking - :param rect1: int, smooth time radius to velocity scan - :param rect2: int, smooth space radius to velocity scan - :param rect3: int, smooth time radius post picking - ''' - - scan = 'scan'+dataCube - nmo = 'nmo'+dataCube - - # Velocity scan - Flow(scan,[dataCube], - ''' - vscan semblance=y - v0=%g dv=%g nv=%d - ''' % (v0,dv,nv)) - - # Velocity picking - Flow(pick,scan, - ''' - pick rect1=%i rect2=%i vel0=%g smooth=y | - smooth rect1=%i - ''' %(rect1,rect2,vel0,rect3)) - - # NMO and stack - Flow(nmo,[dataCube,pick], - 'nmo half=n velocity=${SOURCES[1]}') - - Flow(stack,nmo,'stack') - - # generate RMS velocity : - # sfmul multiply input data by itself to make v**2 - # sfcausint integrate v**2 (this does not include scaling by dt) - # sfmath scale by dt and divide by t - Flow(vrms,pick, - ''' - sfmul $SOURCE | - sfcausint | - sfmath output="sqrt(input*%g/(x1+%g))" | - sfput n3=1 d3=1 o3=0 - '''%(dt,dt)) diff --git a/experiments/test_experiments.py b/experiments/test_experiments.py new file mode 100644 index 0000000..270cfc3 --- /dev/null +++ b/experiments/test_experiments.py @@ -0,0 +1,54 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- +# +# test_experiments.py (Python) +# +# Purpose: Test model building with SCons. +# +# Rum with: python -m unittest test_experiments +# +# Site: https://dirack.github.io +# +# Version 2.0 +# +# Programer: Rodolfo A. C. Neves (Dirack) 04/03/2020 +# +# Email: rodolfo_profissional@hotmail.com +# +# License: GPL-3.0 . + +import unittest +import subprocess + +class TestScons(unittest.TestCase): + ''' + Test of the experiments + ''' + + def test_fullInterpolationAndStack(self): + ''' + Unit test of the fullInterpolationAndStack experiment + ''' + result = subprocess.call('cd fullInterpolationAndStack && scons -j3',shell=True) + self.assertEqual(result,0,'FAILED: experiment fullInterpolationAndStack failed!') + + def test_fullInterpolationAndStackCDS(self): + ''' + Unit test of the fullInterpolationAndStackCDS experiment + ''' + result = subprocess.call('cd fullInterpolationAndStackCDS && scons -j3',shell=True) + self.assertEqual(result,0,'FAILED: experiment fullInterpolationAndStackCDS failed!') + + def test_multiLayerModel(self): + ''' + Unit test of the multiLayerModel experiment + ''' + result = subprocess.call('cd multiLayerModel/cre && scons -j3',shell=True) + self.assertEqual(result,0,'FAILED: experiment multiLayerModel failed!') + + def test_multiLayerModelCDS(self): + ''' + Unit test of the multiLayerModel CDS experiment + ''' + result = subprocess.call('cd multiLayerModel/cds && scons -j3',shell=True) + self.assertEqual(result,0,'FAILED: experiment multiLayerModel CDS failed!')