From d1bdd106f2c145bdd800eb3f6915bf489fa0ba0d Mon Sep 17 00:00:00 2001
From: GeorgeGayno-NOAA <52789452+GeorgeGayno-NOAA@users.noreply.github.com>
Date: Fri, 11 Feb 2022 14:14:16 -0500
Subject: [PATCH] chgres_cube - Run routine 'convert_omega' on all tasks.
 (#627)

When using certain GRIB2 data as input, the vertical velocity must be converted from omega to dzdt. This conversion is controlled by the logical 'conv_omega'. Ensure that logical is set on all MPI tasks.

Fixes #626
---
 sorc/chgres_cube.fd/input_data.F90 | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/sorc/chgres_cube.fd/input_data.F90 b/sorc/chgres_cube.fd/input_data.F90
index dab617f35..965c37df1 100644
--- a/sorc/chgres_cube.fd/input_data.F90
+++ b/sorc/chgres_cube.fd/input_data.F90
@@ -2457,6 +2457,7 @@ end subroutine read_input_atm_tiled_history_file
 !! @author George Gayno NCEP/EMC   
  subroutine read_input_atm_grib2_file(localpet)
 
+ use mpi
  use wgrib2api
  
  use grib2_util, only                   : rh2spfh, rh2spfh_gfs, convert_omega
@@ -2909,6 +2910,8 @@ subroutine read_input_atm_grib2_file(localpet)
    enddo
  endif
 
+ call mpi_bcast(conv_omega,1,MPI_LOGICAL,0,MPI_COMM_WORLD,rc)
+
  if (localpet == 0) print*,"- CALL FieldScatter FOR INPUT DZDT."
  call ESMF_FieldScatter(dzdt_input_grid, dummy3d, rootpet=0, rc=rc)
  if(ESMF_logFoundError(rcToCheck=rc,msg=ESMF_LOGERR_PASSTHRU,line=__LINE__,file=__FILE__)) &
@@ -3081,7 +3084,7 @@ subroutine read_input_atm_grib2_file(localpet)
                     farrayPtr=presptr, rc=rc)
   if(ESMF_logFoundError(rcToCheck=rc,msg=ESMF_LOGERR_PASSTHRU,line=__LINE__,file=__FILE__)) &
     call error_handler("IN FieldGet", rc)
-    
+
   call convert_omega(wptr,presptr,tptr,qptr,clb,cub)
   
  endif