UCVM etree for Garner Valley

From SCECpedia
Jump to navigationJump to search

Garner Valley

Generate etree with different parameters for Garner Valley

ROI: Garner Valley

with SVM

For mpi_Elnaz2_garnervalley_cvms5_svm_5hz_10pts_200ms,

config file etree(size:12G, compressed 3.5G)

For mpi_Elnaz2_garnervalley_cvms5_gtl_5hz_10pts_200ms,

config file etree(size:9.5G, compressed 4.1G)

For mpi_Elnaz2_garnervalley_cvms5_ely_5hz_10pts_200ms,

config file etree(size:9.5G, compressed 4.1G)

For mpi_Elnaz2_garnervalley_cvms5_nogtl_5hz_10pts_200ms,

config file etree(size:3.2G, compressed 1.5G)


For test3_garnervalley_cvms5_svm_0.5hz_10.0pts_1000ms,

config file etree(size:3.6M, compressed 1.8M)

For mpi_test3_garnervalley_cvms5_svm_0.5hz_10.0pts_1000ms,

config file etree(size:3.6M, compressed 1.8M)

For mpi_test3_garnervalley_cvms5_nogtl_0.5hz_10.0pts_1000ms,

config file etree(size:3.6M, compressed 1.8M)

Result

etree_input.txt used to setup the etree generation

##
## etree_input.txt
#
#label,MODEL,GTL,MAX_FREQ, PPWL, VS_MIN
test1,cvmh,gtl,0.5,10.0,1000
test1,cvmh,nogtl,0.5,10.0,1000
test2,cvmh,gtl,5,10.0,1000
test2,cvmh,nogtl,5,10.0,1000
test3,cvms5,gtl,0.5,10.0,1000
test3,cvms5,nogtl,0.5,10.0,1000
test4,cvms5,gtl,5,10.0,1000
test4,cvms5,nogtl,5,10.0,1000
test5,cvmsi,gtl,0.5,10.0,1000
test6,cvmsi,gtl,5,10.0,1000
Elnaz1,cvmh,gtl,5,10,200
Elnaz1,cvmh,nogtl,5,10,200
Elnaz2,cvm5,gtl,5,10,200
Elnaz2,cvm5,nogtl,5,10,200


For test1_garnervalley_cvmh_nogtl_0.5hz_10.0pts_1000ms,

config file etree(size:2.5M, compressed 372K)

For test2_garnervalley_cvmh_nogtl_5hz_10.0pts_1000ms,

config file etree(size:3.8G, compressed 362M)

For test1_garnervalley_cvmh_gtl_0.5hz_10.0pts_1000ms,

config file etree(size:5.4M, compressed 822K)

For test2_garnervalley_cvmh_gtl_5hz_10.0pts_1000ms,

config file etree(size:4.6G, compressed 456M)

For test3_garnervalley_cvms5_nogtl_0.5hz_10.0pts_1000ms,

config file etree(size:5.5M, compressed 2.3M)

For test4_garnervalley_cvms5_nogtl_5hz_10.0pts_1000ms,

config file etree(size:4.7G, compressed 2.2G) ecompact etree(size:3.1G, compressed 1.5G) ecoalesce etree(size:3.1G, compressed 1.5G)

For test3_garnervalley_cvms5_gtl_0.5hz_10.0pts_1000ms,

config file etree(size:6.1M, compressed 2.6M)

For test4_garnervalley_cvms5_gtl_5hz_10.0pts_1000ms,

config file etree(size:4.9G, compressed 2.3G)

For test5_garnervalley_cvmsi_gtl_0.5hz_10.0pts_1000ms,

config file etree(size:2.7M, compressed 886K)


For test6_garnervalley_cvmsi_gtl_5hz_10.0pts_1000ms,

config file etree(size:3.4M, compressed 1.5M)


For Elnaz2_garnervalley_cvms5_gtl_5hz_10pts_200ms,

config file etree(size:31G, compressed 12G) ecompact-etree(size:20G, compressed 7.3G)

For Elnaz2_garnervalley_cvms5_nogtl_5hz_10pts_200ms,

config file etree(size:4.9G, compressed 2.3G)

For Elnaz1_garnervalley_cvmh_gtl_5hz_10pts_200ms,

config file etree(size:9.2G, compressed 937M)

For Elnaz1_garnervalley_cvmh_nogtl_5hz_10pts_200ms,

config file etree(size:24G, compressed 2.0G)

Validation

Using UCVM, plot vertical profiles around

-s 34.15,-117.165

cvm5, baseline
cvm5, svm1d
cvm5, elygtl


