compute_opersvarsfiles(values, varinfo)
  [prevalues]: values from the previous operations
  NOTE: to start with a multiple operations one can use as operation of the first variable 'addc,0'
  values= '%' separated list of dimension ranges and file,operation,variable to compute
    [dimranges1]@[operfile1var]%[dimranges2]@[operfile2var]%[...[dimrangesM]@[operfileMvar]]
    dimranges: ';' separated list of dimension, dimension-variable names and their ranges ('|' separated)
      [dim1]|[dimv1]|[range];[dim2]|[dimv2]|[range];[...|[dimN]|[dimvN]|[range]]
      with [range] as:
        -1, all the values
        -9, last values
        int, a single value
        [beg,end,frq], from a beginning to an end with a given frequency
        &nsbp;&bnsp;'WRFtime': as time variable computed from WRF 'Times' variable
        NOTE: use 'contoperation' to continue an operation with the result of previous
          operations [prevalues] without using a variable from a file
    operfilevar: [oper]|[file]|[var] '|' separated list of triplets of [operation], [file], [variable name to use]
      * [oper]: operations (being [prevalues] the values computed until the operation)
        'add': adding [var] ([prevalues] + [var])
        'addc',[modval1]: add [modval1] to [prevalues]
        'arctan', arctangent [-p1, pi]
        'centerderiv',[N],[ord],[dim]: un-scaled center [N]-derivative of order [ord] along dimension [dim] of [var]
        'div': dividing by [var] ([prevalues] / [var])
        'divc',[modval1]: [prevalues] divide by [modval1]
        'forwrdderiv',[N],[ord],[dim]: un-scaled forward [N]-derivative of order [ord] along dimension [dim] of [var]
        'ifreq_anom',[stepdimn],[stepvardimn],[istep]: computing anomalies by substracting sub-means at each
            step by averaging from there all values taken every [istep] along dimension [stepdim].
            mean(j) = sum(matA[j+k*istep]_k=0,Nstep)/Nstep; Nstep = len(stepdimn)/istep; j=[0,istep]
          [stepdimn]= name of the dimension along which to sample
          [stepvardimn]= name of the variable-dimension with the values for [stepdimn]
          [istep]= frequency to sample
        'ifreq_mean',[stepdimn],[stepvardimn],[istep]: computing sub-means at each step by averaging from
            there all values taken every [istep] along dimension [stepdim].
            mean(j) = sum(matA[j+k*istep]_k=0,Nstep)/Nstep; Nstep = len(stepdimn)/istep; j=[0,istep]
          [stepdimn]= name of the dimension along which to sample
          [stepvardimn]= name of the variable-dimension with the values for [stepdimn]
          [istep]= frequency to sample
        'ifreq_normmeanstd',[stepdimn],[stepvardimn],[istep]: nbormalizing anomalies by substracting mean(j)/dtsv(j)
            sub-stats at each step by averaging from there all values taken every [istep] along dimension[stepdim].
            mean(j) = sum(matA[j+k*istep]_k=0,Nstep)/Nstep; std(j) = std(matA[j+k*istep]_k=0,Nstep);
            Nstep = len(stepdimn)/istep; j=[0,istep]
          [stepdimn]= name of the dimension along which to sample
          [stepvardimn]= name of the variable-dimension with the values for [stepdimn]
          [istep]= frequency to sample
        'inv': inverting [prevalues] (1/[prevalues])
        'lowthres',[modval1],[modval2]: if [prevalues] < [modval1]; prevalues = [modval2]
        'lowthres@oper',[modval1],[oper],[modval2]: if [prevalues] < [modval1]; prevalues = [oper] (operation as [modval2])
        'mul': multiplying by [var] ([prevvalues] * [var])
        'mulc',[modval1]: [prevalues] multiplied by [modval1]
        'norm_meanminmax',[NOTnormdims]: normalization of data as: (val-)/(max(val)-min(val)) except along
          dimensions [NOTnormdims] (':' list of dimension names, or 'any' for using all dimensions)
        'norm_minmax',[NOTnormdims]: normalization of data as: [2val-(max(val)-min(val)]/[max(val)+min(val)]
          except along dimensions [NOTnormdims] (':' list of dimension names, or 'any' for using all dimensions)
        'norm_meanstd',[NOTnormdims]: normalization of data as: (val-)/stdev(val) except along
          dimensions [NOTnormdims] (':' list of dimension names, or 'any' for using all dimensions)
          dimensions [NOTnormdims] (':' list of dimension names)
        'pot': powering with [var] ([prevalues] ** [var])
        'potc',[modval1]: [prevalues] ** [modval1]
        'repl',[modval1]: replace values of [prevalues] with values from [operfile2var] except
          values in second file = [modval1]
        'replbig',[modval1]: replace values of [prevalues] with values from [operfile2var] except
          values in second file > [modval1]
        'repless',[modval1]: replace values of [prevalues] with values from [operfile2var] except
          values in second file < [modval1]
        'sub': substracting [var] ([prevalues] - [var])
        'subc',[modval1]: remove [modval1] of [prevalues]
        'upthres',[modval1],[modval2]: if [prevalues] > [modval1]; prevalues = [modval2]
        'upthres@oper',[modval1],[oper],[modval2]: if [prevalues] > [modval1]; prevalues = [oper] (operation with [modval2])
      * [file]: name of the file
      * [var]: variable to use ('None' for the case with constant operations: 'addc', 'divc', 'inv', 'lowthres',
        'lowthres@oper', 'mulc', 'potc', 'subc', 'upthres', 'upthres@oper'
  varinfo= [varname],[Lvarname],[varunits]
    varname: name of the final variable
    Lvarname: long name of the final variable ('!' for spaces)
    varunits: units of the final variable ('!' for spaces)