Difference between revisions of "UCVM etree for Garner Valley"

From SCECpedia
Jump to navigationJump to search
Line 105: Line 105:
  
 
{|
 
{|
| [[File:cvmsi_gtl_garnervalley_depth_cross.png|thumb|300px|cvmsi, gtl]]
+
| [[File:cvmsi_x_garnervalley_depth_cross.png|thumb|300px|cvmsi, gtl??]]
| [[File:cvmsi_gtl_garnervalley_elevation_cross.png|thumb|300px|cvmsi, gtl]]
+
| [[File:cvmsi_x_garnervalley_elevation_cross.png|thumb|300px|cvmsi, gtl??]]
| [[File:cvmsi_gtl_garnervalley_horizontal.png|thumb|300px|cvmsi, gtl]]
+
| [[File:cvmsi_x_garnervalley_horizontal.png|thumb|300px|cvmsi, gtl??]]
 
|}
 
|}
  

Revision as of 16:31, 15 March 2019

Garner Valley

Generate etree with different parameters for Garner Valley

ROI: Garner Valley

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
Elnaz1,cvmh,gtl,5,10,200
Elnaz1,cvmh,nogtl,5,10,200
Elnaz2,cvmh,gtl,5,10,150
Elnaz2,cvmh,nogtl,5,10,150
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

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 XM)

For test5_garnervalley_cvmsi_gtl_0.5hz_10.0pts_1000ms,

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


Validation

Using UCVM, plot horizontal 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
cvm5, nogtl
cvm5, nogtl
cvm5, nogtl
cvmh, gtl
cvmh, gtl
cvmh, gtl
cvm5, gtl
cvm5, gtl
cvm5, 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