Driver_MOD_PR06OD.f90

Go to the documentation of this file.

      subroutine driver_mod_pr06od( )
 
!-------------------------------------------------------------------
!
!!F77
!
!!Description:
! Driver program for MODIS cloud property algorithm.
!
!!Input Parameters:
! None
!
!!Output Parameters:
! None
!
!!Revision History:
! Revision 3.0  2004/08/24 mag
! the new driver
!
!!Team-unique Header:
! Cloud Retrieval Group, NASA/GSFC, Code 913, Greenbelt, Maryland, USA
!
!!Credit and Reference:
! Programmed by:
! Mark Gray (mag)
! L-3 GSI
! NASA/GSFC, code 913
! Greenbelt, Maryland, USA
!
!  W. Robinson, SAIC -re-worked to remove work file
!
!!End
!---------------------------------------------------------------------
  use names
  use mod06_run_settings
  use ch_xfr, only: cm_from_l2, c2_g_year, c2_g_day, c2_sensor_id
  use cld_tbl_names ! virtually use them all
 
  implicit none
 
  !.....Parameter Declarations
 
 
  !*****************************************************
  !     Variable Declarations
  !*****************************************************
 
  real :: statistics(7)
  integer :: ErrorFlag
  integer :: i
 
  integer :: status
  !
  !
  !--------------------------------------------------------------------------
  !   now, this mainly sets the table names wheras the rest of the information 
  !   comes through l2gen
  !---------------------------------------------------------------------------
  mytime = 0
  ! WDR we have made time the input so...
  myyear = c2_g_year
  myday = c2_g_day
  !  derive the month - it is used - just get a rough value
  mymonth = c2_g_day /30 + 1
  if ( mymonth > 12 ) mymonth = 12
  !
  !  WDR set the cloud mask to be from the l2
  cm_from_l2 = 1
    
  ! water and ice (over land) no-wind tables
  call get_cld_tbl( c2_sensor_id, scat_ice, aice_library, status )
 
  call get_cld_tbl( c2_sensor_id, scat_water, awater_library, status )
    
  ! ice wind-dependent libs
  call get_cld_tbl( c2_sensor_id, scat_ice_03_m_s, alibnames_ice(1), status )
  call get_cld_tbl( c2_sensor_id, scat_ice_sd_03_m_s, &
    alibnames_ice_sdev(1), status )
  call get_cld_tbl( c2_sensor_id, scat_ice_07_m_s, alibnames_ice(2), status )
  call get_cld_tbl( c2_sensor_id, scat_ice_sd_07_m_s, &
    alibnames_ice_sdev(2), status )
  call get_cld_tbl( c2_sensor_id, scat_ice_15_m_s, alibnames_ice(3), status )
  call get_cld_tbl( c2_sensor_id, scat_ice_sd_15_m_s, &
    alibnames_ice_sdev(3), status )
 
  ! water wind-dependent libs
  call get_cld_tbl( c2_sensor_id, scat_water_03_m_s, alibnames_water(1), &
    status )
  call get_cld_tbl( c2_sensor_id, scat_water_sd_03_m_s, &
   alibnames_water_sdev(1), status )
  call get_cld_tbl( c2_sensor_id, scat_water_07_m_s, alibnames_water(2), &
   status )
  call get_cld_tbl( c2_sensor_id, scat_water_sd_07_m_s, &
   alibnames_water_sdev(2), status )
  call get_cld_tbl( c2_sensor_id, scat_water_15_m_s, alibnames_water(3), &
    status )
  call get_cld_tbl( c2_sensor_id, scat_water_sd_15_m_s, &
    alibnames_water_sdev(3), status )
 
  ! single scattering / phase functions
  call get_cld_tbl( c2_sensor_id, single_scat, aphase_library, status )
  
  ! transmittance coefficients
  call get_cld_tbl( c2_sensor_id, transmittance_coeffs, &
    atransmittance_library, status )
    
  ! sfc, snow albedo map
  call get_cld_tbl( c2_sensor_id, ecology_map, aecosystem_data_name, status )
  call get_cld_tbl( c2_sensor_id, snow_albedo, asnowicealbedo_data_name, &
    status )
    
  ! emissivity map
  call get_cld_tbl( c2_sensor_id, global_emiss, aemissivity_name, status )
 
  !...........Perform cloud optical properties retrieval
  call mod_pr06od(statistics,errorflag)
 
 
 
End subroutine driver_mod_pr06od