[Dart-dev] DART/branches Revision: 12190

[dart at ucar.edu]

hendric at ucar.edu
2017-12-08 13:54:05 -0700 (Fri, 08 Dec 2017)
72

adding in N2, O, O2, H molar masses to compute the 
forward operator.




Added: DART/branches/recam/observations/forward_operators/obs_def_oxygen_ion_density_mod.f90
===================================================================
--- DART/branches/recam/observations/forward_operators/obs_def_oxygen_ion_density_mod.f90	                        (rev 0)
+++ DART/branches/recam/observations/forward_operators/obs_def_oxygen_ion_density_mod.f90	2017-12-08 20:54:05 UTC (rev 12190)
@@ -0,0 +1,293 @@
+! 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: obs_def_ion_density_mod.f90 11692 2017-06-02 21:00:44Z nancy at ucar.edu $
+
+!-----------------------------------------------------------------------------
+!   DART Code:  Johnny Hendricks , hendric at ucar.edu
+!   Original DART/Radar work: Nancy Collins  
+!-----------------------------------------------------------------------------
+
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS KIND LIST
+! ELECTRON_DENSITY,            QTY_ELECTRON_DENSITY
+! DENSITY_ION_E,               QTY_DENSITY_ION_E
+! MOLEC_OXYGEN_MIXING_RATIO,   QTY_MOLEC_OXYGEN_MIXING_RATIO
+! ATOMIC_OXYGEN_MIXING_RATIO,  QTY_ATOMIC_OXYGEN_MIXING_RATIO
+! DENSITY_ION_OP,              QTY_DENSITY_ION_OP
+! ION_O_MIXING_RATIO,          QTY_ION_O_MIXING_RATIO
+! ATOMIC_H_MIXING_RATIO,       QTY_ATOMIC_H_MIXING_RATIO
+! END DART PREPROCESS KIND LIST
+!-----------------------------------------------------------------------------
+                             
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS USE OF SPECIAL OBS_DEF MODULE
+!    use obs_def_ion_density_mod, only : get_expected_oxygen_ion_val
+! END DART PREPROCESS USE OF SPECIAL OBS_DEF MODULE
+!-----------------------------------------------------------------------------
+
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
+!  case(DENSITY_ION_OP)
+!     call get_expected_oxygen_ion_val(state_handle, ens_size, location, expected_obs, istatus)
+! END DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
+!-----------------------------------------------------------------------------
+
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS READ_OBS_DEF
+!   case(ELECTRON_DENSITY, &
+!        DENSITY_ION_E, &
+!        MOLEC_OXYGEN_MIXING_RATIO, &
+!        ATOMIC_OXYGEN_MIXING_RATIO, &
+!        DENSITY_ION_OP, &
+!        ION_O_MIXING_RATIO, &
+!        ATOMIC_H_MIXING_RATIO)
+!      continue
+! END DART PREPROCESS READ_OBS_DEF
+!-----------------------------------------------------------------------------
+
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS WRITE_OBS_DEF
+!   case(ELECTRON_DENSITY, &
+!        DENSITY_ION_E, &
+!        MOLEC_OXYGEN_MIXING_RATIO, &
+!        ATOMIC_OXYGEN_MIXING_RATIO, &
+!        DENSITY_ION_OP, &
+!        ION_O_MIXING_RATIO, &
+!        ATOMIC_H_MIXING_RATIO)
+!      continue
+! END DART PREPROCESS WRITE_OBS_DEF
+!-----------------------------------------------------------------------------
+
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS INTERACTIVE_OBS_DEF
+!   case(ELECTRON_DENSITY, &
+!        DENSITY_ION_E, &
+!        MOLEC_OXYGEN_MIXING_RATIO, &
+!        ATOMIC_OXYGEN_MIXING_RATIO, &
+!        DENSITY_ION_OP, &
+!        ION_O_MIXING_RATIO, &
+!        ATOMIC_H_MIXING_RATIO)
+!      continue
+! END DART PREPROCESS INTERACTIVE_OBS_DEF
+!-----------------------------------------------------------------------------
+
+!-----------------------------------------------------------------------------
+! BEGIN DART PREPROCESS MODULE CODE
+module obs_def_ion_density_mod
+
+use        types_mod, only : r8, missing_r8, PI, deg2rad
+use    utilities_mod, only : register_module, error_handler, E_ERR, E_MSG, &
+                             check_namelist_read, find_namelist_in_file,   &
+                             nmlfileunit, do_output, do_nml_file, do_nml_term, &
+                             ascii_file_format
+use     location_mod, only : location_type, write_location, read_location, &
+                             interactive_location, get_location
+use  assim_model_mod, only : interpolate
+use     obs_kind_mod, only : QTY_ELECTRON_DENSITY,           & ! Right QTY to use?
+                             QTY_DENSITY_ION_E,              & ! Right QTY to use?
+                             QTY_MOLEC_OXYGEN_MIXING_RATIO,  &
+                             QTY_ATOMIC_OXYGEN_MIXING_RATIO, &
+                             QTY_DENSITY_ION_OP,             &                  
+                             QTY_ION_O_MIXING_RATIO,         & ! newly defined
+                             QTY_ATOMIC_H_MIXING_RATIO,      & ! newly defined
+                             QTY_TEMPERATURE,                &
+                             QTY_PRESSURE