如何准备毕业生的.ctl文件将接地的文件转换成.nc文件以批处理方式?- 江南体育网页版- - - - -地球科学堆江南电子竞技平台栈交换 最近30从www.hoelymoley.com 2023 - 07 - 08 - t20:26:15z //www.hoelymoley.com/feeds/question/12135 https://creativecommons.org/licenses/by-sa/4.0/rdf //www.hoelymoley.com/q/12135 2 如何准备毕业生的.ctl文件将接地的文件转换成.nc文件以批处理方式? 阿米特·戈什 //www.hoelymoley.com/users/10815 2017 - 08 - 23 - t09:24:43z 2020 - 03 - 10 - t18:13:15z < p >我有很多研磨(毕业生)文件,我想把这些数据转换成使用cdo netCDF文件。操作需要所有相应的ctl和我需要准备这些文件。所以,我的问题有什么解决方法使用bash, python或r . ctl文件的一个例子:< / p > < pre > <代码> * * _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 < /代码> < / pre > //www.hoelymoley.com/questions/12135/-/12140 # 12140 2 Amit Ghosh回答的如何准备毕业生的.ctl文件将接地的文件转换成.nc文件以批处理方式? 阿米特·戈什 //www.hoelymoley.com/users/10815 2017 - 08 - 23 - t21:14:54z 2017 - 08 - 23 - t21:14:54z < p >经过一些试验n错误我想出来使用下面的shell脚本< / p > < pre > <代码> # !/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。grd TITLE Gridded rainfall UNDEF -999.0 XDEF 135 LINEAR 66.5 0.25 YDEF 129 LINEAR 6.5 0.25 ZDEF 1 linear 1 1 * CHANGE TDEF TO 366 FOR LEAP YEARS TDEF ${days} LINEAR 1jan${j} 1DY VARS 1 rf 0 99 GRIDDED RAINFALL ENDVARS">year${j}_rfp25.ctl done
Baidu
map