### ASCII file to manage `model_graphics.bash' from different models and experiments 
#  from each one
#
#  1: Variables choice using its CF name
#  2: Kind of statisitcs to apply to the variable
#  3: Kind of plot for the variable
#  4: Kind of difference among models/experiments
#
# Different assumptions are made:
#  - model outputs are organized as ${ifold}/${model}/${exp}
#  - script outputs will be organized as ${ofold}/${model}/${exp}
#  - when diagnostic might be required a new file will be created
#  - WRF outputs will have added 'WRFtime' variable
#
# NOTE: in this file all the values are ':' separated lists without '

# python HOME
pyHOME = /home/lluis/etudes/WRF_LMDZ/svn/LMDZ_WRF/tools
# cdo HOME
cdoHOME = /home/lluis/bin/gcc_cdo/1.6.4rc6/bin/cdo

# Srcatch: should be any previous work removed and start from the scratch?
scratch = false

# Scratch files: should be any files previous work removed and start from the scratch?
filescratch = false

# Scratch figures: should be any figures previous work removed and start from the scratch?
figscratch = true

# Have new files been added and do they need to be processed?
addfiles = false

# Have new figures been added and do they need to be drawn?
addfigures = false

# Debug
debug = false

# Folder with the files
##ifold = /bdd/PCER/workspace/lfita/etudes/WRF_LMDZ/WaquaL/WRF_LMDZ
ifold = /home/lluis/etudes/WRF_LMDZ/WaquaL_highres/tests/model_graphics

# Output folder
##ofold = /bdd/PCER/workspace/lfita/etudes/WRF_LMDZ/WaquaL/WRF_LMDZ
ofold = /home/lluis/etudes/WRF_LMDZ/WaquaL_highres/tests/model_graphics

# Experiments (model runs)
#   WRF: WRF
#   LMDZ: LMDZ model
#   WRF_LMDZ: LMDZ physics coupled to WRF
#models=WRF:WRF_LMDZ:LMDZ
models=WRF

# Experiments (labels of each experiment)
WRFexps = control:pbl1
LMDZexps = AR40:NPv31
WRF_LMDZexps = AR40:NPv31

# Headers of the files to use
WRFheaders = wrfout
WRF_LMDZheaders = wrfout
LMDZheaders = histins

# CF times characteristics
#   Reference time and units of the CF time for the output files
CFreftime = 19491201000000
CFunitstime = minutes

# Kind of calculations (adding the equivalent `surname' to each variable)
#   acc: temporal accumulated values
#   diff: differences between models
#   direct: no statistics
#   last: last temporal value
#   Lmean: latitudinal mean values
#   Lsec: latitudinal section (latitudinal value must be given, [var]@[lat]) 
#   lmean: longitudinal mean values
#   lsec: longitudinal section (longitudinal value must be given, [var]@[lat])
#   pinterp: pressure interpolation (to the given $plevels, see below)
#   tmean: temporal mean values
#   turb: Taylor's turbulence decomposition value
#   xmean: x-axis mean values
#   ymean: y-axis mean values
#   zsum: vertical aggregated values
# variables (providing CF-name, script must provide equivalent variable for each model)
#   in foudre values from: ${HOME}/UNSW-CCRC-WRF/tools/postprocess/GMS-UC/WRF4G/util/postprocess/wrfncxnj/wrfncxnj.table
#   must exist on the table $pyHOME/variables_values.dat
#
# Variables to compute (as possible multiple consecutive combination of operations must start by 'VAR_' separated by '+'
#   VAR_[calc1]+[calc2] = tas:wss
#     will compute first [calc1] and then [calc2] for 'tas' and 'wss'
VAR_last = tas:uas:vas:hfss:hfls:wss
VAR_tmean = tas:uas:vas:hfss:hfls:wss
VAR_xmean = tas:uas:vas:hfss:hfls:wss
#VAR_last+pinterp+xmean = ta:ua:va:hur:hus:turb
VAR_pinterp+tmean+xmean = ta:ua:va:hur:hus

