subroutine get_scene_type(debug, landsea, andn, land_frc, scene_type) c----------------------------------------------------------------------------------- c Description: Set scene type that is used as an index in output CSRB output arrays. c c Input parameters: c debug debug level index c landsea land/water flag from NWP ancillary data (0=water, 1=land) c andn satellite node (ascending=1, descending=2) c land_frc land fraction in current HIRS FOV from collocated PATMOS-X c AVHRR (GAC) data. c c Output parameters: c scene_type scene type index (1=water, 2=AN land, 3=DN land) c c Revision history: c 04/11 R. Frey Original version c c Calls c None c c----------------------------------------------------------------------------------- IMPLICIT NONE c----------------------------------------------------------------------------------- c Arguments. integer debug, landsea, andn, scene_type real land_frc c Local variables. logical water c----------------------------------------------------------------------------------- c Set scene type. AVHRR-derived 'land_frc' is used whenever available. water = .false. if(land_frc .ne. -99.99) then if(land_frc .lt. 0.5) then water = .true. end if else if(landsea .eq. 0) then water = .true. end if end if if(water) then scene_type = 1 else if(andn .eq. 1) then scene_type = 2 else scene_type= 3 end if end if c----------------------------------------------------------------------------------- if(debug .gt. 0) then write(*,'(/, ''scene type determination'')') write(*,'(''land_frc, water, landsea, andn, scene_type: '', f7.2, l5, 3i5)') * land_frc, water, landsea, andn, scene_type end if c----------------------------------------------------------------------------------- return end c-----------------------------------------------------------------------------------