Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

convert %cs_i_indg and %cs_j_indg to simple arrays in Assim GridComp #111

Closed
gmao-rreichle opened this issue Feb 14, 2020 · 4 comments · Fixed by #295
Closed

convert %cs_i_indg and %cs_j_indg to simple arrays in Assim GridComp #111

gmao-rreichle opened this issue Feb 14, 2020 · 4 comments · Fixed by #295
Assignees

Comments

@gmao-rreichle
Copy link
Contributor

The fields %i_indg and %j_indg in TileCoord are w.r.t. to the global grid that cuts tiles (regardless of whether that is a cube-sphere, EASE, or lat-long grid).

The second set of fields %cs_i_indg and %cs_j_indg should be w.r.t. to the (lat/long) "hash" grid that is used to quickly identify tiles that are within the lat/long grid cell(s) of the "hash" grid. The name "cs_[x]_indg" is confusing, because the grid in question is NOT a cs grid. Also, the current implementation actually overwrites %i_indg and %j_indg and, consequently changes their meaning.

Solution: remove %cs_i_indg and %cs_j_indg from the TileCoord structure. In the Assim Grid Comp, use simple hash_i/j_ind arrays for "hash" functionality and leave the original %i_indg and %j_indg unchanged.

@gmao-rreichle
Copy link
Contributor Author

gmao-rreichle commented Apr 29, 2020

@saraqzhang : For @weiyuan-jiang to work on this issue, he would need a test case for the assimilation on the cube-sphere tiles. Please provide him with a 1-day verified GEOSldas run that is doing assimilation on the cube-sphere tiles (otherwise similar to the GLOBAL/assim test case that mimics the SMAP L4_SM system). The test case doesn't have to be scientifically perfect, just some result that @weiyuan-jiang can use for zero-diff testing. We may want to make this test case (assim on cs tiles) part of the nightly test suite.

@saraqzhang
Copy link
Contributor

@gmao-rreichle @weiyuan-jiang
/discover/nobackup/qzhang/C180ldas/ can be used for this test. it is assimilation on cs, the forcing is FP latlon. The built is v1790beta3, on sles11. It has 10day assimilation results.

I am also running v1790beta4-sles12 innovation-ens24 . I can turn assim. on for a couple of days so to get assimilation results if you prefer to use results from this system.

@weiyuan-jiang
Copy link
Contributor

After we create a lat-lon grid for landpert, we need a tile_coord being with the tile grid. However, the values (like frac_cell, frac_pfap, area) in tile_coord are not usable anymore. Why not create a hash_tile_coord instead of creating hash_i_indg? we set all those values no_data so it won't cause any trouble in the future. It avoid making the situation worse.

@gmao-rreichle gmao-rreichle linked a pull request Aug 18, 2020 that will close this issue
@gmao-rreichle
Copy link
Contributor Author

addressed with pull request #295

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
3 participants