[Dart-dev] DART/branches Revision: 13181

[dart at ucar.edu]

thoar at ucar.edu
2019-06-12 19:26:58 -0600 (Wed, 12 Jun 2019)
402
Converter extended to add 'surface_temperature', which then turned out to
not be the temperature of the soil surface but some 'temperature' that is
useful for the RTM. bugger. The convert_many.sh script is modified to 
faithfully convert daily files using the real year of collection. The previous
version shifted the data to some arbitrary year (for synthetic experiments
using a realistic network).




Modified: DART/branches/cesm_clm/observations/obs_converters/NSIDC/SMAP_L2_to_obs.f90
===================================================================
--- DART/branches/cesm_clm/observations/obs_converters/NSIDC/SMAP_L2_to_obs.f90	2019-06-06 20:05:12 UTC (rev 13180)
+++ DART/branches/cesm_clm/observations/obs_converters/NSIDC/SMAP_L2_to_obs.f90	2019-06-13 01:26:58 UTC (rev 13181)
@@ -1,8 +1,8 @@
-! DART software - Copyright 2004 - 2013 UCAR. This open source software is
-! provided by UCAR, "as is", without charge, subject to all terms of use at
+! DART software - Copyright UCAR. This open source software is provided
+! by UCAR, "as is", without charge, subject to all terms of use at
 ! http://www.image.ucar.edu/DAReS/DART/DART_download
 !
-! $Id$
+! DART $Id$
 
 program SMAP_L2_to_obs
 
@@ -14,9 +14,9 @@
 ! https://nsidc.org/data/ease/tools
 ! https://nsidc.org/data/SPL2SMP/versions/4
 !
-! "This Level-2 (L2) soil moisture product provides estimates of global land 
-! surface conditions retrieved by the Soil Moisture Active Passive (SMAP) 
-! passive microwave radiometer during 6:00 a.m. descending and 6:00 p.m. 
+! "This Level-2 (L2) soil moisture product provides estimates of global land
+! surface conditions retrieved by the Soil Moisture Active Passive (SMAP)
+! passive microwave radiometer during 6:00 a.m. descending and 6:00 p.m.
 ! ascending half-orbit passes. SMAP L-band brightness temperatures are used  to
 ! derive soil moisture data, which are then resampled to an Earth-fixed, global,
 ! cylindrical 36 km Equal-Area Scalable Earth Grid, Version 2.0 (EASE-Grid 2.0).
@@ -25,9 +25,9 @@
 ! SMAP L2 Radiometer Half-Orbit 36 km EASE-Grid Soil Moisture, Version 4
 !
 ! "Surface soil moisture (0-5 cm) in m3/m3 derived from brightness temperatures
-! (TBs) is output on a fixed global 36 km EASE-Grid 2.0. Also included are 
-! brightness temperatures in kelvin representing the weighted average of 
-! Level-1B brightness temperatures whose boresights fall within a 36 km 
+! (TBs) is output on a fixed global 36 km EASE-Grid 2.0. Also included are
+! brightness temperatures in kelvin representing the weighted average of
+! Level-1B brightness temperatures whose boresights fall within a 36 km
 ! EASE-Grid 2.0 cell."
 !
 !=======================================================================
@@ -49,7 +49,7 @@
                                print_time, print_date
 
 use   obs_sequence_mod, only : obs_sequence_type, obs_type, read_obs_seq, &
-                               static_init_obs_sequence, init_obs, write_obs_seq, & 
+                               static_init_obs_sequence, init_obs, write_obs_seq, &
                                init_obs_sequence, get_num_obs, &
                                set_copy_meta_data, set_qc_meta_data
 
@@ -64,7 +64,7 @@
                                h5_get_dimensions, &
                                h5_get_dset_dspace
 
-use       obs_kind_mod, only : SOIL_MOISTURE
+use       obs_kind_mod, only : SOIL_MOISTURE, SOIL_TEMPERATURE
 
 use HDF5
 use H5LT
@@ -102,28 +102,43 @@
 
 integer :: oday, osec, iocode, iunit
 integer :: num_copies, num_qc, max_obs
-           
+
 logical :: first_obs
+logical :: moisture_good, temperature_good
 
-! The EASE grid 
+! The EASE grid
 real(r4),        allocatable, dimension(:) :: longitude
 real(r4),        allocatable, dimension(:) :: latitude
-real(r4),        allocatable, dimension(:) :: observation
-real(r4),        allocatable, dimension(:) :: soil_moisture_error_std
+real(r4),        allocatable, dimension(:) :: moisture
+real(r4),        allocatable, dimension(:) :: moisture_error_std
+real(r4),        allocatable, dimension(:) :: temperature
 type(time_type), allocatable, dimension(:) :: obs_time
 integer,         allocatable, dimension(:) :: retrieval_flag
 
 integer :: icount
-integer :: counts = 20000
+integer :: counts = 40000
 
+integer :: moisture_fill_count
+integer :: moisture_min_count
+integer :: moisture_max_count
+integer :: moisture_count
+integer :: temperature_fill_count
+integer :: temperature_min_count
+integer :: temperature_max_count
+integer :: temperature_count
+integer :: retrieval_flag_count
+integer :: badlat_count
+integer :: badlon_count
+
 character(len=256) :: varname ! HDF variable names can be long
 integer  :: ndims
 
 real(r8) :: qc, obs_val, err_std