Attribute not found while real.exe running

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

A

antonleontiev

New member
Hello! When I run a real.exe, the next error occurs at the rsl.error.0000:

tail rsl.error.0000
d01 2018-11-19_00:00:00 NetCDF error in ext_ncd_get_dom_ti.code INTEGER, line 83 Element JULYR
d01 2018-11-19_00:00:00 NetCDF error: NetCDF: Attribute not found
d01 2018-11-19_00:00:00 NetCDF error in ext_ncd_get_dom_ti.code INTEGER, line 83 Element JULDAY
d01 2018-11-19_00:00:00 mminlu = 'MODIFIED_IGBP_MODIS_NOAH'
NOTE: 1 namelist vs input data inconsistencies found.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1142
NOTE: Please check and reset these options
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

As I understand, there is something wrong with dates, but I don't have any idea, how to fix it. Is anyone know a solution? I use WRF3.9.
Thank you in advance!
 
Please post your namelist.wps and namelist.input for me to take a look.
It is also helpful to let me know which version of WRF and WPS you are running.
 
Thank you for your reply! I solved this issue, but the next issue rose. I successfully passed the WPS stage of the model, but after running the real.exe, I get an error message:

tail rsl.error.0000
d01 2015-01-15_00:00:00 Turning off use of TROPOPAUSE level data in vertical interpolation
Max map factor in domain 1 = 1.03. Scale the dt in the model accordingly.
d01 2015-01-15_00:00:00 checking boundary conditions for grid
d01 2015-01-15_00:00:00 boundary conditions OK for grid
ERROR in psfc: flag_psfc = 1, flag_soilhgt = 1, flag_slp = 0, sfcp_to_sfcp = F
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1367
not enough info for a p sfc computation
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

As I understand, the dt should be: dt = dx * 6 / max_map_factor, in my case it should be 30*6/1.03 = 174 or 175. By default the value of dt was 180. I changed it to 174, then to 175, I tried also another values but nothing was changed, I still get the same issue about the dt.
I use WRFV3.9, and my namelist.wps and namelist.input are in attachment.
Can you help me to solve this new issue?
 

Attachments

When running REAL,time_step doesn't really matter. Usually time_step= 6 x grid interval should work.
Please set sfcp_to_sfcp = .true. in &domain of your nameless.input, then try again.

If REAL still doesn't work, please tell me what input data you are using (e.g., GFS, ERA-I, etc). I will try to repeat this case.
 
Hello! Thank you for your prompt reply!
I added the sfcp_to_sfcp = .true, a new error without notification from my previous post, but with advise to changed e_vert from 28 to 35 and p_top_requsted from 5000 to 3000. I made these changes and at the end get almost the same error:

tail rsl.error.0000
d01 2015-01-15_00:00:00 Turning off use of TROPOPAUSE level data in vertical interpolation
Max map factor in domain 1 = 1.00. Scale the dt in the model accordingly.
Using sfcprs2 to compute psfc
d01 2015-01-15_00:00:00 flag value of flag_qnwfa is 0
d01 2015-01-15_00:00:00 No average surface temperature for use with inland lakes
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 2626
Not enough soil temperature data for Noah LSM scheme.
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

My namelists input and wps are in the attachment, I used GFS data. Do you have any suggestions? And could you please use my set of data while checking? May be I downloaded something wrong because I am a new user of these data and model. I downloaded data from rda.ucar.edu.
Thank you in advance!
 

Attachments

The error message indicates that there might have no soil information in your met_em files.
Please type ncdump -h met_em..d01.yyyy-mm-dd_hh:00:00.nc, and see whether you can find the variables below:
num_st_layers
num_sm_layers
SM
ST

I am suspicious that the above variables don't exist in your met_em files.
if so, please download GFS data from https://rda.ucar.edu/datasets/ds084.1/, and ungrib these data using newer version of WPS like WPSv4.0.

Let me know if you still have problems.
 
Hello!
Thank you for your prompt reply. It seems, that I have all data, you have mentioned. Here is an output of ncdump for variables:

num_st_layers = 4 ;
float SOIL_LAYERS(Time, num_st_layers, south_north, west_east) ;
float ST(Time, num_st_layers, south_north, west_east) ;

num_sm_layers = 4 ;
float SM(Time, num_sm_layers, south_north, west_east) ;

