how to obtain the indices of natural PFTs of each grid of surface data?

[dayon]

Dear scientists,
I attempt to look for the index of PFT of each grid from surface data, but I only find the "natpft" with the long name "indices of natural PFTs". It is just a one-dimensional vector of size 15 instead of gridcell.

 

[oleson]

There should be a variable called "PCT_NAT_PFT" on the primary (h0) CLM history file:

float PCT_NAT_PFT(time, natpft, lat, lon) ;
PCT_NAT_PFT:long_name = "% of each PFT on the natural vegetation (i.e., soil) landunit" ;
PCT_NAT_PFT:units = "%" ;
PCT_NAT_PFT:cell_methods = "time: mean" ;
PCT_NAT_PFT:_FillValue = 1.e+36f ;
PCT_NAT_PFT:missing_value = 1.e+36f ;
PCT_NAT_PFT:landunit_mask = "unknown" ;

 

[dayon]

Oleson, Thanks for your reply, it is so helpful!

 

[adrienD]

There should be a variable called "PCT_NAT_PFT" on the primary (h0) CLM history file:

float PCT_NAT_PFT(time, natpft, lat, lon) ;
PCT_NAT_PFT:long_name = "% of each PFT on the natural vegetation (i.e., soil) landunit" ;
PCT_NAT_PFT:units = "%" ;
PCT_NAT_PFT:cell_methods = "time: mean" ;
PCT_NAT_PFT:_FillValue = 1.e+36f ;
PCT_NAT_PFT:missing_value = 1.e+36f ;
PCT_NAT_PFT:landunit_mask = "unknown" ;

Hi Keith,

How can I know the name of each of those 15 PFT indexes?

I can only find a list with 18 PFTs 2.2. Surface Characterization, Vertical Discretization, and Model Input Requirements — ctsm CTSM master documentation

 

[oleson]

Index 0 of the natpft dimension of PCT_NAT_PFT is bare soil.
Then indices 1-14 correspond to the plant functional types listed in order in Table 2.2.1, e.g., index 1 is Needleleaf evergreen tree – temperate
Those are on the natural vegetation landunit.

 

[strandwg]

Index 0 of the natpft dimension of PCT_NAT_PFT is bare soil.
Then indices 1-14 correspond to the plant functional types listed in order in Table 2.2.1, e.g., index 1 is Needleleaf evergreen tree – temperate
Those are on the natural vegetation landunit.

Keith, could that table have an additional column that specifies the PFT index for each row? That would be ideal. Thanks!

 

[oleson]

Yeah, it gets a bit complicated when crops are included, but it's a good idea, I'll file an issue on it, thanks.

 

[Yuan Sun]

Hi Keith,

How can I know the name of each of those 15 PFT indexes?

I can only find a list with 18 PFTs 2.2. Surface Characterization, Vertical Discretization, and Model Input Requirements — ctsm CTSM master documentation

Hi Andrien,

Do you know how to index land use data with 15 PFT now?

Thanks, Yuan