source: lmdz_wrf/trunk/tools/model_graphics_template.dat @ 1333

### 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/lfita/etudes/WRF_LMDZ/svn/LMDZ_WRF/tools
# cdo HOME
cdoHOME = /usr/bin

# Srcatch: should be any previous work removed and start from the scratch? (overwrites all 'scratch' values)
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 = false

# Scratch differences: should be any experiment differences previous work removed and start from the scratch?
diffscratch = false

# Scratch difference figures: should be any experiment difference figures previous work removed and start from the scratch?
figdiffscratch = false

# Scratch model differences: should be any model differences previous work removed and start from the scratch?
moddiffscratch = false

# Scratch model difference figures: should be any model difference figures previous work removed and start from the scratch?
figmoddiffscratch = false

# Scratch all model-experiments figures: should be any all model-experiments figures previous work removed and start from the scratch?
figallmodexpscratch = false

# Have new files, figures and differences been added and do they need to be processed? (overwrites all 'add' values)
addall = false

# 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 = true

# Have new experiment differences been added and do they need to be processed and drawn?
adddiffs = true

# Have new experiment difference figures been added and do they need to be drawn?
adddifffigures = true

# Have new model differences been added and do they need to be processed and drawn?
addmoddiffs = true

# Have new model difference figures been added and do they need to be drawn?
addmoddifffigures = true

# Have new all model-experiment figures been added and do they need to be drawn?
addallmodexpfigures = true

# Debug
debug = true

# Folder with the files
##ifold = /bdd/PCER/workspace/lfita/etudes/WRF_LMDZ/WaquaL/WRF_LMDZ
ifold = /homedata/lfita/etudes/WRF_LMDZ/WaquaL_HighRes

# Output folder
##ofold = /bdd/PCER/workspace/lfita/etudes/WRF_LMDZ/WaquaL/WRF_LMDZ
# ofold = /homedata/lfita/etudes/WRF_LMDZ/WaquaL_HighRes/analysis
ofold = /data/lfita/etudes/WRF_LMDZ/Waqual_HighRes/analysis

# Models
#   WRF: WRF
#   LMDZ: LMDZ model
#   WRF_LMDZ: LMDZ physics coupled to WRF
models=WRF:WRF_LMDZ:LMDZ

# Graphic characteristics for models
#   Values to be used later in the graphics but not in '2lines'
#    ':' separated list of [model]|[labmodel]|[color]|[linetype]|[marker]|[sizes]
#     [model]: name of the model (equivalent to value in `models')
#     [labmodel]: label of the model as it should appear in the graphics
#     [color]: specific color of the model (and might related to all its experiments in lines plots)
#     [linetype]: specific type of line of the model (and might related to all its experiments in lines plots)
#     [marker]: specific marker of the model (and might related to all its experiments in lines plots)
#     [sizes]: line width and point size for the model (and might related to all its experiments in lines plots)
#     [tmodif]: time modification to apply to the model files ('None' for nothing)
#       'setorigin',[YYYYMMDDHHMISS]: re-set origin of times at [YYYYMMDDHHMISS]
#       
modgraphchar = WRF|wrf|red|-|,|2.|None:WRF_LMDZ|wlmdz|green|-|,|2.|None:LMDZ|lmdz|blue|-|,|2.|setorigin,19800119000000

# Experiments (labels of each experiment from each model)
WRFexps = current:micro1:micro2
LMDZexps = AR40:NPv31
WRF_LMDZexps = AR40:NPv31

# Graphic characteristics for models
#   Values to be used later in the graphics but not in '2lines'
#    ':' separated list of [model]|[exp]|[labexp]|[color]|[linetype]|[marker]|[sizes]
#     [model]: name of the model (equivalent to value in `models')
#     [exp]: name of the experiment (equivalent to value in `[model]exps')
#     [labexp]: label of the experiment as it should appear in the graphics
#   If any of the following values has the value 'asmodel', the equivalent value for [model] will be taken
#     [color]: specific color of the model
#     [linetype]: specific type of line of the model
#     [marker]: specific marker of the model
#     [sizes]: line width and point size for the model
expgraphchar = WRF|current|wcurr|asmodel|asmodel|.|asmodel:WRF|micro1|wmp1|asmodel|asmodel|x|asmodel:WRF|micro2|wmp2|asmodel|asmodel|+|asmodel:LMDZ|AR40|lmdza|asmodel|asmodel|.|asmodel:LMDZ|NPv31|lmdzb|asmodel|asmodel|x|asmodel:WRF_LMDZ|AR40|wlmdza|asmodel|asmodel|.|asmodel:WRF_LMDZ|NPv31|wlmdzb|asmodel|asmodel|x|asmodel

