ERROR: GLACIER_REGION NOT on surfdata fileERROR in glcBehaviorMod.F90

[yinjiyuan]

I'm running the latest version of CESM with compset F2000climo. The model has an error:
ERROR:
ERROR: GLACIER_REGION NOT on surfdata fileERROR in glcBehaviorMod.F90 at line
662
I can't find the variable GLACIER_REGION. What should I do?

 

[katec]

Hi, there is a lot of work going on right now to develop the needed surface data sets for all of the compsets. Could you please be specific about exactly which version of CESM you are running, and which grid resolution? Are you trying to use your own surface data file or the one that is specified by the compset? Could you please also send the path specified for "surfdat" in your "lnd_in" file?

 

[yinjiyuan]

嗨，现在正在进行大量工作，以开发所有竞争对手所需的表面数据集。您能否具体说明您正在运行的是哪个版本的CESM，以及哪种网格分辨率？您是尝试使用自己的曲面数据文件还是竞争对手指定的曲面数据文件？您能否也请在您的“lnd_in”文件中发送为“surfdat”指定的路径？

The version I am running is CESM2.1.3. The grid resolution is f09g17. I specified the surface file with :
SURFDATA_DIR=${CESM_ROOT}/input_data/lnd/clm2/surfdata
forest_filename=surfdata_0.9x1.25_simyr2000_c110921.nc
"echo "fsurdat = '"${SURFDATA_DIR}"/"${forest_filename}"'" >> user_nl_clm" and "echo "fsurdat = '"${SURFDATA_DIR}"/"${forest_filename}"'" >> lnd_in".

 

[yinjiyuan]

Hi, there is a lot of work going on right now to develop the needed surface data sets for all of the compsets. Could you please be specific about exactly which version of CESM you are running, and which grid resolution? Are you trying to use your own surface data file or the one that is specified by the compset? Could you please also send the path specified for "surfdat" in your "lnd_in" file?

I find that the file named "surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr2000_c190214.nc" contains "GLACIER_REGION". So if I use my own surface data file called "surfdata_0.9x1.25_simyr2000_c110921.nc", it doesn't contain "GLACIER_REGION". Can I put the "GLACIER_REGION" from "surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr2000_c190214.nc" into "surfdata_0.9x1.25_simyr2000_c110921.nc" for use?

 

[katec]

Hi there. Your modified file DOES need to contain a field for GLACIER_REGION. The GLACIER_REGION field tells CTSM which areas will contain CISM-run glacier models or expect data that could be used in a glacier model and to handle water fluxes appropriately. For more information on Glacier regions, this website might be helpful:
https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Glacier/CLM50_Tech_Note_Glacier.html

I'll move this thread over to the CTSM forum so you can get help from the experts there about exactly how to set that up. Since I don't know what else is in your current surface data file, it is difficult to know if the field on the default file will make sense. It might be helpful for the CTSM scientists if you could describe what changes you have made in your surface dataset.

 

[yinjiyuan]

Hi there. Your modified file DOES need to contain a field for GLACIER_REGION. The GLACIER_REGION field tells CTSM which areas will contain CISM-run glacier models or expect data that could be used in a glacier model and to handle water fluxes appropriately. For more information on Glacier regions, this website might be helpful:
2.13. Glaciers — ctsm CTSM master documentation

I'll move this thread over to the CTSM forum so you can get help from the experts there about exactly how to set that up. Since I don't know what else is in your current surface data file, it is difficult to know if the field on the default file will make sense. It might be helpful for the CTSM scientists if you could describe what changes you have made in your surface dataset.

Hi！I put the "GLACIER_REGION" from "surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr2000_c190214.nc" into "surfdata_0.9x1.25_simyr2000_c110921.nc" for use. Finally, there was no error. But a new problem has arisen. The following:
pio_support::pio_die:: myrank= -1 : ERROR: nf_mod.F90: 1288 :
NetCDF: Invalid dimension ID or name
Do you know how to solve this problem? Thanks!

 

[oleson]

"Since I don't know what else is in your current surface data file, it is difficult to know if the field on the default file will make sense. It might be helpful for the CTSM scientists if you could describe what changes you have made in your surface dataset."

 

[yinjiyuan]

"Since I don't know what else is in your current surface data file, it is difficult to know if the field on the default file will make sense. It might be helpful for the CTSM scientists if you could describe what changes you have made in your surface dataset."

I didn't make any changes to the surface data. I just used "/trunk/inputdata/lnd/clm2/surfdata/surfdata_0.9x1.25_simyr2000_c110921.nc".
If I do not execute the following command:
SURFDATA_DIR=${CESM_ROOT}/input_data/lnd/clm2/surfdata
forest_filename=surfdata_0.9x1.25_simyr2000_c110921.nc
"echo "fsurdat = '"${SURFDATA_DIR}"/"${forest_filename}"'" >> user_nl_clm"
Mode is used by default. That is "/trunk/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr2000_c190214.nc"
I changed the surface data to use "/trunk/inputdata/lnd/clm2/surfdata/surfdata_0.9x1.25_simyr2000_c110921.nc". Then it made a mistake.

 

[oleson]

I'm not sure why you used surfdata_0.9x1.25_simyr2000_c110921.nc. That is a very old surface dataset (created in 2011) to use in the most recent released version of the model (cesm2.1.3). Generally, you should use the surface dataset that is specified when you create your case with your desired compset and resolution.

 

[yinjiyuan]

Yes. Previously when I used CESM 1.2.2, the default surface data was "/trunk/inputdata/lnd/clm2/surfdata/surfdata_0.9x1.25_simyr2000_c110921.nc". But when I run CESM 2.1.3 now, the default surface data is "/trunk/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr2000_c190214.nc". So I can just use the default data?

 

[oleson]

Yes, exactly, the surface datasets, among other things, change with different model versions. Previous versions of the surface datasets will not likely work with newer versions of the code.

 

[yinjiyuan]

Yes, exactly, the surface datasets, among other things, change with different model versions. Previous versions of the surface datasets will not likely work with newer versions of the code.

Okay, I got it. Thank you very much!!!