ERROR with CRTM_K_Matrix

Post Reply
sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

ERROR with CRTM_K_Matrix

Post by sekluzia » Fri Feb 22, 2019 12:58 pm

Dear Colleges,

I am failing to assimilate the AMSU-A radiance data.

Could you please look the errors in my log.wrfda file and check my namelist.input file to identify the problem(s)?
You can see from my log.wrfda file that 595 radiance observations are available over my simulation domain.
I am using WRFDA v4.0.


It seems my errors are related to the crtm_coeffs files, but I uploaded those from http://www2.mmm.ucar.edu/wrf/users/wrfd ... 2.3.tar.gz as recommended.

In my simulations I use the LU_INDEX/LANDUSEF fields interpolated from 21-class MODIS.

I would much appreciate more helpful comments on my namelist file to improve the radiance data assimilation.

Artur
Attachments
log.wrfda.txt
(7.21 KiB) Downloaded 61 times
namelist.input.txt
(1.7 KiB) Downloaded 60 times

JakeLiu
Posts: 14
Joined: Fri Jun 08, 2018 8:29 am

Re: ERROR with CRTM_K_Matrix

Post by JakeLiu » Mon Feb 25, 2019 6:18 am

This looks like related to your use of 21-class MODIS Landuse.
WRFDA may need some code update to use it. Did you try the older one like USGS?

sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

Re: ERROR with CRTM_K_Matrix

Post by sekluzia » Mon Feb 25, 2019 6:29 am

Dear Jake,

Thanks for your reply! Indeed, I also suspected that the use of 21-class MODIS Landuse could be the reason for this issue. However, my experience with the WRF simulations over our region shows that it is better to use 21-class MODIS Landuse rather than the USGS. On the other hand, the older version of WRFDA (Versions 3.8.1) successfully assimilated satellite radiance with 21-class MODIS Landuse.
Could this issue be solved without changing the landuse scheme?

Artur

JakeLiu
Posts: 14
Joined: Fri Jun 08, 2018 8:29 am

Re: ERROR with CRTM_K_Matrix

Post by JakeLiu » Mon Feb 25, 2019 6:40 am

Ok. Will check.

JJGuerrette
Posts: 19
Joined: Mon Feb 25, 2019 8:14 pm

Re: ERROR with CRTM_K_Matrix

Post by JJGuerrette » Mon Feb 25, 2019 8:28 pm

Dear Artur,

