CLM 4.0 Runtime error, cannot compute rain and snow

[umair@skku_edu]

Hy everyone,


I got problem while running CLM using own forcing data. The error in ccsm log file shows as:



(shr_sys_abort) ERROR: (datm_comp_run) ERROR: cannot compute rain and snow

(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping

(shr_sys_abort) ERROR: (datm_comp_run) ERROR: cannot compute rain and snow

(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping

(shr_sys_abort) ERROR: (datm_comp_run) ERROR: cannot compute rain and snow

(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping



It originates from datm_comp_mod.F90 file, line 1152.

The code in this file at that point which is causing problem is:

!--- rain and snow ---

if (sprecc > 0 .and. sprecl > 0) then

a2x%rAttr(krc,n) = avstrm%rAttr(sprecc,n)

a2x%rAttr(krl,n) = avstrm%rAttr(sprecl,n)

elseif (sprecn > 0) then

a2x%rAttr(krc,n) = avstrm%rAttr(sprecn,n)*0.1_R8

a2x%rAttr(krl,n) = avstrm%rAttr(sprecn,n)*0.9_R8

else

call shr_sys_abort(subName//'ERROR: cannot compute rain and snow')

endif


!--- split precip between rain & snow ---

!--- note: aribitrarily small negative values cause CLM to crash ---

frac = (tbot - tkFrz)*0.5_R8 ! ramp near freezing

frac = min(1.0_R8,max(0.0_R8,frac)) ! bound in [0,1]

a2x%rAttr(ksc,n) = max(0.0_R8, a2x%rAttr(krc,n)*(1.0_R8 - frac) )

a2x%rAttr(ksl,n) = max(0.0_R8, a2x%rAttr(krl,n)*(1.0_R8 - frac) )

a2x%rAttr(krc,n) = max(0.0_R8, a2x%rAttr(krc,n)*( frac) )

a2x%rAttr(krl,n) = max(0.0_R8, a2x%rAttr(krl,n)*( frac) )



enddo



end select





I think there must be some problem in creating input data for precipitation. But i really don't understand the exact reason.


Anyone please, if have any idea behind what could be the cause of this error.





Regards

Umi