chgres.F90

Go to the documentation of this file.

 
 
 program chgres
 
 use mpi_f08
 use esmf
 
 use atmosphere, only          : atmosphere_driver
 
 use program_setup, only       : read_setup_namelist, &
                                 read_varmap, &
                                 convert_atm, &
                                 convert_sfc 
 
 use model_grid, only          : define_target_grid,  &
                                 define_input_grid, &
                                 cleanup_input_target_grid_data
 
 use surface, only             : surface_driver
 
 use utilities, only           : error_handler
 implicit none
 
 integer                      :: ierr, localpet, npets
 
 type(esmf_vm)                :: vm
 
!-------------------------------------------------------------------------
! Initialize mpi and esmf environment.
!-------------------------------------------------------------------------
 
 call mpi_init(ierr)
 
 print*,"- INITIALIZE ESMF"
 call esmf_initialize(rc=ierr)
 if(esmf_logfounderror(rctocheck=ierr,msg=esmf_logerr_passthru,line=__line__,file=__file__)) &
    call error_handler("INITIALIZING ESMF", ierr)
 
 print*,"- CALL VMGetGlobal"
 call esmf_vmgetglobal(vm, rc=ierr)
 if(esmf_logfounderror(rctocheck=ierr,msg=esmf_logerr_passthru,line=__line__,file=__file__)) &
    call error_handler("IN VMGetGlobal", ierr)
 
 print*,"- CALL VMGet"
 call esmf_vmget(vm, localpet=localpet, petcount=npets, rc=ierr)
 if(esmf_logfounderror(rctocheck=ierr,msg=esmf_logerr_passthru,line=__line__,file=__file__)) &
    call error_handler("IN VMGet", ierr)
 
 print*,'- NPETS IS  ',npets
 print*,'- LOCAL PET ',localpet
 
!-------------------------------------------------------------------------
! Read program configuration namelist.
!-------------------------------------------------------------------------
 
 call read_setup_namelist
 
!-------------------------------------------------------------------------
! Read variable mapping file (used for grib2 input data only).
!-------------------------------------------------------------------------
 
 call read_varmap
 
!-------------------------------------------------------------------------
! Create esmf grid objects for input and target grids.
!-------------------------------------------------------------------------
 
 call define_target_grid(localpet, npets)
 
 call define_input_grid(localpet, npets)
 
!-------------------------------------------------------------------------
! Convert atmospheric fields
!-------------------------------------------------------------------------
 
 if (convert_atm) then
 
   call atmosphere_driver(localpet)
 
 end if
 
!-------------------------------------------------------------------------
! Convert surface/nsst fields
!-------------------------------------------------------------------------
 
 if (convert_sfc) then
 
   call surface_driver(localpet)
 
 end if
 
 call cleanup_input_target_grid_data
 
 print*,"- CALL ESMF_finalize"
 call esmf_finalize(endflag=esmf_end_keepmpi, rc=ierr)
 
 call mpi_finalize(ierr)
 
 print*,"- DONE."
 
 end program chgres