Skip to content

Feature/add_CrISobserror#1

Open
zhumingying wants to merge 12 commits intoibanos90:mainfrom
NCAR:feature/add_CrISobserror
Open

Feature/add_CrISobserror#1
zhumingying wants to merge 12 commits intoibanos90:mainfrom
NCAR:feature/add_CrISobserror

Conversation

@zhumingying
Copy link

No description provided.

ibanos90 and others added 12 commits July 12, 2022 13:21
**Description**
This PR fixes some bugs found to convert raw AHI observations from the AWS cloud in Himawari Standard Data (HSD) format to the IODA-v2 needed in `mpas-jedi`. It was added the capability to read files with segment number 01, which corresponds to the no segment division of HSD data. The reading of the variables was modified to read each one explicitly since otherwise it was not working properly (not sure why). Specifically the way the variable `idata`, which is the count (total number of pixels that results from multiplying the number of columns by the number of lines) needed for the radiance calculation is also modified. This modification results in the correct calculation of the radiance and then the brightness temperature. In addition, the variable for the `dateTime` in `epochtime` type was allocated since that was missing. The observation error was also added following the values used in GSI (https://www.nco.ncep.noaa.gov/pmb/codes/nwprod/gfs.v16.2.2/fix/fix_gsi/global_satinfo.txt) for clear-sky and a bug in the function that sets the errors was fixed as well.

**Closes** 
Issue #2
This PR adds the superobbing capability for AHI observations. Giving the half width of the superobbing grid as an entry parameter, the code computes the average brightness temperature in the box defined using that parameter. This code is adapted from https://github.com/wrf-model/WRF/blob/develop/var/da/da_radiance/da_read_obs_netcdf4ahi_jaxa.inc and follows the same implementation as for ABI observations from GOES-16. The debugging build flags for the goes_abi code are removed. The following figures show the brightness temperature results for AHI observations after using this converter for 2018-04-15T00 channel 8 (a) without superobbing and  (b) with superobbing on a grid of 59x59 pixels:
**a)**
![RAW_2018041500_b8coverage_test_20220721_ch8_08_12](https://user-images.githubusercontent.com/33061897/184459639-fa6a39bb-9d83-4458-8724-3b42922898b5.png)
**b)**
![RAW_2018041500_b8coverage_test_20220721_superobbing59X59](https://user-images.githubusercontent.com/33061897/184459648-9ac5d155-76a5-435c-be7a-aa6e5255743e.png)

**Issue closed**
Closes #1
@zhumingying
Copy link
Author

Modified define_mod.f90: Added obserror for CrIS-FSR

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants