|
Last change
on this file since 158 was
101,
checked in by htune, 7 years ago
|
|
New syntax for csv metrics
Fredho
##############################################################################
# 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 ( nebp z dz ) / int ( nebp 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
#
##############################################################################
|
|
File size:
678 bytes
|
| Line | |
|---|
| 1 | ############################################################################# |
|---|
| 2 | # Purpose : define the setup of the experiments, i.e.the case, subcase, |
|---|
| 3 | # target variables and wave number. |
|---|
| 4 | ############################################################################# |
|---|
| 5 | case_name="csvmm" |
|---|
| 6 | subcase_name="thBL" |
|---|
| 7 | case_name="ARMCU" |
|---|
| 8 | subcase_name="REF" |
|---|
| 9 | #if use of compute_metric*sh |
|---|
| 10 | #need to define case_name='csvAR79RIBOAY' & subcase_name="thBL" |
|---|
| 11 | targetvar=c("theta500") |
|---|
| 12 | plotvar="theta" # qv ou rneb |
|---|
| 13 | |
|---|
| 14 | # Exemple for multi-target experiment : |
|---|
| 15 | #targetvar= c("theta500","qv500","nebzave","nebmax","nebzmin","nebzmax") |
|---|
| 16 | nmetrique=length(targetvar) |
|---|
| 17 | |
|---|
| 18 | # Wave number |
|---|
| 19 | WAVEN=1 |
|---|
| 20 | |
|---|
| 21 | tau=1 |
|---|
| 22 | |
|---|
| 23 | Cutoff=3 |
|---|
Note: See
TracBrowser
for help on using the repository browser.