# 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
#   tstd: temporal standard deviation values
#   tturb: Taylor's turbulence decomposition value (x - <x>) for time
#   tvar: temporal variance values
#   xmean: x-axis mean values
#   xvar: x-axis variance 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 = uas:vas
#VAR_tmean = hfss:hfls
#VAR_tstd = tas:wss
#VAR_tvar+xmean = tas:wss:hurs
#VAR_xvar+ymean = tas:wss:hurs
VAR_tturb+xmean+last = tas:wss:hurs
VAR_tturb+xmean+tmean = tas:wss:hurs
VAR_xmean = uas:vas:wss:tas:pr:clt:cll:clm:clh:hurs:prw:hfls:hfss:rsds:zmla
VAR_tmean = uas:vas:wss:tas:pr:clt:cll:clm:clh:hurs:prw:hfls:hfss:rsds:zmla
VAR_xmean+tmean = uas:vas:wss:tas:pr:clt:cll:clm:clh:hurs:prw:hfls:hfss:rsds:zmla
##VAR_last+pinterp+xmean = ua:va:hur:hus:turb
VAR_pinterp+tmean+xmean = ua:va:wa:hur:hus:ta

# ':' seaprated list of operations which can not be re-projected (due to the lack of both coordinates: 'lon' and 'lat')
NOreprojops = lmean:Lmean:xmean:ymean:turb

# 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
#      e.g.: once 'tmean', or 'xmean' are computed, variable [varn] will become [varn]mean
opsur_mean = Lmean:lmean:tmean:xmean:ymean
opsur_std = tstd
opsur_turb = tturb
opsur_var = tvar:xvar

# 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:tstd|tstd:tturb|tturb:tvar|tvar:xmean|xmean:xvar|xvar: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.

# Variables reprojection to compute differences among models (which might use different projections) 
# ':' separated list of models to use as reference for the common projection
RefProj = WRF:WRF_LMDZ

# Reprojection methodologies using `RefProj' as common projection
#
# CDO remapping option to be used for each variable
#  ':' separated list for each of CDO's operators (version depending) to use for each variable
#  List of available remapping operators from CDO v 1.6.4rc6
#    remapbil   Bilinear interpolation
#               Performs a bilinear interpolation on all input fields.
#               This interpolation method only works on quadrilateral curvilinear grids.
#    remapbic   Bicubic interpolation
#               Performs a bicubic interpolation on all input fields.
#               This interpolation method only works on quadrilateral curvilinear grids.
#    remapdis   Distance-weighted average remapping
#               Performs a distance-weighted average remapping of the four nearest neighbor 
#               values on all input fields.
#    remapnn    Nearest neighbor remapping
#               Performs a nearest neighbor remapping on all input fields.
# Operators only available if projections have the corners of the grid points.
#    remapcon   First order conservative remapping
#               Performs a first order conservative remapping on all input fields.
#    remapcon2  Second order conservative remapping
#               Performs a second order conservative remapping on all input fields.
#    remaplaf   Largest area fraction remapping
#               Performs a largest area fraction remapping on all input fields.
####### ####### ####### or additionally ####### ####### #######
# python remapping
#   remapping using functions created on purpose
#    dis        Distance-weighted average remapping
#               Performs a distance-weighted average remapping of the four nearest neighbor 
#               values on all input fields.
#    pnn        Nearest neighbor remapping
#               Performs a nearest neighbor remapping on all input fields.
reprojectvar_dis = uas:vas:wss:pr:ua:va:wa:hur:hus:ta:hfls:hfss:zmla
reprojectvar_pnn = tas:hfls:hfss:hus:hurs:clt:cll:clm:clh:prw:rsds

# Differences to compute
#   These differences will be computed among experiments and among models (and model-experiments)
#  ':' separated list of [var]|[op]
#  DIFFOP_[op]: differences of the [op] of [var]. NOTE: [op] must exist for each experiment
#    Differences are computed directly from the last stage of the operation
#  DIFFVAR_[op]: [op] of the [var] differences [var]. NOTE: [var] must exist for each experiment
#    First are computed the differences from the initial variable file and then operations are made
# NOTE: remember that: meanvar2 - meanvar1 = mean(var2 - var1)
# NOTE2: For that variables with an opertion from which coordinates 'lon', 'lat' are reduced
#   One cannot compute model differences as 'DIFFOP' since there is no way to re-project 
#   data which has already lost one of the coordinates variable
#
DIFFOP_last = uas:vas
DIFFOP_tmean = uas:vas:wss:tas:pr:clt:cll:clm:clh:hurs:prw:hfls:hfss:rsds:zmla

