Changeset 243 in lmdz_wrf


Ignore:
Timestamp:
Feb 2, 2015, 6:44:09 PM (10 years ago)
Author:
lfita
Message:

Adding pressure table eqjuivalencies with the level number

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/tools/lmdz_check.bash

    r242 r243  
    1 #!/bin/bash
     1#!/bin/bash -x
    22## e.g. # lmdz_check.bash -1:6:3:26 wrfout_d01_1979-12-01_00:00:00 /home/lluis/PY /home/lluis/etudes/WRF_LMDZ/WaquaL/LaTeX/NPv31_hgardfou/figs 0 None 0
    33## e.g. # lmdz_check.bash -1:9:9:4 iphysiq1_time_step150/wrfout_d01_check.nc /home/lfita/PY iphysiq1_time_step150/check 1 None 0
    … …  
    1818  echo "  [visualize]: 0, no; 1: yes, whether to visualize results"
    1919else
     20
     21function vals_textab() {
     22# Function to provide an LaTeX like table from a series of ':' separated values
     23#   values= list of values
     24#   ncols= number of columns
     25#   labels= which labels
     26#     'inum': integer ascending numbers
     27#     'labs@'[labels]: as a list of given ':' separated labels
     28#   caption= text as caption '*' for spaces
     29#   tablabel= text as label of the table
     30
     31# Output in `vals_textab.tex' as:
     32#   A LaTeX table as \verb+[ilabel]:+ & [ivalue] & ... & [Ncol]
     33  errmsg='ERROR -- error -- ERROR -- error'
     34  fname='vals_textab'
     35
     36  values=$1
     37  ncols=$2
     38  labels=$3
     39  caption=$4
     40  tablabel=$5
     41
     42  vals=`echo ${values} | tr ':' ' '`
     43  cap=`echo ${caption} | tr '*' ' '`
     44
     45  Nvals=`echo ${vals} | wc -w | awk '{print $1}'`
     46
     47  ofile='vals_textab.tex'
     48
     49  icol=1
     50  tabfmt=''
     51  while test ${icol} -le ${ncols}; do
     52    tabfmt=${tabfmt}'lc'
     53    icol=`expr ${icol} + 1`
     54  done
     55
     56  labk=${labels:0:4}
     57  case ${labk} in
     58    'inum')
     59      in=1
     60      while test ${in} -le ${Nvals}; do
     61        labvals=${labvals}${in}':'
     62        in=`expr ${in} + 1`
     63      done
     64    ;;
     65    'labs')
     66      labvals=`echo ${labels} | tr '@' ' ' | awk '{print $2}'`
     67    ;;
     68    *)
     69      echo ${errmsg}
     70      echo "  "${fname}": kind of labels '"${labk}"' not ready!!"
     71      exit
     72    ;;
     73  esac
     74
     75  cat << EOF > ${ofile}
     76\\begin{table}
     77\\caption{${cap}}
     78\\label{tab:${tablabel}}
     79\\begin{center}
     80\\begin{tabular}{${tabfmt}}
     81EOF
     82  Nmultivalscols=`expr ${Nvals} % ${ncols}`
     83  if test ${Nmultivalscols} -ne 0; then
     84    Nvals0=`expr ${Nvals} - ${Nmultivalscols}`
     85  else
     86    Nvals0=${Nvals}
     87  fi
     88
     89  ival=1
     90  line=''
     91
     92  while test ${ival} -le ${Nvals}; do
     93    modival=`expr ${ival} % ${ncols}`
     94    vali=`echo ${values} | tr ':' '\n' | head -n ${ival} | tail -n 1`
     95    labi=`echo ${labvals} | tr ':' '\n' | head -n ${ival} | tail -n 1`
     96
     97    if test ${modival} -eq 0; then
     98      line=${line}'\verb+'${labi}':+ & '${vali}' \\ '
     99      echo ${line} >> ${ofile}
     100      line=''
     101    else
     102      line=${line}'\verb+'${labi}':+ & '${vali}' & '
     103    fi
     104
     105    ival=`expr ${ival} + 1`
     106  done
     107
     108  if test ${Nmultivalscols} -ne 0; then
     109    Nemptycols=`expr ${ncols} - ${Nmultivalscols}`
     110    if test ${Nemptycols} -ne 0; then
     111      iecol=1
     112      while test ${iecol} -lt ${Nemptycols}; do
     113        line=${line}' & & '
     114        iecol=`expr ${iecol} + 1`
     115      done
     116    fi
     117    line=${line}' & \\ '
     118    echo ${line} >> ${ofile}
     119  fi
     120
     121  cat << EOF >> ${ofile}
     122\\end{tabular}
     123\\end{center}
     124\\end{table}
     125EOF
     126
     127  return
     128}
    20129
    21130function to_LaTeX() {
    … …  
    451560  if test $dolatex -eq 1; then
    452561  echo "  LaTeX file '"${otex}"' generation..."
     562
     563# Vertical pressure values
     564  dvals='Time:-9|bottom_top:-1|south_north:'${yval}'|west_east:'${xval}
     565  pvals=`python ${PYHOME}/nc_var.py -o varout -f ${infile} -S ${dvals} -v LPRES | awk '{print $2}'`
     566  presvals=`echo ${pvals} | tr ' ' ':'`
     567  vals_textab ${presvals} 3 inum 'level*pressure*(Pa)*equivalencies' pres_vals
     568
    453569  if test ${tval} -eq -1; then tvalS='last'; fi
    454570
    … …  
    456572\\documentclass{article}
    457573\\usepackage{graphicx}
    458 \usepackage[colorlinks=true,urlcolor=blue]{hyperref}
     574\\usepackage[colorlinks=true,urlcolor=blue]{hyperref}
    459575
    460576\\textwidth=18cm
    … …  
    473589\\listoffigures
    474590\\newpage
    475 EOF
     591
     592EOF
     593  cat vals_textab.tex >> ${otex}
     594
     595  cat << EOF >> ${otex}
     596
     597\\clearpage
     598
     599EOF
     600
    476601  for vt in ${vtypes}; do
    477602    echo "variable type: "${vt}
    … …  
    688813    esac
    689814  done
     815
    690816  cat << EOF >> ${otex}
    691817\\end{document}
Note: See TracChangeset for help on using the changeset viewer.