Problem in CLM4.0 with irrigation part using CESM1.2.0

[xingyhuang@gmail_com]

Hi there, I have tried to run CESM using CLM4.0 with the irrigation model on. I have created a surface dataset containing irrigation data and also added  "-irrig=.true." to CLM_BLDNML_OPTS followned by the user guide. It seems that the irrigaition model has been turned on.However, I have checked the first three years' results and the QIRRIG (i.e. water added through irrigation) is all zero. I know that irrigation part only works for C3 generic crop, which  corresponds to pft type 16. So, why the irrigation model didnot work? Have I missed something here? How can I add irrigation water over cropland then? Thanks,Huang

 

[rosh_me91@gmail_com]

Hi  Huang,Im getting the same problem for clm4. So I tried working with clm 4.5 with crop also turned on. I still cant get a result for QIRRIG. Please let me know if you come to know the problem.  Thank you

 

[rosh_me91@gmail_com]

Hi  Huang,Im getting the same problem for clm4. So I tried working with clm 4.5 with crop also turned on. I still cant get a result for QIRRIG. Please let me know if you come to know the problem.  Thank you

 

[xingyhuang@gmail_com]

Not yet. However, if turning on the crop model, we also need to turn on BGC. That's really not what I want.Anyway, I am still expecting to get responses about this issue.Thanks for sharing the info.

 

[xingyhuang@gmail_com]

Not yet. However, if turning on the crop model, we also need to turn on BGC. That's really not what I want.Anyway, I am still expecting to get responses about this issue.Thanks for sharing the info.

 

[oleson]

After talking with Xingying about this at the CESM tutorial, I agree that this would seem to be a bug in CLM4.0 for
creating surface datasets related to irrigation.
The code in mklaiMod.F90 reads:

! if irrigation dataset present, copy LAI,SAI,Heights from non-irrigated
! into irrigated
if (firrig /= ' ') then
write(6,*) 'Irrigation dataset present; Copying crop ', &
' LAI, SAI, and heights into irrigated crop '
mlai_o(:,nonIrrigIdx) = mlai_o(:,IrrigIdx)
msai_o(:,nonIrrigIdx) = msai_o(:,IrrigIdx)
mhgtt_o(:,nonIrrigIdx) = mhgtt_o(:,IrrigIdx)
mhgtb_o(:,nonIrrigIdx) = mhgtb_o(:,IrrigIdx)
endif

Contrary to the comment, LAI, SAI, Heights are being copied from irrigated into non-irrigated. So the
recommended fix is, e.g.,

mlai_o(:,IrrigIdx) = mlai_o(:,nonIrrigIdx)

This seems to be OK for CLM4.5

I filed a CLM4.0 bug report (#2200)

Keith 

 

[oleson]

After talking with Xingying about this at the CESM tutorial, I agree that this would seem to be a bug in CLM4.0 for
creating surface datasets related to irrigation.
The code in mklaiMod.F90 reads:

! if irrigation dataset present, copy LAI,SAI,Heights from non-irrigated
! into irrigated
if (firrig /= ' ') then
write(6,*) 'Irrigation dataset present; Copying crop ', &
' LAI, SAI, and heights into irrigated crop '
mlai_o(:,nonIrrigIdx) = mlai_o(:,IrrigIdx)
msai_o(:,nonIrrigIdx) = msai_o(:,IrrigIdx)
mhgtt_o(:,nonIrrigIdx) = mhgtt_o(:,IrrigIdx)
mhgtb_o(:,nonIrrigIdx) = mhgtb_o(:,IrrigIdx)
endif

Contrary to the comment, LAI, SAI, Heights are being copied from irrigated into non-irrigated. So the
recommended fix is, e.g.,

mlai_o(:,IrrigIdx) = mlai_o(:,nonIrrigIdx)

This seems to be OK for CLM4.5

I filed a CLM4.0 bug report (#2200)

Keith 

 

[oleson]

Regarding irrigation in CLM4.5, there are two places where irrigation needs to be set to true:1. -irrig .true.  (in env_run.xml: CLM_BLDNML_OPTS)2. irrigate = .true. (in user_nl_clm) Keith

 

[oleson]

Regarding irrigation in CLM4.5, there are two places where irrigation needs to be set to true:1. -irrig .true.  (in env_run.xml: CLM_BLDNML_OPTS)2. irrigate = .true. (in user_nl_clm) Keith