Difference between revisions of "UCVM cvmsi tapering for CyberShake Study 22.2"

From SCECpedia
Jump to navigationJump to search
 
(8 intermediate revisions by the same user not shown)
Line 27: Line 27:
 
the .firstOrSecond files are the ones that is needed for CyberShake study
 
the .firstOrSecond files are the ones that is needed for CyberShake study
  
=== Plot the binary file ===
+
== Plot the binary file ==
 
To plot these files with ucvm_plotting, data ingest mode of plot_z10_map.py and plot_z25_map.py from ucvm_plotting is used  
 
To plot these files with ucvm_plotting, data ingest mode of plot_z10_map.py and plot_z25_map.py from ucvm_plotting is used  
 
but first needs to change those file name to have .binary file name so plotting scripts know how to parse the data
 
but first needs to change those file name to have .binary file name so plotting scripts know how to parse the data
Line 60: Line 60:
 
</pre>
 
</pre>
  
 +
{|
 +
| [[File:cvmsi_taper_z1.0.png|thumb|300px|CVMSI, taper, z1.0]]
 +
| [[File:cvmsi_taper_z2.5.png|thumb|300px|CVMSI, taper, z2.5]]
 +
|}
  
 +
{|
 +
| [[File:cvmsi_z1.0.png|thumb|300px|CVMSI, no taper, z1.0]]
 +
| [[File:cvmsi_z2.5.png|thumb|300px|CVMSI, no taper, z2.5]]
 +
|}
  
=== Create difference map between 2 binary files ===
+
== Create difference map between 2 binary files ==
  
 
<pre>
 
<pre>
plot_vs30_etree_difference_map.py -s 0.005 -c cvmsi -a d -i /var/www/html/UCVM_web_test/web/model/UCVM_TARGET -b 33.25,-119.38 -u 35.08,-116.49 -o diff_taper_none_z1.0.png -f basin_query_mpi_complete_cvmsi_taper_z1.0x.binary,basin_query_mpi_complete_cvmsi_z1.0x.binary -x 578 -y 366 -t "diff between cvmsi and cvmsi+taper Z1.0" -A -0.52,0.52
+
plot_vs30_etree_difference_map.py -s 0.005 -c cvmsi -a d -i ${UCVM_INSTALL_PATH} -b 33.25,-119.38 -u 35.08,-116.49 -o diff_taper_none_z1.0.png -f cvmsi_taper_z1.0.binary,cvmsi_z1.0.binary -x 578 -y 366 -t "diff between cvmsi and cvmsi+taper Z1.0" -A -0.52,0.52
 +
 
 +
plot_vs30_etree_difference_map.py -s 0.005 -c cvmsi -a d -i ${UCVM_INSTALL_PATH} -b 33.25,-119.38 -u 35.08,-116.49 -o diff_taper_none_z2.5.png -f cvmsi_taper_z2.5.binary,cvmsi_z2.5.binary -x 578 -y 366 -t "diff between cvmsi and cvmsi+taper Z2.5" -A -0.52,0.52
 +
</pre>
 +
 
 +
An example of the output :
 +
 
 +
<pre>
 +
Using parameters:
 +
 
 +
ny = 366
 +
lat1 = 33.25
 +
datafile2 = cvmsi_z1.0.binary
 +
datafile1 = cvmsi_taper_z1.0.binary
 +
lat2 = 35.08
 +
color = d
 +
title = diff between cvmsi and cvmsi+taper Z1.0
 +
spacing = 0.005
 +
nx = 578
 +
outfile = diff_taper_none_z1.0.png
 +
cvm = cvmsi
 +
lon1 = -119.38
 +
scalemax = 0.52
 +
scalemin = -0.52
 +
lon2 = -116.49
 +
installdir = /var/www/html/UCVM_web/web/model/UCVM_TARGET
 +
 
 +
Retrieving data. Please wait...
 +
 
 +
Using --> cvmsi_taper_z1.0.binary
 +
TOTAL number of binary data read:211548
 +
 
 +
 
 +
Using --> cvmsi_z1.0.binary
 +
TOTAL number of binary data read:211548
 
</pre>
 
</pre>
 +
 +
 +
{|
 +
| [[File:diff_taper_none_z1.0.png|thumb|300px|difference map, CVMSI, taper/none, z1.0]]
 +
| [[File:diff_taper_none_z2.5.png|thumb|300px|difference map, CVMSI, taper/none, z2.5]]
 +
|}

Latest revision as of 01:31, 2 March 2023

Generate Z1.0 and Z2.5

on Discovery, commands used to create Z1.0 and Z2.5 basin query results,

