search

amansnama / 2022_Yao_etal_JAMES

Scripts used for study: Implementation and evaluation of irrigation techniques in the Community Land Model
0 stars 0 forks source link

Is nlevgrnd == nlevmaxurbgrnd? #1

Open amansnama opened 2 weeks ago

amansnama commented 2 weeks ago

Replaced nlevgrnd by nlevmaxurbgrnd and got Yi Yao's irrigation running. These are two different variables. Which is the correct variable to use?

amansnama commented 1 week ago

In main/clm_varpar.F90,

public    :: nlevsoi               ! number of hydrologically active soil layers
  integer, public    :: nlevsoifl             ! number of soil layers on input file
  integer, public    :: nlevgrnd              ! number of ground layers 
                                              ! (includes lower layers that are hydrologically inactive)
  integer, public    :: nlevurb               ! number of urban layers
  integer, public    :: nlevmaxurbgrnd        ! maximum of the number of ground and urban layers

nlevgrnd should be number of soil layers + bedrock layers. (20+5 by default). nlevmaxurbgrnd is max of ground and urban layers main/clm_varpar.F90: nlevmaxurbgrnd = max0(nlevurb,nlevgrnd)