Using UCVM, plot slices and cross sections around ROI,

plotting Vs with these 2 points:

 lat1  =  33.4022
 lon1  =  -116.8593
 lat2  =  33.7958
 lon2  =  -116.3984


cvmh, nogtl
cvmh, nogtl
cvmh, nogtl
cvms5, nogtl
cvms5, nogtl
cvms5, nogtl
cvmh, gtl
cvmh, gtl
cvmh, gtl
cvms5, gtl
cvms5, gtl
cvms5, gtl
cvmsi, gtl
cvmsi, gtl
cvmsi, gtl

Utility script

A sh script used to generate the config file and running ucvm2etree

#!/bin/sh  

#Request for etree
#From:   Esmaeilzadeh Seylabi, Elnaz <elnaz@caltech.edu>
#CC:     Ricardo Taborda Rios <rtaborda@eafit.edu.co>, Asimaki,
#        Domniki (Domniki) <domniki@caltech.edu>,
#        Doriam Leidin Restrepo Sánchez <drestre6@eafit.edu.co>


etree_go () {
  ${UCVM_INSTALL}/bin/ucvm2etree -f ${ETREE_CONFIG}
}

## generate the etree.conf
etree_foo() {

configfile=./${ETREE_CONFIG}

/bin/cat >$configfile << END 

#
# ${DATE}
#
# ${ETREE_CONFIG}
#

proj=geo-bilinear 
lon_0=-116.8593
lat_0=33.6403

lon_1=-116.6457
lat_1=33.7958

lon_2=-116.3984
lat_2=33.5572

lon_3=-116.6119
lat_3=33.4022

# Domain dimensions (meters):
x-size=26250.0000
y-size=35000.0000
z-size=4375.0000

# Max freq, points per wavelength, Vs min
max_freq=${MAX_FREQ}
ppwl=${PPWL}
vs_min=${VS_MIN}

# Blocks partition parameters:
nx=6
ny=8
 
# Max allowed size of octants in meters
max_octsize=10000.0

# Etree parameters and info
title=GarnerValley_${MAX_FREQ}hz_${PPWL}pts_${VS_MIN}ms
author=Mei_Hui_Su4Elnza
date=02/27/2019
outputfile=./${LABEL}_garnervalley_${MODEL}_${GTL}_${MAX_FREQ}hz_${PPWL}pts_${VS_MIN}ms.e
format=etree 

# UCVM parameters
ucvmstr=${MODEL}
ucvm_interp_zrange=0.0,350.0 
ucvmconf=${UCVM_INSTALL}/conf/ucvm.conf

# Scratch
scratch=./scratch

#
# Buffering parameters
#
# Etree buffer size in MB
buf_etree_cache=128
# Max octants to buffer for flat file during extraction
buf_extract_mem_max_oct=4194304
# Max octants to save in flat file before reporting full during extraction
buf_extract_ffile_max_oct=16000000
# Max octants to read from input flat file during sorting
buf_sort_ffile_max_oct=20000000
# Minimum number of octants between reports during merging

buf_merge_report_min_oct=10000000
# MPI send/recv octant buffer size during merging
buf_merge_sendrecv_buf_oct=4096
# Etree read/write octant buffer size during merging
buf_merge_io_buf_oct=4194304


END
}

##
## etree_input.txt
#
#label,MODEL,GTL,MAX_FREQ, PPWL, VS_MIN
#test1,cvmh,gtl,0.5,10.0,1000
#test1,cvmh,nogtl,0.5,10.0,1000
#test2,cvmh,gtl,5,10.0,1000
#test2,cvmh,nogtl,5,10.0,1000
#Elnaz1,cvmh,gtl,5,10,200
#Elnaz1,cvmh,nogtl,5,10,200
#Elnaz2,cvmh,gtl,5,10,150
#Elnaz2,cvmh,nogtl,5,10,150

UCVM_INSTALL="/home/scec-00/mei/ucvmc_18_5/install"
DATE=`date`

input="etree_input.txt"
while IFS= read -r file; do
  arrData=(${file//,/ })
  [[ $arrData =~ ^#.* ]] && continue
  LABEL=${arrData[0]}
  MODEL=${arrData[1]}
  GTL=${arrData[2]}
  MAX_FREQ=${arrData[3]}
  PPWL=${arrData[4]}
  VS_MIN=${arrData[5]}
  ETREE_CONFIG=${LABEL}_garnervalley_${MODEL}_${GTL}_${MAX_FREQ}hz_${PPWL}pts_${VS_MIN}ms.conf
  
  etree_foo
  etree_go

done < "$input"

Related Links