Update defaults for ne16 grids.

[mvdebolskiy]

Goes together with NorESMhub/fates#51.

• Make LUH files for fates.
• Add tests for new LUH files
• Check if LUH does not cause NorESM3 emission-driven restart issue NorESM#770

[mvdebolskiy]

@rosiealice @JessicaNeedham are we still in need for 16pft datasets for ne16? Or now running with 78 ones is fine, do you know by any chance?

[JessicaNeedham]

@mvdebolskiy , I'm actually not sure. I will try to figure this out. It is possible that with land use on the pft distribution in the surface datafile is not used but we need to confirm.

[JessicaNeedham]

@mvdebolskiy my read of the code is that the FATES pft distributions in fixed biogeography mode with land use are coming entirely from the fates_landuse_pft_map and not from the surface datafile.

GetLandusePFTData reads in the landuse_pft_file and fills in a landuse_pft_map:

CTSM/src/utils/clmfates_interfaceMod.F90

Line 4135 in 5acb791

real(r8), allocatable, intent(inout) :: landuse_pft_map(:,:,:)

which is then transferred to pft_areafrac_lu :

CTSM/src/utils/clmfates_interfaceMod.F90

Line 1007 in 5acb791

this%fates(nc)%bc_in(s)%pft_areafrac_lu(:,1:num_landuse_pft_vars) = landuse_pft_map(g,:,1:num_landuse_pft_vars)

and this is used in EDInitMod to initialise sites(s)%area_pft

https://github.com/NorESMhub/fates/blob/133dd3379e25dfa59421d5bbf63416a01d836761/main/EDInitMod.F90#L571

If fixed_biogeography is true but land use is not on, I think the surface data file might be used - see here:

CTSM/src/utils/clmfates_interfaceMod.F90

Line 1015 in 5acb791

this%fates(nc)%bc_in(s)%pft_areafrac(ft)=wt_nat_patch(g,m)

For CMIP and other MIPs we will only be running with land use so I don't think it is a high priority to have surface datafiles for non-landuse runs. But I'm not aware of everything that people have planned.

@rosiealice and @kjetilaas - does this match your understanding?

[JessicaNeedham]

I ran two ne16 runs with either the 78 or 16 HLM pft surface datafile. Here is FATES_NOCOMP_PATCHAREA_PF in the 78 run - 16 run.

They look identical.

[mvdebolskiy]

Ok, great. I am closing this then.