Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Outputing roughness length

X

xtian15

Member
Hi,
I am trying to output the roughness length in either the diag or history files and have two question.
1. I am seeing a variable `gz1oz0` with a description of "log(z/z0) where z0 is roughness length". The `z0` named "roughness length" here is what I am after. But further down the file, there are two variables `z0` named "roughness height" and `znt` named "roughness length". Which one among `z0` and `znt` is the `z0` in "log(z/z0) where z0 is roughness length"?
2. I tried adding either z0 or znt into the list of output within Registry.xml, recompile, and also in stream_list.atmosphere.output in actual jobs. But the complaint of "ERROR: MPAS IO Error: Undefined field" keeps showing up and no `z0/znt` is written.
Please advise on which and how the roughness length can be output in history or diag files.
Thank you!
 
Hi,
Please see my answers below:
xtian15 said:
Hi,
I am trying to output the roughness length in either the diag or history files and have two question.
1. I am seeing a variable `gz1oz0` with a description of "log(z/z0) where z0 is roughness length". The `z0` named "roughness length" here is what I am after. But further down the file, there are two variables `z0` named "roughness height" and `znt` named "roughness length". Which one among `z0` and `znt` is the `z0` in "log(z/z0) where z0 is roughness length"?
Click to expand...
z0 is roughness height, which is the one used in log(z/z0)
xtian15 said:
2. I tried adding either z0 or znt into the list of output within Registry.xml, recompile, and also in stream_list.atmosphere.output in actual jobs. But the complaint of "ERROR: MPAS IO Error: Undefined field" keeps showing up and no `z0/znt` is written.
Click to expand...
Please keep Registry.xml unchanged, and simply add 'z0' and 'znt' to stream_list.atmosphere.output. Then rerun the case. Let me know whether the model works and whether you have the two variable sin your history output files.
xtian15 said:
Please advise on which and how the roughness length can be output in history or diag files.
Thank you!
Click to expand...
 
Ming Chen said:
Hi,
Please see my answers below:

z0 is roughness height, which is the one used in log(z/z0)

Please keep Registry.xml unchanged, and simply add 'z0' and 'znt' to stream_list.atmosphere.output. Then rerun the case. Let me know whether the model works and whether you have the two variable sin your history output files.
Click to expand...

Thanks for your response! "simply add 'z0' and 'znt' to stream_list.atmosphere.output" will cause an error message written out in error logs as below:

Code: 
----------------------------------------------------------------------
Beginning MPAS-atmosphere Error Log File for task     159 of     160
    Opened at 2024/09/03 13:21:35
----------------------------------------------------------------------

ERROR: MPAS IO Error: Undefined field
The simulation will continue and no z0 or znt variables will be written out as instructed in the history files.
 
I just did a quick test and I can successfully output znt and z0 to my history files. Below is the variables included in my
"stream_list.atmosphere.output"
Code: 
rthratensw
rthratenlw
rthcuten
znt
z0

The MPAS version I used is V8.1.

Can you double check to make sure you follow the correct approach?
 
You must log in or register to reply here.
Top