Shell Programming and Scripting

awk '{
  split(sprintf(fmt,$1),N,/E/)            # Convert input number to scientific format and split it into array N, using "E" as field separator
  ofmt=fmt                                # Set output format to input format 
  sub(/E/,"fE",ofmt)                      # Replace "E" by "fE" in the output format
  printf ofmt "%+03d\n",N[1]/10,N[2]+1    # Print the part before the E /10 and after E+1
}' fmt='%.11E'                            # Set f to the precision

awk '{
  i = index($0, ".");
  split(sprintf(fmt,$1),N,/E/)
  ofmt=fmt
  sub(/E/,"fE",ofmt)
  printf ofmt "%+03d\n",N[1]/10,N[2]+(i-1)
}' fmt='%.11E'

Input_value     Scrutinizer_output      Don_Cragun_output
123456789555    0.12345678955E+12       0.12345678956E+12
0               0.00000000000E+01       0.00000000000E+01
1e400           infE+01                 0.IE-9223372036854775808
1e-312          1.00000000000E-312      0.10000000000E-311

awk '{ does some search, replace, algebra and etc.
q1 = some algebraic calculation
q2 = some algebraic calculation
q3 = some algebraic calculation
q4 = some algebraic calculation

printf ("%s  %.11E  %.11E  %.11E  %.11E\n",block,q1,q2,q3,q4)
}' data

awk '
  function pre( fmt, nr,	ofmt, N) {
    split(sprintf(fmt, nr),N,/E/)                  # Convert input number to scientific format and split it into array N, using "E" as field separator
    ofmt=fmt                                       # Set output format to input format 
    sub(/E/,"fE",ofmt)                             # Replace "E" by "fE" in the output format
    return sprintf(ofmt "%+03d",N[1]/10,N[2]+1)    # return the part before the E /10 and after E+1
  }
  { 
    q1=3.01234567891E+03
    q2=123.45678901
    q3=12345e200
    q4=1234e-200 
    printf "%s %s %s %s\n", pre( "%.11E", q1), pre( "%.11E", q2),pre( "%.11E", q3),pre( "%.11E", q4) 
  }
'

#!/usr/bin/env bash

# @(#) s1	Demonstrate Fortran formatting.

# Utility functions: print-as-echo, print-line-with-visual-space, debug.
# export PATH="/usr/local/bin:/usr/bin:/bin"
pe() { for _i;do printf "%s" "$_i";done; printf "\n"; }
pl() { pe;pe "-----" ;pe "$*"; }
db() { ( printf " db, ";for _i;do printf "%s" "$_i";done;printf "\n" ) >&2 ; }
db() { : ; }
C=$HOME/bin/context && [ -f $C ] && $C gfortran

FILE=${1-data1}

s=f1.f90
pl " Fortran source, $s:"
cat $s

pl " Input data file $FILE:"
cat $FILE

pl " Results:"
rm -f a.out
gfortran $s
./a.out < data1

exit 0

% ./s1

Environment: LC_ALL = C, LANG = C
(Versions displayed with local utility "version")
OS, ker|rel, machine: Linux, 2.6.26-2-amd64, x86_64
Distribution        : Debian GNU/Linux 5.0.8 (lenny) 
bash GNU bash 3.2.39
gfortran GNU Fortran (Debian 4.3.2-1.1) 4.3.2

-----
 Fortran source, f1.f90:
     program f1
! @(#) f1	Demonstrate Fortran-90 formatting.
       do
         read(*,*,end=100) x
         write(*,fmt="(e18.11)") x
       end do
100    continue
     end program f1

-----
 Input data file data1:
3.01234567891E+03
123.45678901
12345e200
1234e-200
-123456
1e400
1e-312
0

-----
 Results:
 0.30123457031E+04
 0.12345678711E+03
         +Infinity
 0.00000000000E+00
-0.12345600000E+06
         +Infinity
 0.00000000000E+00
 0.00000000000E+00