Commit 85650b70 authored by Evans, Katherine's avatar Evans, Katherine
Browse files

add structure to acct for glc percent not turned on yet

parent 95ef9500
......@@ -46,8 +46,10 @@ ncwa -O -a lat,lon -w gris_area -B "gris_mask == 1" $outpath/lnd/climos/${case}_
elif [ "$version" == "e3sm" ]; then
# sum% glacier coverage over all elevations to create total % coverage
ncwa -N -a glc_nec -v PCT_GLC_MEC $outpath/lnd/climos/${case}_ANN_climo.nc $outpath/lnd/climos/${case}_ANN_climo_pct_glc.nc
echo "sum pct glacier coverage over all elevations to create total pct coverage"
ncwa -O -N -a glc_nec -v PCT_GLC_MEC $outpath/lnd/climos/${case}_ANN_climo.nc $outpath/lnd/climos/${case}_ANN_climo_pct_glc.nc
ncrename -v PCT_GLC_MEC,PCT_GLC_SUM $outpath/lnd/climos/${case}_ANN_climo_pct_glc.nc
climos_path=$outpath/lnd/climos
......
......@@ -46,13 +46,14 @@ fi
elif [ "$version" == "e3sm" ]; then
#for fl in $outpath/lnd/tseries/*clm2.h0.yrly*; do
#ncks -A -C -v gris_mask $mask_file $fl
#done
#echo "add pct coverage of glacier to file before averaging"
#ncks -A -v PCT_GLC_SUM $outpath/lnd/climos/${case}_ANN_climo_pct_glc.nc $outpath/lnd/tseries/$case.$lnd.${freq}.yrly.$ivar.nc
echo "average over ice sheet region for yearly averaged lnd data for $ivar"
ncwa -O -a lat,lon -w gris_area -B "gris_mask == 1" $outpath/lnd/tseries/$case.$lnd.${freq}.yrly.$ivar.nc $outpath/lnd/tseries/$case.$lnd.aavg.${freq}.yrly.$ivar.nc
#ncap2 -O -s 'areaice=PCT_GLC_SUM*0.01*gris_area' $outpath/lnd/tseries/$case.$lnd.${freq}.yrly.$ivar.nc -o $outpath/lnd/tseries/$case.$lnd.${freq}.yrly_ai.$ivar.nc
if test -a $outpath/lnd/tseries/${case}.$lnd.aavg.${freq}.yrly.$ivar.nc; then
echo "Global avg yearly masked $lnd averages of $ivar made"
else
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment