NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Topics specifically related to the wrf.exe program
Post Reply
ashreeva
Posts: 5
Joined: Fri May 07, 2021 5:18 am

NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by ashreeva » Tue Jun 15, 2021 11:31 pm

Hello,

I am trying to simulate a heatwave over California using WRF (NLCD Landuse). In order to make it work with NLCD, I added the line: geog_data_res = 'nlcd2011_30m+5m', 'nlcd2011_30m+2m', 'nlcd2011_30m+urbfrac_nlcd2011+30s', to my namelist.wps. I then did the following remapping procedure to make sure it works with the Urban Canopy Model (UCM, sf_urban_physics=1):
31 (Dwarf scrub) --> 7 (Open shrublands)
32 (Shrubs) --> 7 (Open shrublands)
33 (Grasslands) --> 10 (Grasslands)
23 (Developed open space) --> 31
24 (Low intensity) --> 31
25 (Mid intensity) --> 32
26 (High intensity) --> 33
With these changes, I was able to run real.exe successfully, but now my wrf.exe is crashing with the following message:

"INPUT LandUse = "NLCD40"
LANDUSE TYPE = "NLCD40" FOUND 40 CATEGORIES 2 SEASONS WATER CATEGORY = 17 SNOW CATEGORY = 15
Using NLCD40 for Noah, redefine urban categories
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
Skipping over LUTYPE = MODIFIED_IGBP_MODIS_NOAH
LANDUSE TYPE = NLCD40 FOUND 40 CATEGORIES
At line 1981 of file module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')
Fortran runtime error: Bad integer for item 1 in list input "

I have attached my rsl.error and namelist.input files to this post. Please let me know how can I proceed and what else can I share to make the debugging easier?

Thanks, Anamika
Attachments
Archive.zip
(14.89 KiB) Downloaded 3 times

Ming Chen
Posts: 1588
Joined: Mon Apr 23, 2018 9:42 pm

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by Ming Chen » Mon Jun 21, 2021 8:53 pm

Anamika,
Would you please run this case using WRFv4.2 and let me know whether it works?
I am trying to find a solution for this issue and hope to get back to you soon. Thanks for your patience.
WRF Help Desk

Ming Chen
Posts: 1588
Joined: Mon Apr 23, 2018 9:42 pm

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by Ming Chen » Thu Jun 24, 2021 2:00 pm

Anamika,
Would you please use the VEGPRM.TBL I attached here to replace the one in WRFV4.3, then rerun your case? Please let me know whether your case can work. Thanks.
Attachments
VEGPARM.TBL
(22.5 KiB) Downloaded 6 times
WRF Help Desk

ashreeva
Posts: 5
Joined: Fri May 07, 2021 5:18 am

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by ashreeva » Thu Jun 24, 2021 6:42 pm

Thank you, Ming Chen, for your response.
I executed my WRF run with the attached VEGPARM.TBL. With this, I was able to get past my initial error and I got a new error this time.

This is the last line of my rsl.error.0000 file:
-------------- FATAL CALLED ---------------
USING 10 WUDAPT LCZ WITHOUT URBPARM_LCZ.TBL. SET USE_WUDAPT_LCZ=1
------------------------------------------------
rsl files.zip
(2.68 KiB) Downloaded 3 times
I have the URBPARM_LCZ.TBL in the same WRF/run/ folder though. So I am not sure what the error is.I am using NLCD urban classification system so could you please tell me if it is compatible with the LCZ urban morphology parameters (URBPARM_LCZ.TBL)? Also, I had remapped the urban classes as discussed in my original post. Could you please tell me if that is correct?

Ming Chen
Posts: 1588
Joined: Mon Apr 23, 2018 9:42 pm

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by Ming Chen » Fri Jun 25, 2021 3:10 am

Please look at the code phys/module_physics_init.F, find the following piece of code:

Code: Select all

               IF (config_flags%use_wudapt_lcz.eq.0 .and. max_utype_urb2d.gt.3.0) THEN  !new LCZ
                 CALL wrf_error_fatal &
                 ('USING 10 WUDAPT LCZ WITHOUT URBPARM_LCZ.TBL. SET USE_WUDAPT_LCZ=1')
               ENDIF
You can comment out the above lines, then recompile WRF and run your case.
Let me know whether it works.
WRF Help Desk

