From 8a29cf11037248fe7c0a75cdd7fa1e8dda242c9d Mon Sep 17 00:00:00 2001 From: Tracy Date: Wed, 8 Apr 2020 08:25:56 -0600 Subject: [PATCH 1/3] Add GOES 16/17 to GRIB2 table, update release tag in UG --- docs/UPP_GRIB2_Table.lyx | 1150 ++++++++++++++++++++++++++++++++- docs/UPP_Users_guide_V4.1.lyx | 2 +- 2 files changed, 1146 insertions(+), 6 deletions(-) diff --git a/docs/UPP_GRIB2_Table.lyx b/docs/UPP_GRIB2_Table.lyx index 448f90809..01b9308c0 100644 --- a/docs/UPP_GRIB2_Table.lyx +++ b/docs/UPP_GRIB2_Table.lyx @@ -116,7 +116,7 @@ postcntrl.xml \begin_layout Standard \begin_inset Tabular - + @@ -25169,7 +25169,7 @@ NA - + \begin_inset Text \begin_layout Plain Layout @@ -25180,7 +25180,7 @@ Insat 3d brightness temperature IR channels 1-4 \end_inset - + \begin_inset Text \begin_layout Plain Layout @@ -25191,7 +25191,7 @@ Not currently available for grib2 \end_inset - + \begin_inset Text \begin_layout Plain Layout @@ -25202,7 +25202,7 @@ NA \end_inset - + \begin_inset Text \begin_layout Plain Layout @@ -25222,6 +25222,1146 @@ NA 865 \end_layout +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 7 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA167_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA167 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +927 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 8 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA168_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA168 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +928 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 9 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA169_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA169 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +929 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 10 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1610_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1610 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +930 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 11 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1611_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1611 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +931 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 12 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1612_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1612 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +932 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 13 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1613_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1613 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +933 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 14 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1614_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1614 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +934 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 15 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1615_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1615 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +935 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-16 brightness temperature channel 16 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1616_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1616 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +936 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 7 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA177_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA177 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +937 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 8 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA178_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA178 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +938 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 9 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA179_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA179 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +939 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 10 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1710_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1710 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +940 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 11 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1711_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1711 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +941 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 12 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1712_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1712 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +942 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 13 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1713_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1713 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +943 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 14 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1714_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1714 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +944 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 15 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1715_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1715 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +945 +\end_layout + +\end_inset + + + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +GOES-17 brightness temperature channel 16 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1716_ON_TOP_OF_ATMOS +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +SBTA1716 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +109 +\end_layout + +\end_inset + + +\begin_inset Text + +\begin_layout Plain Layout + +\size footnotesize +946 +\end_layout + \end_inset diff --git a/docs/UPP_Users_guide_V4.1.lyx b/docs/UPP_Users_guide_V4.1.lyx index b304c5157..eb13a5f39 100644 --- a/docs/UPP_Users_guide_V4.1.lyx +++ b/docs/UPP_Users_guide_V4.1.lyx @@ -394,7 +394,7 @@ release-tag-name \end_inset is the release tag you wish to clone (e.g. - for UPPV4.1 use the release tag v4.1.0). + for UPPV4.1 use the release tag dtc_post_v4.1.0). \end_layout \end_deeper From 74be3967671935a886ccb54397dd88a3f70cf4a8 Mon Sep 17 00:00:00 2001 From: Tracy Date: Tue, 14 Apr 2020 16:29:25 -0600 Subject: [PATCH 2/3] Add Himawari-8 for grib2 --- parm/params_grib2_tbl_new | 10 ++ parm/post_avblflds.xml | 70 ++++++++++++++ sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f | 109 ++++++++++++++++------ sorc/ncep_post.fd/RQSTFLD.F | 33 ++++++- 4 files changed, 191 insertions(+), 31 deletions(-) diff --git a/parm/params_grib2_tbl_new b/parm/params_grib2_tbl_new index 81b86e105..6ea5fffd6 100755 --- a/parm/params_grib2_tbl_new +++ b/parm/params_grib2_tbl_new @@ -809,6 +809,16 @@ 3 192 52 0 SBTAGR7 3 192 53 0 SBTAGR8 3 192 54 0 SBTAGR9 + 3 192 65 0 SBTAHI10 + 3 192 66 0 SBTAHI11 + 3 192 67 0 SBTAHI12 + 3 192 68 0 SBTAHI13 + 3 192 69 0 SBTAHI14 + 3 192 70 0 SBTAHI15 + 3 192 71 0 SBTAHI16 + 3 192 62 0 SBTAHI7 + 3 192 63 0 SBTAHI8 + 3 192 64 0 SBTAHI9 3 0 2 0 SBTMP 4 2 10 0 SCINT 0 1 84 0 SCLIWC diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml index 339af0ac1..c50b7b723 100755 --- a/parm/post_avblflds.xml +++ b/parm/post_avblflds.xml @@ -6573,6 +6573,76 @@ 4.0 + + 912 + SBTAHI7_ON_TOP_OF_ATMOS + SBTAHI7 + top_of_atmos + 4.0 + + + 913 + SBTAHI8_ON_TOP_OF_ATMOS + SBTAHI8 + top_of_atmos + 4.0 + + + 914 + SBTAHI9_ON_TOP_OF_ATMOS + SBTAHI9 + top_of_atmos + 4.0 + + + 915 + SBTAHI10_ON_TOP_OF_ATMOS + SBTAHI10 + top_of_atmos + 4.0 + + + 916 + SBTAHI11_ON_TOP_OF_ATMOS + SBTAHI11 + top_of_atmos + 4.0 + + + 917 + SBTAHI12_ON_TOP_OF_ATMOS + SBTAHI12 + top_of_atmos + 4.0 + + + 918 + SBTAHI13_ON_TOP_OF_ATMOS + SBTAHI13 + top_of_atmos + 4.0 + + + 919 + SBTAHI14_ON_TOP_OF_ATMOS + SBTAHI14 + top_of_atmos + 4.0 + + + 920 + SBTAHI15_ON_TOP_OF_ATMOS + SBTAHI15 + top_of_atmos + 4.0 + + + 921 + SBTAHI16_ON_TOP_OF_ATMOS + SBTAHI16 + top_of_atmos + 4.0 + 927 SBTA167_ON_TOP_OF_ATMOS diff --git a/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f b/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f index 82b86813e..e01ba1b42 100644 --- a/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f +++ b/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f @@ -13,6 +13,7 @@ SUBROUTINE CALRAD_WCLOUD ! FUNCTION EFFR TO COMPUTE EFFECTIVE PARTICLE RADII ! CHANNEL SELECTION USING LVLS FROM WRF_CNTRL.PARM ! 19-04-01 Sharon NEBUDA - Added output option for GOES-16 & GOES-17 ABI IR Channels 7-16 + ! 20-04-09 Tracy Hertneky - Added Himawari-8 AHI CH7-CH16 ! ! USAGE: CALL MDLFLD ! INPUT ARGUMENT LIST: @@ -121,7 +122,7 @@ SUBROUTINE CALRAD_WCLOUD ! integer,parameter:: n_clouds = 4 integer,parameter:: n_aerosols = 0 ! Add your sensors here - integer(i_kind),parameter:: n_sensors=21 + integer(i_kind),parameter:: n_sensors=22 character(len=20),parameter,dimension(1:n_sensors):: sensorlist= & (/'imgr_g15 ', & 'imgr_g13 ', & @@ -143,29 +144,31 @@ SUBROUTINE CALRAD_WCLOUD 'imgr_insat3d ', & 'abi_gr ', & 'abi_g16 ', & - 'abi_g17 '/) - character(len=12),parameter,dimension(1:n_sensors):: obslist= & - (/'goes_img ', & - 'goes_img ', & - 'goes_img ', & - 'goes_img ', & - 'amsre ', & - 'tmi ', & - 'ssmi ', & - 'ssmi ', & - 'ssmi ', & - 'ssmis ', & - 'ssmis ', & - 'ssmis ', & - 'ssmis ', & - 'ssmis ', & - 'seviri ', & - 'imgr_mt2 ', & - 'imgr_mt1r ', & - 'imgr_insat3d', & - 'abi ', & - 'abi ', & - 'abi '/) + 'abi_g17 ', & + 'ahi_himawari8 '/) + character(len=13),parameter,dimension(1:n_sensors):: obslist= & + (/'goes_img ', & + 'goes_img ', & + 'goes_img ', & + 'goes_img ', & + 'amsre ', & + 'tmi ', & + 'ssmi ', & + 'ssmi ', & + 'ssmi ', & + 'ssmis ', & + 'ssmis ', & + 'ssmis ', & + 'ssmis ', & + 'ssmis ', & + 'seviri ', & + 'imgr_mt2 ', & + 'imgr_mt1r ', & + 'imgr_insat3d ', & + 'abi ', & + 'abi ', & + 'abi ', & + 'ahi_himawari8'/) character(len=20),dimension(1:n_sensors):: sensorlist_local ! integer(i_kind) sensorindex @@ -196,7 +199,7 @@ SUBROUTINE CALRAD_WCLOUD real,parameter:: constoz = 604229.0_r_kind real sublat,sublon real RHO,RHOX - character(12)::obstype + character(13)::obstype character(20)::isis character(20)::isis_local @@ -209,6 +212,7 @@ SUBROUTINE CALRAD_WCLOUD logical micrim,microwave logical post_abig16, post_abig17, post_abigr ! if true, user requested at least one abi channel logical fix_abig16, fix_abig17 ! if true, abi_g16, abi_g17 fix files are available + logical post_ahi8 ! if true, user requested at least on ahi channel (himawari8) ! logical,dimension(nobs):: luse logical, parameter :: debugprint = .false. type(crtm_atmosphere_type),dimension(1):: atmosphere @@ -286,13 +290,17 @@ SUBROUTINE CALRAD_WCLOUD do n = 958, 958+9 ! 958 set in RQSTFLD.f if (iget(n) > 0) post_abigr=.true. enddo + post_ahi8=.false. + do n = 912, 912+9 ! 912 set in RQSTFLD.f + if (iget(n) > 0) post_ahi8=.true. + enddo ! DO NOT FORGET TO ADD YOUR NEW IGET HERE (IF YOU'VE ADDED ONE) ! START SUBROUTINE CALRAD. ifactive: if (iget(327) > 0 .or. iget(328) > 0 .or. iget(329) > 0 & .or. iget(330) > 0 .or. iget(446) > 0 .or. iget(447) > 0 & - .or. iget(448) > 0 .or. iget(449) > 0 .or. iget(456) > 0 & + .or. iget(448) > 0 .or. iget(449) > 0 .or. iget(456) > 0 & .or. iget(457) > 0 .or. iget(458) > 0 .or. iget(459) > 0 & .or. iget(460) > 0 .or. iget(461) > 0 .or. iget(462) > 0 & .or. iget(463) > 0 .or. iget(483) > 0 .or. iget(484) > 0 & @@ -326,7 +334,8 @@ SUBROUTINE CALRAD_WCLOUD .or. iget(871) > 0 .or. iget(872) > 0 .or. iget(873) > 0 & .or. iget(874) > 0 .or. iget(875) > 0 .or. iget(876) > 0 & .or. iget(877) > 0 .or. iget(878) > 0 .or. iget(879) > 0 & - .or. iget(880) > 0 .or. iget(881) > 0 .or. iget(882) > 0 & + .or. iget(880) > 0 .or. iget(881) > 0 .or. iget(882) > 0 & + .or. post_ahi8 & .or. post_abig16 .or. post_abig17 .or. post_abigr ) then ! specify numbers of cloud species @@ -339,7 +348,6 @@ SUBROUTINE CALRAD_WCLOUD .or. imp_physics==28 .or. imp_physics==11)then n_clouds=5 else - n_clouds=0 print*,'Warning: number of cloud species (n_clouds) being set to zero for imp_physics=',imp_physics end if @@ -433,6 +441,20 @@ SUBROUTINE CALRAD_WCLOUD endif endif + ! Himawari-8 AHI infrared + if(post_ahi8)then + nchanl=0 + do n = 912, 912+9 ! 912 set in RQSTFLD.f + if (iget(n) > 0) then + nchanl = nchanl+1 + endif + enddo + if (nchanl > 0 .and. nchanl <10) then + do n = 912, 912+9 ! 912 set in RQSTFLD.f + if (iget(n) == 0) channelinfo(22)%Process_Channel(n-912+1)=.False. ! turn off channel processing + enddo + endif + endif ! SSMI, F13-F15 (19H,19V,??H,37H,37V,85H,85V) if(iget(800)>0)then @@ -473,10 +495,11 @@ SUBROUTINE CALRAD_WCLOUD call select_channels_L(channelinfo(17),4,(/ 1,2,3,4 /),lvls(1:4,iget(864)),iget(864)) endif ! INSAT 3D (Kalpana) - if(iget(864)>0)then + if(iget(865)>0)then call select_channels_L(channelinfo(18),4,(/ 1,2,3,4 /),lvls(1:4,iget(865)),iget(865)) endif + ! Loop over data types to process sensordo: do isat=1,n_sensors @@ -512,7 +535,8 @@ SUBROUTINE CALRAD_WCLOUD (isis=='abi_g16' .and. post_abig16) .OR. & (isis=='abi_g17' .and. post_abig17) .OR. & (isis=='abi_gr' .and. post_abigr) .OR. & - (isis=='seviri_m10' .and. iget(876)>0) )then + (isis=='seviri_m10' .and. iget(876)>0) .OR. & + (isis=='ahi_himawari8' .and. post_ahi8) )then print*,'obstype, isis= ',obstype,isis ! isis='amsua_n15' @@ -1231,6 +1255,7 @@ SUBROUTINE CALRAD_WCLOUD (isis=='abi_g16' .and. post_abig16) .OR. & (isis=='abi_g17' .and. post_abig17) .OR. & (isis=='seviri_m10' .and. iget(876)>0) .OR. & + (isis=='ahi_himawari8' .and. post_ahi8) .OR. & (isis=='imgr_g12' .and. (iget(456)>0 .or. & iget(457)>0 .or. iget(458)>0 .or. iget(459)>0)) .or. & (isis=='imgr_g11' .and. (iget(460)>0 .or. & @@ -1271,6 +1296,9 @@ SUBROUTINE CALRAD_WCLOUD else if(isis=='imgr_insat3d') then sublat=0.0 sublon=74.0 + else if(isis=='ahi_himawari8') then + sublat=0.0 + sublon=140.7 end if ! use zenith angle = 53.1 for SSMI and SSMIS: @@ -2092,6 +2120,27 @@ SUBROUTINE CALRAD_WCLOUD endif enddo ! channel loop end if ! end of outputting goes 17 + if(isis=='ahi_himawari8') then ! writing Himawari-8 AHI to grib + do ichan=1,10 + igot=iget(912+ichan) + if(igot>0)then + do j=jsta,jend + do i=1,im + grid1(i,j)=tb(i,j,ichan) + enddo + enddo + id(1:25) = 0 + id(02) = 2 + id(08) = 118 + id(09) = 109 + if(grib=="grib2" )then + cfld=cfld+1 + fld_info(cfld)%ifld=IAVBLFLD(igot) + datapd(1:im,1:jend-jsta+1,cfld)=grid1(1:im,jsta:jend) + endif + endif + enddo + endif ! end of outputting himawari-8 ahi end if nonnadir ! end if for computing simulated radiance with zenith angle correction diff --git a/sorc/ncep_post.fd/RQSTFLD.F b/sorc/ncep_post.fd/RQSTFLD.F index 50410d134..8cb6fc668 100644 --- a/sorc/ncep_post.fd/RQSTFLD.F +++ b/sorc/ncep_post.fd/RQSTFLD.F @@ -2465,7 +2465,38 @@ module RQSTFLD_mod DATA IFILV(911),AVBL(911),IQ(911),IS(911),AVBLGRB2(911) & & /1,'VTEMP AT FD HEIGHTS ',012,103, & & 'VTMP ON spec_alt_above_mean_sea_lvl'/ - +! LVLS-DRIVEN CHANNEL SELECTION (L(1)-L(10)> CRTM CH1-CH10 > AHI CH7-16) +! HWRF satellite additions: Himawari-8, replacement for MTSAT-2 + DATA IFILV(912),AVBL(912),IQ(912),IS(912),AVBLGRB2(912) & + & /1,'HIM-8 AHI CH7 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH7 NON-NADIR BRTEMP'/ + DATA IFILV(913),AVBL(913),IQ(913),IS(913),AVBLGRB2(913) & + & /1,'HIM-8 AHI CH8 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH8 NON-NADIR BRTEMP'/ + DATA IFILV(914),AVBL(914),IQ(914),IS(914),AVBLGRB2(914) & + & /1,'HIM-8 AHI CH9 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH9 NON-NADIR BRTEMP'/ + DATA IFILV(915),AVBL(915),IQ(915),IS(915),AVBLGRB2(915) & + & /1,'HIM-8 AHI CH10 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH10 NON-NADIR BRTEMP'/ + DATA IFILV(916),AVBL(916),IQ(916),IS(916),AVBLGRB2(916) & + & /1,'HIM-8 AHI CH11 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH11 NON-NADIR BRTEMP'/ + DATA IFILV(917),AVBL(917),IQ(917),IS(917),AVBLGRB2(917) & + & /1,'HIM-8 AHI CH12 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH12 NON-NADIR BRTEMP'/ + DATA IFILV(918),AVBL(918),IQ(918),IS(918),AVBLGRB2(918) & + & /1,'HIM-8 AHI CH13 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH13 NON-NADIR BRTEMP'/ + DATA IFILV(919),AVBL(919),IQ(919),IS(919),AVBLGRB2(919) & + & /1,'HIM-8 AHI CH14 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH14 NON-NADIR BRTEMP'/ + DATA IFILV(920),AVBL(920),IQ(920),IS(920),AVBLGRB2(920) & + & /1,'HIM-8 AHI CH15 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH15 NON-NADIR BRTEMP'/ + DATA IFILV(921),AVBL(921),IQ(921),IS(921),AVBLGRB2(921) & + & /1,'HIM-8 AHI CH16 NNADR',118,109, & !table 130 + & 'HIM-8 AHI CH16 NON-NADIR BRTEMP'/ ! Add GOES-16 & GOES-17 ABI IR Channels 7-16 ! Grib2 defines each channel as a separate output field DATA IFILV(927),AVBL(927),IQ(927),IS(927),AVBLGRB2(927) & From a5dbe9f57cf948819f259bef038eaf602292df6a Mon Sep 17 00:00:00 2001 From: Tracy Date: Tue, 14 Apr 2020 16:46:05 -0600 Subject: [PATCH 3/3] Link Himawari-8 coefficient files in run_unipost script --- scripts/run_unipost | 2 ++ 1 file changed, 2 insertions(+) diff --git a/scripts/run_unipost b/scripts/run_unipost index 765cacccb..a755e2b59 100755 --- a/scripts/run_unipost +++ b/scripts/run_unipost @@ -279,6 +279,8 @@ ln -fs $CRTMDIR/TauCoeff/ODPS/Big_Endian/seviri_m10.TauCoeff.bin ./ ln -fs $CRTMDIR/SpcCoeff/Big_Endian/v.seviri_m10.SpcCoeff.bin ./ ln -fs $CRTMDIR/TauCoeff/ODPS/Big_Endian/abi_gr.TauCoeff.bin ./ ln -fs $CRTMDIR/SpcCoeff/Big_Endian/abi_gr.SpcCoeff.bin ./ +ln -fs $CRTMDIR/TauCoeff/ODPS/Big_Endian/ahi_himawari8.TauCoeff.bin ./ +ln -fs $CRTMDIR/SpcCoeff/Big_Endian/ahi_himawari8.SpcCoeff.bin ./ ####################################################### # 1. Run Unipost