Changeset 926 in lmdz_wrf for trunk/tools


Ignore:
Timestamp:
Jun 22, 2016, 8:00:03 PM (8 years ago)
Author:
lfita
Message:

Working version with first figure! `shadcount2Dsfc'

Location:
trunk/tools
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/tools/model_graphics.bash

    r924 r926  
    33
    44main='model_graphics.bash'
     5
     6# All purpouse
     7######## ###### ##### #### ### ## #
    58
    69function uploadvars() {
    … …  
    8083  fi
    8184}
     85
     86function list_add_label(){
     87# Function to add the label to a list
     88#   list= ':' separated list to wchich add a ticket
     89#   label= label to add to the list
     90#   ch= character to separate list values and new label
     91#   pos= position of the label ('beg': beginning, 'end': ending)
     92# $ list_add_label 1:2:3:4:5 abc # beg
     93# abc#1:abc#2:abc#3:abc#4:abc#5
     94
     95  fname='list_add_label'
     96
     97  list=$1
     98  label=$2
     99  ch=$3
     100  pos=$4
     101
     102  ival=1
     103  lvs=`echo ${list} | tr ':' ' '`
     104 
     105  newlist=''
     106  for listv in ${lvs}; do
     107    if test ${pos} = 'beg'; then
     108      nlv=${label}${ch}${listv}
     109    else
     110      nlv=${listv}${ch}${label}
     111    fi
     112    if test ${ival} -eq 1; then
     113      newlist=${nlv}
     114    else
     115      newlist=${newlist}':'${nlv}
     116    fi
     117    ival=`expr ${ival} + 1`
     118  done
     119
     120  echo ${newlist}
     121}
     122
     123function stats_filename() {
     124# Function to provide the name of a given statisitcs file
     125#   CFvarn= CF variable name
     126#   fkind= kind of the file
     127#   headf= head of the file
     128# $ stats_filename ta xmean wrfout
     129# ta_wrfout_xmean.nc
     130# $ stats_filename ta pinterp@last@xmean wrfout
     131# ta_wrfout_pinterp_last_xmean.nc
     132  fname='stats_filename'
     133
     134  CFvarn=$1
     135  fkind=$2
     136  headf=$3
     137
     138  # Number of operations in a combo file they are '@' separated
     139  nopers=`echo ${fkind} | tr '@' ' ' | wc -w | awk '{print $1}'`
     140
     141  if test $(index_string ${fkind} pinterp) -ne -1; then
     142    fhead=${headf}'p'
     143  else
     144    fhead=${headf}
     145  fi
     146
     147  if test ${nopers} -eq 1; then
     148    filename=${CFvarn}'_'${fhead}'_'${fkind}
     149  else
     150    opers=`echo ${fkind} | tr '@' ' '`
     151    filename=${CFvarn}'_'${fhead}
     152    for op in ${opers}; do
     153      filename=${filename}'_'${op}
     154    done
     155  fi
     156
     157  echo ${filename}.nc
     158}
     159
     160function index_string(){
     161# Function to provide the index of a given word inside a string
     162#   string= string to find in
     163#   word= word to find
     164#  $ index_string 123abc456ab789bca0 bca
     165#  14
     166  fname='index_string'
     167  string=$1
     168  word=$2
     169
     170  Lstring=`expr length ${string}`
     171  Lword=`expr length ${word}`
     172 
     173  Lfinalstring=`expr ${Lstring} - ${Lword}`
     174
     175  iw=0
     176  index=-1
     177  while test ${iw} -le ${Lfinalstring}; do
     178    if test ${string:${iw}:${Lword}} = ${word}; then
     179      index=${iw}
     180      break
     181    fi
     182    iw=`expr ${iw} + 1`
     183  done
     184
     185  echo ${index}
     186  return
     187}
     188
     189function list_filter() {
     190# Function to filter by a value a list
     191#   list: list to filter
     192#   value: value to filter with
     193#  $ list_filter '1|i:2|ii:3|iii:4|iv:5|v:1|I:2|II:3|III:4|iv:5|v' 3
     194#  3|iii:3|III
     195  fname='list_filter'
     196
     197  list=$1
     198  varn=$2
     199
     200  lvs=`echo ${list} | tr ':' ' '`
     201
     202  newlist=''
     203  il=1
     204  for lv in ${lvs}; do
     205    if test $(index_string ${lv} ${varn}) -ne -1; then
     206      if test ${il} -eq 1; then
     207        newlist=${lv}
     208      else
     209        newlist=${newlist}':'${lv}
     210      fi
     211      il=`expr ${il} + 1`
     212    fi
     213  done
     214
     215  echo ${newlist}
     216
     217  return
     218}
     219
     220
     221# Specific
     222######## ###### ##### #### ### ## #
    82223
    83224function variable_compute(){
    … …  
    313454        mv ${pyout} ${filen}
    314455
     456        # Keeping the operations
     457        pyins=${pyHOME}"/nc_var.py -f "${cf}" -o DataSetSection_multivars -v "${vs}
     458        pyins=${pyins}" -S "${values}
     459        echo " " >> ${odir}/all_computevars.inf
     460        echo "# ${CFvarn}" "${modvar}" >> ${odir}/all_computevars.inf
     461        echo ${pyins} >> ${odir}/all_computevars.inf
     462
    315463        pyout=`${pyHOME}/nc_var.py -f ${filen} -o chvarname -v ${modvar} -S ${CFvarn}`
    316464        pyn=$?
    … …  
    338486        ferrmsg ${pyn} ${fname} "python!'diagnostics'!failed#"${Spyout}
    339487        mv diagnostics.nc ${filen}
     488
     489        # Keeping the operations
     490        pyins=${pyHOME}"/diagnostics.py -f "${cf}" -d "${dims}" -v '"${diagn}"|"
     491        pyins=${pyins}${diagc}"'"
     492        echo " " >> ${odir}/all_computevars.inf
     493        echo "# ${CFvarn}" "${diagvar}" >> ${odir}/all_computevars.inf
     494        echo ${pyins} >> ${odir}/all_computevars.inf
    340495      fi
    341496
    … …  
    483638        ferrmsg ${pyn} ${fname} "python!'DataSetSection'!'last'!failed#"${Spyout}
    484639        mv ${cfiles:0:${Lcfiles3}}_time_B-9-E0-I0.nc ${filen}
     640
     641        # Keeping the operations
     642        pyins=${pyHOME}"/nc_var.py -o DataSetSection -S "${vals}" -f "${cfiles}
     643        pyins=${pyins}"-v "${CFvarn}
     644        echo " " >> ${odir}/all_statsvars.inf
     645        echo "# ${CFvarn}" "${vark}" >> ${odir}/all_statsvars.inf
     646        echo ${pyins} >> ${odir}/all_statsvars.inf
    485647      ;;
    486648      # latitudinal mean values
    … …  
    514676        cp pinterp.nc ${filen}
    515677
     678        # Keeping the operations
     679        pyins=${pyHOME}"/nc_var.py -o pinterp -S "${vals}" -f "${cfiles}
     680        pyins=${pyins}"-v "${CFvarn}
     681        echo " " >> ${odir}/all_statsvars.inf
     682        echo "# ${CFvarn}" "${vark}" >> ${odir}/all_statsvars.inf
     683        echo ${pyins} >> ${odir}/all_statsvars.inf
     684
    516685        # adding CF lon,lat,time in WRF files
    517686        if test ${headerf:0:3} = 'wrf'; then
    … …  
    530699        ferrmsg ${pyn} ${fname} "python!'file_oper_alongdims'!'tmean'!failed#"${Spyout}
    531700        mv file_oper_alongdims_mean.nc ${filen}
     701
     702        # Keeping the operations
     703        pyins=${pyHOME}"/nc_var.py -o file_oper_alongdims -S '"${vals}"' -f "${cfiles}
     704        pyins=${pyins}" -v "${CFvarn}
     705        echo " " >> ${odir}/all_statsvars.inf
     706        echo "# ${CFvarn}" "${vark}" >> ${odir}/all_statsvars.inf
     707        echo ${pyins} >> ${odir}/all_statsvars.inf
     708
    532709        varkeep=':'${CFvarn}'mean:timestats'
    533710      ;;
    … …  
    543720        ferrmsg ${pyn} ${fname} "python!'file_oper_alongdims'!'tmean'!failed#"${Spyout}
    544721        mv file_oper_alongdims_mean.nc ${filen}
     722
     723        # Keeping the operations
     724        pyins=${pyHOME}"/nc_var.py -o file_oper_alongdims -S '"${vals}"' -f "${cfiles}
     725        pyins=${pyins}" -v "${CFvarn}
     726        echo " " >> ${odir}/all_statsvars.inf
     727        echo "# ${CFvarn}" "${vark}" >> ${odir}/all_statsvars.inf
     728        echo ${pyins} >> ${odir}/all_statsvars.inf
     729
    545730        varkeep=':'${CFvarn}'mean:lonstats'
    546731      ;;
    … …  
    556741        ferrmsg ${pyn} ${fname} "python!'file_oper_alongdims'!'tmean'!failed#"${Spyout}
    557742        mv file_oper_alongdims_mean.nc ${filen}
     743
     744        # Keeping the operations
     745        pyins=${pyHOME}"/nc_var.py -o file_oper_alongdims -S '"${vals}"' -f "${cfiles}
     746        pyins=${pyins}" -v "${CFvarn}
     747        echo " " >> ${odir}/all_statsvars.inf
     748        echo "# ${CFvarn}" "${vark}" >> ${odir}/all_statsvars.inf
     749        echo ${pyins} >> ${odir}/all_statsvars.inf
     750
    558751        varkeep=':'${CFvarn}'mean:latstats'
    559752      ;;
    … …  
    569762      ;;
    570763    esac
     764    cleaning_varsfile ${filen} ${CFvarn}':lon:lat:pres:time'${varkeep}
    571765  fi
    572766
    573   cleaning_varsfile ${filen} ${CFvarn}':lon:lat:pres:time'${varkeep}
     767}
     768
     769function draw_plot(){
     770# Function to carry out the drawing of the plots
     771#   plot= [kplot]+[varn1]|[kind1]#[varn2]|[kind2]#....
     772#     diffmap2Dsfc: 2D map of surface differences values of 1 variable
     773#     diffmap2Dz: 2D map of 3D differences values of 1 variable
     774#     hovmsfc: Hovmoeller diagrams of 1 variable at the surface
     775#     map2Dsfc: 2D map of surface values of 1 variable
     776#     map3D: 2D map of 3D values of 1 variable
     777#     shadcount2D: 2D map of shadow (1st variable) and countour (2nd variable) [stvar1]#[stvar2]
     778#
     779#     ':' separated list of statitsics variable values are given as: [var]|[kind(including combo values)]
     780#     in figures with more than 1 variable, use '#' to separate them
     781#   odir= working directory
     782#   headf= head of the file
     783  fname='draw_plot'
     784
     785  plot=$1
     786  odir=$2
     787  headf=$3
     788  scratch=$4
     789
     790  cd ${odir}
     791
     792  kplot=`echo ${plot} | tr '+' ' ' | awk '{print $1}'`
     793  plotvals=`echo ${plot} | tr '+' ' ' | awk '{print $2}'`
     794
     795  # Assuming that multiple variables will be separated by '#'
     796  nvars=`echo ${plotvals} | tr '#' ' ' | wc -w | awk '{print $1}'`
     797
     798  figfiles=''
     799  varinfilen=''
     800  CFvars=''
     801  figvarparam=''
     802  varkinds=''
     803  if test ${nvars} -ge 1; then
     804    figname=${kplot}_${headf}
     805    iv=1
     806    while test $iv -le ${nvars}; do
     807      varkind=`echo ${plotvals} | tr '#' '\n' | head -n ${iv} | tail -n 1`
     808      varn=`echo ${varkind} | tr '|' ' ' | awk '{print $1}'`
     809      kind=`echo ${varkind} | tr '|' ' ' | awk '{print $2}'`
     810      # Do we have 'pinterp' in the plot?
     811
     812      if $(isInlist ${varmeanname} ${kind}); then
     813        vn=${varn}'mean'
     814      else
     815        vn=${varn}
     816      fi
     817      filen=`stats_filename ${varn} ${kind} ${headf}`
     818      if test ${iv} -eq 1; then
     819        varkinds=${varkind}
     820        CFvars=${varn}
     821        varinfilen=${vn}
     822        figfiles=${filen}
     823      else
     824        varkinds=${varkinds}':'${varkind}
     825        CFvars=${CFvars}':'${varn}
     826        varinfilen=${varinfilen}':'${vn}
     827        figfiles=${figfiles}':'${filen}
     828      fi
     829
     830      figname=${figname}'_'${vn}'-'${kind}
     831
     832      # Getting plot variable configuration
     833      specific=false
     834      if $(isInlist ${specificvarplot} ${varn}); then
     835        vals=`listfilter ${specificvarplot} ${varn}`
     836        kvals=`listfilter ${specificvarplot} ${kind}`
     837        Lkvals=`expr length ${kvals}'0'`
     838        if test ${Lkvals} -gt 1; then
     839          fkvals=`listfilter ${kvals} ${kplot}`
     840          Lfkvals=`expr length ${fkvals}'0'`
     841          if test ${Lfkvals} -gt 1; then
     842            specific=true
     843            minval=`echo ${fkvals} | tr '|' ' ' | awk '{print $4}'`
     844            maxval=`echo ${fkvals} | tr '|' ' ' | awk '{print $5}'`
     845            colorbar=`echo ${fkvals} | tr '|' ' ' | awk '{print $6}'`
     846            cntformat=`echo ${fkvals} | tr '|' ' ' | awk '{print $7}'`
     847            colorcnt=`echo ${fkvals} | tr '|' ' ' | awk '{print $8}'`
     848          fi
     849        fi
     850      fi
     851      if ! ${specific}; then
     852        allvals=`python ${pyHOME}/drawing.py -o variable_values -S ${varn} |         \
     853          grep all_values | awk '{print $3}'`
     854        pyn=$?
     855        ferrmsg ${pyn} ${fname} "python!'variable_values'!'${varn}'!failed#"
     856        minval=`echo ${allvals} | tr ',' ' ' | awk '{print $3}'`
     857        maxval=`echo ${allvals} | tr ',' ' ' | awk '{print $4}'`
     858        colorbar=`echo ${allvals} | tr ',' ' ' | awk '{print $7}'`
     859        cntformat='%g'
     860        colorcnt='black'
     861      fi
     862
     863      # Graphical paramters for each variable in the plot
     864      if test ${iv} -eq 1; then
     865        figvarparam=${minval}'|'${maxval}'|'${colorbar}'|'${cntformat}'|'${colorcnt}
     866      else
     867        figvarparam=${figvarparam}':'${minval}'|'${maxval}'|'${colorbar}'|'
     868        figvarparam=${figvarparam}${cntformat}'|'${colorcnt}
     869      fi
     870
     871      iv=`expr ${iv} + 1`
     872    # End of number of variables of the plot
     873    done
     874  fi
     875  figname=${figname}.${kindfig}
     876
     877  if ${scratch}; then rm ${figname} >& /dev/null; fi
     878  if test ! -f ${figname}; then
     879    case ${kplot} in
     880      'diffmap2Dsfc')
     881        echo "  "${fname}": kind of plot '"${kplot}"' not ready !!"
     882        exit
     883      ;;
     884      'diffmap2Dz')
     885        echo "  "${fname}": kind of plot '"${kplot}"' not ready !!"
     886        exit
     887      ;;
     888      'hovmsfc')
     889        echo "  "${fname}": kind of plot '"${kplot}"' not ready !!"
     890        exit
     891      ;;
     892      'map2Dsfc')
     893        echo "  "${fname}": kind of plot '"${kplot}"' not ready !!"
     894        exit
     895      ;;
     896      'map3D')
     897        echo "  "${fname}": kind of plot '"${kplot}"' not ready !!"
     898        exit
     899      ;;
     900      'shadcount2Dsfc')
     901        figtit=`echo ${varkinds} | tr ':' '|'`
     902        shdstdn=`echo ${CFvars} | tr ':' ' ' | awk '{print $1}'`
     903        cntstdn=`echo ${CFvars} | tr ':' ' ' | awk '{print $2}'`
     904        figfs=`echo ${figfiles} | tr ':' ','`
     905        srange=`echo ${figvarparam} | tr ':' ' ' | awk '{print $1}' | tr '|' ' ' |   \
     906          awk '{print $1","$2}'`
     907        cbar=`echo ${figvarparam} | tr ':' ' ' | awk '{print $1}' | tr '|' ' ' |     \
     908          awk '{print $3}'`
     909        crange=`echo ${figvarparam} | tr ':' ' ' | awk '{print $2}' | tr '|' ' ' |   \
     910          awk '{print $1","$2}'`
     911        cline=`echo ${figvarparam} | tr ':' ' ' | awk '{print $2}' | tr '|' ' ' |    \
     912          awk '{print $5}'`
     913        cfmt=`echo ${figvarparam} | tr ':' ' ' | awk '{print $2}' | tr '|' ' ' |     \
     914          awk '{print $4}'`
     915
     916        graphvals=$(echo ${CFvars} | tr ':' ',')
     917        graphvals=${graphvals}':lon|-1,lat|-1:lon|-1,lat|-1:lon:lat:'${cbar}':fixc,'
     918        graphvals=${graphvals}${cline}':'${cfmt}':'${srange}':'${crange}',9:'
     919        graphvals=${graphvals}${figtit}':'${kindfig}':False:'${mapval}
     920 
     921        plotins="python "${pyHOME}"/drawing.py -f "${figfs}" -o draw_2D_shad_cont "
     922        plotins=${plotins}"-S '"${graphvals}"' -v "${shdstdn}"mean,"${cntstdn}"mean"
     923        echo " " >> ${odir}/all_figures.inf
     924        echo "#"$(echo $f{igtit} | tr '|' ' ') >> ${odir}/all_figures.inf
     925        echo ${plotins} >> ${odir}/all_figures.inf
     926        pyout=`python ${pyHOME}/drawing.py -f ${figfs} -o draw_2D_shad_cont          \
     927          -S ${graphvals} -v ${shdstdn}'mean,'${cntstdn}'mean'`
     928        pyn=$?
     929        Spyout=`echo ${pyout} | tr '\n' '#' | tr ' ' '!'`
     930        ferrmsg ${pyn} ${fname} "python!'draw_2D_shad_cont'!failed#"${Spyout}
     931        mv 2Dfields_shadow-contour.${kindfig} ${figname}
     932#        exit
     933      ;;
     934      *)
     935        echo ${errormsg}
     936        echo "  "${fname}": plot kind '"${kplot}"' not ready !!"
     937        exit
     938      ;;
     939    esac
     940  fi
     941
     942#LLUIS
     943
    574944}
    575945
    … …  
    587957  echo "  "${main}": starting from the SCRATCH !!"
    588958  echo "    10 seconds left!!"
     959  filescratch=true
     960  figscratch=true
    589961  sleep 10
    590962else
    591963  scratch=false
     964  if test ${filescratch} = 'true'; then
     965    filescratch=true
     966    echo ${warnmsg}
     967    echo "  "${main}": files starting from the SCRATCH !!"
     968    echo "    5 seconds left!!"
     969    sleep 5
     970  else
     971    filescratch=false
     972  fi
     973  if test ${figscratch} = 'true'; then
     974    figscratch=true
     975    echo ${warnmsg}
     976    echo "  "${main}": figures starting from the SCRATCH !!"
     977    echo "    5 seconds left!!"
     978    sleep 5
     979  else
     980    figscratch=false
     981  fi
    592982fi
    593983
    … …  
    6731063
    6741064    # Need to pass to analyze all the data?
    675     if ${scratch}; then rm ${owdir}/varcompute.inf >& /dev/null; fi
     1065    if ${filescratch}; then
     1066      rm ${owdir}/varcompute.inf >& /dev/null
     1067      rm ${owdir}/all_computevars.inf >& /dev/null
     1068      rm ${owdir}/all_statsvars.inf >& /dev/null
     1069
     1070      echo "## Computation of variables " > ${owdir}/all_computevars.inf
     1071      echo "## Computation of statistics " > ${owdir}/all_statsvars.inf
     1072
     1073    fi
     1074
     1075
    6761076    if test ! -f ${owdir}/varcompute.inf; then
    6771077
    … …  
    6821082      ih=1
    6831083      for fh in ${fheaders}; do
    684         if ${scratch}; then rm ${owdir}/*_${fh}*.nc >& /dev/null; fi
     1084        if ${filescratch}; then rm ${owdir}/*_${fh}*.nc >& /dev/null; fi
    6851085
    6861086        files1h=`ls -1 ${fh}* | tr '\n' '@'`
    … …  
    8081208      done
    8091209
    810       # Outwritting the varcompute to avoid next time (if it is not scratch!)
     1210      # Outwritting the varcompute to avoid next time (if it is not filescratch!)
    8111211      cat << EOF > ${owdir}/varcompute.inf
    8121212files: ${files}
    … …  
    8431243      if test ! ${modv} = 'None' || test ! ${diagv} = 'None'; then
    8441244        compute_variable ${iwdir} ${files} ${owdir} ${cvar} ${moddims} ${modvdims}   \
    845           ${scratch}
     1245          ${filescratch}
    8461246        ic=`expr ${ic} + 1`
    8471247
    … …  
    8531253          fi
    8541254          compute_statistics ${iwdir} ${CFv}_${fhead}.nc ${owdir} ${cvar}            \
    855             ${moddims} ${modvdims} ${scratch}
     1255            ${moddims} ${modvdims} ${filescratch}
    8561256          isc=`expr ${isc} + 1`
    8571257        else
    … …  
    8921292          combcvar=${CFv}'|'${comb}'|'${fileh}'|'${modv}'|'${diagv}
    8931293          compute_statistics ${iwdir} ${ifile}.nc ${owdir} ${combcvar} ${moddims}    \
    894             ${modvdims} ${scratch} ${ifile}.nc
     1294            ${modvdims} ${filescratch} ${ifile}.nc
    8951295          ifile=${ifile}_${comb}
    8961296          isc=`expr ${isc} + 1`
    … …  
    9041304    echo "  "${main}": "${isc}" statistics has been computed"
    9051305
    906     exit
    9071306# Direct plotting
    9081307##
    909     dirplots=`echo ${directplots} | tr ':' ' '`
     1308    echo "  "${main}": Plotting direct figures ..."
     1309    if ${figscratch}; then rm directplotsdraw.inf; fi
    9101310    if test ! -f directplotsdraw.inf; then
     1311      drwplts=`echo ${drawplots} | tr ':' ' '`
     1312      ikp=1
    9111313      plots=''
    912       ikp=1
    913       for kplot in ${dirplots}; do
    914         case ${kplot} in
    915           'hovmsfc')
    916             kplots=${plthovmsfc}
    917           ;;
    918           'map2Dsfc')
    919             kplots=${pltmap2Dsfc}
    920           ;;
    921           'map3D')
    922             kplots=${pltmap3D}
    923           ;;
    924           'shadcount2Dsfc')
    925             kplots=${pltshadcount2Dsfc}
    926           ;;
    927           *)
    928             echo ${errmsg}
    929             echo "  "${main}": plot kind '"${kplot}"' not ready !!"
    930             exit
    931         esac
    932         if test ${ikp} -eq 0; then
    933           plots=${kplots}
    934         else
    935           plots=${plots}':'${kplots}
    936         fi
    937         ikp=`expr ${ikp} + 1`
     1314      for drw in ${drwplts}; do
     1315        if $(isInlist ${directplots} ${drw}); then
     1316          case ${drw} in
     1317            'hovmsfc')
     1318              plts=${plothovmsfc}
     1319            ;;
     1320            'map2Dsfc')
     1321              plts=${pltmap2Dsfc}
     1322            ;;
     1323            'map3D')
     1324              plts=${pltmap3D}
     1325            ;;
     1326            'shadcount2Dsfc')
     1327              plts=${pltshadcount2Dsfc}
     1328            ;;
     1329            *)
     1330              echo ${errmsg}
     1331              echo "  "${main}": plot kind '"${drw}"' not ready !!"
     1332              exit
     1333          esac
     1334          # Do we have plots for this kind?
     1335          Lkplots=`expr length ${plts}'0'`
     1336          if test ${Lkplots} -lt 2; then
     1337            ferrmsg 1 ${main} "plot!kind!"${drw}"!without!variables"
     1338          fi
     1339          kplots=`list_add_label ${plts} ${drw} + beg`
     1340          if test ${ikp} -eq 1; then
     1341            plots=${kplots}
     1342          else
     1343            plots=${plots}':'${kplots}
     1344          fi
     1345          ikp=`expr ${ikp} + 1`
     1346        fi
     1347# End of coincident direct plots
    9381348      done
     1349      cat << EOF > directplotsdraw.inf
     1350plots= ${plots}
     1351EOF
    9391352    else
    9401353      echo ${warnmsg}
    9411354      echo "  "${main}": retrieving direct plots from existing file 'directplotsdraw.inf' !!"
    942     fi
     1355      plots=`cat directplotsdraw.inf | grep plots | awk '{print $2}'`
     1356    fi
     1357
     1358    if ${figscratch}; then rm ${owdir}/all_figures.inf; fi
     1359    if test ! -f ${owdir}/all_figures.inf; then
     1360      echo "###### List of all figures" > ${owdir}/all_figures.inf
     1361    fi
     1362
     1363    pts=`echo ${plots} | tr ':' ' '`
     1364    for pt in ${pts}; do
     1365      echo "        "${pt}
     1366      draw_plot ${pt} ${owdir} ${fileh} ${figscratch}
     1367    done
    9431368
    9441369# LLUIS

  • trunk/tools/model_graphics_template.dat

    r924 r926  
    2020# Srcatch
    2121scratch=false
     22
     23# Scratch files
     24filescratch=false
     25
     26# Scratch figures
     27figscratch=false
    2228
    2329# Debug
    … …  
    6975#   must exist on the table $pyHOME/variables_values.dat
    7076#   varcombo = [varn];[vark1]@...@[varkN]: variable as consecutive kinds of operation
    71 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
    72 vartmean = prc:prls:zmla:hfss:hfls:evspsbl:bils:uas:vas:tas:ps:clt:cll:clm:clh:prw:huss:pr:wakes:wakeh:stherm:zmaxth:wss
     77varlast = 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
     78vartmean = 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
    7379varxmean = ta:ua:va:hur
    7480vardiff = ualmean:valmean:talmean:huslmean:uas:vas:ps:pr
    7581varpinterp = ta:ua:va:hur
    7682varcombo = ta;pinterp@last@xmean:ua;pinterp@last@xmean:va;pinterp@last@xmean:hur;pinterp@last@xmean
     83
     84# Operations which attach 'mean' at the end of the variable name
     85varmeanname='Lmean:lmean:tmaen:xmean:ymean'
    7786
    7887# Pressure levels
    … …  
    8998#
    9099#   ':' separated list of statitsics variable values are given as: [var]|[kind(including combo values)]
     100#     in figures with more than 1 variable, use '#' to separate them
     101drawplots = shadcount2Dsfc
    91102
    92103# Figures which can be made directly without another model/experiment output
    93 directplots='hovmsfc:map2Dsfc:map3D:shadcount2Dsfc'
     104directplots = hovmsfc:map2Dsfc:map3D:shadcount2Dsfc
    94105
    95 pltshadcount2Dsfc = huss|tmean#tas|tmean:vas|tmean#uas|tmean:ps|tmean#wss|mean: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|mean: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
     106pltshadcount2Dsfc = 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|mean: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
    96107
    97108pltdiffmap2Dsfc = uas:vas:ps:pr
    98109pltdiffmap2Dz = uapxmean:vapxmean:tapxmean:huspxmean
    99110
    100 # Figures output
     111# Specific variable-statistisc graphic parameters (if no value is given, they will be retrieved from 'variables_values.dat')
     112#  ':' separated list of specific parameters for plotting as:
     113#      [varname]|[statistics]|[minval]|[maxval]|[colorbar]|[cntformat]|[colorcnt]
     114#        [varname]: name of the variable
     115#        [statistics]: statistics (also combo)
     116#        [figure]: figure to be used
     117#        [minval]: minimum value
     118#        [maxval]: minimum value
     119#        [colorbar]: name of the colorbar (from matplotlib) to use
     120#        [cntformat]: format of the contour labels
     121#        [colorcnt]: color for the countor lines
     122specificvarplot = cll|last|shadcount2Dsfc|0.|1.|gist_gray|black:clh|last|shadcount2Dsfc|0.|1.|gist_gray|black
     123
     124# Figures output kind
    101125kindfig = pdf
     126
     127# Map to use
     128mapval = lcc,l
    102129
    103130# Generic
Note: See TracChangeset for help on using the changeset viewer.