ashreeva
Posts: 5
Joined: Fri May 07, 2021 5:18 am

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by ashreeva » Fri Jul 02, 2021 4:49 am

Hello, Ming Chen,

I believe that I am past the previous error that we discussed here. But I am now getting a Segmentation error.


Here's the rsl.error.0000 file:

USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
ThompMP: read qr_acr_qgV3.dat instead of computing
ThompMP: read qr_acr_qsV2.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
mediation_integrate.G 1944 DATASET=HISTORY
mediation_integrate.G 1945 grid%id 1 grid%oid 1
Timing for Writing /scratch_lg/uhi-modeling/WRF_Home/OUTPUTS/LA2020/wrfout_d01_2020-08-01_00:00:00 for domain 1: 0.65176 elapsed seconds
d01 2020-08-01_00:00:00 Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain 1: 0.03110 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 234 JS 1 JE 180
WRF NUMBER OF TILES = 1

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0 0x2aaaab11021d in ???
#1 0x2aaaab10f463 in ???
#2 0x2aaaabcf03ff in ???
#3 0x2aaaabd3af3b in ???

Could it be a time_step issue? I am attaching my namelist.input file here. Please advice.

Anamika
Attachments
namelist.input
(3.91 KiB) Downloaded 3 times

ashreeva
Posts: 5
Joined: Fri May 07, 2021 5:18 am

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by ashreeva » Fri Jul 02, 2021 4:49 am

Hello, Ming Chen,

I believe that I am past the previous error that we discussed here. But I am now getting a Segmentation error.


Here's the rsl.error.0000 file:

USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
USING DEFAULT URBAN MORPHOLOGY
ThompMP: read qr_acr_qgV3.dat instead of computing
ThompMP: read qr_acr_qsV2.dat instead of computing
ThompMP: read freezeH2O.dat instead of computing
mediation_integrate.G 1944 DATASET=HISTORY
mediation_integrate.G 1945 grid%id 1 grid%oid 1
Timing for Writing /scratch_lg/uhi-modeling/WRF_Home/OUTPUTS/LA2020/wrfout_d01_2020-08-01_00:00:00 for domain 1: 0.65176 elapsed seconds
d01 2020-08-01_00:00:00 Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain 1: 0.03110 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 234 JS 1 JE 180
WRF NUMBER OF TILES = 1

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0 0x2aaaab11021d in ???
#1 0x2aaaab10f463 in ???
#2 0x2aaaabcf03ff in ???
#3 0x2aaaabd3af3b in ???

Could it be a time_step issue? I am attaching my namelist.input file here. Please advice.

Anamika

Ming Chen
Posts: 1588
Joined: Mon Apr 23, 2018 9:42 pm

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by Ming Chen » Fri Jul 02, 2021 9:03 pm

Hi, Anamika,
I looked at your namelist.input and everything looks just fine.
In WRFV4.3, a new option (LCZ) is introduced and some codes are modified. Despite the tests with Modis and USGS as input, I guess there is no test is conducted using NLCD as input. This is why we have troubles in your case.
With the changes in VEGPARM.TBl and phys/module_physics_init.F, I suppose the first two issues you reported have been resolved. Please let me know if I am wrong.
I also would like to know whether the same case works when using WRRFV4.2? We need to confirm that this problem is specific for WRFV4.3, which will give us some hints what is wrong.
WRF Help Desk

ashreeva
Posts: 5
Joined: Fri May 07, 2021 5:18 am

Re: NLCD40 WRF run crashes because of module_sf_noahdrv.f90 (unit = 19, file = 'VEGPARM.TBL')

Post by ashreeva » Thu Jul 15, 2021 5:58 pm

UPDATE: I found a solution to the SEGMENTATION FAULT issue.

I noticed that my model ran just fine if I only un the largest domain, but only when I add the smallest sub-domain does it break down. I then found the following resource: viewtopic.php?t=8508 which suggests that it might be a wrong choice of processors. I then followed the guidelines provided here on the post "How many processors should I use to run WRF?": viewtopic.php?f=73&t=5082 which helped me figure out that a square number of processors is optimal for me.

Running the job on 36 processors on 4 nodes seems to be working just fine.

All this to say that the records kept here at the WRF forum are pretty cool and I hope that this thread helps someone too. :)

Cheers,
Anamika

Post Reply

Return to “wrf.exe”