source: src/compute_metrics_csv.sh @ 108

#!/bin/bash

##############################################################################
# Computing metrics from LES and SCM
# Result in csv format
# Auteur: F Couvreux, R Honnert, F Hourdin, N Villefranque, C Rio, n'co
##############################################################################
#
# metrics are specified in list_case
# metrics names follow the syntax
#
# CASE_SUBCASE_METRICS_T1_T2
# T1 and T2 are initial and final time step for time averages
#
# available METRICS :
# ===================
#
# 1/ zav-400-600-var  -> variable "var" averaged between 400 and 600 m
#
# 2/ Ay-var -> integral ( min ( var -var(t=0) ) dz ) / integral ( dz )
#           integral taken from 0 to  zmax
#
# 3/ nebzave, neb2zave, neb4zave : Effective cloud height 
#           = int ( neb^p z dz ) / int ( neb^p dz ) with  p=1, 2 or 4
#
# 4/ nebmax : maximum cloud fraction on the column
#
# 5/ nebzmin, nebzmax : minimum/maximum cloud height
#
# 6/ lwp : liquid water path
#
# TBD :
# =====
# 1/ integrals are computed assuming rho=1 because rho is not systematically avalble
# 2/ the time average is coded for zav metrics only
#
##############################################################################

nWAVE=1


list_case=""
while (($# > 0)) ; do
   case $1 in
      -wave) nWAVE=$2 ; echo WAVE $nWAVE ;  shift ; shift ; echo OPTION $* ;;
      -help|-h) echo Usage "$0 [-wave N] [metrics1] [metrics2] ..." ; exit ;;
      *) list_case=( ${list_case[*]} $1 ) ; shift ;;
   esac
done

# Default metrics
if [ $list_case = "" ] ; then list_case=( ARMCU_REF_zav-400-600-theta_9_9 ) ; fi
        
wavedir=WAVE$nWAVE
NAME_OUT=$nWAVE
starting=1
rm -f metrics_$NAME_OUT.csv metrics_LES_$NAME_OUT.csv 

for name in ${list_case[@]};do
    for dir in LES ${wavedir} ; do
         echo ./extract_onemetric_csv.sh $name $dir
         ./extract_onemetric_csv.sh $name $dir
         #========================
         # Concatenation
         #========================
         if [ $starting == 1 ];then
            cp -f ${dir}_$name.csv metrics_${dir}_$NAME_OUT.csv
         else
            cut -d, -f2 ${dir}_$name.csv | paste -d, metrics_${dir}_$NAME_OUT.csv - > temp
            mv temp metrics_${dir}_$NAME_OUT.csv
         fi
    done
    starting=0
done # end loop $name

Rscript htune_csv2Rdata.R