[Dart-dev] DART/branches Revision: 11139

[dart at ucar.edu]

mizzi at ucar.edu
2017-02-22 13:53:02 -0700 (Wed, 22 Feb 2017)
34
Update WRF-Chem/DART repository.




Modified: DART/branches/mizzi/obs_def/obs_def_IASI_CO_mod.f90
===================================================================
--- DART/branches/mizzi/obs_def/obs_def_IASI_CO_mod.f90	2017-02-22 20:51:31 UTC (rev 11138)
+++ DART/branches/mizzi/obs_def/obs_def_IASI_CO_mod.f90	2017-02-22 20:53:02 UTC (rev 11139)
@@ -39,12 +39,11 @@
 !
 ! 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 mpi_utilities_mod, only : my_task_id
+use        types_mod, only : 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
 implicit none
 private
 public :: write_iasi_co, read_iasi_co, interactive_iasi_co, &
@@ -79,7 +78,6 @@
 !----------------------------------------------------------------------------
 ! subroutine initialize_module
 call register_module(source, revision, revdate)
-!print *, 'APM: Initialize IASI CO - Task ', my_task_id()
 module_initialized = .true.
 end subroutine initialize_module
 !
@@ -98,7 +96,6 @@
 real(r8), dimension(iasi_dimp)	:: pressure_1
 integer 			:: keyin
 if ( .not. module_initialized ) call initialize_module
-!print *, 'APM: Read IASI CO - Task ', my_task_id()
 fileformat = "ascii"   ! supply default
 if(present(fform)) fileformat = trim(adjustl(fform))
 
@@ -141,7 +138,6 @@
 real(r8), dimension(iasi_dim)   :: avg_kernels_temp
 real(r8), dimension(iasi_dimp)  :: pressure_temp
 if ( .not. module_initialized ) call initialize_module
-!print *, 'APM: Write IASI CO - Task ', my_task_id()
 fileformat = "ascii"   ! supply default
 if(present(fform)) fileformat = trim(adjustl(fform))
 !
@@ -177,7 +173,6 @@
 integer, intent(out) :: key
 character(len=129) :: msgstring
 if ( .not. module_initialized ) call initialize_module
-!print *, 'APM: Interactive IASI CO - Task ', my_task_id()
 !
 ! Make sure there's enough space, if not die for now (clean later)
 if(num_iasi_co_obs >= max_iasi_co_obs) then
@@ -207,149 +202,149 @@
 subroutine get_expected_iasi_co(state, location, key, val, istatus)
 !----------------------------------------------------------------------
 !subroutine get_expected_iasi_co(state, location, key, val, istatus)
-   real(r8), intent(in)            :: state(:)
-   type(location_type), intent(in) :: location
-   integer, intent(in)             :: key
-   real(r8), intent(out)           :: val
-   integer, intent(out)            :: istatus
-   integer :: i,kstr
-   type(location_type) :: loc2
-   real(r8)            :: mloc(3)
-   real(r8)            :: obs_val,wrf_psf,level
-   real(r8)            :: co_min,obs_val_min,iasi_prs_mid,iasi_psf,iasi_psf_save,iasi_prs
-   integer             :: nlevels,nlevelsp,nnlevels
-   integer             :: iflg
-   character(len=129)  :: msgstring
-   real(r8)            :: vert_mode_filt
+real(r8), intent(in)            :: state(:)
+type(location_type), intent(in) :: location
+integer, intent(in)             :: key
+real(r8), intent(out)           :: val
+integer, intent(out)            :: istatus
+integer :: i,kstr
+type(location_type) :: loc2
+real(r8)            :: mloc(3)
+real(r8)	    :: obs_val,wrf_psf,level,missing
+real(r8)            :: co_min,iasi_prs_mid,iasi_psf,iasi_psf_save
+integer             :: nlevels,nlevelsp,nnlevels
+integer             :: iflg
+character(len=129)  :: msgstring
+real(r8)            :: vert_mode_filt
 !
 ! Initialize DART
-   if ( .not. module_initialized ) call initialize_module
-!   write(msgstring, *)'APM DEBUG: Enter IASI CO obs def'
-!   call error_handler(E_MSG,'set_obs_def_iasi_co',msgstring,source,revision,revdate)
-!   print *, 'APM: Enter Obs Def - Task ', my_task_id()
+if ( .not. module_initialized ) call initialize_module
 !
 ! Initialize variables
-   co_min=1.e-2
-   obs_val_min=4.e-2
+co_min=1.e-4
+missing=-888888.0_r8
 !
 ! Get iasi data