Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

ARWpost: WARNING --- I do not recognize this data.

L

loly_

New member
I have a problem with post processing. I am using ARWpost_V3 with gfs 0.5 data. This is the error that appeared when trying to run /ARWpost.exe :
START PROCESSING DATA

WARNING --- I do not recognize this data.

Will make an attempt to read it.

At line 150 of file process_domain_module.f (unit = 11)
Fortran runtime error: RECL parameter is non-positive in OPEN statement

Error termination. Backtrace:
#0 0x7f0954bf6d21 in ???
#1 0x7f0954bf7869 in ???
#2 0x7f0954bf854f in ???
#3 0x7f0954e37a03 in ???
#4 0x7f0954e37ccc in ???
#5 0x5591778514bf in ???
#6 0x559177841f8a in ???
#7 0x7f09548b9082 in ???
#8 0x5591778405fd in ???
#9 0xffffffffffffffff in ???

I'm working with data gfs: gfs.t00z.pgrb2.0p50.f000
 

Attachments

  • namelist.input
    3.2 KB · Views: 7
  • rsl.out.0000
    267.6 KB · Views: 1
I know that the file to be used in ARWpost is the one generated after executing ./wrf.exe that generated me wrfout_d01_2023-06-27_00:00:00. But when trying to run the command ./ARWpost.exe. it returns an error message.

namelist.ARWpost
&datetime
start_date = '2023-06-27_00:00:00',
end_date = '2023-06-28_00:00:00',
interval_seconds = 10800,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/home/clima/Build_WRF/WRF-4.5-ARW/test/em_real/wrfout_d01_2023-06-27_00:00:00'
output_root_name = '/home/clima/Build_WRF/ARWpost/wrfd1_2023062700'
plot = 'all_list'
fields = 'height, geopt, theta, tc, td2, td, rh, rh2, umet, vmet,
slp, pressure, u10m, v10m, wdir, wspd, wd10, ws10, cape, mcin,
lcl, lfc, cape, cin, clfr'
output_type = 'grads'
mercator_defs = .true.
/

split_output = .true.
frames_per_outfile = 2

plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10, slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'

&interp
interp_method = 1,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.

interp_method = 1, ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,
 
loly_ said:
I know that the file to be used in ARWpost is the one generated after executing ./wrf.exe that generated me wrfout_d01_2023-06-27_00:00:00. But when trying to run the command ./ARWpost.exe. it returns an error message.

namelist.ARWpost
&datetime
start_date = '2023-06-27_00:00:00',
end_date = '2023-06-28_00:00:00',
interval_seconds = 10800,
tacc = 0,
debug_level = 0,
/

&io
input_root_name = '/home/clima/Build_WRF/WRF-4.5-ARW/test/em_real/wrfout_d01_2023-06-27_00:00:00'
output_root_name = '/home/clima/Build_WRF/ARWpost/wrfd1_2023062700'
plot = 'all_list'
fields = 'height, geopt, theta, tc, td2, td, rh, rh2, umet, vmet,
slp, pressure, u10m, v10m, wdir, wspd, wd10, ws10, cape, mcin,
lcl, lfc, cape, cin, clfr'
output_type = 'grads'
mercator_defs = .true.
/

split_output = .true.
frames_per_outfile = 2

plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10, slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'

&interp
interp_method = 1,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.

interp_method = 1, ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,
Click to expand...
found any solution?
 
hello all
i too have the same problem
START PROCESSING DATA

WARNING --- I do not recognize this data.

Will make an attempt to read it.

At line 152 of file process_domain_module.f (unit = 11)
Fortran runtime error: RECL parameter is non-positive in OPEN statement

Error termination. Backtrace:
#0 0x7f98799102ed in ???
#1 0x7f9879910ed5 in ???
#2 0x7f987991169d in ???
#3 0x7f9879a83f53 in ???
#4 0x7f9879a84214 in ???
#5 0x5628fb60e126 in ???
#6 0x5628fb5ff5d7 in ???
#7 0x7f9879187c86 in ???
#8 0x5628fb5fdfc9 in ???
#9 0xffffffffffffffff in ???
can you please help me to solve this?
 
/Models/WRF-Chem/ARW$ ./ARWpost.exe

!!!!!!!!!!!!!!!!
ARWpost v3.1
!!!!!!!!!!!!!!!!

FOUND the following input files:
./wrfout_d01_2015-09-07_00:00:00

START PROCESSING DATA

WARNING --- I do not recognize this data.

Will make an attempt to read it.

At line 150 of file process_domain_module.f (unit = 11)
Fortran runtime error: RECL parameter is non-positive in OPEN statement

Error termination. Backtrace:
#0 0x7ffa98823960 in ???
#1 0x7ffa988244d9 in ???
#2 0x7ffa9882510f in ???
#3 0x7ffa98a73a34 in ???
#4 0x7ffa98a741dd in ???
#5 0x55f11aa9e4af in ???
#6 0x55f11aa8ef7a in ???
#7 0x7ffa98429d8f in __libc_start_call_main
at ../sysdeps/nptl/libc_start_call_main.h:58
#8 0x7ffa98429e3f in __libc_start_main_impl
at ../csu/libc-start.c:392
#9 0x55f11aa8d5ed in ???
#10 0xffffffffffffffff in ???


