[Dart-dev] DART/branches Revision: 10905

[dart at ucar.edu]

mizzi at ucar.edu
2017-01-17 12:01:53 -0700 (Tue, 17 Jan 2017)
49
Saving in preparation to move to the RMA_TRUNK.




Modified: DART/branches/mizzi/obs_def/obs_def_IASI_CO_mod.f90
===================================================================
--- DART/branches/mizzi/obs_def/obs_def_IASI_CO_mod.f90	2017-01-17 17:20:34 UTC (rev 10904)
+++ DART/branches/mizzi/obs_def/obs_def_IASI_CO_mod.f90	2017-01-17 19:01:53 UTC (rev 10905)
@@ -1,123 +1,104 @@
-! 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
-! http://www.image.ucar.edu/DAReS/DART/DART_download
+! Data Assimilation Research Testbed -- DART
+! Copyright 2004, 2005, Data Assimilation Initiative, University Corporation for Atmospheric Research
+! Licensed under the GPL -- www.gpl.org/licenses/gpl.html
 !
-! $Id$
-
 ! BEGIN DART PREPROCESS KIND LIST
 ! IASI_CO_RETRIEVAL, KIND_CO
 ! END DART PREPROCESS KIND LIST
-
+!
 ! BEGIN DART PREPROCESS USE OF SPECIAL OBS_DEF MODULE
 !   use obs_def_iasi_co_mod, only : write_iasi_co, read_iasi_co, &
 !                                  interactive_iasi_co, get_expected_iasi_co, &
 !                                  set_obs_def_iasi_co
 ! END DART PREPROCESS USE OF SPECIAL OBS_DEF MODULE
-
+!
 ! BEGIN DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
-!         case(IASI_CO_RETRIEVAL)
-!            call get_expected_iasi_co(state, location, obs_def%key, obs_val, istatus)
+!         case(IASI_CO_RETRIEVAL)                                                           
+!            call get_expected_iasi_co(state, location, obs_def%key, obs_val, istatus)  
 ! END DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
-
+!
 ! BEGIN DART PREPROCESS READ_OBS_DEF
 !      case(IASI_CO_RETRIEVAL)
 !         call read_iasi_co(obs_def%key, ifile, fileformat)
 ! END DART PREPROCESS READ_OBS_DEF
-
+!
 ! BEGIN DART PREPROCESS WRITE_OBS_DEF
 !      case(IASI_CO_RETRIEVAL)
 !         call write_iasi_co(obs_def%key, ifile, fileformat)
 ! END DART PREPROCESS WRITE_OBS_DEF
-
+!
 ! BEGIN DART PREPROCESS INTERACTIVE_OBS_DEF
 !      case(IASI_CO_RETRIEVAL)
 !         call interactive_iasi_co(obs_def%key)
 ! END DART PREPROCESS INTERACTIVE_OBS_DEF
-
+!
 ! BEGIN DART PREPROCESS SET_OBS_DEF_IASI_CO
 !      case(IASI_CO_RETRIEVAL)
 !         call set_obs_def_iasi_co(obs_def%key)
 ! END DART PREPROCESS SET_OBS_DEF_IASI_CO
-
+!
 ! BEGIN DART PREPROCESS MODULE CODE
-
 module obs_def_iasi_CO_mod
-
-use        types_mod, only : r8, MISSING_R8
-use    utilities_mod, only : register_module, error_handler, E_ERR, E_MSG
-use     location_mod, only : location_type, set_location, get_location, VERTISPRESSURE, VERTISSURFACE
-use  assim_model_mod, only : interpolate
-use     obs_kind_mod, only : KIND_CO, KIND_SURFACE_PRESSURE
-
+use         types_mod, only : r8, missing_r8
+use     utilities_mod, only : register_module, error_handler, E_ERR, E_MSG
+use      location_mod, only : location_type, set_location, get_location, VERTISPRESSURE, VERTISSURFACE
+use   assim_model_mod, only : interpolate
+use      obs_kind_mod, only : KIND_CO, KIND_SURFACE_PRESSURE
+use mpi_utilities_mod, only : my_task_id
 implicit none
 private
+public :: write_iasi_co, read_iasi_co, interactive_iasi_co, &
+          get_expected_iasi_co, set_obs_def_iasi_co
 
-public :: write_iasi_co, &
-          read_iasi_co, &
-          interactive_iasi_co, &
-          get_expected_iasi_co, &
-          set_obs_def_iasi_co
-
 ! Storage for the special information required for observations of this type
-integer, parameter :: MAX_IASI_CO_OBS = 10000000
-integer, parameter :: IASI_DIM = 19
-integer, parameter :: IASI_DIMP = 20
-integer            :: num_iasi_co_obs = 0
-integer            :: counts1 = 0
-
-real(r8), dimension(MAX_IASI_CO_OBS,IASI_DIM)  :: avg_kernel
-real(r8), dimension(MAX_IASI_CO_OBS,IASI_DIMP) :: pressure
-real(r8), dimension(MAX_IASI_CO_OBS)           :: iasi_prior
-real(r8), dimension(MAX_IASI_CO_OBS)           :: iasi_psurf
-integer,  dimension(MAX_IASI_CO_OBS)           :: iasi_nlevels
-integer,  dimension(MAX_IASI_CO_OBS)           :: iasi_nlevelsp
-
+integer, parameter               :: max_iasi_co_obs = 10000000
+integer, parameter               :: iasi_dim = 19