Floating-Point Exception Error

[mehmoodmz]

Hi,
I am getting a floating-point exception error while running CROPSIM-Wheat. What would be a better way to resolve this issue? Do I need to remove RSWT on Line 1022 or add RSWT on Line 1023?

dssat-csm-os/Plant/CROPSIM/CRP_Integrate.for

Line 1021 in cb4043a

! Originally included retained dead matter and reserves

dssat-csm-os/Plant/CROPSIM/CRP_Integrate.for

Line 1022 in cb4043a

IF (((LFWT+STWT+RSWT)*PLTPOP*10.0).GT.0.0)

dssat-csm-os/Plant/CROPSIM/CRP_Integrate.for

Line 1023 in cb4043a

& VANC = VNAD/((LFWT+STWT)*PLTPOP*10.0)

[palderman]

It seems like these comments indicate a desire to remove reserves from calculation of VNPCM:

dssat-csm-os/Plant/CROPSIM/CSCRP.for

Lines 102 to 106 in cb4043a

	! Changes for KSAS 02/12/2013
	! 1. Established flag set when LAI>0.0,not when shoot/root>2
	! (But ? this as allows re-use in the fall)
	! 2. May need to take reserve CH2O out of calculation of VNPCM

VNPCM is itself calculated based on VANC:

dssat-csm-os/Plant/CROPSIM/CRP_Integrate.for

Line 1556 in cb4043a

VNPCM = VANC*100.0

So presumably the calculation on line 1023 is as intended and it is the if statement on the preceding line that should be corrected:

dssat-csm-os/Plant/CROPSIM/CRP_Integrate.for

Lines 1021 to 1023 in cb4043a

	! Originally included retained dead matter and reserves
	IF (((LFWT+STWT+RSWT)*PLTPOP*10.0).GT.0.0)
	& VANC = VNAD/((LFWT+STWT)*PLTPOP*10.0)

If so, then line 1022 would become:

            IF (((LFWT+STWT)*PLTPOP*10.0).GT.0.0)

@lpmorenoc , @GerritHoogenboom, @chporter , what do y'all think?

[chporter]

@palderman Nice detective work! Yes, I agree that this should solve the problem.