float ST(Time, num_st_layers, south_north, west_east) ;
ST:FieldType = 104 ;
ST:MemoryOrder = "XYZ" ;
ST:units = "" ;
ST:description = "" ;
ST:stagger = "M" ;
ST:sr_x = 1 ;
ST:sr_y = 1 ;
float ST100200(Time, south_north, west_east) ;
ST100200:FieldType = 104 ;
ST100200:MemoryOrder = "XY " ;
ST100200:units = "K" ;
ST100200:description = "T 100-200 cm below ground layer (Bottom)" ;
ST100200:stagger = "M" ;
ST100200:sr_x = 1 ;
ST100200:sr_y = 1 ;
float ST040100(Time, south_north, west_east) ;
ST040100:FieldType = 104 ;
ST040100:MemoryOrder = "XY " ;
ST040100:units = "K" ;
ST040100:description = "T 40-100 cm below ground layer (Upper)" ;
ST040100:stagger = "M" ;
ST040100:sr_x = 1 ;
ST040100:sr_y = 1 ;
float ST010040(Time, south_north, west_east) ;
ST010040:FieldType = 104 ;
ST010040:MemoryOrder = "XY " ;
ST010040:units = "K" ;
ST010040:description = "T 10-40 cm below ground layer (Upper)" ;
ST010040:stagger = "M" ;
ST010040:sr_x = 1 ;
ST010040:sr_y = 1 ;
float ST000010(Time, south_north, west_east) ;
ST000010:FieldType = 104 ;
ST000010:MemoryOrder = "XY " ;
ST000010:units = "K" ;
ST000010:description = "T 0-10 cm below ground layer (Upper)" ;
ST000010:stagger = "M" ;
ST000010:sr_x = 1 ;
ST000010:sr_y = 1 ;
:SIMULATION_START_DATE = "2018-11-10_00:00:00" ;
:WEST-EAST_GRID_DIMENSION = 40 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 39 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 40 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 68 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 69 ;
:STAND_LON = 35.f ;
:FLAG_ST000010 = 1 ;
:FLAG_ST010040 = 1 ;
:FLAG_ST040100 = 1 ;
:FLAG_ST100200 = 1 ;

float SM(Time, num_sm_layers, south_north, west_east) ;
SM:FieldType = 104 ;
SM:MemoryOrder = "XYZ" ;
SM:units = "" ;
SM:description = "" ;
SM:stagger = "M" ;
SM:sr_x = 1 ;
SM:sr_y = 1 ;
float SM100200(Time, south_north, west_east) ;
SM100200:FieldType = 104 ;
SM100200:MemoryOrder = "XY " ;
SM100200:units = "fraction" ;
SM100200:description = "Soil Moist 100-200 cm below gr layer" ;
SM100200:stagger = "M" ;
SM100200:sr_x = 1 ;
SM100200:sr_y = 1 ;
float SM040100(Time, south_north, west_east) ;
SM040100:FieldType = 104 ;
SM040100:MemoryOrder = "XY " ;
SM040100:units = "fraction" ;
SM040100:description = "Soil Moist 40-100 cm below grn layer" ;
SM040100:stagger = "M" ;
SM040100:sr_x = 1 ;
SM040100:sr_y = 1 ;
float SM010040(Time, south_north, west_east) ;
SM010040:FieldType = 104 ;
SM010040:MemoryOrder = "XY " ;
SM010040:units = "fraction" ;
SM010040:description = "Soil Moist 10-40 cm below grn layer" ;
SM010040:stagger = "M" ;
SM010040:sr_x = 1 ;
SM010040:sr_y = 1 ;
float SM000010(Time, south_north, west_east) ;
SM000010:FieldType = 104 ;
SM000010:MemoryOrder = "XY " ;
SM000010:units = "fraction" ;
SM000010:description = "Soil Moist 0-10 cm below grn layer (Up)" ;
SM000010:stagger = "M" ;
SM000010:sr_x = 1 ;
SM000010:sr_y = 1 ;
:FLAG_SM000010 = 1 ;
:FLAG_SM010040 = 1 ;
:FLAG_SM040100 = 1 ;
:FLAG_SM100200 = 1 ;

Anyway, I will try to use WPS4 due to your advice. I will get back to you shortly with the update.
Thank you!
 
Hello! I tried to change my version of WPS to WPS4, but there is almost the same problem. I downloaded data from link, that you mentioned, and in my opinion I did exactly what you told me to do. My new error is:

tail rsl.error.0000
d01 2015-01-15_00:00:00 Turning off use of TROPOPAUSE level data in vertical interpolation
Max map factor in domain 1 = 1.00. Scale the dt in the model accordingly.
Using sfcprs2 to compute psfc
d01 2015-01-15_00:00:00 flag value of flag_qnwfa is 0
d01 2015-01-15_00:00:00 No average surface temperature for use with inland lakes
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 2626
Not enough soil temperature data for Noah LSM scheme.
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

I really don't have any idea, why is it happens. May be you have another solution?
Thank you in advance!
 
Hello! Finally, I solved the problem. My data were corrupted, after downloading new data from here https://rda.ucar.edu/datasets/ds084.1/index.html#sfol-wl-/data/ds084.1 the problem has gone.
 
Hi, I had the same error message today switching to WPS V4.1
Code: 
ERROR in psfc: flag_psfc =  1, flag_soilhgt =  1, flag_slp =  0, sfcp_to_sfcp = F
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1404
not enough info for a p sfc computation

The solution was to switch from GFS downloaded variable PRMSL to MSLET as it is said in the VTable.GFS file.
 
You must log in or register to reply here.
Back
Top