Pecan: sipnet.out 中的空格问题

当然，其他人也遇到过这种情况？ 到目前为止，我只在 US-Me2 遇到了这个问题......还有其他人吗？

我目前的“修复”是从使用 PFT“boreal.conifrous”切换到“
tempate.needleleaf.evergreen”无论如何都更准确，但这可能是一个SIPNET错误

我记得遇到过这个错误。 已经有一段时间了，所以我不完全记得细节，但我认为我的解决方案是更改 sipnet 代码本身，允许fprinf更多数字

@istfer谢谢！ @para2x尝试运行 SDA 代码，这仅在我设置 sipnet 运行时发生，特别是在使用 boreal.coniferous PFT 时。 你能用那个 PFT 运行网站吗？ 似乎每当我尝试使用该 PFT 时，它都会因这个问题而失败。 也许您正在使用@istfer在更新她提到的代码后创建的模型二进制文件。

看起来这个问题就像你说的粗根C，特别是价值5.0010000.62 。 这是 5.00 和 10000.62 的串联。 粗根C可以有多大？ 现在它假定最大值为 9999.99。 就像@serbinsh所说的那样，我们可能只想为其添加一些额外的空间以使其正确适应。 这不应破坏任何现有文件。

@robkooper我怀疑问题是山核桃对该 PFT 的参数化，因为我一直在 SIPNET 中找到非常大的结果参数，用于具有一些 pfts 的根，然后可能会在 sipnet 中创建大型根池。 我确实喜欢只增加最大值的想法......尽管这可能只会让更多无意义的数字吐出。

也许您正在使用@istfer在更新她提到的代码后创建的模型二进制文件。

是的，我愿意。 我在 Istems 中创建了我的 xml，我们在 xml 中明确地有一个指向 SIPNET 版本的标记。

是的，这通常是一个问题，因为现在我在运行 SDA 测试时遇到了这个问题....是 BC pft，@para2x？ 想法？

> require (PEcAn.SIPNET)
Loading required package: PEcAn.SIPNET
Loading required package: PEcAn.data.atmosphere
>     model2netcdf.SIPNET('/data/sserbin/Modeling/sipnet/NASA_CMS/SDA/out/ENS-00013-1000025731', 40.6658, -77.904, '1980/01/01', '2010/12/31', FALSE, 'r136')
Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  :
  line 11126 did not have 28 elements
Calls: model2netcdf.SIPNET -> read.table -> scan
Execution halted

我真的不太了解 SIPNET 代码本身，无法在这里提出建议。

不用担心， @para2x我的问题更多是关于你如何使用文件https://github.com/PecanProject/pecan/blob/develop/modules/assim.sequential/inst/MultiSite-Exs/SDA/Bartlett.param可能有一个参数值影响 rootC 导致非常大的值

我使用Bartlett.param的唯一原因是指定第一行的植物初始木材。 您可以在这里尝试两件事：一是不使用该文件或仅保留第一行。

BU 的人，您能告诉我您的 sipnet.c 文件在第 782-785 和 813 行上是什么样子吗？ 那是针对 BETY 型号 ID 编号 1000000014。我想看看这是不是已编辑输出的版本

我觉得这就是

LAI 之后标题中的最后一个变量名称被某种方式切割，但它应该是木头创建的（请注意，这可能不在其他 sipnet.c 文件中，这是我创建的另一个跟踪器（我应该将它移到另一个 sipnet 下）版本）

void outputHeader(FILE *out) {
  fprintf(out, "Notes: (PlantWoodC, PlantLeafC, Soil and Litter in g C/m^2; Water and Snow in cm; SoilWetness is fraction of WHC;\n");
  fprintf(out, "loc year day time plantWoodC plantLeafC ");

        #if SOIL_MULTIPOOL
            int counter;

            for(counter=0; counter<NUMBER_SOIL_CARBON_POOLS; counter++) {
                fprintf(out, "soil(%8.2f) ",envi.soil[counter]);
            }
                fprintf(out,"totSoilC ");

        #else
            fprintf(out, "soil ");
        #endif

            fprintf(out, "microbeC ");
            fprintf(out, "coarseRootC fineRootC ");
  fprintf(out, "litter litterWater soilWater soilWetnessFrac snow ");
  fprintf(out, "npp nee cumNEE gpp rAboveground rSoil rRoot ra rh rtot evapotranspiration fluxestranspiration fPAR LAI ation\n");
}

