source: src/compute_metrics_csv.sh @ 176

#!/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=1) ) 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
wave=$nWAVE
starting=1

# Remove old files
rm -f metrics_$wave.csv metrics_LES_$wave.csv 

for name in ${list_case[@]};do
  # Where are the reference (target) nc files?
  if [ ${name:0:3} == RAD ] ; then 
    REF=RAD
  else 
    REF=LES 
  fi
  for dir in $REF ${wavedir} ; do
    echo ./extract_onemetric_csv.sh $name $dir
    ./extract_onemetric_csv.sh $name $dir
    #========================
    # Concatenation
    #========================
    if [ $dir = $REF ] ; then dirname=REF 
    else dirname=$dir ; fi
    if [ $starting == 1 ];then
      cp -f ${dir}_$name.csv metrics_${dirname}_$wave.csv
    else
      cut -d, -f2 ${dir}_$name.csv | paste -d, metrics_${dirname}_$wave.csv - > temp
      mv temp metrics_${dirname}_$wave.csv
    fi
  done
  starting=0
done # end loop $name

Rscript --vanilla htune_csv2Rdata.R $wave