Have you tried using the WRF/WRFDA code in the develop branch on github (https://github.com/wrf-model/WRF/tree/develop)? WRFDA V4.0 still uses CRTM v2.2.3. It looks like you are using CRTM coefficient files for CRTM V2.3.0. The latest WRFDA code about to be released, and available in that develop branch, is compatible with CRTM V2.3.0. You could try the soon to be released version 4.1 branch instead (https://github.com/wrf-model/WRF/tree/release-v4.1), which may be more stable.

EDIT: Although looking at the CRTM code itself, the problem appears to be with the Land_Type used as input to CRTM. The error you report is only given when Land_Type (derived from IVGTYP in the fg file) is set to 0 in the nearest grid cell to a particular pixel. Therefore does your fg file have any zeros in IVGTYP?

JJ

sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

Re: ERROR with CRTM_K_Matrix

Post by sekluzia » Tue Feb 26, 2019 6:07 am

Hi,

Thanks for your reply and help!
You are right I edited my Registry.EM_COMMON file trying to reduce the wrfout* files size (viewtopic.php?f=40&t=629), I also excluded the IVGTYP parameter. I will recompile the WRF model with the default Registry.EM_COMMON file and try again to assimilate satellite radiance.

Artur

sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

Re: ERROR with CRTM_K_Matrix

Post by sekluzia » Tue Feb 26, 2019 5:10 pm

Thank you so much for your help! After the presence of the IVGTYP data in my fg file it works.

Artur

JJGuerrette
Posts: 19
Joined: Mon Feb 25, 2019 8:14 pm

Re: ERROR with CRTM_K_Matrix

Post by JJGuerrette » Fri Mar 01, 2019 8:46 pm

Glad we could help!

sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

Re: ERROR with CRTM_K_Matrix

Post by sekluzia » Tue Mar 05, 2019 12:46 pm

Hi,

I have the same issue but now with assimilation of HIRS4, NOAA 18 data.

I have already updated the WRFDA version to v4.1 and I am using the CRTM V2.3.0 as you recommended. Now, I have all default fields in my "fg" file as I mentioned in my previous post (the file is attached). However, the error message is the following:
---------------------------- FATAL ERROR -----------------------
Fatal error in file: da_get_innov_vector_crtm.inc LINE: 210
modify crtm_irland_coef to use the coef file that is consistent with your model land use type
----------------------------------------------------------------
Note: The following floating-point exceptions are signalling: IEEE_DENORMAL
STOP wrf_abort
You can find attached the entire information of WRFDA run in log.out file. The hirs4.bufr, LANDUSE.TBL and namelist.input files are also attached.

I think I need to make some modifications in the da_get_innov_vector_crtm.inc file. I am not sure, but it seems that this issue is due to use of Modis landuse instead of USGS.
Could you help me?

Artur
Attachments
LANDUSE.TBL.txt
(29.12 KiB) Downloaded 50 times
namelist.input.txt
(1.94 KiB) Downloaded 56 times
hirs4.tar.gz
(24.72 MiB) Downloaded 56 times
fg.tar.gz
(91.18 MiB) Downloaded 55 times
log.out.txt
(6.52 KiB) Downloaded 55 times

JJGuerrette
Posts: 19
Joined: Mon Feb 25, 2019 8:14 pm

Re: ERROR with CRTM_K_Matrix

Post by JJGuerrette » Thu Mar 07, 2019 12:15 am

Which coefficient file are you using for IR over land in CRTM? This should be linked/copied into your run directory. The landuse categories associated with that file need to match with the num_land_cat setting in the namelist and with the fg file (and corresponding met_em* files from WPS).

USGS
WPS namelist:
geog_data_res = 'usgs_30s+default','usgs_30s+default',

WRF namelist:
num_land_cat = 24

CRTM file:
USGS.IRland.EmisCoeff.bin

MODIS/IGBP

WPS namelist:
geog_data_res = 'default','default',

WRF namelist:
num_land_cat = 21

CRTM file:
IGBP.IRland.EmisCoeff.bin

JJGuerrette
Posts: 19
Joined: Mon Feb 25, 2019 8:14 pm

Re: ERROR with CRTM_K_Matrix

Post by JJGuerrette » Thu Mar 07, 2019 12:19 am

I forgot to mention that the namelist option "crtm_irland_coef" also needs to match the coefficient file you are using. That appears to be the cause for your error.

sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

Re: ERROR with CRTM_K_Matrix

Post by sekluzia » Thu Mar 07, 2019 5:24 am

Thank you so much!

The
&wrfvar14

crtm_irland_coef='IGBP.IRland.EmisCoeff.bin'
works!

Artur

Rekha_Bharali1
Posts: 4
Joined: Tue Nov 27, 2018 12:09 pm

Re: ERROR with CRTM_K_Matrix

Post by Rekha_Bharali1 » Thu Mar 14, 2019 10:32 am

Hi,

I am also facing the same "CRTM_K_Matrix(FAILURE) : Input data check failed for profile #1" issue. The error massage is:

CRTM_LandSurface_IsValid(INFORMATION) : Invalid Land Surface data
CRTM_K_Matrix(FAILURE) : Input data check failed for profile #1
---------------------------- FATAL ERROR -----------------------
Fatal error in file: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/ ... m_k.inc</a> LINE: 38
Error in calling CRTM_K_Matrix
----------------------------------------------------------------
wrf_abort


I have assimilated AMSUA data successfully. But when I tried to assimilate MHS data I got the said issue. I have tried with both USGS and IGBP Landuse data. However, the errors are same. I am using CRTM_2.1.3 available with WRFDA_V3.8. I have attached the namelist and log files for your reference. Kindly help me to solve the issue.
Attachments
namelist.input_wrfda.txt
(5.18 KiB) Downloaded 61 times
log.txt
(10.83 KiB) Downloaded 64 times

Rekha_Bharali1
Posts: 4
Joined: Tue Nov 27, 2018 12:09 pm

Re: ERROR with CRTM_K_Matrix

Post by Rekha_Bharali1 » Mon Mar 18, 2019 11:55 am

As of now , there is no reply. I hope somebody will help me. I have further checked and found that the error that I have mentioned in my previous post is with only NOAA-19 whether its is MHS or AMSUA, it really doesn't matter.

sekluzia
Posts: 187
Joined: Mon Oct 15, 2018 12:42 pm

Re: ERROR with CRTM_K_Matrix

Post by sekluzia » Mon Mar 18, 2019 12:26 pm

Hi,

The following record is not present in your namelist.input file:
&wrfvar14

crtm_irland_coef='IGBP.IRland.EmisCoeff.bin'

If you are using MODIS landuse in your WRF simulations, add this in your namelist and try to run again.

Rekha_Bharali1
Posts: 4
Joined: Tue Nov 27, 2018 12:09 pm

Re: ERROR with CRTM_K_Matrix

Post by Rekha_Bharali1 » Wed Mar 20, 2019 4:21 am

Thank you for your reply. I have tried with both MODIS and USGS LULC data. As per your suggestion, I have put the line in my namelist file and run again. However, the error is the same. I could run for NOAA-15-18 successfully for both AMSU and MHS. But the error I got is only with NOAA-19. Whenever I tried to use options such as (1-19-3, 1-19-15), the error arises.

JJGuerrette
Posts: 19
Joined: Mon Feb 25, 2019 8:14 pm

Re: ERROR with CRTM_K_Matrix

Post by JJGuerrette » Mon Jun 03, 2019 4:37 pm

Please consider the following lines in var/external/crtm_2.3.0/libsrc/CRTM_Surface_Define.f90:

Code: Select all

    IF ( Sfc%Land_Temperature      < ZERO .OR. &
         Sfc%Soil_Moisture_Content < ZERO .OR. &
         Sfc%Canopy_Water_Content  < ZERO .OR. &
         Sfc%Vegetation_Fraction   < ZERO .OR. &
         Sfc%Soil_Temperature      < ZERO .OR. &
         Sfc%LAI                   < ZERO      ) THEN
      msg = 'Invalid Land Surface data'
      CALL Display_Message( ROUTINE_NAME, TRIM(msg), INFORMATION )
      IsValid = .FALSE.
    ENDIF
This is a slightly different error than the one described above, which was caused by the land type only. Here there are several possibilities. All of these quantities are set in var/da/da_radiance/da_get_innov_vector_crtm.inc, and can then be traced back to WRF variables stored in the "grid" as described in the Registry directory. These need to be checked one-by-one to determine the cause of the error.

hclin2018
Posts: 31
Joined: Wed Jun 13, 2018 5:38 pm

Re: ERROR with CRTM_K_Matrix

Post by hclin2018 » Thu Jul 25, 2019 10:13 pm

WRFDA only handles the land type conversion from WRF to CRTM for
num_land_cat=24 (USGS no lake) and
num_land_cat=20 (MODIS no lake)

Post Reply

Return to “Radiance DA”