#varlast = tas:uas:vas:ps:pr:pracc:rsds:prw:evspsbl:hfss:hfls:hurs:huss:zmla:hufs:wakes:lwakeh:stherm:zmaxth:clt:cll:clm:clh:prc:prls:bils:wss
#vartmean = tas:uas:vas:ps:pr:pracc:rsds:prw:evspsbl:hfss:hfls:hurs:huss:zmla:hufs:wakes:lwakeh:stherm:zmaxth:clt:cll:clm:clh:prc:prls:bils:wss
#varxmean = ta:ua:va:hur:hurs:tas:wss:ps
#varpinterp = ta:ua:va:hur:ws:hus
#varcombo = ta;pinterp@last@xmean:ua;pinterp@last@xmean:va;pinterp@last@xmean:hur;pinterp@last@xmean:ws;pinterp@turb@xmean:hus;pinterp@turb@xmean

vardiff = tas|last:tas|xmean:tas|tmean:wss|last:wss|xmean:wss|tmean:ua|pinterp@last@xmean:va|pinterp@last@xmean

# Operations which attach a 'surname' at the end of the variable name
#  Values are given as opsur_[surname] = [operation1]:[operation2][:...[operationN]]
#    opsur_mean = tmean:xmean
#      once 'tmean', or 'xmean' are computed, variable [varn] will become [varn]mean
opsur_mean = Lmean:lmean:tmean:xmean:ymean

# Text as title for each operation
#   as ':' separated list of [op]|[explanation '!' as spaces]
titleoperations = acc|tacc:diff|diff:direct|!:last|last:Lmean|latmean:Lsec|latsec:lmean|lonmean:lsec|lonsec:pinterp|pinterp:tmean|tmean:turb|turbulence:xmean|xmean:ymean|ymean:zsum|zsum

# Pressure levels
#   Common vertical pressure levels (in Pa) to vertically interpolate
#39Lowp#plevels = 100500.:100000.:98500.:97500.:95000.:92500.:90000.:87500.:85000.:80000.:75000.:70000.:65000.:60000.:55000.:50000.:45000.:40000.:35000.:30000.:25000.:20000.:15000.:10000.:900.:800.:700.:600.:500.:400.:300.:200.::1000.:85.:50.:40.:30.:20.:10.
#26#plevels = 100000.:97500.:95000.:92500.:90000.:85000.:80000.:75000.:70000.:65000.:60000.:55000.:50000.:45000.:40000.:35000.:30000.:25000.:20000.:15000.:10000.:5000.:2500.:1000.:500.:250.
plevels = 100000.:97500.:95000.:92500.:90000.:85000.:80000.:75000.:70000.:65000.:60000.:55000.:50000.:45000.:40000.:35000.:30000.:25000.:20000.:15000.:10000.:5000.:2500.:1000.:500.:250.

# Common projection
#   File with the information of the projection to be used for inter-model differences
#   File obtained from 'cdo griddes'
commonproj = WaqualCommon.griddes

# Kind of plots
#   diffmap2Dsfc: 2D map of surface differences values of 1 variable
#   diffmap2Dz: 2D map of 3D differences values of 1 variable
#   map2Dsfc: 2D map of surface values of 1 variable
#   shadconthovmsfc: Hovmoeller diagrams of 2 variables at the surface one in shadow and the other in contourn
#   shadcont2Dsfc: 2D map of shadow (1st variable) and countour (2nd variable) [stvar1]#[stvar2]
#   shadcont2Dzsec: 2D map of vertical section of 2 variables one in shadow and the other in contourn
#
#   ':' separated list of statitsics variable values are given as: [var]|['+' separated list of operations]
#     in figures with more than 1 variable, use '#' to separate them
# DIRPLT_: Figures which can be made directly without another model/experiment output

