How do I get the LCZ part of CGLC_MODIS_LCZ_global to work?

[Rtsquared]

I am not sure why, but when I try using the above dataset and run geogrid.exe the result is pure MODIS- no LCZ categories at all. The paper on the dataset says that the GEOGRID.TBL entry needs to have higher priority than default MODIS, but if the priority is different in the GEOGRID.TBL.ARW_LCZ file I get the following error:
Parsed 51 entries in GEOGRID.TBL
Processing domain 1 of 2
ERROR: Could not find interpolator sequence for requested resolution for LANDUSEF. The sources specified in namelist.wps is not listed in GEOGRID.TBL.

Adjusting the file slightly (see attached) so everything is priority 1 does let it run, but the result is pure MODIS (there also doesn't seem to be a difference between using "CGLC_MODIS_LCZ_global" and "CLGL_MODIS_LCZ_global+default" for my domain). I don't expect the LCZ data to actually work in my WRF run (it doesn't seem to be compatible with NoahMP), but I would like to copy over the building variables that are created after running real.exe as the default NUDAPT data for my cities is quite lacking. What am I missing here that would fix this problem?

As a secondary question, this one about NCAR's HPC systems in general (I am running things on Derecho), what software is installed for comparing model data to observations? When I was working at the EPA’s ORD center in Durham last year I used a tool they had developed called AMET and it looks like NCAR’s MELODIES MONET might be similar, though it’s focus appears to be atmospheric chemistry (I just need surface temperature, humidity, and winds). Is that what I should be trying to use, and if so how?

 

[kwerner]

Hi,
I first would like to apologize for the long delay in response. We have been occupied with preparing for our upcoming code release and have gotten behind with forum posts. Thank you so much for your patience. Since it's been so long, can you let me know if this is still an issue for you? Thanks.

 

[Rtsquared]

It is, I've been trying to fix a different issue in the meantime which requires making some edits to WRF's code which is tricky as I don't really know Fortran and it's not all well commented.

 

[kwerner]

I just ran a test and was able to create new landuse (LU_INDEX) fields when I use the LCZ option. The differences in my files and yours are that I did not modify my GEOGRID.TBL.ARW_LCZ file, and my geog_data_res setting in namelist.wps is listed as:

geog_data_res = 'cglc_modis_lcz+default','cglc_modis_lcz+default',

Yours is:

geog_data_res = 'CGLC_MODIS_LCZ_global','CGLC_MODIS_LCZ_global',

I'm sure you are, but make sure you are also changing the link in the WPS/geogrid directory for the GEOGRID.TBL. If you do a full listing in that directory, you should see

GEOGRID.TBL -> GEOGRID.TBL.ARW_LCZ

Just to note, I'm using WPSV4.5. When I run with the geog_data_res just set to 'default' for each domain, I only get 21 landuse types, but when I use 'cglc_modis_lcz+default' I get 58 types.

Can you try that and see if that makes any difference?

 

[Rtsquared]

I changed things to match what you have- I hadn't noticed the rel_path = cglc_modis_lcz:CGLC_MODIS_LCZ_global/ line in GEOGRID.TBL.ARW_LCZ which was the issue. I just took the name from the directory where the tables are located and didn't realize it had been renamed in the geogrid table which is why the priority 1 part had been failing. It works now.

Do you have any suggestions about the second half of my question?

 

[kwerner]

I'm so glad to hear that you were able to get past the error! Thanks for the update.

Regarding your second question - since this is a different question than the original one posted here, can you create a new thread to ask about it? We like to keep threads separate with only a single issue per thread - this helps with readability and searching in the future. When you post it, please provide some additional details, like what exactly you're trying to do, why you need to do it, the version of the code, and attach the files you've modified. We aren't really able to help with code modifications due to a lack of resources, but if it's something simple, we may be able to guide you in the right direction.

 

[likeyou511]

I just ran a test and was able to create new landuse (LU_INDEX) fields when I use the LCZ option. The differences in my files and yours are that I did not modify my GEOGRID.TBL.ARW_LCZ file, and my geog_data_res setting in namelist.wps is listed as:

geog_data_res = 'cglc_modis_lcz+default','cglc_modis_lcz+default',

Yours is:

geog_data_res = 'CGLC_MODIS_LCZ_global','CGLC_MODIS_LCZ_global',

I'm sure you are, but make sure you are also changing the link in the WPS/geogrid directory for the GEOGRID.TBL. If you do a full listing in that directory, you should see

GEOGRID.TBL -> GEOGRID.TBL.ARW_LCZ

Just to note, I'm using WPSV4.5. When I run with the geog_data_res just set to 'default' for each domain, I only get 21 landuse types, but when I use 'cglc_modis_lcz+default' I get 58 types.

Can you try that and see if that makes any difference?

how to change the link to GEOGRID.TBL.ARW_LCZ

 

[kwerner]

how to change the link to GEOGRID.TBL.ARW_LCZ

Try this:

1. Go to your WPS/geogrid directory
2. Issue this command:

ln -sf GEOGRID.TBL.ARW_LCZ GEOGRID.TBL

Then if you issue "ls -ls" you should see that the GEOGRID.TBL is now linked to the LCZ table.