How can I get LAI of each PFT using CLM5

[Sean Chen]

fsurdat='/scratch01/csy/surfdata/surfdata_360x720cru_16pfts_Irrig_CMIP6_simyr2000_c170824.nc'(my new surfdata)
hist_avgflag_pertape = 'A'
hist_nhtfrq = -192
hist_mfilt = 1
hist_empty_htapes = .true
hist_fincl1='TLAI'
The 'TLAI' is a 2D array for the global LAI. But I need to output the LAI of each PFT, what should I do with the user_nl_clm file?
I'm currently using CESM2.2, CLM5.0.

Any suggestions will be appreciated!
Thanks!
Sean

 

[slevis]

@Sean Chen
Please see namelist instructions here:

1.2.3. Customizing CLM’s namelist — ctsm release-clm5.0 documentation

and even more specifically this section:

1.2.3. Customizing CLM’s namelist — ctsm release-clm5.0 documentation

 

[Sean Chen]

@Sean Chen
Please see namelist instructions here:

1.2.3. Customizing CLM’s namelist — ctsm release-clm5.0 documentation

and even more specifically this section:

1.2.3. Customizing CLM’s namelist — ctsm release-clm5.0 documentation

Hi Slevis,
Thanks a lot! I run the CLM5 according to the tech documents from your linkages. But I only got the 2D (pft, time) LAI outputs. But we need 3D (pft, lon, lat) outputs. Could you please help me with this problem?
This is the setup in the user_nl_clm file:
hist_fincl2 = 'TLAI'
hist_fincl3 = 'TLAI'
hist_fincl4 = 'TLAI'
hist_dov2xy = .true., .false., .false., .false.
hist_type1d_pertape = ' ', 'GRID', 'COLS', 'PFTS'
hist_nhtfrq = 0, -24, -24, -24
hist_mfilt = 12, 30, 28,24

 

[slevis]

The model cannot generate (pft, lon, lat) output. It can generate (lon, lat) output at the grid-cell level OR pft or col or land-unit output as vectors. To convert the vector output to (lon, lat), use the various index variables provided in the history files.

 

[slevis]

Additional information that you may find helpful is here:

Output the clm results at the pft level

Dear official staff, I want to get the result on the PFT level, and I set the following parameters: hist_fincl2 = 'TG' hist_fincl3 = 'TG' hist_fincl4 = 'TG' hist_dov2xy = .true., .false., .false., .false. hist_type1d_pertape = ' ', 'GRID', 'COLS', 'PFTS' hist_nhtfrq = 0, -24, -24, -24...

bb.cgd.ucar.edu