DIRPLT_shadcont2Dsfc = uas|last#vas|last:hfls|last#hfss|last:wss|last#tas|last
DIRPLT_shadconthovmsfc = wss|xmean#tas|xmean
DIRPLT_shadcont2Dzsec = ua|pinterp+tmean+xmean#va|pinterp+tmean+xmean:hus|pinterp+tmean+xmean#ta|pinterp+tmean+xmean

#pltshadcont2Dsfc = huss|tmean#tas|tmean:vas|tmean#uas|tmean:ps|tmean#wss|tmean:pr|tmean#rsds|tmean:prw|tmean#clt|tmean:hfls|tmean#hfss|tmean:evspsbl|tmean#bils|tmean:zmaxth|tmean#stherm|tmean:wakeh|tmean#wakes|tmean:zmaxth|tmean#stherm|tmean:prls|tmean#prc|tmean:clt|tmean#cll|tmean:clm|tmean#clh|tmean:hurs|tmean#zmla|tmean:huss|last#tas|last:vas|last#uas|last:ps|last#wss|tmean:pr|last#rsds|last:prw|last#clt|last:hfls|last#hfss|last:evspsbl|last#bils|last:zmaxth|last#stherm|last:wakeh|last#wakes|last:zmaxth|last#stherm|last:prls|last#prc|last:clt|last#cll|last:clm|last#clh|last:hurs|last#zmla|last
#pltshadconthovmsfc = hurs|xmean#tas|xmean:wss|xmean#ps|xmean
#pltshadcont2Dzsec = ua|pinterp@last@xmean#va|pinterp@last@xmean:hur|pinterp@last@xmean#ta|pinterp@last@xmean:ws|pinterp@turb@xmean#hus|pinterp@turb@xmean

pltdiffshadcont2Dsfc = tas|tmean#wss|tmean:wss|last#tas|last
pltdiffshadconthovmsfc = wss|xmean#tas|xmean
pltdiffshadcont2Dzsec = ua|pinterp@last@xmean#va|pinterp@last@xmean

# Specific variable-statistisc graphic parameters (if no value is given, they will be retrieved from 'variables_values.dat')
#  ':' separated list of specific parameters for plotting as:
#      [varname]|[operation]|[minval]|[maxval]|[colorbar]|[cntformat]|[colorcnt]
#        [varname]: name of the variable
#        [operation]: '+' separated list of operations
#        [figure]: figure to be used
#        [minval]: minimum value
#        [maxval]: minimum value
#        [colorbar]: name of the colorbar (from matplotlib) to use
#        [cntformat]: format of the contour labels
#        [colorcnt]: color for the countor lines
specificvarplot = cll|last|shadcont2Dsfc|0.|1.|gist_gray|%g|black:clh|last|shadcont2Dsfc|0.|1.|gist_gray|%g|black:husturbmean|pinterp@turb@xmean|shadcont2Dzsec|0.|1.e-6|Blues|%g|black
diffspecificvarplot = tas|tmean|shadcont2Dsfc|-20.|20.|seismic|%g|black:wss|tmean|shadcont2Dsfc|-30.|30.|seismic|%g|black:tas|last|shadcont2Dsfc|-20.|20.|seismic|%g|black:wss|last|shadcont2Dsfc|-30.|30.|seismic|%g|black

# Figures output kind
kindfig = pdf

# Map to use
mapval = lcc,l

# Time parameters
#  [timekind]; kind of output for the time ticks
#     'Nval': according to a given number of values as 'Nval',[Nval]
#     'exct': according to an exact number [Nunits] of time unit as 'exct',[tunit]; 
#       tunit= [Nunits],[tu]
#         [tu]= 'c': centuries, 'y': year, 'm': month, 'w': week, 'd': day, 'h': hour, 'i': minute, 's': second, 'l': milisecond
#  [timefmt]; desired format for the time ticks (combination of LaTeX and 'linux date formats')
#  [timelabel]; label at the graph ('!' for spaces)
timekind = exct,6,h
timefmt = $%d^{%H}$
timelabel = date!($[DD]^{[HH]}$)

# Generic
errmsg = ERROR--error--ERROR--error
warnmsg = WARNING--warning--WARNING--warning