[Dart-dev] DART/branches Revision: 12747

[dart at ucar.edu]

nancy at ucar.edu
2018-07-13 10:18:44 -0600 (Fri, 13 Jul 2018)
1017
interface to cloud cleared level 2 radiances - granule info
from example data reader.  obs converter now making obs_seq
files  with metadata for satellite/solar location, and 
platform/sat_id/sensor/channel.  can read/write them now.
created a test input obs_seq from a downloaded CC_RAD hdf file.

forward operator code is building temp/pres/moisture profiles
on model levels.

now trying to fold the example program code into an rttov_interface
module so it can be called by the forward operator and keep the
rttov calls isolated as much as possible.  still lots of questions
about surface fields and some of the other inputs the rt model
code needs.  this file needs to be preprocessed.  that's tricky
on the macs where the filesystems ignore case (which is insane).
so bob.f90 and bob.F90 can't be distinguished.

these .f90 filenames need to be revamped eventually to be clear.
we still have the existing airs retrieval converter and right
now the names are confusing.  do this later when more code is
working.




Modified: DART/branches/radiance/observations/forward_operators/obs_def_rttov_mod.f90
===================================================================
--- DART/branches/radiance/observations/forward_operators/obs_def_rttov_mod.f90	2018-07-13 16:07:06 UTC (rev 12746)
+++ DART/branches/radiance/observations/forward_operators/obs_def_rttov_mod.f90	2018-07-13 16:18:44 UTC (rev 12747)
@@ -15,6 +15,7 @@
 !    sat az/el
 !    sun az/el
 !    platform
+!    sat_id
 !    instrument
 !    channel
 !    <anything else useful>
@@ -21,7 +22,7 @@
 !----------------------------------------------------------------------
 
 ! BEGIN DART PREPROCESS KIND LIST
-! AIRS_AMSU_RADIANCE,    QTY_RADIANCE
+! AQUA_AIRS_AMSU_RADIANCE,    QTY_INFRARED_RADIANCE
 ! END DART PREPROCESS KIND LIST
 
 
@@ -34,25 +35,25 @@
 
 
 ! BEGIN DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
-!      case(AIRS_AMSU_RADIANCE)
+!      case(AQUA_AIRS_AMSU_RADIANCE)
 !         call get_expected_radiance(state_handle, ens_size, location, obs_def%key, expected_obs, istatus)
 ! END DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
 
 
 ! BEGIN DART PREPROCESS READ_OBS_DEF
-!   case(AIRS_AMSU_RADIANCE)
+!   case(AQUA_AIRS_AMSU_RADIANCE)
 !      call read_rttov_metadata(obs_def%key, key, ifile, fform)
 ! END DART PREPROCESS READ_OBS_DEF
 
 
 ! BEGIN DART PREPROCESS WRITE_OBS_DEF
-!   case(AIRS_AMSU_RADIANCE)
+!   case(AQUA_AIRS_AMSU_RADIANCE)
 !      call write_rttov_metadata(obs_def%key, ifile, fform)
 ! END DART PREPROCESS WRITE_OBS_DEF
 
 
 ! BEGIN DART PREPROCESS INTERACTIVE_OBS_DEF
-!   case(AIRS_AMSU_RADIANCE)
+!   case(AQUA_AIRS_AMSU_RADIANCE)
 !      call interactive_rttov_metadata(obs_def%key)
 ! END DART PREPROCESS INTERACTIVE_OBS_DEF
 
@@ -60,20 +61,21 @@
 ! BEGIN DART PREPROCESS MODULE CODE
 module obs_def_rttov_mod
 
-use        types_mod, only : r8, PI, metadatalength, MISSING_R8
+use        types_mod, only : r8, PI, metadatalength, MISSING_R8, MISSING_I
 use    utilities_mod, only : register_module, error_handler, E_ERR, E_WARN, E_MSG, &
                              logfileunit, get_unit, open_file, close_file, nc_check, &
                              file_exist, ascii_file_format
-use     location_mod, only : location_type, set_location, get_location, &
+use     location_mod, only : location_type, set_location, get_location, VERTISUNDEF, &
                              VERTISHEIGHT, VERTISLEVEL, set_location_missing
-use     obs_kind_mod, only : QTY_GEOPOTENTIAL_HEIGHT, QTY_SOIL_MOISTURE
+use     obs_kind_mod, only : QTY_GEOPOTENTIAL_HEIGHT, QTY_PRESSURE, QTY_TEMPERATURE, &
+                             QTY_SPECIFIC_HUMIDITY
 use  assim_model_mod, only : interpolate
 
 use obs_def_utilities_mod, only : track_status
 use ensemble_manager_mod,  only : ensemble_type
 
-use typesizes
-use netcdf
+!FIXME
+use rttov_interface_mod,   only : a, b, c
 
 implicit none
 private
@@ -96,23 +98,22 @@
 
 ! Metadata for rttov observations.
 
+! AQUA and TERRA are platform 9 (EOS)
+! AQUA may be satellite id 2  ! FIXME, check this
+
 ! AIRS is sensor 11 w/ 1-2378 channels (visible/near infrared/infrared)
 ! AMSU-A is sensor 3 with 1-15 channels (infrared/microwave)
 
-!FIXME
+!FIXME - add additional fields here as needed
 type obs_metadata
    private
-!    sat az/el
-!    sun az/el
-!    platform
-!    instrument
-!    channel
    real(r8)            :: sat_az     ! azimuth of satellite position
    real(r8)            :: sat_ze     ! azimuth of satellite position
    real(r8)            :: sun_az     ! zenith of solar position
    real(r8)            :: sun_ze     ! zenith of solar position