// pre: out is open for writing
// print current state to output file
void outputState(FILE *out, int loc, int year, int day, double time) {


        fprintf(out,"%8d %4d %3d %5.2f %8.2f %8.2f ",loc,year,day,time,envi.plantWoodC,envi.plantLeafC);


        #if SOIL_MULTIPOOL
            int counter;

            for(counter=0; counter<NUMBER_SOIL_CARBON_POOLS; counter++) {
                fprintf(out, "%8.2f ",envi.soil[counter]);
            }
                fprintf(out,"%8.2f ",trackers.totSoilC);



        #else
            fprintf(out, "%8.2f ",envi.soil);
        #endif



            fprintf(out, "%8.2f ", envi.microbeC);

            fprintf(out, "%8.2f %8.2f", envi.coarseRootC,envi.fineRootC);

    fprintf(out, " %8.2f %8.3f %8.2f %8.3f %8.2f ",
         envi.litter, envi.litterWater, envi.soilWater, trackers.soilWetnessFrac, envi.snow);
    fprintf(out,"%8.2f %8.2f %8.2f %8.8f %8.3f %8.3f %8.3f %8.3f %8.3f %8.3f %8.8f %8.4f %8.4f %8.4f %8.4f\n", trackers.npp, trackers.nee, trackers.totNee, trackers.gpp, trackers.rAboveground,
           trackers.rSoil, trackers.rRoot, trackers.ra, trackers.rh, trackers.rtot, trackers.evapotranspiration, fluxes.transpiration, trackers.fpar, trackers.LAI, trackers.wcr);

//note without modeling root dynamics

//trackers.fa, trackers.fr, fluxes.rLeaf*climate->length,trackers.evapotranspiration

}

是的，在我的 sipnet.c 中不一样，它是基于 r136 的

第 787-840 行

// pre: out is open for writing
// print current state to output file
void outputState(FILE *out, int loc, int year, int day, double time) {


                fprintf(out,"%8d %4d %3d %5.2f %8.2f %8.2f ",loc,year,day,time,envi.plantWoodC,envi.plantLeafC);


                #if SOIL_MULTIPOOL
                        int counter;

                        for(counter=0; counter<NUMBER_SOIL_CARBON_POOLS; counter++) {
                                fprintf(out, "%8.2f ",envi.soil[counter]);
                        }
                                fprintf(out,"%8.2f ",trackers.totSoilC);



                #else
                        fprintf(out, "%8.2f ",envi.soil);
                #endif



                        fprintf(out, "%8.2f", envi.microbeC);

                        fprintf(out, "%8.2f %8.2f", envi.coarseRootC,envi.fineRootC);

        fprintf(out, " %8.2f %8.3f %8.2f %8.3f %8.2f ",
                 envi.litter, envi.litterWater, envi.soilWater, trackers.soilWetnessFrac, envi.snow);
        fprintf(out,"%8.2f %8.2f %8.2f %8.2f %8.3f %8.3f %8.3f %8.3f %8.3f %8.3f %8.8f %8.4f %8.4f\n", trackers.npp, trackers.nee, trackers.totNee, trackers.gpp, trackers.rAboveground,
           trackers.rSoil, trackers.rRoot, trackers.ra, trackers.rh, trackers.rtot, trackers.evapotranspiration, fluxes.transpiration, trackers.fpar);

//note without modeling root dynamics

//trackers.fa, trackers.fr, fluxes.rLeaf*climate->length,trackers.evapotranspiration
}

void outputStatecsv(FILE *out, int loc, int year, int day, double time) {
  fprintf(out, "%8d , %4d , %3d , %5.2f , %8.2f , %8.2f , ", loc, year, day, time,
                envi.plantWoodC, envi.plantLeafC);

        #if SOIL_MULTIPOOL
                fprintf(out, "%8.2f ,", trackers.totSoilC);
        #else
                fprintf(out, "%8.2f ,", envi.soil);
        #endif

        fprintf(out, "%8.2f , %8.3f, %8.2f , %8.3f , %8.2f , %8.2f , %8.2f , %8.2f , %8.2f , %8.3f , %8.3f , %8.3f, %8.3f , %8.3f , %8.3f %8.8f %8.4f %8.4f\n",
          envi.litter, envi.litterWater, envi.soilWater, trackers.soilWetnessFrac, envi.snow,
          trackers.npp, trackers.nee, trackers.totNee, trackers.gpp, trackers.rAboveground, trackers.rSoil, trackers.rRoot, trackers.ra, trackers.rh, trackers.rtot, trackers.evapotranspiration, fluxes.transpiration, trackers.fpar);

}