same problem, maybe someone can help?
 
This error is caused by update of WRF codes. To make ARWPost work, please modify the code input_module.f90, in which you will find the following piece of code:
istatus = NF_GET_ATT_TEXT(ncid, nf_global, 'TITLE', title)


IF ( INDEX(title,'OUTPUT FROM GEOGRID') /= 0 ) iprogram = 1 !! geogrid output


IF ( INDEX(title,'OUTPUT FROM GRIDGEN') /= 0 ) iprogram = 1 !! old geogrid output


IF ( INDEX(title,'OUTPUT FROM METGRID') /= 0 ) iprogram = 3 !! metgrid output


IF ( INDEX(title,'OUTPUT FROM OBSGRID') /= 0 ) iprogram = 3 !! obsgrid output


IF ( INDEX(title,'OUTPUT FROM REAL_EM') /= 0 ) iprogram = 6 !! real.exe output


IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 ) iprogram = 8 !! wrf.exe output


IF ( iprogram == 0 ) THEN


print*," "


print*," WARNING --- I do not recognize this data."


print*," ",trim(title)


print*," Will make an attempt to read it."


print*," "


END IF

Please manually set the the value of iprogram and comment out the original code. Then recompile ARWPOst and ty again.
Note that ARWpost is an obsolete package and we recommend you to switch to python or NCL.
 
Ming Chen said:
This error is caused by update of WRF codes. To make ARWPost work, please modify the code input_module.f90, in which you will find the following piece of code:
istatus = NF_GET_ATT_TEXT(ncid, nf_global, 'TITLE', title)


IF ( INDEX(title,'OUTPUT FROM GEOGRID') /= 0 ) iprogram = 1 !! geogrid output


IF ( INDEX(title,'OUTPUT FROM GRIDGEN') /= 0 ) iprogram = 1 !! old geogrid output


IF ( INDEX(title,'OUTPUT FROM METGRID') /= 0 ) iprogram = 3 !! metgrid output


IF ( INDEX(title,'OUTPUT FROM OBSGRID') /= 0 ) iprogram = 3 !! obsgrid output


IF ( INDEX(title,'OUTPUT FROM REAL_EM') /= 0 ) iprogram = 6 !! real.exe output


IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 ) iprogram = 8 !! wrf.exe output


IF ( iprogram == 0 ) THEN


print*," "


print*," WARNING --- I do not recognize this data."


print*," ",trim(title)


print*," Will make an attempt to read it."


print*," "


END IF

Please manually set the the value of iprogram and comment out the original code. Then recompile ARWPOst and ty again.
Note that ARWpost is an obsolete package and we recommend you to switch to python or NCL.
Click to expand...
thank you Ming Chen, its my request you that, can you please provide any GitHub repositories for post processing?
 
tharundola said:
thank you Ming Chen, its my request you that, can you please provide any GitHub repositories for post processing?
Click to expand...
post processing tools

github.com

GitHub - NCAR/wrf-python: A collection of diagnostic and interpolation routines for use with output from the Weather Research and Forecasting (WRF-ARW) Model.

A collection of diagnostic and interpolation routines for use with output from the Weather Research and Forecasting (WRF-ARW) Model. - NCAR/wrf-python
github.com github.com
 
Dear Minch Chen,
i have two doubts to ask you
1. you have suggested the below modifications in input_module.f90. i have checked the input_module.f90 file, the script is similar to the code what you have given, i didn't properly understood what you suggested us to change manually. could you please explain us properly.
2. is it possible to update NetCDF version in WRF? if so can you tell me?

istatus = NF_GET_ATT_TEXT(ncid, nf_global, 'TITLE', title)


IF ( INDEX(title,'OUTPUT FROM GEOGRID') /= 0 ) iprogram = 1 !! geogrid output


IF ( INDEX(title,'OUTPUT FROM GRIDGEN') /= 0 ) iprogram = 1 !! old geogrid output


IF ( INDEX(title,'OUTPUT FROM METGRID') /= 0 ) iprogram = 3 !! metgrid output


IF ( INDEX(title,'OUTPUT FROM OBSGRID') /= 0 ) iprogram = 3 !! obsgrid output


IF ( INDEX(title,'OUTPUT FROM REAL_EM') /= 0 ) iprogram = 6 !! real.exe output


IF ( INDEX(title,'OUTPUT FROM WRF') /= 0 ) iprogram = 8 !! wrf.exe output


IF ( iprogram == 0 ) THEN


print*," "


print*," WARNING --- I do not recognize this data."


print*," ",trim(title)


print*," Will make an attempt to read it."


print*," "


END IF

Please manually set the the value of iprogram and comment out the original code.
 
You must log in or register to reply here.
Top