var/da/da_radiance/da_det_crtm_climat.inc

   1 subroutine da_det_crtm_climat (lat, month, crtm_climat)
   2
   3 ! determine the CRTM climatology model according to the
   4 ! latitude and month
   5 ! CRTM climatology model is categorized in CRTM_Atmosphere_Define.f90 as
   6 ! INTEGER, PARAMETER :: INVALID_MODEL          = 0
   7 ! INTEGER, PARAMETER :: TROPICAL               = 1
   8 ! INTEGER, PARAMETER :: MIDLATITUDE_SUMMER     = 2
   9 ! INTEGER, PARAMETER :: MIDLATITUDE_WINTER     = 3
  10 ! INTEGER, PARAMETER :: SUBARCTIC_SUMMER       = 4
  11 ! INTEGER, PARAMETER :: SUBARCTIC_WINTER       = 5
  12 ! INTEGER, PARAMETER :: US_STANDARD_ATMOSPHERE = 6
  13
  14    implicit none
  15
  16    real,    intent(in)  :: lat
  17    integer, intent(in)  :: month
  18    integer, intent(out) :: crtm_climat
  19
  20    crtm_climat = 6  ! initialized to be us_standard_atmosphere
  21
  22    ! 1: Tropical           (-23.4378 ~ 23.4378)
  23    if ( lat >= -23.4378 .and. lat <= 23.4378  ) crtm_climat = 1
  24
  25    ! 2: Midlatitude summer (23.4378 ~ 66.561)
  26    ! 3: Midlatitude winter
  27    if ( lat > 23.4378 .and. lat <= 66.561  ) then ! North Mid-Lat
  28       if ( month >= 3 .and. month <= 8 ) then
  29          crtm_climat = 2  ! Summer
  30       else
  31          crtm_climat = 3  ! Winter
  32       end if
  33    end if
  34
  35    if ( lat < -23.4378 .and. lat >= -66.561  ) then ! South Mid-Lat
  36       if ( month >= 3 .and. month <= 8 ) then
  37          crtm_climat = 3  ! Winter
  38       else
  39          crtm_climat = 2  ! Summer
  40       end if
  41    end if
  42
  43    ! 4: Subarctic summer   ( > 66.561)
  44    ! 5: Subarctic winter
  45    if ( lat > 66.561  ) then ! Subarctic
  46       if ( month >= 3 .and. month <= 8 ) then
  47          crtm_climat = 4  ! Summer
  48       else
  49          crtm_climat = 5  ! Winter
  50       end if
  51    end if
  52
  53    if ( lat < -66.561  ) then ! Subantarctic
  54       if ( month >= 3 .and. month <= 8 ) then
  55          crtm_climat = 5  ! Winter
  56       else
  57          crtm_climat = 4  ! Summer
  58       end if
  59    end if
  60
  61 end subroutine da_det_crtm_climat