chgres.F90

Go to the documentation of this file.

 program chgres

 use mpi
 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

 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