Importing Planetary Boundary Layer Height (PBLH) from Coupling

[hga007]

Description

This PR adds the Planetary Boundary Layer Height (PBLH), also known as the atmosphere boundary layer thickness (m), to the Import and Export state of ROMS and WRF ESMF-based NUOPC coupling modules. The PBLH is defined as the height of the turbulent mixed layer of the atmosphere in contact with the Earth's surface. The PBLH is crucial in computing turbulent heat, moisture, and momentum fluxes.

In bulk_fluxes.F is used to estimate wind gustiness:

!
!  Compute gustiness in wind speed.
!
            Bf=-g/TairK(i)*                                             &
     &         Wstar(i)*(Tstar(i)+0.61_r8*TairK(i)*Qstar(i))
            IF (Bf.gt.0.0_r8) THEN
# if defined FRC_COUPLING && defined PBLH
              Wgus(i)=blk_beta*(Bf*PBLh(i,j))**r3
# else
              Wgus(i)=blk_beta*(Bf*blk_Zabl)**r3
# endif
            ELSE
              Wgus(i)
           END IF
           delW(i)=SQRT(Wmag(i)*Wmag(i)+Wgus(i)*Wgus(i))

This option is possible when BULK_FLUXES and PBLH are activated during ESMF coupling using FRC_COUPLING and WRF_COUPLING.