Skip to content

Commit

Permalink
Fix meteosat 11 error messages in DA log files and addition of Meteos…
Browse files Browse the repository at this point in the history 
…at 11 into crtm. Correction for NOAA 19 for hirs4 and amsua in crtm descriptor (#1762)

TYPE: Bug fix

DESCRIPTION OF CHANGES:
Problem:
Meteosat 11 was added to the GFS bufr files in Feb 2018 meaning when the seviri bufr read was turned on in the namelist it would throw out an error message for each data point in the log file making it very large. The change to the error loop removes this message but also prepares the bufr read file to ingest meteosat 11 data into wrf. Dont think a crtm currently supports meteosat 11 but i have changed the VARBC.in file to include data for this sat number and added an msg-4-seviri.info file to the radiance_info folder so people can use data from this satellite if desired if crtm supports it. 
Crtm descriptor file was missing noaa 19 in hirs4 and amsua.
Added noaa 19 hirs info file as it seemed to be missing.
Added templates for noaa 15-17 HIRS3 and noaa 18-19 HIRS4 to VARBC.in as were missing.

FIles Changed:
var/da/da_radiance/da_read_obs_bufrseviri.inc
var/da/da_radiance/da_crtm_sensor_descriptor.inc
var/run/VARBC.in
var/run/radiance_info/msg-4-seviri.info
var/run/radiance_info/noaa-19-hirs.info

Solution:
Added identifier for meteosat 11 and then changed the upper limit of the identifier error if statement with the new identifier significantly higher than previous. Added data MSG4, NOAA 15-17 HIRS3 and NOAA 18-19 HIRS4 to VARBC.in. Added MSG-seviri.info file. Added m11 descriptor to crtm sensor descriptor so crtm recognizes it.
Added an "or" statement to the crtm descriptor file to recognize noaa 19 and not just noaa 18 in hirs4 crtm. Added hirs noaa 19 info file in the radiance info folder.

RELEASE NOTE: Updated bufrseviri.inc so that if seviri is turned on in WRFDA after Feb 2018 it suppresses many error messages in log file and meteosat 11 data is ingested into WRFDA to be used by the crtm.
Corrected crtm sensor info and added hirs 4 info file to recognize noaa 19 and not just noaa 18 in hirs 4 and amsua files. Updated VARBC.in to include templates of NOAA 15-17 HIRS3 and NOAA 18-19 HIRS4.
  • Loading branch information
@WQadam
WQadam committed Sep 13, 2022
1 parent f34b11d commit e288837
Show file tree
Hide file tree
Showing 5 changed files with 194 additions and 4 deletions.
5 changes: 3 additions & 2 deletions var/da/da_radiance/da_crtm_sensor_descriptor.inc
View file
Expand Up @@ -33,6 +33,7 @@ subroutine da_crtm_sensor_descriptor(nsensor,sensor_descriptor)
if (satellite_id == 1) crtm_sat='m08'
if (satellite_id == 2) crtm_sat='m09'
if (satellite_id == 3) crtm_sat='m10'
if (satellite_id == 4) crtm_sat='m11'
else if (trim(crtm_platform_name(platform_id)) == 'gcom-w') then
if (satellite_id == 1) crtm_sat='gcom-w1'
else
Expand All @@ -47,11 +48,11 @@ subroutine da_crtm_sensor_descriptor(nsensor,sensor_descriptor)
elseif ( trim(crtm_sensor_name(sensor_id)) == 'hirs' ) then
if (satellite_id <= 14) crtm_sensor='hirs2'
if (satellite_id >= 15 .and. satellite_id <= 17) crtm_sensor='hirs3'
if (satellite_id == 18) crtm_sensor='hirs4'
if (satellite_id == 18 .or. satellite_id == 19) crtm_sensor='hirs4'
elseif ( trim(crtm_sensor_name(sensor_id)) == 'avhrr' ) then
if (satellite_id <= 14) crtm_sensor='avhrr2'
if (satellite_id >= 15 .and. satellite_id <= 17) crtm_sensor='avhrr3'
if (satellite_id == 18) crtm_sensor='avhrr4'
if (satellite_id == 18 .or. satellite_id == 19) crtm_sensor='avhrr4'
else
crtm_sensor=crtm_sensor_name(sensor_id)
end if
Expand Down
5 changes: 4 additions & 1 deletion var/da/da_radiance/da_read_obs_bufrseviri.inc
View file
Expand Up @@ -294,7 +294,8 @@ subroutine da_read_obs_bufrseviri (obstype,iv,infile)
! SAID 55 is meteosat-8 or msg-1
! SAID 56 is meteosat-9 or msg-2
! SAID 57 is meteosat-10 or msg-3
if ( ( kidsat > 57) .or. ( kidsat < 55) ) then
! SAID 70 is meteosat-11 or msg-4
if ( ( kidsat > 70) .or. ( kidsat < 55) ) then
write(unit=message(1),fmt='(A,I6)') 'Unknown platform: ', kidsat
call da_warning(__FILE__,__LINE__,message(1:1))
end if
Expand All @@ -305,6 +306,8 @@ subroutine da_read_obs_bufrseviri (obstype,iv,infile)
satellite_id = 2
else if ( kidsat == 57 ) then
satellite_id = 3
else if ( kidsat == 70 ) then
satellite_id = 4
end if

if (clrsky) then ! asr bufr has no sza
Expand Down
159 changes: 158 additions & 1 deletion var/run/VARBC.in
View file
@@ -1,5 +1,5 @@
VARBC version 1.0 - Number of instruments:
42
48
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
Expand Down Expand Up @@ -870,6 +870,146 @@
17 17 0 0 0 0 0 0 0 0
18 18 0 0 0 0 0 0 0 0
19 19 0 0 0 0 0 0 0 0
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
1 15 0 19 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
1.0 9281.9 11619.0 294.4 30.5 7.2 80.0 1107.7
0.0 147.2 172.9 6.7 11.9 5.4 104.1 1955.8
10000 10000 10000 10000 10000 10000 10000 10000
-----> Chanl_id Chanl_nb Pred_use(-1/0/1) Param
1 1 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0
3 3 0 0 0 0 0 0 0 0
4 4 0 0 0 0 0 0 0 0
5 5 0 0 0 0 0 0 0 0
6 6 0 0 0 0 0 0 0 0
7 7 0 0 0 0 0 0 0 0
8 8 0 0 0 0 0 0 0 0
9 9 0 0 0 0 0 0 0 0
10 10 0 0 0 0 0 0 0 0
11 11 0 0 0 0 0 0 0 0
12 12 0 0 0 0 0 0 0 0
13 13 0 0 0 0 0 0 0 0
14 14 0 0 0 0 0 0 0 0
15 15 0 0 0 0 0 0 0 0
16 16 0 0 0 0 0 0 0 0
17 17 0 0 0 0 0 0 0 0
18 18 0 0 0 0 0 0 0 0
19 19 0 0 0 0 0 0 0 0
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
1 16 0 19 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
1.0 9281.9 11619.0 294.4 30.5 7.2 80.0 1107.7
0.0 147.2 172.9 6.7 11.9 5.4 104.1 1955.8
10000 10000 10000 10000 10000 10000 10000 10000
-----> Chanl_id Chanl_nb Pred_use(-1/0/1) Param
1 1 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0
3 3 0 0 0 0 0 0 0 0
4 4 0 0 0 0 0 0 0 0
5 5 0 0 0 0 0 0 0 0
6 6 0 0 0 0 0 0 0 0
7 7 0 0 0 0 0 0 0 0
8 8 0 0 0 0 0 0 0 0
9 9 0 0 0 0 0 0 0 0
10 10 0 0 0 0 0 0 0 0
11 11 0 0 0 0 0 0 0 0
12 12 0 0 0 0 0 0 0 0
13 13 0 0 0 0 0 0 0 0
14 14 0 0 0 0 0 0 0 0
15 15 0 0 0 0 0 0 0 0
16 16 0 0 0 0 0 0 0 0
17 17 0 0 0 0 0 0 0 0
18 18 0 0 0 0 0 0 0 0
19 19 0 0 0 0 0 0 0 0
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
1 17 0 19 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
1.0 9281.9 11619.0 294.4 30.5 7.2 80.0 1107.7
0.0 147.2 172.9 6.7 11.9 5.4 104.1 1955.8
10000 10000 10000 10000 10000 10000 10000 10000
-----> Chanl_id Chanl_nb Pred_use(-1/0/1) Param
1 1 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0
3 3 0 0 0 0 0 0 0 0
4 4 0 0 0 0 0 0 0 0
5 5 0 0 0 0 0 0 0 0
6 6 0 0 0 0 0 0 0 0
7 7 0 0 0 0 0 0 0 0
8 8 0 0 0 0 0 0 0 0
9 9 0 0 0 0 0 0 0 0
10 10 0 0 0 0 0 0 0 0
11 11 0 0 0 0 0 0 0 0
12 12 0 0 0 0 0 0 0 0
13 13 0 0 0 0 0 0 0 0
14 14 0 0 0 0 0 0 0 0
15 15 0 0 0 0 0 0 0 0
16 16 0 0 0 0 0 0 0 0
17 17 0 0 0 0 0 0 0 0
18 18 0 0 0 0 0 0 0 0
19 19 0 0 0 0 0 0 0 0
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
1 18 0 19 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
1.0 9281.9 11619.0 294.4 30.5 7.2 80.0 1107.7
0.0 147.2 172.9 6.7 11.9 5.4 104.1 1955.8
10000 10000 10000 10000 10000 10000 10000 10000
-----> Chanl_id Chanl_nb Pred_use(-1/0/1) Param
1 1 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0
3 3 0 0 0 0 0 0 0 0
4 4 0 0 0 0 0 0 0 0
5 5 0 0 0 0 0 0 0 0
6 6 0 0 0 0 0 0 0 0
7 7 0 0 0 0 0 0 0 0
8 8 0 0 0 0 0 0 0 0
9 9 0 0 0 0 0 0 0 0
10 10 0 0 0 0 0 0 0 0
11 11 0 0 0 0 0 0 0 0
12 12 0 0 0 0 0 0 0 0
13 13 0 0 0 0 0 0 0 0
14 14 0 0 0 0 0 0 0 0
15 15 0 0 0 0 0 0 0 0
16 16 0 0 0 0 0 0 0 0
17 17 0 0 0 0 0 0 0 0
18 18 0 0 0 0 0 0 0 0
19 19 0 0 0 0 0 0 0 0
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
1 19 0 19 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
1.0 9281.9 11619.0 294.4 30.5 7.2 80.0 1107.7
0.0 147.2 172.9 6.7 11.9 5.4 104.1 1955.8
10000 10000 10000 10000 10000 10000 10000 10000
-----> Chanl_id Chanl_nb Pred_use(-1/0/1) Param
1 1 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0
3 3 0 0 0 0 0 0 0 0
4 4 0 0 0 0 0 0 0 0
5 5 0 0 0 0 0 0 0 0
6 6 0 0 0 0 0 0 0 0
7 7 0 0 0 0 0 0 0 0
8 8 0 0 0 0 0 0 0 0
9 9 0 0 0 0 0 0 0 0
10 10 0 0 0 0 0 0 0 0
11 11 0 0 0 0 0 0 0 0
12 12 0 0 0 0 0 0 0 0
13 13 0 0 0 0 0 0 0 0
14 14 0 0 0 0 0 0 0 0
15 15 0 0 0 0 0 0 0 0
16 16 0 0 0 0 0 0 0 0
17 17 0 0 0 0 0 0 0 0
18 18 0 0 0 0 0 0 0 0
19 19 0 0 0 0 0 0 0 0
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
Expand Down Expand Up @@ -2209,6 +2349,23 @@
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
12 4 21 8 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
10000 10000 10000 10000 10000 10000 10000 10000 10000
-----> Chanl_id Chanl_nb Pred_use(-1/0/1) Param
1 1 0 0 0 0 0 0 0 0 4.457 -0.211 -0.054 0.213 -0.101 1.770 -0.283 -0.794
2 2 0 0 0 0 0 0 0 0 4.457 -0.211 -0.054 0.213 -0.101 1.770 -0.283 -0.794
3 3 0 0 0 0 0 0 0 0 4.457 -0.211 -0.054 0.213 -0.101 1.770 -0.283 -0.794
4 4 0 0 0 0 0 0 0 0 4.457 -0.211 -0.054 0.213 -0.101 1.770 -0.283 -0.794
5 5 0 0 0 0 0 0 0 0 4.457 -0.211 -0.054 0.213 -0.101 1.770 -0.283 -0.794
6 6 0 0 0 0 0 0 0 0 1.195 0.009 -0.048 0.082 -0.055 -1.055 0.894 -0.481
7 7 0 0 0 0 0 0 0 0 1.653 0.037 -0.072 0.066 -0.097 -1.610 0.345 0.302
8 8 0 0 0 0 0 0 0 0 0.582 0.253 -0.041 -0.042 -0.031 -0.803 0.227 0.076
------------------------------------------------
Platform_id Sat_id Sensor_id Nchanl Npredmax
------------------------------------------------
29 1 63 14 8
-----> Bias predictor statistics: Mean & Std & Nbgerr
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
Expand Down
9 changes: 9 additions & 0 deletions var/run/radiance_info/msg-4-seviri.info
View file
@@ -0,0 +1,9 @@
sensor channel IR/MW use idum varch polarisation(0:vertical;1:horizontal)
207 1 1 -1 0 1.8000000000E+00 2.0000000000E+00 msg-4-seviri
207 2 1 1 0 2.5000000000E+00 2.0000000000E+00
207 3 1 1 0 2.2500000000E+00 2.0000000000E+00
207 4 1 1 0 1.2500000000E+00 2.0000000000E+00
207 5 1 -1 0 1.2500000000E+00 2.0000000000E+00
207 6 1 -1 0 1.2500000000E+00 2.0000000000E+00
207 7 1 1 0 1.4500000000E+00 2.0000000000E+00
207 8 1 -1 0 1.2500000000E+00 3.0000000000E+00
20 changes: 20 additions & 0 deletions var/run/radiance_info/noaa-19-hirs.info
View file
@@ -0,0 +1,20 @@
sensor channel IR/MW use idum varch polarisation(0:vertical;1:horizontal)
18 1 0 -1 19 0.1000000000E+05 0.0000000000E+00 noaa-19-hirs
18 2 0 1 19 0.6000000000E+00 0.0000000000E+00
18 3 0 1 19 0.5300000000E+00 0.0000000000E+00
18 4 0 1 19 0.4000000000E+00 0.0000000000E+00
18 5 0 1 19 0.3600000000E+00 0.0000000000E+00
18 6 0 1 19 0.4600000000E+00 0.0000000000E+00
18 7 0 1 19 0.5700000000E+00 0.0000000000E+00
18 8 0 1 15 0.1000000000E+01 0.0000000000E+00
18 9 0 1 12 0.1100000000E+01 0.0000000000E+00
18 10 0 1 18 0.6000000000E+00 0.0000000000E+00
18 11 0 1 8 0.1200000000E+01 0.0000000000E+00
18 12 0 1 8 0.1600000000E+01 0.0000000000E+00
18 13 0 1 7 0.3640000000E+00 0.0000000000E+00
18 14 0 1 7 0.2600000000E+00 0.0000000000E+00
18 15 0 1 7 0.2600000000E+00 0.0000000000E+00
18 16 0 -1 7 0.1000000000E+05 0.0000000000E+00
18 17 0 -1 7 0.1000000000E+05 0.0000000000E+00
18 18 0 -1 6 0.1000000000E+05 0.0000000000E+00
18 19 0 -1 4 0.1000000000E+05 0.0000000000E+00

0 comments on commit e288837

Please sign in to comment.