DIFFVAR_tturb+xmean = tas:wss:hurs
DIFFVAR_tturb+xmean+last = tas:wss:hurs
DIFFVAR_tturb+xmean+tmean = tas:wss:hurs
DIFFVAR_xmean = uas:vas:tas:wss:pr:clt:cll:clm:clh:hurs:prw:hfls:hfss:rsds:zmla
DIFFVAR_xmean+tmean = uas:vas:wss:tas:pr:clt:cll:clm:clh:hurs:prw:hfls:hfss:rsds:zmla
DIFFVAR_pinterp+tmean+xmean = ua:va:wa:hur:hus:ta

# 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
#   2lines: 1D plot with 2 lines with values at each x/x2 or y/y2 axes
#   2lines_time: 1D plot with 2 time-lines with values at each x/x2 or y/y2 axes
#   Nlines: 1D plot with N lines (only for all models/experimenst figures)
#   Nlines_time: 1D plot with N time-lines (only for all models/experimenst figures)
#   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

# tmean
DIRPLT_shadcont2Dsfc = uas|tmean#vas|tmean:hfls|tmean#hfss|tmean:wss|tmean#tas|tmean:pr|tmean#clt|tmean:hurs|tmean#cll|tmean:clm|tmean#clh|tmean:prw|tmean#zmla|tmean
# xmean
DIRPLT_shadconthovmsfc = uas|xmean#vas|xmean:hfls|xmean#hfss|xmean:wss|xmean#tas|xmean:pr|xmean#clt|xmean:hurs|xmean#cll|xmean:clm|xmean#clh|xmean:prw|xmean#zmla|xmean
# xmean+tmean
DIRPLT_2lines = uas|xmean+tmean#vas|xmean+tmean:hfls|xmean+tmean#hfss|xmean+tmean:wss|xmean+tmean#tas|xmean+tmean:pr|xmean+tmean#clt|xmean+tmean:hurs|xmean+tmean#cll|xmean+tmean:clm|xmean+tmean#clh|xmean+tmean:prw|xmean+tmean#zmla|xmean+tmean
# pinterp+tmean+xmean
DIRPLT_shadcont2Dzsec = hus|pinterp+tmean+xmean#ta|pinterp+tmean+xmean:va|pinterp+tmean+xmean#ua|pinterp+tmean+xmean:wa|pinterp+tmean+xmean#hur|pinterp+tmean+xmean

# tmean
PLOTDIFFOP_shadcont2Dsfc = uas|tmean#vas|tmean:hfls|tmean#hfss|tmean:wss|tmean#tas|tmean:pr|tmean#clt|tmean:hurs|tmean#cll|tmean:clm|tmean#clh|tmean:prw|tmean#zmla|tmean

# xmean
DIFFPLOTVAR_shadconthovsfc = uas|xmean#vas|xmean:hfls|xmean#hfss|xmean:wss|xmean#tas|xmean:pr|xmean#clt|xmean:hurs|xmean#cll|xmean:clm|xmean#clh|xmean:prw|xmean#zmla|xmean
# xmean+tmean
DIFFPLOVAR_2lines = uas|xmean+tmean#vas|xmean+tmean:hfls|xmean+tmean#hfss|xmean+tmean:wss|xmean+tmean#tas|xmean+tmean:pr|xmean+tmean#clt|xmean+tmean:hurs|xmean+tmean#cll|xmean+tmean:clm|xmean+tmean#clh|xmean+tmean:prw|xmean+tmean#zmla|xmean+tmean
# pinterp+tmean+xmean
DIFFPLOTVAR_shadcont2Dzsec = hus|pinterp+tmean+xmean#ta|pinterp+tmean+xmean:va|pinterp+tmean+xmean#ua|pinterp+tmean+xmean:wa|pinterp+tmean+xmean#hur|pinterp+tmean+xmean

# All mod/exp
PLOTALLMODEXP_Nlines = uas|xmean+tmean:vas|xmean+tmean:hfls|xmean+tmean:hfss|xmean+tmean:wss|xmean+tmean:tas|xmean+tmean:pr|xmean+tmean:clt|xmean+tmean:hurs|xmean+tmean:cll|xmean+tmean:clm|xmean+tmean:clh|xmean+tmean:prw|xmean+tmean:zmla|xmean+tmean

