MOZBC报错trying to increment file 已解决

klingel

使用CAM-Chem数据，使用MOZBC来制作边界场和初始场，如果你遇到了和我一样的问题，如下：我估计你下的时间和你模拟的时间不一致，因为这个报错的本根原因是，mozbc觉得你的CAM-Chem数据时间和你的wrfbdy时间不匹配，然后增加新的文件去看看有没有时间匹配，最后发现没有这个文件。其实我最开始知道是时间出问题了，才有这个报错，但是时间到底怎么样才能对呢，通过对比朋友成功的案例，我把time这个维度的数据明明改成了我想要的时间了，并且在前后多加了5天时间，还是报错。通过不断的查找，发现很多人与我遇到了一样的问题，他们可能重新下载了对应模拟时间的数据，就好了，可是我的模拟时间是2022年，没有这个时间的数据啊，好像CAM-Chem目前最新只有2020年12月底之前的，我只能改时间了，最后发现，这个文件里面有个叫date和time_bnds的变量。这两个变量的值也必须得改掉。终于试错，成功了！！！看到bc_wrfchem completed successfully.

janelle713

请问怎么批量修改文件的时间呀

上帝专用小马甲

谢谢，亲测可用。

Enthusiasm

谢谢楼主的提醒，我想来更新一下这问题。

首先，我也遇到了同样的问题，时间不匹配。我用的是WACCM的数据，我根据楼主的提示首先修改了date、time_bnds、time的变量。但运行时发现仍报错。后来突然间看到了datesec的秒的变量，突发奇想也做了一下修改，就成功运行了。
对比了一下CAM和WACCM的时间格式，我感觉WACCM的格式可能更加容易修改。
供大家参考！

nuist_zhangk

f=addfile("camchem-20240419054837203305.nc", "w")

; 读取变量数据
date_data = f->date;
time_bnds_data = f->time_bnds;
time_data = f->time;

print(date_data)
print(time_bnds_data)
print(time_data)

time_data := (/7914, 7914.25, 7914.5, 7914.75, 7915, 7915.25, 7915.5, 7915.75, 7916, 7916.25, 7916.5, 7916.75, 7917, 7917.25, 7917.5, 7917.75, 7918, 7918.25, \
              7918.5, 7918.75, 7919, 7919.25, 7919.5, 7919.75, 7920, 7920.25, 7920.5, 7920.75, 7921/)
time_bnds_data = (/ (/7913.75, 7914/), \
                   (/7914, 7914.25/), \
                   (/7914.25, 7914.5/), \
                   (/7914.5, 7914.75/), \
                   (/7914.75, 7915/), \
                   (/7915, 7915.25/), \
                   (/7915.25, 7915.5/), \
                   (/7915.5, 7915.75/), \
                   (/7915.75, 7916/), \
                   (/7916, 7916.25/), \
                   (/7916.25, 7916.5/), \
                   (/7916.5, 7916.75/), \
                   (/7916.75, 7917/), \
                   (/7917, 7917.25/), \
                   (/7917.25, 7917.5/), \
                   (/7917.5, 7917.75/), \
                   (/7917.75, 7918/), \
                   (/7918, 7918.25/), \
                   (/7918.25, 7918.5/), \
                   (/7918.5, 7918.75/), \
                   (/7918.75, 7919/), \
                   (/7919, 7919.25/), \
                   (/7919.25, 7919.5/), \
                   (/7919.5, 7919.75/), \
                   (/7919.75, 7920/), \
                   (/7920, 7920.25/), \
                   (/7920.25, 7920.5/), \
                   (/7920.5, 7920.75/), \
                   (/7920.75, 7921/) /)
date_data := (/20231030,20231030,20231030,20231030,20231031,20231031,20231031,20231031,20231101,20231101,20231101,20231101,20231102,20231102,20231102,20231102,20231103,20231103,20231103,20231103,20231104,20231104,20231104,20231104,20231105,20231105,20231105,20231105,20231106/)

f->date = date_data;
f->time_bnds = time_bnds_data;
f->time = time_data;

print(date_data)
print(time_bnds_data)
print(time_data)
~
接着上面的话，这是如何修改变量的脚本，我今天也遇见这个问题了，希望对大家有帮助
time是指从2003年3.1开始到当前时间的天数