Printing an Array Slice

The following macro prints a Fortran array slice. This macro, like others shown in this chapter, relies heavily on Tcl and uses unique CLI commands sparingly.

proc pf2Dslice {anArray i1 i2 j1 j2 {i3 1} {j3 1} \

{width 20}} {

for {set i $i1} {$i <= $i2} {incr i $i3} {

set row_out ""

for {set j $j1} {$j <= $j2} {incr j $j3} {

set ij [capture dprint $anArray\($i,$j\)]

set ij [string range $ij \

[expr [string first "=" $ij] + 1] end]

set ij [string trimright $ij]

if {[string first "-" $ij] == 1} {

set ij [string range $ij 1 end]}

append ij " "

append row_out " " \

[string range $ij 0 $width] " "

}

puts $row_out

}

}

The CLI’s dprint command lets you specify a slice. For example, you can type: dprint a(1:4,1:4).

After invoking this macro, the CLI prints a two-dimensional slice (i1:i2:i3, j1:j2:j3) of a Fortran array to a numeric field whose width is specified by the width argument. This width doesn’t include a leading minus sign (-).

All but one line is standard Tcl. This line uses the dprint command to obtain the value of one array element. This element’s value is then captured into a variable. The CLI capture command allows a value that is normally printed to be sent to a variable. For information on the difference between values being displayed and values being returned, see “About CLI Output”.

The following shows how this macro is used:

d1.<> pf2Dslice a 1 4 1 4

0.841470956802 0.909297406673 0.141120001673-0.756802499294

0.909297406673-0.756802499294-0.279415488243 0.989358246326

0.141120001673-0.279415488243 0.412118494510-0.536572933197

-0.756802499294 0.989358246326-0.536572933197-0.287903308868

d1.<> pf2Dslice a 1 4 1 4 1 1 17

0.841470956802 0.909297406673 0.141120001673-0.756802499294

0.909297406673-0.756802499294-0.279415488243 0.989358246326

0.141120001673-0.279415488243 0.412118494510-0.536572933197

-0.756802499294 0.989358246326-0.536572933197-0.287903308868

d1.<> pf2Dslice a 1 4 1 4 2 2 10

0.84147095 0.14112000

0.14112000 0.41211849

d1.<> pf2Dslice a 2 4 2 4 2 2 10

-0.75680249 0.98935824

0.98935824 -0.28790330

d1.<>