New Soil Color Field on CESM

[Jairo]

Hi,
I am conducting an experiment to examine the impact of artificially altering the albedo of the ground in a specific area on the CESM fully coupled. The soil cover is irrelevant on this experiment. However, I’m struggling to create a new SOIL_COLOR field.

Currently, the model uses 20 prescribed soil colors. I want to create a new soil color and modify its albedos. I know how to modify the surfdata input, but I still need to modify the SOIL_COLOR properties. I was expecting to find its properties in the paramfile, but they are not there.

I was wondering if someone can orientate me on how to create the new soil color type. I am assuming this information is in some input file.

Thank you.

 

[oleson]

In release-cesm2.1.5, soil color translates into albedo in biogeophys/SurfaceAlbedoMod.F90:

if (mxsoil_color == 8) then
albsat(1:8,1) = (/0.12_r8,0.11_r8,0.10_r8,0.09_r8,0.08_r8,0.07_r8,0.06_r8,0.05_r8/)
albsat(1:8,2) = (/0.24_r8,0.22_r8,0.20_r8,0.18_r8,0.16_r8,0.14_r8,0.12_r8,0.10_r8/)
albdry(1:8,1) = (/0.24_r8,0.22_r8,0.20_r8,0.18_r8,0.16_r8,0.14_r8,0.12_r8,0.10_r8/)
albdry(1:8,2) = (/0.48_r8,0.44_r8,0.40_r8,0.36_r8,0.32_r8,0.28_r8,0.24_r8,0.20_r8/)
else if (mxsoil_color == 20) then
albsat(1:20,1) = (/0.25_r8,0.23_r8,0.21_r8,0.20_r8,0.19_r8,0.18_r8,0.17_r8,0.16_r8,&
0.15_r8,0.14_r8,0.13_r8,0.12_r8,0.11_r8,0.10_r8,0.09_r8,0.08_r8,0.07_r8,0.06_r8,0.05_r8,0.04_r8/)
albsat(1:20,2) = (/0.50_r8,0.46_r8,0.42_r8,0.40_r8,0.38_r8,0.36_r8,0.34_r8,0.32_r8,&
0.30_r8,0.28_r8,0.26_r8,0.24_r8,0.22_r8,0.20_r8,0.18_r8,0.16_r8,0.14_r8,0.12_r8,0.10_r8,0.08_r8/)
albdry(1:20,1) = (/0.36_r8,0.34_r8,0.32_r8,0.31_r8,0.30_r8,0.29_r8,0.28_r8,0.27_r8,&
0.26_r8,0.25_r8,0.24_r8,0.23_r8,0.22_r8,0.20_r8,0.18_r8,0.16_r8,0.14_r8,0.12_r8,0.10_r8,0.08_r8/)
albdry(1:20,2) = (/0.61_r8,0.57_r8,0.53_r8,0.51_r8,0.49_r8,0.48_r8,0.45_r8,0.43_r8,&
0.41_r8,0.39_r8,0.37_r8,0.35_r8,0.33_r8,0.31_r8,0.29_r8,0.27_r8,0.25_r8,0.23_r8,0.21_r8,0.16_r8/)
else
write(iulog,*)'maximum color class = ',mxsoil_color,' is not supported'
call endrun(msg=errMsg(sourcefile, __LINE__))
end if

 

[Jairo]

Thank you so much! That is working for me.

Best,