我将首先尝试将 sipnet.C 中的粗根 C 的输出大小增加到 9.2f 之类的值； 不确定我们是否应该将其增加得更大

我的直觉是 write.table 有一个特定的浮点数格式（4.2%f 是它的接缝），它导致一个无效的文件（仍在试验这个）。

@robkooper不确定我理解.....你的意思是 read.table 吗？ 我更改为 9.4f 并且到目前为止还没有遇到相同站点/pft 组合的错误

也就是说，在 sipnet.C 文件中为粗根 C 输出变量

是的 。 没有帮助。 我想我们只需要增加空格

试图弄清楚您粘贴的文件是如何创建的。 这是由 PEcAn 创建的吗？ 我尝试运行相同的运行并生成输入文件（https://pecan-docker.ncsa.illinois.edu/minio/dbfiles/CRUNCEP_SIPNET_site_0-763/）我在那里的文件中没有看到这个问题。 数字应该用制表符分隔，而不是空格。

@bailsofhay仅供参考 - 这可能是您神秘的 SIPNET 崩溃的原因。 我们应该进一步调查这个

尝试通过更新 sipnet.C 来解决此问题

        fprintf(out,"%8.2f %8.2f %8.2f %8.8f %8.3f %8.3f %8.3f %8.3f %8.3f %8.3f %8.8f %8.4f %8.4f %8.4f\n", trackers.npp, trackers.nee, trackers.totNee, trackers.gpp, trackers.rAboveground,
           trackers.rSoil, trackers.rRoot, trackers.ra, trackers.rh, trackers.rtot, trackers.evapotranspiration, fluxes.transpiration, trackers.fpar, trackers.LAI);

4 改为 8.8

还添加了基于@istfer示例的

我更新了标题：

void outputHeader(FILE *out) {
  fprintf(out, "Notes: (PlantWoodC, PlantLeafC, Soil and Litter in g C/m^2; Water and Snow in cm; SoilWetness is fraction of WHC);\n");
  fprintf(out, "loc year day time plantWoodC plantLeafC ");

        #if SOIL_MULTIPOOL
                        int counter;

                        for(counter=0; counter<NUMBER_SOIL_CARBON_POOLS; counter++) {
                                fprintf(out, "soil(%8.2f) ",envi.soil[counter]);
                        }
                                fprintf(out,"totSoilC ");

                #else
                        fprintf(out, "soil ");
                #endif

                        fprintf(out, "microbeC ");
                        fprintf(out, "coarseRootC fineRootC ");
  fprintf(out, "litter litterWater soilWater soilWetnessFrac snow ");
  fprintf(out, "npp nee cumNEE gpp rAboveground rSoil rRoot ra rh rtot evapotranspiration fluxestranspiration fPAR LAI\n");
}

好的，这已经奏效了，除了我发现 PEcAn 在 nc 文件中为 LAI 提供的内容与 SIPNET 输出的 LAI 列中的内容之间存在脱节。 我知道这是@istfer和@mdietze之前指出的一个问题，但我忘记了我们应该如何解决这个问题。 所以 1) PEcAn 在哪里获得 LAI，它是不是在事后使用其他输出数据进行计算，而不是假设 LAI 是从 SIPNET 输出的？ 2）我们应该使用计算出的LAI还是在添加trackers.LAI后从SIPNET吐出的LAI？ 3）对于SDA，我们应该使用它，因为这可能是SIPNET运行中某些奇怪的原因

例如，这里是 SIPET 中的 2010 LAI 输出（来自 sipnet.out），介于 0.2 和 ~23 m2/m2 之间，但这里是 PEcAN 输出 LAI

@istfer您是否使用过 SIPNET LAI 或 PEcAn 计算的 LAI？

好吧，嗯...还有 SIPNET LAI 不会到 0...经常在 0.1 或 0.2 左右徘徊，但有时在冬天高达 0.5？！ 对于阔叶网站

我今天会和@istem谈谈，看看她是否记得对

这是一个我认为与此相关的新问题，因为我们更改了 signet.c 文件

