Interpolate multiple maps and get some raster statistics

#!/bin/sh
#first array
monthNamesAbr=(jan feb mar apr mai jun jul aug sep okt nov dec)
#second array
modelNames=(m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 m13 m14 m15 m16 m17 m18 m19)

#helper string
periods1="fut"

for i in "${monthNamesAbr[@]}"  #outer loop over months
do
    for j in "${modelNames[@]}" #inner loop over models
    do
        v.surf.rst input=${i}_${periods1}_norm elev=${i}_${periods1}_norm_${j} --overwrite zcolumn=${j}
        #line number 9 containts maximum of raster value.
        r.univar map=${i}_${periods1}_norm_${j} |awk -v month=$i -v model=$j 'NR==9{
        print month",",model,",",$0
        }'>>maximum.txt
    done    #end of inner loop
done        #end of outer loop