I am trying to assimilate IASI radiance data obtained from NCEP GDAS Satellite Data server (NCAR RDA Dataset ds735.0 that can be directly assimilated with WRFDA)
into WRF-ARW 4.6.0, using RTTOV 12.3.
However, the 3DVAR process stopped at the module "rttov_check_profiles.F90", with the reason "invalid zenith angle".
In the subroutine that reads the IASI data (da_read_obs_bufriasi.inc) there is a check that transforms the satellite zenith angle in negative;
infact, if (ifov <= 60) sat_zenang = -sat_zenang.
Thus, positive angles are transformed in negative ones and these values do not pass the zenith angle check in the subroutine rttov_check_profiles.F90, which works only with positive angles.
I attached the log files (rs.error.0000, rsl.out.0000) and the namelist.input file.
I think it's a bug that needs to be fixed between the subroutine that handles the IASI data and RTTOV.
I wonder if someone can give me some ideas to solve this problem or some person contacts.
Thank you !
into WRF-ARW 4.6.0, using RTTOV 12.3.
However, the 3DVAR process stopped at the module "rttov_check_profiles.F90", with the reason "invalid zenith angle".
In the subroutine that reads the IASI data (da_read_obs_bufriasi.inc) there is a check that transforms the satellite zenith angle in negative;
infact, if (ifov <= 60) sat_zenang = -sat_zenang.
Thus, positive angles are transformed in negative ones and these values do not pass the zenith angle check in the subroutine rttov_check_profiles.F90, which works only with positive angles.
I attached the log files (rs.error.0000, rsl.out.0000) and the namelist.input file.
I think it's a bug that needs to be fixed between the subroutine that handles the IASI data and RTTOV.
I wonder if someone can give me some ideas to solve this problem or some person contacts.
Thank you !