loc year day time plantWoodC plantLeafC soil microbeC coarseRootC fineRootC litter litterWater soilWater soilWetnessFrac snow npp nee cumNEE gpp rAboveground rSoil rRoot ra rh rtot evapotranspiration fluxestranspiration fPAR LAI
       0 2002 197  0.00 15311.11     0.00 143942.45    71.97 911070.92 23093.15   280.00    0.500     6.89    0.537     0.08   -17.83    64.48    64.48 0.00000000    0.128   64.356   17.699   17.827   46.657   64.484 0.00697267   0.0000   0.0000   0.0000
       0 2002 197  6.00 15310.95     0.00 143943.92    71.97 911018.37 23072.85   280.00    0.500     6.96    0.577     0.00   -17.85    71.53   136.02 0.00000000    0.126   71.408   17.721   17.846   53.687   71.533 0.00000001   0.0000   0.0000   0.0000
       0 2002 197 12.00 15310.74     0.00 143944.66    71.97 910965.80 23052.57   280.00    0.500     6.92    0.578     0.00   -17.92    72.31   208.33 -0.00000000    0.165   72.148   17.752   17.917   54.396   72.313 0.03588458   0.0000   0.0000   0.0000
       0 2002 197 18.00 15310.55     0.00 143946.14    71.97 910913.31 23032.32   280.00    0.500     6.89    0.575     0.00   -17.82    71.45   279.78 -0.00000000    0.156   71.297   17.661   17.816   53.637   71.453 0.03566278   0.0000   0.0000   0.0000
       0 2002 198  0.00 15310.39     0.00 143948.22    71.97 910860.90 23012.08   280.00    0.500     6.85    0.572     0.00   -17.72    70.74   350.53 0.00000000    0.128   70.616   17.596   17.723   53.020   70.743 0.03615626   0.0000   0.0000   0.0000
       0 2002 198  6.00 15310.23     0.00 143950.43    71.97 910808.46 22991.85   280.00    0.500     6.81    0.569     0.00   -17.74    70.61   421.14 0.00000000    0.125   70.483   17.618   17.743   52.866   70.609 0.03677418   0.0000   0.0000   0.0000
       0 2002 198 12.00 15310.01     0.00 143952.75    71.97 910755.99 22971.65   280.00    0.500     6.79    0.567     0.00   -17.82    70.57   491.71 -0.00000000    0.176   70.393   17.647   17.824   52.746   70.569 0.02752184   0.0000   0.0000   0.0000
       0 2002 198 18.00 15309.82     0.00 143955.88    71.97 910703.65 22951.46   280.00    0.500     6.76    0.564     0.00   -17.69    69.60   561.30 -0.00000000    0.160   69.436   17.527   17.687   51.909   69.596 0.02520775   0.0000   0.0000   0.0000
       0 2002 199  0.00 15309.65     0.00 143959.58    71.97 910651.39 22931.30   280.00    0.500     6.73    0.562     0.00   -17.58    68.90   630.20 0.00000000    0.135   68.762   17.449   17.584   51.314   68.897 0.02564393   0.0000   0.0000   0.0000
       0 2002 199  6.00 15309.48     0.00 143963.49    71.97 910599.13 22911.15   280.00    0.500     6.71    0.560     0.00   -17.58    68.67   698.87 0.00000000    0.134   68.540   17.443   17.577   51.097   68.674 0.02672652   0.0000   0.0000   0.0000
       0 2002 199 12.00 15309.28     0.00 143967.57    71.97 910546.87 22891.02   280.00    0.500     6.87    0.566     0.00   -17.60    68.50   767.37 0.00000000    0.159   68.341   17.441   17.600   50.900   68.501 0.03973389   0.0000   0.0000   0.0000
       0 2002 199 18.00 15309.10     0.00 143970.78    71.97 910494.70 22870.90   280.00    0.500     6.99    0.578     0.00   -17.51    69.26   836.63 0.00000000    0.145   69.112   17.363   17.508   51.748   69.256 0.03744034   0.0000   0.0000   0.0000

为什么我们得到 -0.000 的 AGB 值？？？？

就这样结束了。 实际值为负但 > -0.0000000005。 并且列标题是分隔的，没有对齐，所以我认为这实际上不是 AGB

现在将其更改为 sci 符号。 修复后会更新

将第 19 列（gpp）更改为：%8.3e

好的，所以实际上我之前犯了一个错误。 我实际上需要做的是修复两个不同的列。

第一个解决最初的问题，第 8 列和第 9 列有时会合并，因为没有足够的空间，第二个是 GPP 的表示法（第 19 列）。

@serbinsh看起来我们还需要在科学记数法问题之上更多地考虑 GPP。 我正在查看我的 signet.out 文件以查找失败的集合，这次运行中没有负数，但是当 GPP 值从 1 位变为 10 位时，列间距在今年晚些时候搞砸了。 这是一个示例（列与我上面发布的示例相同：
0 2000 221 18.00 14575.71 113.85 56661.11 17.39 0.00 285788.90 280.00 0.500 10.78 0.898 0.00 4.43 18.77 1947.13 7.93938489 2.561 24.148 0.951 3.512 23.197 26.709 0.17197340 0.6069 0.5578 3.8823 0 2000 222 0.00 14579.71 114.58 56875.78 17.39 0.00 285551.70 280.00 0.500 10.69 0.894 0.00 5.23 17.92 1965.05 8.50201393 2.330 24.090 0.941 3.271 23.149 26.420 0.09273126 0.3709 0.5594 3.9070 0 2000 222 6.00 14583.60 115.44 57090.45 17.39 0.00 285314.87 280.00 0.500 10.59 0.887 0.00 5.41 17.54 1982.59 8.63300350 2.283 23.893 0.935 3.218 22.957 26.176 0.09346984 0.3566 0.5613 3.9363 0 2000 222 12.00 14585.45 116.43 57305.10 17.39 0.00 285078.41 280.00 0.500 10.68 0.886 0.00 3.68 19.10 2001.69 7.51178589 2.897 23.716 0.931 3.828 22.785 26.612 0.20875203 0.7065 0.5635 3.9701 0 2000 222 18.00 14588.05 117.52 57519.55 17.39 0.00 284842.29 280.00 0.500 10.78 0.894 0.00 4.67 18.13 2019.82 8.21516699 2.635 23.712 0.915 3.550 22.797 26.347 0.17083126 0.5641 0.5660 4.0074 0 2000 223 0.00 14591.06 118.73 57733.69 17.39 0.00 284606.53 280.00 0.500 10.68 0.894 0.00 5.36 17.57 2037.39 8.70971254 2.447 23.830 0.904 3.351 22.926 26.277 0.10070652 0.4028 0.5688 4.0487 0 2000 223 6.00 14594.25 120.08 57947.87 17.39 0.00 284371.14 280.00 0.500 10.59 0.886 0.00 5.84 16.85 2054.24 9.03531674 2.297 23.590 0.898 3.195 22.693 25.888 0.08617136 0.3222 0.5718 4.0945 0 2000 223 12.00 14595.80 121.57 58162.01 17.39 0.00 284136.13 280.00 0.500 10.74 0.889 0.00 4.53 18.02 2072.27 8.22739039 2.807 23.443 0.894 3.702 22.549 26.251 0.20809168 0.6766 0.5752 4.1453 0 2000 223 18.00 14597.85 123.17 58375.78 17.39 0.00 283901.48 280.00 0.500 10.96 0.904 0.00 5.29 17.45 2089.72 8.77340379 2.601 23.622 0.882 3.483 22.741 26.223 0.19588011 0.5974 0.5788 4.2000 0 2000 224 0.00 14600.46 124.91 58588.93 17.39 0.00 283667.19 280.00 0.500 11.04 0.917 0.00 6.16 17.01 2106.73 9.39612713 2.366 24.045 0.873 3.239 23.171 26.411 0.10046154 0.3110 0.5827 4.2592 0 2000 224 6.00 14603.09 126.79 58801.70 17.39 0.00 283433.30 280.00 0.500 10.96 0.916 0.00 6.52 16.85 2123.58 9.65203856 2.265 24.233 0.870 3.135 23.363 26.497 0.08316609 0.2768 0.5868 4.3235 0 2000 224 12.00 14604.52 128.70 59014.36 17.39 0.00 283199.63 280.00 0.500 11.10 0.919 0.00 5.37 17.92 2141.49 9.03953546 2.803 24.153 0.869 3.672 23.284 26.956 0.18632810 0.6005 0.5912 4.3884 0 2000 224 18.00 14606.58 130.63 59226.61 17.39 0.00 282966.20 280.00 0.500 11.23 0.930 0.00 6.07 17.43 2158.92 9.52919823 2.598 24.361 0.859 3.457 23.502 26.959 0.17821633 0.5664 0.5956 4.4545 0 2000 225 0.00 14609.44 132.61 59438.41 17.39 0.00 282733.02 280.00 0.500 11.21 0.935 0.00 6.95 16.80 2175.73 10.11249382 2.305 24.612 0.853 3.158 23.759 26.917 0.12228206 0.3314 0.6000 4.5217 0 2000 225 6.00 14612.64 134.62 59649.97 17.39 0.00 282500.08 280.00 0.500 11.09 0.929 0.00 7.40 16.42 2192.15 10.36770066 2.118 24.668 0.854 2.972 23.814 26.786 0.11599917 0.2490 0.6044 4.5903 0 2000 225 12.00 14614.56 136.67 59861.39 17.39 0.00 282267.39 280.00 0.500 10.89 0.916 0.00 6.20 17.55 2209.70 9.75227886 2.691 24.614 0.859 3.549 23.756 27.305 0.23872529 0.8320 0.6088 4.6604 0 2000 225 18.00 14617.15 138.80 60073.06 17.39 0.00 282035.00 280.00 0.500 10.74 0.901 0.00 7.05 16.27 2225.97 10.32246739 2.422 24.175 0.853 3.275 23.322 26.597 0.23678459 0.7050 0.6132 4.7330 0 2000 226 0.00 14620.34 140.99 60284.78 17.39 0.00 281802.87 280.00 0.500 10.61 0.889 0.00 7.78 15.30 2241.28 10.81901992 2.186 23.938 0.853 3.039 23.085 26.124 0.14153777 0.3429 0.6176 4.8077 0 2000 226 6.00 14623.66 143.24 60496.39 17.39 0.00 281571.01 280.00 0.500 10.47 0.879 0.00 8.03 14.98 2256.26 10.99722667 2.109 23.865 0.858 2.968 23.007 25.975 0.13726935 0.2748 0.6221 4.8842 0 2000 226 12.00 14623.64 145.53 60707.87 17.39 0.00 281339.40 280.00 0.500 10.35 0.868 0.00 4.81 18.13 2274.39 8.49493616 2.818 23.811 0.866 3.685 22.944 26.629 0.28877915 0.9217 0.6267 4.9626 0 2000 226 18.00 14625.66 147.86 60919.42 17.39 0.00 281108.02 280.00 0.500 10.26 0.859 0.00 6.91 15.78 2290.17 10.38533933 2.612 23.551 0.862 3.474 22.690 26.164 0.29844127 0.9112 0.6312 5.0419 0 2000 227 0.00 14628.99 150.24 61130.90 17.39 0.00 280876.90 280.00 0.500 10.19 0.852 0.00 8.34 14.22 2304.39 11.49566120 2.293 23.421 0.861 3.154 22.560 25.714 0.16840103 0.4209 0.6357 5.1229 0 2000 227 6.00 14632.40 152.68 61342.16 17.39 0.00 280646.05 280.00 0.500 10.03 0.842 0.00 8.56 14.04 2318.43 11.70045859 2.271 23.474 0.867 3.139 22.607 25.745 0.15826026 0.3612 0.6403 5.2061 0 2000 227 12.00 14629.28 155.18 61553.35 17.39 0.00 280415.47 280.00 0.500 9.89 0.830 0.00 2.16 20.32 2338.75 6.26260256 3.230 23.352 0.875 4.105 22.477 26.582 0.27769493 0.8826 0.6448 5.2914 0 2000 227 18.00 14626.96 157.65 61764.71 17.39 0.00 280185.05 280.00 0.500 9.65 0.815 0.00 2.90 19.22 2357.97 6.84262945 3.074 22.989 0.867 3.941 22.122 26.063 0.26190036 0.8707 0.6494 5.3757 0 2000 228 0.00 14630.43 160.08 61976.42 17.39 0.00 279954.77 280.00 0.500 9.43 0.795 0.00 8.61 12.97 2370.94 12.07581213 2.600 22.446 0.862 3.462 21.585 25.046 0.22249240 0.7077 0.6538 5.4586 0 2000 228 6.00 14634.17 162.58 62188.30 17.39 0.00 279724.77 280.00 0.500 9.27 0.779 0.00 9.03 12.19 2383.14 12.42523436 2.528 22.092 0.865 3.393 21.227 24.619 0.16641807 0.5000 0.6580 5.5438 0 2000 228 12.00 14627.63 165.14 62400.19 17.39 0.00 279495.04 280.00 0.500 9.05 0.763 0.00 -1.10 22.13 2405.26 3.72972367 3.961 21.894 0.870 4.831 21.024 25.855 0.22022551 0.8153 0.6623 5.6311 0 2000 228 18.00 14622.48 167.59 62612.55 17.39 0.00 279265.37 280.00 0.500 8.82 0.744 0.00 0.03 20.34 2425.60 4.51935059 3.632 21.223 0.856 4.488 20.367 24.855 0.22887585 0.7960 0.6666 5.7145 0 2000 229 0.00 14625.24 169.92 62825.31 17.39 0.00 279035.75 280.00 0.500 8.60 0.725 0.00 7.72 12.05 2437.65 11.61133388 3.047 20.616 0.846 3.893 19.771 23.663 0.22073243 0.7758 0.6706 5.7941 0 2000 229 6.00 14629.30 172.29 63038.32 17.39 0.00 278806.37 280.00 0.500 8.40 0.708 0.00 9.10 10.23 2447.88 12.96732463 3.023 20.172 0.844 3.867 19.328 23.195 0.19535639 0.6865 0.6743 5.8748 0 2000 229 12.00 14621.47 174.70 63251.53 17.39 0.00 278577.25 280.00 0.500 8.20 0.692 0.00 -2.68 21.64 2469.51 2.74794197 4.588 19.798 0.843 5.432 18.954 24.386 0.20192009 0.7392 0.6781 5.9572 0 2000 229 18.00 14616.93 176.95 63465.20 17.39 0.00 278348.13 280.00 0.500 8.38 0.691 0.00 0.26 18.03 2487.54 4.97890595 3.889 19.115 0.825 4.714 18.290 23.004 0.23649558 0.7214 0.6818 6.0340 0 2000 230 0.00 14622.07 179.08 63678.37 17.39 0.00 278119.05 280.00 0.500 8.22 0.691 0.00 9.69 8.91 2496.45 13.60916949 3.107 19.411 0.814 3.921 18.598 22.519 0.16278224 0.5916 0.6852 6.1066 0 2000 230 6.00 14627.67 181.27 63891.65 17.39 0.00 277890.24 280.00 0.500 8.13 0.681 0.00 10.26 8.04 2504.49 13.96623695 2.889 19.113 0.813 3.703 18.300 22.003 0.12493134 0.4129 0.6884 6.1811 0 2000 230 12.00 14620.98 183.51 64104.82 17.39 0.00 277661.68 280.00 0.500 7.94 0.670 0.00 -1.88 20.11 2524.60 3.22887897 4.295 19.047 0.816 5.112 18.230 23.342 0.18911457 0.7155 0.6916 6.2575

@bailsofhay这实际上不是问题，而且经常发生。 我发布的问题更严重，1) 列之间没有足够的间距，以至于两列合并在一起，然后 PEcAN 找不到正确的列数。 另一个问题是 GPP 可能存在的舍入误差。 您上面显示的内容对于 sipnet 非常常见，到目前为止我已经忽略了它，但也许这是一个问题？ 我以为没有。 其他的？

@serbinsh我认为这是错误的。 如果您查看小数点前位数增加的其他一些列，小数点保持对齐，小数点左侧的数字移动。 GPP 则相反，小数点左边的数字排成一排，而小数点右边的数字不再排成一排。 您可能是对的，这实际上不是问题，但有可能是导致某些运行崩溃的问题。 我会看看一些成功的运行，看看它们是否有数字排​​列的变化，如果它们存在于其他文件中，那么这可能不是问题

@bailsofhay我认为这不是问题的原因是，如果列数没有达到 28 或对齐，PEcAn 将失败。 所以在日志中寻找那个错误，如果你没有发现它不是问题

另外我只是更新 sipnet 以使用 GPP 的 sci 表示法，这也可能解决这个问题，请尝试重新运行您的 SDA

kk会的。 我还检查了其他一些运行，并且有像失败一样的数字变化，但是所有的空格都在那里，所以你是对的。 那不是问题。

这个问题是陈旧的，因为它已经开放了 365 天，没有任何活动。