mauro_morichetti
New member
Hi all,
I am facing a problem when I run real.exe using the chemistry (chem_opt= 112), the code stops after looking for "Subroutine input_chem_profile" with a segmentation fault, while if I do not use chemistry (chem_opt=0) or other chemistry (i.e. chem_opt =300) not related to MOZART gas phase, I have not problem.
I have read something about, and it seems to be a memory limit problem, but I have already set "limit stacksize unlimited' on my enviroment and, I used the code successfully before (with the same namelist and domain), so I can't understand what is happening.
Someone adviced me to set "chem_in_opt =1", but the model stops since I do not have the "wrf_chem_input_d01" (actually I have nevere used before).
WRF-Chem is compiled with the traceback option. If you open the rsl.error.0009 attached you can see the backtrace error.
I opened the module_input_chem_data, and I noticed that he problem is that the "is" indice is equal 0, since "iref" is 0
in the module_input_chem_data.F:
I looked trough the code and I found that "iref" is defined and associated to a precise value, for different subroutines (i.e., setup_gasprofile_map_saprcnov, setup_gasprofile_map_radm_racm...). For the "SUBROUTINE make_chem_profile" it is defined, but it has not a value associated. So I copied the lines 411-413 (the iref for SUBROUTINE setup_gasprofile_map_radm_racm)
from:
to:
Now the real.exe (with chemistry) is running, but I would like to understand why it is happing. As I said before I have already used the same simulation (setting, physics, chemistry, domain) without problem, and I have read different people with the same problem, so I believe could be a bug or something.
Any Ideas?
Thanks
I am facing a problem when I run real.exe using the chemistry (chem_opt= 112), the code stops after looking for "Subroutine input_chem_profile" with a segmentation fault, while if I do not use chemistry (chem_opt=0) or other chemistry (i.e. chem_opt =300) not related to MOZART gas phase, I have not problem.
I have read something about, and it seems to be a memory limit problem, but I have already set "limit stacksize unlimited' on my enviroment and, I used the code successfully before (with the same namelist and domain), so I can't understand what is happening.
Someone adviced me to set "chem_in_opt =1", but the model stops since I do not have the "wrf_chem_input_d01" (actually I have nevere used before).
WRF-Chem is compiled with the traceback option. If you open the rsl.error.0009 attached you can see the backtrace error.
I opened the module_input_chem_data, and I noticed that he problem is that the "is" indice is equal 0, since "iref" is 0
in the module_input_chem_data.F:
Code:
847 !UMM = Update Mauro Morichetti
848 !print*,'iref',iref
849 DO l=2, numgas
850 is=iref(l-1)
851 DO j=ny1,ny2
852 DO k=nz1,nz2
853 DO i=nx1,nx2
854 chem(i,k,j,l)=fracref(l-1)*stor(i,k,j,is) * 1.E6
855 ENDDO
856 ENDDO
857 ENDDO
858 ENDDO
I looked trough the code and I found that "iref" is defined and associated to a precise value, for different subroutines (i.e., setup_gasprofile_map_saprcnov, setup_gasprofile_map_radm_racm...). For the "SUBROUTINE make_chem_profile" it is defined, but it has not a value associated. So I copied the lines 411-413 (the iref for SUBROUTINE setup_gasprofile_map_radm_racm)
from:
Code:
411 iref(:) = 7 !default value
412 iref(1:41) = (/12,19,2,2,1,3,4,9,8,5,5,32,6,6,6,30,30,10,26,13,11,6,6, &
413 14,15,15,23,23,32,16,23,31,17,23,23,23,23,23,7,28,29/)
to:
Code:
797 !UMM = Update Mauro Morichetti
798 iref(:) = 7 !default value
799 iref(1:41) = (/12,19,2,2,1,3,4,9,8,5,5,32,6,6,6,30,30,10,26,13,11,6,6, &
800 14,15,15,23,23,32,16,23,31,17,23,23,23,23,23,7,28,29/)
Now the real.exe (with chemistry) is running, but I would like to understand why it is happing. As I said before I have already used the same simulation (setting, physics, chemistry, domain) without problem, and I have read different people with the same problem, so I believe could be a bug or something.
Any Ideas?
Thanks