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

WRFChem-3.9.1 Compilation Error

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.

Hi,

It would be extremely helpful if you could get the size of the seasonal_wes array to be printed (are you missing a closing parentheses?) or at least the other print statements that I suggested (i.e., the value of i, j, iland, and current month). Currently the modeel is failing on one of the index and therefore will not print out a value; i.e., you are now just printing the values associated with indexes that are valid. We need to know what the indexes are before it tries to access that array so we know which one may be out of bounds.

Jordan
 
Hi Jordan,

Thanks for your response.
Yes, I missed the closing parentheses.

I have put these lines before L212 in module_dep_simple.F.

write(*,*) 'i = ', i
write(*,*) 'j = ', j
write(*,*) 'id = ', id
write(*,*) 'iland = ' , iland
write(*,*) 'current_month = ', current_month
write(*,*) 'size of seasonal_wes = ', size(seasonal_pft(id)%seasonal_wes:),:,:,:))

Now the compilation is successful.

The rsl.out file gives following size of seasonal_wes :

d01 2009-01-01_00:01:30 calling dry_deposition_driver
d01 2009-01-01_00:01:30 DOING DRY DEP VELOCITIES WITH WESELY METHOD
d01 2009-01-01_00:01:30 in dry_dep_wesely
i = 1
j = 193
id = 1
iland = 16
current_month = 1
size of seasonal_wes = 534528
i = 2
j = 193
id = 1
iland = 16
current_month = 1
size of seasonal_wes = 534528
i = 3
j = 193
id = 1
iland = 16
current_month = 1
size of seasonal_wes = 534528
......................................................
.....................................................

You can check these values at the end of rsl.out.0072,rsl.out.0100 files.
I have attached the rsl.out.0072 and rsl.error.0100 files for your reference.

Kindly suggest me further.


Regards,

Saurabh kumar
 

Attachments

  • rsl.error.0100.txt
    3.5 MB · Views: 12
  • rsl.out.0072.txt
    3.7 MB · Views: 12
  • rsl.out.0100.txt
    3.7 MB · Views: 13
Hi Saurabhm,

Ok, your error is error is occurring because iland = 0.

iland comes from ivgtyp. Please check your wrfinput file, it likely has 0's where they shouldnt be for this field.

Jordan
 
Hi, Jordan

Thanks for your email.

I have checked IVGTYP in wrfinput_d01.
em_real]$ ncdump -v IVGTYP wrfinput_d01 > & ivgtyp.txt &

Its value ranges between 1 to 31. It seems IVGTYP has not zero value in wrfinput file.

I have attached the ivgtyp.txt file and ncview plot for the same.

Kindly help me.


Regards,

Saurabh Kumar
 

Attachments

  • IVGTYP.pdf
    515.2 KB · Views: 13
  • ivgtyp.txt
    1.5 MB · Views: 11
Hi Saurabh,

I think your namelist setting num_land_cat is inconsistent with your use of MODIFIED_IGBP_MODIS_NOAH land use types. I think setting num_land_cat = 20 or 21 depending on if you have included the lake category should fix this.

Jordan
 
Top