#!/bin/bash
#
# -- Make the variable mesh definitions for the AM files 
# -- 1 parameter - the file name of the text file (mesh.txt)
#
 PROG=`basename $0`
# echo Start $PROG
 if [ $# -lt 1 ]; then
   echo Error - no filename specified
   exit 1
 fi

 awk '
   $1~/X/ {for(i=2;i<=NF;i++) xl[i-1]=$i; nx=NF-1; xmn=$2; xmx=$NF}
   $1~/K/ {for(i=2;i<=NF;i++) xd[i-1]=$i; nk=NF-1}
   $1~/Y/ {for(i=2;i<=NF;i++) yl[i-1]=$i; ny=NF-1; ymn=$2; ymx=$NF}
   $1~/L/ {for(i=2;i<=NF;i++) yd[i-1]=$i; nl=NF-1}
   END{
       k=1;
       for(i=1;i<=nx;i++){
         if(i<nx){
           kk=int((xl[i+1]-xl[i])/xd[i]+0.1);
           k+=kk;
           nds[i]=k;
         }
         if(i>1 && i<nx){
           printf"xreg%d =%7.1f, ",i-1,xl[i];
         }
       }
       printf"\n";
       for(i=1;i<nx-1;i++) printf"kreg%d =%5d  , ",i,nds[i];
       printf"kmax =%5d  , \n",nds[nx-1];
       printf";x step"                  ;
       for(i=1;i<nx;i++) printf" %7.3f        ",xd[i];
       printf"\n";

       k=1;
       for(i=1;i<=ny;i++){
         if(i<ny){
           kk=int((yl[i+1]-yl[i])/yd[i]+0.1);
           k+=kk;
           nds[i]=k;
         }
         if(i>1 && i<ny){
           printf"yreg%d =%7.1f, ",i-1,yl[i];
         }
       }
       printf"\n";
       for(i=1;i<ny-1;i++) printf"lreg%d =%5d  , ",i,nds[i];
       printf"lmax =%5d  , \n",nds[ny-1];
       printf";y step"                  ;
       for(i=1;i<ny;i++) printf" %7.3f        ",yd[i];
       printf"\n";
       print "; the full area";
       printf"&po y= %8.2f , x= %8.2f &\n",ymn,xmn;
       printf"&po y= %8.2f , x= %8.2f &\n",ymx,xmn;
       printf"&po y= %8.2f , x= %8.2f &\n",ymx,xmx;
       printf"&po y= %8.2f , x= %8.2f &\n",ymn,xmx;
       printf"&po y= %8.2f , x= %8.2f &\n",ymn,xmn;

      }
   ' $1