Water field capacity (watfc) as output

[NiCosentino]

Hello,
"watfc" appears in the history field list for CLM5, however I get error "in histFileMod.F90, htapes_fieldlist, watfc in fincl(1) for history tape 1 not found"

How can I output "watfc"?

I am using compset:
1850_CAM60_CLM50%SP_CICE%PRES_DOCN%DOM_RTM_SGLC_SWAV

thanks in advance

 

[oleson]

I see this in SoilStateType.F90:

if (use_cn) then
this%watfc_col(begc:endc,:) = spval
call hist_addfld2d (fname='watfc', units='m^3/m^3', type2d='levgrnd', &
avgflag='A', long_name='water field capacity', &
ptr_col=this%watfc_col, default='inactive')
end if

So it's only available by default if you are running BGC. Your compset is SP mode.
You could try removing the use_cn if statement. The calculation of watfc is done in SoilStateInitTimeConstMod.F90 and I think that calculation is done regardless of BGC or SP.

 

[NiCosentino]

Thanks! This worked perfectly.