my $theta = $T * (( 1013 / $_ ) ** ( 287/1004 ));
    push ( @press, $_ );
    push ( @thetas, $theta );