如何准备毕业生的.ctl文件将接地的文件转换成.nc文件以批处理方式?

Question

我有很多研磨(毕业生)文件,我想把这些数据转换成使用cdo netCDF文件。操作需要所有相应的ctl和我需要准备这些文件。所以,我的问题有什么解决方法使用bash, python或r . ctl文件的一个例子:

* * _rfp25需要更改文件名年。为每个文件研磨DSET /home/amit/climate_data /射频/ year1901_rfp25。接地的标题网格降雨UNDEF -999.0 XDEF 135线性线性6.5 - 0.25 66.5 - 0.25 129年YDEF ZDEF 1线性1 1 * TDEF更改为366 365年闰年TDEF线性jan1901 1 dy var 1射频0 99网格降雨ENDVARS

Answer 1

经过一些试验n错误我想使用下面的shell脚本

# !/bin/bash day_number(){年= 1美元的飞跃=“false”如果(((年% 4))- ne 0美元);然后:elif[$((% 400年))情商0];然后跳= " true " elif[$((% 100年))情商0];其他:飞跃= " true " fi如果[" $飞跃" = = " true ");然后回声还有366与365 j的fi} {1901 . . 2015};做天= $ (day_number j)回声“DSET /home/amit/climate_data /射频/年$ {j} _rfp25。接地的标题网格降雨UNDEF -999.0 XDEF 135线性线性6.5 - 0.25 66.5 - 0.25 129年YDEF ZDEF 1线性1 1 * TDEF更改为366美元闰年TDEF线性1 1月$ {j}{天}1 dy var 1射频0 99网格降雨ENDVARS " > {j} _rfp25美元。ctl做