enforce time range for ASCAT and SMAP obs reader

Author: amfox37

GEOSlandassim_GridComp/clsm_ensupd_read_obs.F90

@@ -1629,6 +1629,7 @@ subroutine read_obs_sm_ASCAT_EUMET(                                         &
     type(date_time_type) :: date_time_tmp
     type(date_time_type) :: date_time_low, date_time_low_fname
     type(date_time_type) :: date_time_up
+    type(date_time_type) :: date_time_first_obs, date_time_last_obs
 
     integer :: ii, ind, N_tmp, N_files, kk, N_obs, N_fnames, N_fnames_tmp, obs_dir_hier
 
@@ -1688,6 +1689,38 @@ subroutine read_obs_sm_ASCAT_EUMET(                                         &
     nullify( tmp_obs, tmp_lat, tmp_lon, tmp_tile_num, tmp_jtime )    
 
     ! ---------------
+
+    ! check if we are within time periods when observations are available for metop-a, metop-b, or metop-c
+    ! if not, then return
+
+    if (trim(this_obs_param%descr)     == 'ASCAT_META_SM') then
+       date_time_first_obs = date_time_type(2007, 6, 1, 1,31, 0, 0, 0)
+       date_time_last_obs  = date_time_type(2021,11,15, 9, 0, 0, 0, 0)
+    elseif (trim(this_obs_param%descr) == 'ASCAT_METB_SM') then
+       date_time_first_obs = date_time_type(2013, 4,24, 8, 0, 0, 0, 0)
+       date_time_last_obs  = date_time_type(2100, 1, 1, 0, 0, 0, 0, 0)
+    elseif (trim(this_obs_param%descr) == 'ASCAT_METC_SM') then
+       date_time_first_obs = date_time_type(2019,11,25,12, 0, 0, 0, 0)
+       date_time_last_obs  = date_time_type(2100, 1, 1, 0, 0, 0, 0, 0)
+    else
+       err_msg = 'Unknown observation time range for this ASCAT observation type'
+       call ldas_abort(LDAS_GENERIC_ERROR, Iam, err_msg)
+    end if
+
+    if (datetime_lt_refdatetime(date_time, date_time_first_obs) .or.                        &
+       datetime_lt_refdatetime(date_time_last_obs, date_time)) then
+       found_obs    = .false.
+       ASCAT_sm     = this_obs_param%nodata
+       ASCAT_lon    = this_obs_param%nodata
+       ASCAT_lat    = this_obs_param%nodata
+       ASCAT_time   = real(this_obs_param%nodata,kind(0.0D0))
+       ASCAT_sm_std = this_obs_param%nodata
+
+       err_msg = 'Looking for ' // trim(this_obs_param%descr) // ' observations before/after they are available'
+       call ldas_warn(LDAS_GENERIC_WARNING, Iam, err_msg)
+
+       return
+    end if
 
     ! initialize
 
@@ -6849,6 +6882,7 @@ subroutine read_obs_SMAP_halforbit_Tb( date_time, N_catd, this_obs_param,  &
 
     type(date_time_type) :: date_time_low,       date_time_upp
     type(date_time_type) :: date_time_low_fname, date_time_tmp
+    type(date_time_type) :: date_time_first_obs
 
     integer              :: ii, jj, kk, nn, mm
     integer              :: N_fnames, N_fnames_tmp, N_obs_tmp
@@ -6900,6 +6934,24 @@ subroutine read_obs_SMAP_halforbit_Tb( date_time, N_catd, this_obs_param,  &
 
     ! -------------------------------------------------------------------
 
+    ! check if we are in the time range of SMAP observations which start on 2015-03-31
+
+    date_time_first_obs = date_time_type(2015, 3,31, 0, 0, 0, 0, 0)
+
+    if (datetime_lt_refdatetime(date_time, date_time_first_obs)) then
+       found_obs = .false.
+       SMAP_data     =      this_obs_param%nodata
+       SMAP_lon      =      this_obs_param%nodata
+       SMAP_lat      =      this_obs_param%nodata
+       SMAP_time     = real(this_obs_param%nodata,kind(0.0D0))
+       std_SMAP_data =      this_obs_param%nodata
+
+       err_msg = 'Looking for SMAP observations before 2015-03-31'
+       call ldas_warn(LDAS_GENERIC_WARNING, Iam, err_msg)
+
+       return
+    end if
+
     ! check inputs
 
     ! the subroutine makes sense only if dtstep_assim <= 3 hours