Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Importing Planetary Boundary Layer Height (PBLH) from Coupling #52

Open
wants to merge 5 commits into
base: develop
Choose a base branch
from
Open

Importing Planetary Boundary Layer Height (PBLH) from Coupling #52

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
b1da466
Import Planetary Boundary Layer Height from Coupling
hga007 Feb 8, 2025
70d24e0
Add missing logic to set Qout(idPBLh)
hga007 Mar 2, 2025
423ba42
Merging latest develop into feature/pblh_coupling
hga007 Mar 9, 2025
3dac379
Merging latest develop into feature/pblh_coupling
hga007 Mar 9, 2025
eaeff2c
Merging latest develop into feature/pblh_coupling
hga007 Mar 19, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions ESM/coupling_esmf.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -253,6 +253,7 @@ WeightsFile(atmos) = meld_weights_atm.nc
! Ustr atmos roms, cice surface U-momentum stress
! Vstr atmos roms, cice surface V-momentum stress
! Wstar atmos roms turbulent frictional wind magnitude
! pblh atmos roms atmosphere boundary layer thickness
!
! zlvl atmos cice atmspheric height lowest level
! rhoa atmos cice surface air density
Expand Down
15 changes: 15 additions & 0 deletions ESM/coupling_esmf.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -349,3 +349,18 @@ metadata:
connected_to: *ATM
regrid_method: bilinear
extrapolate_method: nearest

- standard_name: atmosphere_boundary_layer_thickness
long_name: atmosphere boundary layer thickness
short_name: pblh
data_variables: [pblh, time]
source_units: m
destination_units: m
source_grid: cell_center
destination_grid: cell_center
add_offset: 0.0d0
scale: 1.0d0
debug_write: false
connected_to: *ATM
regrid_method: bilinear
extrapolate_method: nearest
15 changes: 15 additions & 0 deletions ESM/coupling_esmf_wrf.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -469,3 +469,18 @@ metadata:
connected_to: *ATM
regrid_method: bilinear
extrapolate_method: nearest

- standard_name: atmosphere_boundary_layer_thickness
long_name: atmosphere boundary layer thickness
short_name: pblh
data_variables: [pblh, time]
source_units: m
destination_units: m
source_grid: cell_center
destination_grid: cell_center
add_offset: 0.0d0
scale: 1.0d0
debug_write: false
connected_to: *ATM
regrid_method: bilinear
extrapolate_method: nearest
15 changes: 15 additions & 0 deletions ESM/roms_cmeps.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -431,3 +431,18 @@ import:
connected_to: *ATM
map_type: mapbilnr
map_norm: none