#!/bin/bash
if [ -z "$UCVM_INSTALL_PATH" ]; then
  echo "Need to set UCVM_INSTALL_PATH to run >" ${0##*/} 
  exit
fi
source $UCVM_INSTALL_PATH/conf/ucvm_env.sh

BIN_DIR=${UCVM_INSTALL_PATH}/bin
CONF_DIR=${UCVM_INSTALL_PATH}/conf

TEST=basin_query_mpi_complete_cvmsi_taper_z2.5
salloc ${UCVM_SALLOC_ENV} -Q --nodes=2 --ntasks=4 --mem=16GB --time=01:00:00 srun -Q -o ${TEST}.srun.out ${BIN_DIR}/basin_query_mpi_complete -b ${TES
T}.first,${TEST}.firstOrSecond,${TEST}.last,${TEST}.secondOnly,${TEST}.threeLast -o ${TEST}.result,${TEST}.meta.json -f ${CONF_DIR}/ucvm.conf -m cvms
i,elygtl:taper -i 10 -v 2500 -l 33.25,-119.38 -s 0.005 -x 578 -y 366

TEST=basin_query_mpi_complete_cvmsi_taper_z1.0
salloc ${UCVM_SALLOC_ENV} -Q --nodes=2 --ntasks=4 --mem=16GB --time=01:00:00 srun -Q -o ${TEST}.srun.out ${BIN_DIR}/basin_query_mpi_complete -b ${TES
T}.first,${TEST}.firstOrSecond,${TEST}.last,${TEST}.secondOnly,${TEST}.threeLast -o ${TEST}.result,${TEST}.meta.json -f ${CONF_DIR}/ucvm.conf -m cvms
i,elygtl:taper -i 10 -v 1000 -l 33.25,-119.38 -s 0.005 -x 578 -y 366

the .firstOrSecond files are the ones that is needed for CyberShake study

Plot the binary file

To plot these files with ucvm_plotting, data ingest mode of plot_z10_map.py and plot_z25_map.py from ucvm_plotting is used but first needs to change those file name to have .binary file name so plotting scripts know how to parse the data

plot_z10_map.py -a s -s 0.0050 -b 33.25,-119.38 -u 35.08,-116.49 -c cvmsi,elygtl:taper -x 578 -y 366 -f cvmsi_taper_z1.0.binary -o cvmsi_taper_z1.0.png
plot_z25_map.py -a s -s 0.005 -b 33.25,-119.38 -u 35.08,-116.49 -c cvmsi,elygtl:taper -x 578 -y 366 -f cvmsi_taper_z2.5.binary -o cvmsi_taper_z2.5.png 

plot_z10_map.py -a s -s 0.0050 -b 33.25,-119.38 -u 35.08,-116.49 -c cvmsi -x 578 -y 366 -f cvmsi_z1.0.binary -o cvmsi_z1.0.png
plot_z25_map.py -a s -s 0.005 -b 33.25,-119.38 -u 35.08,-116.49 -c cvmsi -x 578 -y 366 -f cvmsi_z2.5.binary -o cvmsi_z2.5.png 

An example of the output,

Using parameters:

datafile = cvmsi_taper_z1.0.binary
lat1 = 33.25
lat2 = 35.08
color = s
outfile = cvmsi_taper_z1.0.png
spacing = 0.0050
nx = 578
ny = 366
cvm = cvmsi,elygtl:taper
lon1 = -119.38
lon2 = -116.49

Retrieving data. Please wait...
TOTAL number of binary data read:211548
CVMSI, taper, z1.0
CVMSI, taper, z2.5
CVMSI, no taper, z1.0
CVMSI, no taper, z2.5

Create difference map between 2 binary files

plot_vs30_etree_difference_map.py -s 0.005 -c cvmsi -a d -i ${UCVM_INSTALL_PATH} -b 33.25,-119.38 -u 35.08,-116.49 -o diff_taper_none_z1.0.png -f cvmsi_taper_z1.0.binary,cvmsi_z1.0.binary -x 578 -y 366 -t "diff between cvmsi and cvmsi+taper Z1.0" -A -0.52,0.52

plot_vs30_etree_difference_map.py -s 0.005 -c cvmsi -a d -i ${UCVM_INSTALL_PATH} -b 33.25,-119.38 -u 35.08,-116.49 -o diff_taper_none_z2.5.png -f cvmsi_taper_z2.5.binary,cvmsi_z2.5.binary -x 578 -y 366 -t "diff between cvmsi and cvmsi+taper Z2.5" -A -0.52,0.52

An example of the output :

Using parameters:

ny = 366
lat1 = 33.25
datafile2 = cvmsi_z1.0.binary
datafile1 = cvmsi_taper_z1.0.binary
lat2 = 35.08
color = d
title = diff between cvmsi and cvmsi+taper Z1.0
spacing = 0.005
nx = 578
outfile = diff_taper_none_z1.0.png
cvm = cvmsi
lon1 = -119.38
scalemax = 0.52
scalemin = -0.52
lon2 = -116.49
installdir = /var/www/html/UCVM_web/web/model/UCVM_TARGET

Retrieving data. Please wait...

Using --> cvmsi_taper_z1.0.binary
TOTAL number of binary data read:211548


Using --> cvmsi_z1.0.binary
TOTAL number of binary data read:211548


difference map, CVMSI, taper/none, z1.0
difference map, CVMSI, taper/none, z2.5