diff --git a/NAERSTADMO.py b/NAERSTADMO.py index 9dcec5512ecc880a13b6e6b202b90738daf32603..035ce094a51538d7c3016aa64d06feec3c0de40f 100644 --- a/NAERSTADMO.py +++ b/NAERSTADMO.py @@ -185,7 +185,7 @@ def create_matrix(): date = datetime.fromtimestamp(timestep[0]).astimezone(local_timezone).strftime('%Y-%m-%d') hour = datetime.fromtimestamp(timestep[0]).astimezone(local_timezone).strftime('%H') - subprocess.run(f'cdo -expr,"TSWH=(WHS==0)?0:sumWHTM;IR=(TSWH>=40?1:0);RISK=(IR==1?(TSWH/40)*VRS:0);VRS=(VRS==0?0:VRS);" {tmpfile_path}bg_final_{str(time)}.nc {tmpfile_path}final_{str(date)}_{str(hour)}.nc' , shell=True) + subprocess.run(f'cdo -expr,"TSWH=(WHS==0)?0:sumWHTM;IR=(TSWH>=40?1:0);RISK=(IR==1?(TSWH/40)*VRS:0);VRS=(VRS==0?0:VRS);" {tmpfile_path}bg_final_{str(time)}.nc {outfile_path}final_{str(date)}_{str(hour)}.nc' , shell=True) if (str(hour) == '00' or time == 1): subprocess.run(f'cp {tmpfile_path}final_{str(date)}_{str(hour)}.nc {outfile_path}final_{str(date)}.nc', shell=True)