# 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
#        [cntkind]: kind of contours 
#          'cmap': as it gets from colorcnt (which in this case must be a valid colorbar)
#          'fixc': fixed color [colname], all stright lines
#          'fixsigc': fixed color [colname], >0 stright, <0 dashed  line
#        [colorcnt]: color for the countor lines
#  NOTE: meaning of values change by lines-plots
#  ':' separated list of specific parameters for plotting as:
#      [varname]|[operation]|[minval]|[maxval]|[lcolor]|[lformat]|[lmark]|[wdthsz]
#        [varname]: name of the variable
#        [operation]: '+' separated list of operations
#        [figure]: figure to be used
#        [minval]: minimum value
#        [maxval]: minimum value
#        [lcolor]: name of the color for the line to use
#        [lformat]: format of the line ('-', straight, '-.', etc...)
#        [lmark]: marker of the line
#        [wdthsz]: width of line and size of marker
specificvarplot = cll|last|shadcont2Dsfc|0.|1.|gist_gray|%g|fixc|black:clh|last|shadcont2Dsfc|0.|1.|gist_gray|%g|fixc|black:husturbmean|pinterp@turb@xmean|shadcont2Dzsec|0.|1.e-6|Blues|%g|fixc|black:tas|tstd|shadcont2Dsfc|0.|10.|Reds|%g|fixc|black:wss|tstd|shadcont2Dsfc|0.|10.|Blues|%g|fixc|black:tas|tstd|shadcont2Dsfc|0.|100.|Reds|%g|fixc|black:wss|tstd|shadcont2Dsfc|0.|100.|Blues|%g|fixc|black:hfss|xmean|shadconthovmsfc|-40.|40.|Reds|%g|fixc|black:tas|tstd|shadcont2Dsfc|0.|1.5|Reds|%g|fixc|black:tas|tvar|shadcont2Dsfc|0.|1.5|Reds|%g|fixc|black:wss|tstd|shadcont2Dsfc|0.|8.|Greens|%g|fixc|black:wss|xvar+ymean|2linesTime|0.|15.|green|-|,|2.:tas|xvar+ymean|2linesTime|0.|2.|blue|-|,|2.
specificdiffopplot = hfls|tmean|shadcont2Dsfc|-800.|800.|seismic|%g|fixsigc|black:hfss|tmean|shadcont2Dsfc|-300.|300.|seismic|%g|fixsigc|black:hfls|last|shadcont2Dsfc|-800.|800.|seismic|%g|fixsigc|black:hfss|last|shadcont2Dsfc|-400.|400.|seismic|%g|fixsigc|black:uas|xmean|shadcont2Dsfc|-20.|20.|seismic|%g|fixsigc|black:vas|xmean|shadcont2Dsfc|-30.|30.|seismic|%g|fixsigc|black:ua|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|fixsigc|black:va|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|fixsigc|black:tas|tstd|shadcont2Dsfc|0.|10.|Reds|%g|fixsigc|black:wss|tstd|shadcont2Dsfc|0.|10.|Blues|%g|fixsigc|black:tas|tvar|shadcont2Dsfc|0.|4.|Reds|%g|fixsigc|black:wss|tvar|shadcont2Dsfc|0.|30.|Blues|%g|fixsigc|black:hfss|xmean|shadconthovmsfc|-40.|40.|Reds|%g|fixsigc|black:hus|pinterp+tmean+xmean|shadcont2Dzsec|-0.001|0.001|seismic|%g|fixsigc|black:ta|pinterp+tmean+xmean|shadcont2Dzsec|-5.|5.|seismic|%g|fixsigc|black:wss|xvar+ymean|2linesTime|0.|15.|green|-|,|2.:tas|xvar+ymean|2linesTime|0.|2.|blue|-|,|2.
specificdiffvarplot = hfls|tmean|shadcont2Dsfc|-800.|800.|seismic|%g|fixsigc|black:hfss|tmean|shadcont2Dsfc|-300.|300.|seismic|%g|fixsigc|black:hfls|last|shadcont2Dsfc|-800.|800.|seismic|%g|fixsigc|black:hfss|last|shadcont2Dsfc|-400.|400.|seismic|%g|fixsigc|black:uas|xmean|shadcont2Dsfc|-20.|20.|seismic|%g|fixsigc|black:vas|xmean|shadcont2Dsfc|-30.|30.|seismic|%g|fixsigc|black:ua|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|fixsigc|black:va|pinterp+tmean+xmean|shadcont2Dsfc|-20.|20.|seismic|%g|fixsigc|black:tas|tstd|shadcont2Dsfc|0.|10.|Reds|%g|fixsigc|black:wss|tvar|shadcont2Dsfc|0.|40.|Blues|%g|fixsigc|black:tas|tvar|shadcont2Dsfc|0.|4.|Reds|%g|fixsigc|black:wss|tstd|shadcont2Dsfc|0.|100.|Blues|%g|fixsigc|black:hfss|xmean|shadconthovmsfc|-40.|40.|Reds|%g|fixsigc|black:hus|pinterp+tmean+xmean|shadcont2Dzsec|-0.003|0.003|seismic|%g|fixsigc|black:ta|pinterp+tmean+xmean|shadcont2Dzsec|-80.|80.|seismic|%g|fixsigc|black:wss|xvar+ymean|2linesTime|0.|25.|green|-|,|2.:tas|xvar+ymean|2linesTime|0.|3.|blue|-|,|2.

# Figures output kind
kindfig = png

# Map to use
mapval = None

# 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,1,d
timefmt = $%d^{%H}$
timelabel = date!($[DD]^{[HH]}$)

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