- standard_name: atmosphere_boundary_layer_thickness
long_name: atmosphere boundary layer thickness
short_name: pblh
data_variables: [pblh, time]
source_units: m
destination_units: m
source_grid: cell_center
destination_grid: cell_center
add_offset: 0.0d0
scale: 1.0d0
debug_write: false
connected_to: *ATM
map_type: mapbilnr
map_norm: none
42 changes: 39 additions & 3 deletions Master/cmeps_roms.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -126,9 +126,9 @@
USE mod_iounits, ONLY : Iname, SourceFile, stdout
USE mod_mixing, ONLY : MIXING
USE mod_ncparam, ONLY : Iinfo, idLdwn, idLrad, idPair, &
& idQair, idrain, idSrad, idTair, &
& idTsur, idUair, idUsms, idVair, &
& idVsms
& idPBLh, idQair, idrain, idSrad, &
& idTair, idTsur, idUair, idUsms, &
& idVair, idVsms
# ifdef TIME_INTERP
USE mod_netcdf, ONLY : netcdf_get_ivar, &
& netcdf_get_svar, &
Expand Down Expand Up @@ -5538,6 +5538,42 @@
& Vstress)
END IF
# endif
# if defined BULK_FLUXES & defined PBLH
!
! Atmosphere boundary layer thickness (m).
!
CASE ('pblh', 'PBLh')
romsScale=scale
ifield=idPBLh
gtype=r2dvar
Tindex=3-Iinfo(8,ifield,ng)
DO j=JstrR,JendR
DO i=IstrR,IendR
IF (ABS(ptr2d(i,j)).lt.TOL_dp) THEN
Fval=scale*ptr2d(i,j)+add_offset
ELSE
Fval=0.0_dp
END IF
MyFmin(1)=MIN(MyFmin(1),ptr2d(i,j))
MyFmax(1)=MAX(MyFmax(1),ptr2d(i,j))
MyFmin(2)=MIN(MyFmin(2),Fval)
MyFmax(2)=MAX(MyFmax(2),Fval)
FORCES(ng)%PBLh(i,j)=Fval
END DO
END DO
IF (localDE.eq.localDEcount-1) THEN
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& FORCES(ng)%PBLh)
END IF
CALL mp_exchange2d (ng, tile, iNLM, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& FORCES(ng)%PBLh)
END IF
# endif
# if defined WIND_MINUS_CURRENT && !defined BULK_FLUXES
!
! Surface air density (kg/m3).
Expand Down
21 changes: 21 additions & 0 deletions Master/esmf_atm_wrf.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4580,6 +4580,27 @@
ptr2d(i,j)=Fval
END DO
END DO
# if defined BULK_FLUXES & defined PBLH
!
! Atmosphere boundary layer thickness (m).
!
CASE ('pblh', 'PBLh')
MyFmin(1)= MISSING_dp
MyFmax(1)=-MISSING_dp
DO j=Jstr,Jend
DO i=Istr,Iend
# ifdef WRF_TIMEAVG
! Fval=REAL(grid%pblhavg(i,j),dp)
Fval=REAL(grid%pblh(i,j),dp)
# else
Fval=REAL(grid%pblh(i,j),dp)
# endif
MyFmin(1)=MIN(MyFmin(1),Fval)
MyFmax(1)=MAX(MyFmax(1),Fval)
ptr2d(i,j)=Fval
END DO
END DO
# endif
!
! Surface (2m) air temperature (K).
!
Expand Down
42 changes: 39 additions & 3 deletions Master/esmf_roms.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -108,9 +108,9 @@
USE mod_iounits, ONLY : Iname, SourceFile, stdout
USE mod_mixing, ONLY : MIXING
USE mod_ncparam, ONLY : Iinfo, idLdwn, idLrad, idPair, &
& idQair, idrain, idSrad, idTair, &
& idTsur, idUair, idUsms, idVair, &
& idVsms
& idPBLh, idQair, idrain, idSrad, &
& idTair, idTsur, idUair, idUsms, &
& idVair, idVsms
# ifdef TIME_INTERP
USE mod_netcdf, ONLY : netcdf_get_ivar, &
& netcdf_get_svar, &
Expand Down Expand Up @@ -3934,6 +3934,42 @@
& Vstress)
END IF
# endif
# if defined BULK_FLUXES & defined PBLH
!
! Atmosphere boundary layer thickness (m).
!
CASE ('pblh', 'PBLh')
romsScale=scale
ifield=idPBLh
gtype=r2dvar
Tindex=3-Iinfo(8,ifield,ng)
DO j=JstrR,JendR
DO i=IstrR,IendR
IF (ABS(ptr2d(i,j)).lt.TOL_dp) THEN
Fval=scale*ptr2d(i,j)+add_offset
ELSE
Fval=0.0_dp
END IF
MyFmin(1)=MIN(MyFmin(1),ptr2d(i,j))
MyFmax(1)=MAX(MyFmax(1),ptr2d(i,j))
MyFmin(2)=MIN(MyFmin(2),Fval)
MyFmax(2)=MAX(MyFmax(2),Fval)
FORCES(ng)%PBLh(i,j)=Fval
END DO
END DO
IF (localDE.eq.localDEcount-1) THEN
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& FORCES(ng)%PBLh)
END IF
CALL mp_exchange2d (ng, tile, iNLM, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& FORCES(ng)%PBLh)
END IF
# endif
# if defined WIND_MINUS_CURRENT && !defined BULK_FLUXES
!
! Surface air density (kg/m3).
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/da_roms.tmp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -767,6 +767,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2886,6 +2887,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/nl_roms.tmp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -767,6 +767,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2886,6 +2887,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_basin.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -766,6 +766,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2885,6 +2886,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_benchmark1.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_benchmark2.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_benchmark3.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_bio_toy.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2894,6 +2895,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_bl_test.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_canyon2d.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_canyon3d.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_channel.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
2 changes: 2 additions & 0 deletions ROMS/External/roms_damee_4.in
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -768,6 +768,7 @@ Qout(idWdis) == F ! Dissip_roller wave roller dissipation
Qout(idWrol) == F ! roller_action wave roller action density

Qout(idPair) == F ! Pair surface air pressure
Qout(idPBLh) == F ! PBLh atmosphere boundary layer thickness
Qout(idTair) == F ! Tair near surface air temperature
Qout(idUair) == F ! Uair near surface U-wind
Qout(idVair) == F ! Vair near surface V-wind
Expand Down Expand Up @@ -2887,6 +2888,7 @@ PIO_I2C_Preq = 65 ! Maximum pending I2C requests
! Qout(idWrol) Write out wave roller action density.
!
! Qout(idPair) Write out surface air pressure.
! Qout(idPBLh) Write out atmosphere boundary layer height.
! Qout(idTair) Write out near surface air temperature.
! Qout(idUair) Write out near surface U-wind component.
! Qout(idVair) Write out near surface V-wind component.
Expand Down
Loading