changing sfc albedo in SurfaceAlbedoMod.F90

[debasish@gi_alaska_edu]

Dear Sam,

According to your advice, I checked in the program SurfaceAlbedoMod.F90 and I found examples that make use of the correspondence between p and gridcell (g).

do p = lbp, ubp
g = pgridcell(p)
coszen_pft(p) = coszen_gcell(g)
end do

The dimension of my albedo change is 8192 (lon=128Xlat=64) gridcells but when I use the above mentioned conversion, I will get NaNS values in place of real values.

do p =lbp, ubp
g1 = pgridcell(p) ! g1 varies from 1 to 8192
alpha_pft(p)=alpha(g1)
end do

Any clue regarding the error?

Thanks

Debasish

 

[slevis]

Given that alpha and alpha_pft are new variables that you added to the code, I suspect that you have not initialized or assigned values to alpha.

My purpose as liaison is to help you get started in your work with the clm, not to help you with your debugging. So you will have to proceed on your own from here.

Sam