haiqingsong
Member
Can someone run WRF-V4.5 with the Thompson scheme (mp_physics=38)？
It will crash when i setup mp_physics=38.
It will crash when i setup mp_physics=38.
You are using an out of date browser. It may not display this or other websites correctly.
You should upgrade or use an alternative browser.
You should upgrade or use an alternative browser.
Can someone run WRF with the Thompson scheme (mp_physics=38)？
- Thread starter haiqingsong
- Start date
-
- Tags
- wrf physics
haiqingsong
Member
This case crashed immediately！Did this case crashed immediately? What error message did you get?
Please upload your namelist.input for me to take a look. Thanks.
When I setup mp_physics=8 or others， wrf runs well.
I use the default namelist.input from WRF.
haiqingsong
Member
The error：Did this case crashed immediately? What error message did you get?
Please upload your namelist.input for me to take a look. Thanks.
DEBUGT k, hgt = 51 20647.57
Apparently there are no initial CCN aerosols.
checked column at point (i,j) = 1 1
Apparently there are no initial IN aerosols.
checked column at point (i,j) = 1 1
Apparently there are no initial BC aerosols.
checked column at point (i,j) = 1 1
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
libc.so.6 0000146DF0372D90 Unknown Unknown Unknown
wrf.exe 00000000021DF35D Unknown Unknown Unknown
wrf.exe 00000000028CC112 Unknown Unknown Unknown
wrf.exe 00000000028B650E Unknown Unknown Unknown
wrf.exe 000000000205A6DE Unknown Unknown Unknown
wrf.exe 0000000001CF5983 Unknown Unknown Unknown
wrf.exe 000000000167EB94 Unknown Unknown Unknown
haiqingsong
Member
The WRF4.6.0 crashes immediately with mp_physics=38 ！Did this case crashed immediately? What error message did you get?
Please upload your namelist.input for me to take a look. Thanks.
Hi,
Could you upload your namelist.input and rsl.error.0000? What are pre-computed microphysics tables in your run directory (*.dat files)? You need a different microphysics table for mp=38, so maybe that's the issue, but without seeing your rsl.error file, I don't know whether the tables have been found or not.
You can download the table from Index of /wrf/src/wrf_files or set it to calculate at initialization by including write_thompson_mp38table = .true. on the namelist.
Could you upload your namelist.input and rsl.error.0000? What are pre-computed microphysics tables in your run directory (*.dat files)? You need a different microphysics table for mp=38, so maybe that's the issue, but without seeing your rsl.error file, I don't know whether the tables have been found or not.
You can download the table from Index of /wrf/src/wrf_files or set it to calculate at initialization by including write_thompson_mp38table = .true. on the namelist.
Last edited:
haiqingsong
Member
Thanks！ I have downloaded qr_acr_qg_mp38V1.dat to my test/em_realHi,
Could you upload your namelist.input and rsl.error.0000? What are pre-computed microphysics tables in your run directory (*.dat files)?
My max_dom=1.
The real.exe runs successfully.
The wrf.exe crashed.
Attachments
It looks like there's an issue when reading the tables... When I compare your rsl file with one I have here, the subsequent step should be to read the tables:
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND 20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
ThompMP: read qr_acr_qg_mp38V1.dat instead of computing
ThompMP: read qr_acr_qsV2.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
Could you try to add this option on your namelist under the &physics section:
write_thompson_mp38table = .true.,
and remove the file qr_acr_qg_mp38V1.dat from your run directory.
If it still crashes, could you pls upload your wrfinput_d01 and wrfbdy_d01 files?
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND 20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
ThompMP: read qr_acr_qg_mp38V1.dat instead of computing
ThompMP: read qr_acr_qsV2.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
Could you try to add this option on your namelist under the &physics section:
write_thompson_mp38table = .true.,
and remove the file qr_acr_qg_mp38V1.dat from your run directory.
If it still crashes, could you pls upload your wrfinput_d01 and wrfbdy_d01 files?
haiqingsong
Member
It does not work.It looks like there's an issue when reading the tables... When I compare your rsl file with one I have here, the subsequent step should be to read the tables:
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND 20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
ThompMP: read qr_acr_qg_mp38V1.dat instead of computing
ThompMP: read qr_acr_qsV2.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
Could you try to add this option on your namelist under the &physics section:
write_thompson_mp38table = .true.,
and remove the file qr_acr_qg_mp38V1.dat from your run directory.
If it still crashes, could you pls upload your wrfinput_d01 and wrfbdy_d01 files?
Attachments
haiqingsong
Member
The wrfinput_d01 and wrfbdy_d01 files are too large for the server to process. I will email you now！It looks like there's an issue when reading the tables... When I compare your rsl file with one I have here, the subsequent step should be to read the tables:
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND 20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES
ThompMP: read qr_acr_qg_mp38V1.dat instead of computing
ThompMP: read qr_acr_qsV2.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
Could you try to add this option on your namelist under the &physics section:
write_thompson_mp38table = .true.,
and remove the file qr_acr_qg_mp38V1.dat from your run directory.
If it still crashes, could you pls upload your wrfinput_d01 and wrfbdy_d01 files?
Could you tell me if you ran real using the same exact namelist you shared above?
I'm thinking the segmentation fault could be related to not having the aerosol profile in the inputs. When I check the variables in your input, the variables QNIFA and QNWFA which should contain the aerosol number concentration, are all empty. The model *shouldn't* crash because of it, but this can be the quickest solution to get your simulation to run. It's hard for me to pinpoint exactly where the model crashes without recompiling it to debug properly.
I was able to confirm that the option mp=28 runs without an issue whereas mp=38 crashes. This is a good start and it looks like there may be a bug in the code. However, investigating this more closely will take me a little time.
Maybe @Ming Chen or someone from wrfhelp could help me narrow it down? I have the user's input files on /glade/derecho/scratch/frediani/test_mp38/files_from_user
I'm thinking the segmentation fault could be related to not having the aerosol profile in the inputs. When I check the variables in your input, the variables QNIFA and QNWFA which should contain the aerosol number concentration, are all empty. The model *shouldn't* crash because of it, but this can be the quickest solution to get your simulation to run. It's hard for me to pinpoint exactly where the model crashes without recompiling it to debug properly.
I was able to confirm that the option mp=28 runs without an issue whereas mp=38 crashes. This is a good start and it looks like there may be a bug in the code. However, investigating this more closely will take me a little time.
Maybe @Ming Chen or someone from wrfhelp could help me narrow it down? I have the user's input files on /glade/derecho/scratch/frediani/test_mp38/files_from_user
haiqingsong
Member
I have downloaded https://www2.mmm.ucar.edu/wrf/src/wps_files/QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat.bz2Could you tell me if you ran real using the same exact namelist you shared above?
I'm thinking the segmentation fault could be related to not having the aerosol profile in the inputs. When I check the variables in your input, the variables QNIFA and QNWFA which should contain the aerosol number concentration, are all empty. The model *shouldn't* crash because of it, but this can be the quickest solution to get your simulation to run. It's hard for me to pinpoint exactly where the model crashes without recompiling it to debug properly.
I was able to confirm that the option mp=28 runs without an issue whereas mp=38 crashes. This is a good start and it looks like there may be a bug in the code. However, investigating this more closely will take me a little time.
Maybe @Ming Chen or someone from wrfhelp could help me narrow it down? I have the user's input files on /glade/derecho/scratch/frediani/test_mp38/files_from_user
And added constants_name = '/share/wrfinputdata/geog/QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat' in my namelist.wps
Should I delete this line in my namelist.wps?
haiqingsong
Member
Thank you so much！Could you tell me if you ran real using the same exact namelist you shared above?
I'm thinking the segmentation fault could be related to not having the aerosol profile in the inputs. When I check the variables in your input, the variables QNIFA and QNWFA which should contain the aerosol number concentration, are all empty. The model *shouldn't* crash because of it, but this can be the quickest solution to get your simulation to run. It's hard for me to pinpoint exactly where the model crashes without recompiling it to debug properly.
I was able to confirm that the option mp=28 runs without an issue whereas mp=38 crashes. This is a good start and it looks like there may be a bug in the code. However, investigating this more closely will take me a little time.
Maybe @Ming Chen or someone from wrfhelp could help me narrow it down? I have the user's input files on /glade/derecho/scratch/frediani/test_mp38/files_from_user
When I delete lines below in namelist.wps，it runs now！
aer_opt = 3,
use_aero_icbc =.true.,
So，the aerosol options affects MP=38 ？There may be a bug in the code.
Great! Thanks for letting me know!
Were you able to get the aerosol fields QNIFA and QNWFA in your wrfinput file once you made the change in WPS? If you have the fields with valid data, I'd expect you'd be able to run with
aer_opt = 3,
use_aero_icbc =.true.
Either way, it does look like there's a bug.
Were you able to get the aerosol fields QNIFA and QNWFA in your wrfinput file once you made the change in WPS? If you have the fields with valid data, I'd expect you'd be able to run with
aer_opt = 3,
use_aero_icbc =.true.
Either way, it does look like there's a bug.
haiqingsong
Member
Thanks!Great! Thanks for letting me know!
Were you able to get the aerosol fields QNIFA and QNWFA in your wrfinput file once you made the change in WPS? If you have the fields with valid data, I'd expect you'd be able to run with
aer_opt = 3,
use_aero_icbc =.true.
Either way, it does look like there's a bug.
I have downloaded https://www2.mmm.ucar.edu/wrf/src/wps_files/QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat.bz2
And added constants_name = '/share/wrfinputdata/geog/QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat' in my namelist.wps
The wrfinput_d01 emailed to you contains wrong aerosol data?
what is wrong in my namelist?
https://www2.mmm.ucar.edu/wrf/src/wps_files/QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat.bz2 covers Asia?
how to resolve this?
haiqingsong
Member
And running wrf with MP=38 needs aerosol data？Great! Thanks for letting me know!
Were you able to get the aerosol fields QNIFA and QNWFA in your wrfinput file once you made the change in WPS? If you have the fields with valid data, I'd expect you'd be able to run with
aer_opt = 3,
use_aero_icbc =.true.
Either way, it does look like there's a bug.
The file QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat contains a global aerosol climatology. You don't necessarily need the aerosol climatology to run the model with mp=38.
You can find more information about the aerosols and how the scheme uses them here: WRF Model Users Site
This description is also valid for the hail-aware version of the scheme (mp=38).
You can find more information about the aerosols and how the scheme uses them here: WRF Model Users Site
This description is also valid for the hail-aware version of the scheme (mp=38).
haiqingsong
Member
Update：The file QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat contains a global aerosol climatology. You don't necessarily need the aerosol climatology to run the model with mp=38.
You can find more information about the aerosols and how the scheme uses them here: WRF Model Users Site
This description is also valid for the hail-aware version of the scheme (mp=38).
We must delete the following 2 lines in namelist.input at the same time, otherwise MP38 will crash.
wif_input_opt = 2,
num_wif_levels = 30,