The surface data produced using the CLM are all 0

[weixinyu]

Dear all
I defined my own resolution and region and successfully followed the process to produce map, domain and surfer data, although the surface data production shows that it has been successful, at the end of the day, when I viewed the surfdata_syh_78pfts_CMIP6_simyr1850_c240617 it showed all zeros, and I put the output logs in a text file.
My procedure is as follows：
1、Create a new case "syh" with a resolution of 0.1°.
The latitude and longitude of the four points are 101,40, 105,40, 101,37, 105,37
And I successfully created map data, domain data and surfdata，as shown in the figure
But I don't know why, when I open the last made suferdata, it shows all 0. Please help me with this, thank you very much!

 

[oleson]

I'm not sure we can help with this since we don't have access to the machine you are using, but please attach your surface dataset namelist for 1850. What version of cesm/clm are you using? I guess you should also check the domain file you created to make sure it has land in it.

 

[weixinyu]

I'm not sure we can help with this since we don't have access to the machine you are using, but please attach your surface dataset namelist for 1850. What version of cesm/clm are you using? I guess you should also check the domain file you created to make sure it has land in it.

Thank you very much for your reply.
I am using CLM version 5.0 and I use the command. /create_newcase --case Hist_syh --res CLM_USRDAT --compset IHistClm50Bgc successfully created Hist_syh,I put the nc data image of my domain in the attachment, please you you check it, thank you very much!

 

[weixinyu]

I'm not sure we can help with this since we don't have access to the machine you are using, but please attach your surface dataset namelist for 1850. What version of cesm/clm are you using? I guess you should also check the domain file you created to make sure it has land in it.

Another thing I thought of is that when I created the MAP data, domain data I didn't let him define it as clm version 5.0, I see that the default generated is version 4.5, could this be the reason?

 

[oleson]

Check the mask variable in your domain.lnd.... file to make sure there is land. Land points should have 1 for the mask.

What is the output of the command "git describe" in your clm5.0 directory?

 

[weixinyu]

Dear oleson
As you said I checked the domain mask and found that the value is 0. Is there a problem in making my domain file? How should I modify it?
Besides that, the output value of git describe is release-clm5.0.03, how to change it, thanks for your reply!

 

[oleson]

If your mask in your domain.lnd..... file is all zeros then yes there is a problem in making your domain file. You'll need to look at the process you used to create that file. Maybe your mapping file is incorrect.

 

[weixinyu]

If your mask in your domain.lnd..... file is all zeros then yes there is a problem in making your domain file. You'll need to look at the process you used to create that file. Maybe your mapping file is incorrect.

Thank you very much for your reply.
How should I see if my mapping file is correct or not?

 

[weixinyu]

如果您的掩码在 domain.lnd 中.....文件都是零，那么是的，制作域文件有问题。您需要查看用于创建该文件的过程。可能您的映射文件不正确。

我可以麻烦你解释一下吗？因为映射数据很复杂，需要大量的内存，我需要用别人的工作站来做，我用别人的工作站做了测绘数据，然后复制到我自己的电脑上，这有什么效果吗？使用相同的 clm 配置。
您能具体告诉我如何查看映射数据是否正确吗？
谢谢！

 

[oleson]

I think I would look at the variables mask_a and mask_b in that file. I think mask_a would have all zeros and mask_b would have all ones.

 

[weixinyu]

Dear oleson
I made a total of 17 mapping data, I opened each one to view it and found that mask_a is not 0 and mask_b is not 1
Can I ask what is the reason for this? I have put my data details in the attachment below, could you please take a look, thank you very much!

 

[oleson]

I would look at the mapping file that was specified as an argument in the gen_domain step.

 

[weixinyu]

I would look at the mapping file that was specified as an argument in the gen_domain step.

Dear oleson
I'm very sorry because I'm a beginner in clm, could you please point out specifically what file you're looking at in that place? I didn't understand a bit of what you were saying, thank you very much for your reply!

 

[weixinyu]

I would look at the mapping file that was specified as an argument in the gen_domain step.

Dear oleson
I looked at mask_a and mask_b in my map file and they were 0 and 1 respectively, other than that I reran it and found that making the mksurfdata_map showed
* number of non-zero elements: 0
(gridmap_map_read) NetCDF: Index exceeds dimension bound
The result of the data generated in this section is 0.
How should this be solved? Could you please answer again?

 

[oleson]

Please confine your posts that are related to this subject to a single thread (this one here). It's confusing to have to go back and forth between threads. So far I don't see the cause of the problem. At this point, I ask that you post all of the steps you took from beginning to end that produced the surface dataset with the zeros. Perhaps then we can see the source of the problem or as a last resort try to replicate the process here on our machine.