Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

NSOIL8 not read properly when running sys_cpl != 1 #571

Open
aubreyd opened this issue Jul 8, 2021 · 2 comments
Open

NSOIL8 not read properly when running sys_cpl != 1 #571

aubreyd opened this issue Jul 8, 2021 · 2 comments
Assignees
Labels
bug Something isn't working community dev Issue open for community development

Comments

@aubreyd
Copy link
Collaborator

aubreyd commented Jul 8, 2021

Current Behavior (issue identified by @TimLahmers... thanks Tim!)

It looks like the ZSOIL8 data field (soil layer discretization) is currently being read into WRF-Hydro through a hrldas routine. When running NOT under hrldas (so sys_cpl != 1, e.g., LIS coupling), these routines are not called so the ZSOIL8 variable is not being set.

The ZSOIL8 variable appears to be read in hrldas_cpl_HYDRO_ini (maybe based on the specs in namelist.hrldas?) and then ALSO in the orchestrator in the hydro.namelist read, but the place where those values are mapped to the nlst object in the orchestrator is commented out (see lines 639-642 in https://github.com/NCAR/wrf_hydro_nwm_public/blob/master/trunk/NDHMS/OrchestratorLayer/config.f90).

Possible Solutions (courtesy of @danrosen25)

1 - Activate those commented out lines in the orchestrator when sys_cpl != 1
2 - Call the hrldas init routines from the NUOPC cap

@aubreyd aubreyd added bug Something isn't working community dev Issue open for community development labels Jul 8, 2021
@danrosen25
Copy link
Collaborator

@rcabell
I plan to make the following change if you agree?

if(sys_cpl .ne. 1) then
    nlst(did)%NSOIL = NSOIL
    allocate(nlst(did)%ZSOIL8(NSOIL))
    nlst(did)%ZSOIL8 = ZSOIL8
endif

https://github.com/NCAR/wrf_hydro_nwm_public/blob/master/trunk/NDHMS/OrchestratorLayer/config.f90#L639

@rcabell
Copy link
Collaborator

rcabell commented Jul 23, 2021

Yep, that looks fine to me. @kafitzgerald or @aubreyd, any concerns?

@rcabell
I plan to make the following change if you agree?

if(sys_cpl .ne. 1) then
    nlst(did)%NSOIL = NSOIL
    allocate(nlst(did)%ZSOIL8(NSOIL))
    nlst(did)%ZSOIL8 = ZSOIL8
endif

https://github.com/NCAR/wrf_hydro_nwm_public/blob/master/trunk/NDHMS/OrchestratorLayer/config.f90#L639

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working community dev Issue open for community development
Projects
None yet
Development

No branches or pull requests

3 participants