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

Missing Data in tslist

K

Khin Myat

Member
Hi,
I run the simulation with 62 vertical levels. When I check the ts list u,v,w and ph, it has only 61 levels. Why are there not 62 levels in the tslist?
Doesn't the u,v,w and ph calculate the wind on the first eta level (eta =1)?
 
Can you check how many levels do you have for U and V, and how many levels do you have for W and PH?
 
Dear Ming Chen,
Thank you for your reply. I checked again. In each pfx.dNN.UU, VV, WW, and PH file, there are 61 levels.
But, when I take the PHB value from the wrfout file, there are 62 levels. When I take the "ua,va;wa" from wrfout file, there are 61 levels because they are on mass point.
As far as I understand, if I activate "unstagger" function, there will be 61 levels. But I did not do "unstagger'. So, I am confused.

Thank you in advance for your time.
Best,
Khin
 
Dear Ming chen,
Thank you very much for your reply.
I still haven't got your point. As far as I understand, u, v points are at the center of left and right, top and bottom of grid cell. So I thought that they were on the same level of mass points. Doesn't it mean that they are on the half eta level?
When I read the tslist,it is described that " by default ,the u, v, and w component winds are output on the staggered grid" . To my understanding, stagger points are on the full eta level and mass points are on the half eta level. PH values are calculated on full level. If so, why are there only 61 levels for u,v,w and ph in tslist vertical profile?
Where do I misunderstand? Would you please explain to me when you have free time?
Thank you in advance for your time.

Best,
Khin
 
Khin,
Please take a look the WRF Technote here, where you can find how staggering grids are employed.
Specifically, Fig. 3.2 illustrates horizontal and vertical structures of grids and variable distributions.

Hope this is helpful for you. let me know if you still have questions.
 
Dear Ming Chen,
Thank you very much for your information.
I still have questions.
If geopotential height is defined at w points, is the geopotential height from ts list stagger or destagger value? As far as I understand, it is destagger value. Is it right?
If the geopotential height is destagger, should the u,v w value be destagger?
Are mass point levels the same as half eta model levels?
Thank you very much in advance for your time.

Best,
Khin
 
geopotential height is defined in full model level, and it should be in full levels in tslist files.
mass points are located at half model levels.
u and v are not destagged if you don't specify the option tslist_unstagger_winds.
 
Dear Ming Chen,
Thank you very much for your reply.
I also checked and looked ph and phb values from wrfout file with ncl command and there are also 62 levels. I also knew that ph values from wrfout file and ts list files were also different at the same grid point and at the same time. It is strange if they are on the same level. I also attached here. I think it should be the same.
If geopotential height should be in full levels in tslist files, all my ts list files are missing one level. Where should I check to solve that problem?

I am sorry if I am bothering you.
Thank you very much for your time.

Best,
Khin
 

Attachments

  • PH+PHB.png
    PH+PHB.png
    13.5 KB · Views: 3
  • PH.png
    PH.png
    18.5 KB · Views: 3
Did you specify max_ts_level in your namelist.input? If so, what is the value you set?

This is the highest model level for time series output.
 
Dear Ming Chen,
I found that "
"max_ts_level : number of model levels for time series vertical profiles (default is 15). The maximum number of max_ts_level is e_vert-1 (the number of half layers in the model run)" in wrf menu.
I think this is the reason why I get only 61 levels.

Thank you very much for your time.

Best ,
Khin
 
Khin,
max_ts_level must be smaller than e_vert. if you set it wrong, WRF will automatically changes it to e_vert -1.
 
Dear Ming Chen,
If the e_vert is 62, the max_ts_level must be 61. If so, the ts list file doesn't produce the variation of the PH value of the topmost full layer. Does " the number of half layers in the model run" not mean that PH values are on the half-model layer?
Are all stagger points of U,V, and W on the full model level ?
Are only stagger points of W on the full model level?

Thank you so much for your time.
Best,
Khin
 
Last edited:
Khin Myat said:
Dear Ming Chen,
If the e_vert is 62, the max_ts_level must be 61. If so, the ts list file doesn't produce the variation of the PH value of the topmost full layer. Does " the number of half layers in the model run" not mean that PH values are on the half-model layer?
Are all stagger points of U,V, and W on the full model level ?
Are only stagger points of W on the full model level?

Thank you so much for your time.
Best,
Khin
Click to expand...
Hi,Khin
Can you tell me if the first layer height inside the ts.UU file corresponds to the layer with eta_level=0, or does it correspond to the layer above eta=0. Also, is each layer of U and V in the ts file a full mode layer.
I did a test and found that UU increases monotonically in height, but I used the ua variable in wrfout and found that the 4th layer ua (lowest being 0) is smaller than the 4th value in ts.UU, which I think is problematic.
 
Ming Chen said:
Khin,
I believe W and PH should be on full model level, while U and V should be on half model level.
Click to expand...
Dear Ming Chen, as you said, uu is in the model half layer, so can I assume that the uu of ts is the same as the uava of wrfout output?But when I extract the ua and ts.uu at the same moment I realize they are not the same.As a comparison, I also compared the u10v10 in the ts file and found that they are the same as the u10v10 in the corresponding moment of wrfout, do you know anything about this please.
In general, I would like to know if each level inside the ts.uu is a model half layer or a whole layer.
 
In ts.UU, the u and v components are located at half model levels.

Note that ts.UU and ts.VV are winds at every time step, while wrfout gives values at specified output times.

When you compare ts and wrfout, please make sure their output times match.
 
Note that U and V in TS files (time series output ) are earth-relative, while U and V in wrfout are grid-relative. This explains why there are differences between winds form TS and wrfout files.
 
You must log in or register to reply here.
Top