171 | | * For sharing through the module 'mod_assimilation', we further initialize an array 'coords_l' containing the coordinates that describe the local domain. These coordinates have to describe one location in space that is used in the OMI observation modules to compute the distance from observations. This requires that the coordinates in 'coords_l' have the same units as those used for the observations. For geographic distance computations, the unit of the coordinates needs to be radian (0, 2*pi) or (-pi,pi). |
172 | | * Any form of local domain is possible as long as it can be describe as a single location. If observations are only horizontally distributed (a common situation with satellite data in the ocean), the local analysis domain can be a single vertical column of the model grid. In this case, the size of the state in the local analysis domain will be just the number of vertical grid points at this location and the horizontal coordinates are used in 'coords_l' |
| 171 | * For sharing through the module `mod_assimilation`, we further initialize an array `coords_l` containing the coordinates that describe the local domain. |
| 172 | * These coordinates have to describe one location in space that is used in the OMI observation modules to compute the distance from observations. |
| 173 | * The coordinates in `coords_l` have the same units as those used for the observations |
| 174 | * For geographic distance computations, the unit of the coordinates needs to be radian, thus (0, 2*pi) or (-pi,pi) for longitude and (-pi/2, pi/2) for latitude. |
| 175 | * Any form of local domain is possible as long as it can be describe as a single location. |
| 176 | * If the local domain is a single grid point, `dim_l` will be the number of model variables at this grid point. |
| 177 | * The local analysis domain can also be a single vertical column of the model grid if observations are only horizontally distributed (a common situation with satellite data in the ocean). |
| 178 | * In this case, `dim_l` will be the number of vertical grid points at this location times the number of model fields that exist in the vertical, plus possible variables at e.g. the surface. |
| 179 | * In this case only the horizontal coordinates are used in `coords_l`. |