WRF-SUEWS momentum and heat coupling

hamidrezaomidvar commented 5 years ago

I opened this issue so we can discuss the momentum coupling of WRF-SUEWS here.

hamidrezaomidvar commented 5 years ago

Here are the variables that we need to make sure are updated in the land surface model for passing to the pbl scheme. @sunt05 would you please help me on this to see which ones are necessary to correct, and how are they calculated:

`in`/`inout` by `MYJ`

[x] AKHS: SFC EXCH COEFF FOR HEAT: Tstar*Ustar/dT, note the unit [m s-1]
[x] AKMS: SFC EXCH COEFF FOR MOMENTUM: Ustar*Ustar/dU, note the unit [m s-1]
[x] USTAR: friction velocity, calculated by SUEWS and can be passed into WRF.
[x] ZNT: Time varying roughness length (is is heat or momentum??): momentum, calculated by SUEWS.
[x] QSFC: already fine and is calculated in SUEWS side
[x] PBLH: output of PBL schemes: we are fine.
[x] THZ0: Potential temperature at ZNT. Do we need to calculate this in SUEWS or it is calculated elsewhere if we pass the correct ZNT
[x] UZ0, VZ0: wind speed at ZNT
[x] QZ0: I think these are being calculated in the surface layer scheme using the ZNT of previous time step. But needs to double check
[x] HFX (the current pbl scheme does not use this): already fine and is calculated in SUEWS side
[x] TSK: already fine and is calculated in SUEWS side

~~- [ ] CT: COUNTERGRADIENT TERM. I am not sure about this variable!!~~

sunt05 commented 5 years ago

I've converted the list into todo list so we can mark the resolved ones.

hamidrezaomidvar commented 5 years ago

flhc: not sure what is the difference between this and akhs

sunt05 commented 5 years ago

CT is an out variable in MYJSFC. So I think we don't have to consider it.

hamidrezaomidvar commented 5 years ago

flhc: not sure what is the difference between this and akhs

The difference is the units and is not needed for the current pbl scheme , but might be needed for other pbl schemes. It is the output of MYJSFC

sunt05 commented 5 years ago

flhc: not sure what is the difference between this and akhs
FLHC=RLOW*CP*AKHS

sunt05 commented 5 years ago

I agree with you that we sort out the MYJ scheme first; then move onto other schemes.

hamidrezaomidvar commented 5 years ago

RMOL: 1./Monin Ob. Length

seems it is a INOUT for MYJSFC but not being calculated or used in this subroutine!

hamidrezaomidvar commented 5 years ago

Z0BASE used by MYJSFC: no definition is given though

sunt05 commented 5 years ago

THZ0 seems not to be calculated by any physics scheme but associated with grid; so I assume it's provided by the core solver.

hamidrezaomidvar commented 5 years ago

We finalised the variables and I start the implementation.

hamidrezaomidvar commented 5 years ago

The implementation is completed. I am testing this, and will close these issue when the test is completed.

sunt05 commented 5 years ago

Do we output z0 in this implementation? It would be good to know how z0 varies across grids, which is also a very demonstrative feature of SUEWS.

hamidrezaomidvar commented 5 years ago

not yet. I can output it to see how it looks like.

hamidrezaomidvar commented 5 years ago

I tested this and the system is working well. I am going to output z0 as well and compare the z0 output with velocities to see if we are capturing the feedback.

sunt05 commented 5 years ago

well done!

Urban-Meteorology-Reading / WRF-SUEWS