- 积分
- 7
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-12-4
- 最后登录
- 1970-1-1
|
发表于 2023-12-3 12:55:50
|
显示全部楼层
各位前辈,请问我运行create_UV.perl的时候,没有报错,但是目标文件中没有生成UV文件
create_UV.perl的内容如下:
#!/usr/bin/perl
$filename = "create_UV";
system ("mpif90 -o $filename $filename.f90 -L/public/home/fdu_qiusy/qiusy/hrldas-release-3.7/HRLDAS/w3lib-2.2.0 -lw3");
@nums = ("00","01","02","03","04","05","06","07","08","09",
"10","11","12","13","14","15","16","17","18","19",
"20","21","22","23","24","25","26","27","28","29",
"30","31");
@yrs = ("16");
$day_start = 279;
$day_end = 281;
@hrs = ("00","03","06","09","12","15","18","21");
@noleap_days = (0,31,59,90,120,151,181,212,243,273,304,334,365);
@leap_days = (0,31,60,91,121,152,182,213,244,274,305,335,366);
$data_dir = "/public/home/fdu_qiusy/qiusy/hrldas-release-3.7/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/extracted_GLADS_data_path";
$results_dir = "/public/home/fdu_qiusy/qiusy/hrldas-release-3.7/HRLDAS/HRLDAS_forcing/run/examples/GLDAS/extracted_GLDAS_data_path";
for $yy (@yrs)
{
# This will be the jday time loop
for($julday=$day_start;$julday<=$day_end;$julday++)
{
# This little section finds the text month and day
@modays = @noleap_days;
if($yy == "00" || $yy == "04" || $yy == "08" || $yy == "12") {@modays = @leap_days}
for($mo=1;$mo<=12;$mo++)
{
if($julday>$modays[$mo-1] && $julday<=$modays[$mo])
{
$mon = $mo;
$day = $julday - $modays[$mo-1];
}
}
for $hr (@hrs)
{
$infile = "$data_dir/Wind/GLDAS_Wind.20$yy$nums[$mon]$nums[$day]$hr.grb";
$outfile1 = " $results_dir/U/GLDAS_U.20$yy$nums[$mon]$nums[$day]$hr.grb";
$outfile2 = " $results_dir/V/GLDAS_V.20$yy$nums[$mon]$nums[$day]$hr.grb";
print ("$infile \n");
system ("$filename $infile $outfile1 $outfile2");
}
}
} # End of outer time loop
system("rm $filename"); |
|