请求帮助：代码查错

凋谢1991

请走过路过的大神们帮忙查找一下为什么不能输出Script.txt文件以及内容，代码如下：附件中是调运子程序castepcalculation计算得到的1.castep文件。

1. #!perl
   
2. 
   
3. use strict;
   
4. use Getopt::Long;
   
5. use MaterialsScript qw(:all);
   
6. 
   
7. my $a     = 1;
   
8. my $myDoc = $Documents{"$a.xsd"};
   
9. my $stressz;
   
10. 
    
11. sub castepcalculation ($stressz) {
    
12.     my ($stressz) = @_;
    
13. 
    
14.     my $castep = Modules->CASTEP;
    
15.     $castep->ChangeSettings(
    
16.         Settings(
    
17.             Task                  => "Geometry Optimization",
    
18.             NonLocalFunctional    => "PBE",
    
19.             UseCustomEnergyCutoff => "Yes",
    
20.             EnergyCutoff          => 350,
    
21.             SCFConvergence        => 5e-002,
    
22.             MaximumSCFCycles      => 200,
    
23.             CalculateStress       => "Yes",
    
24.             OptimizeCell          => "Yes",
    
25.             EnergyConvergence     => 5e-001,
    
26.             Sxx                   => 0,
    
27.             Syy                   => 0,
    
28.             Szz                   => "$stressz",
    
29.             Quality               => "Ultra-fine",
    
30. 
    
31.             # PropertiesKPointQuality => "Fine"
    
32.             KPointDerivation => "CustomGrid",
    
33.             ParameterA       => 1,
    
34.             ParameterB       => 1,
    
35.             ParameterC       => 1
    
36.         )
    
37.     );
    
38.     $castep->GeometryOptimization->Run($myDoc);
    
39. }
    
40. 
    
41. for ( my $stressz = 0 ; $stressz <= 3 ; ++$stressz ) {
    
42.     &castepcalculation($stressz);  #calculate the results files by castep module     
    
43. 
    
44.     #Processing result file
    
45.     open( result, "$a.castep" || die("open 1.castep fasiled") );
    
46.     open( OUTFILE, ">>Script.txt" );
    
47.     my @line  = ();
    
48.     my $count = 1;
    
49.     my $perlscript = "/Documents/script/try2";
    
50.     my $outputName;
    
51.     @line = <result>;
    
52. 
    
53.     while ( $count <= @line )      #$line[count]
    
54.     {
    
55.         if ($line[$count] =~ /" BFGS : Geometry optimization completed successfully."/ )
    
56.         {
    
57.             open (OUTFILE, ">>Script.txt") || die ("write Script.out fasiled");
    
58.             printf OUTFILE ("stressz=%d Lattice=%2.8f \n",$stressz,$line[$count+17]);
    
59.             opendir(filedir,"$perlscript" || die "open dir-try2 fasiled" );
    
60.             last;
    
61.         }
    
62.         else{
    
63.              if ( $count < @line ){
    
64.                 $count++;
    
65.              }
    
66.              redo &castepcalculation if($count = @line);
    
67.         }
    
68.     }
    
69. }
    

复制代码

凋谢1991

目前发现问题一：
if ($line[$count] =~ /" BFGS : Geometry optimization completed successfully."/ )应将字符串的引号去掉；
目前发现问题二：printf输出的Lattice应改为%s，且利用sunstr函数抓取部分；