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

Add ocean temp and rough heat in for shield coupling #143

Merged
merged 3 commits into from
Sep 27, 2024
Merged

Add ocean temp and rough heat in for shield coupling #143

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
2aa025b
add ocean temp and rough heat
JosephMouallem Aug 15, 2024
b8c12ae
Merge branch 'NOAA-GFDL:main' into shiemom_fluxes
JosephMouallem Sep 11, 2024
eadc455
reduce line length
JosephMouallem Sep 26, 2024
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
23 changes: 20 additions & 3 deletions full/atm_land_ice_flux_exchange.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -107,7 +107,7 @@ module atm_land_ice_flux_exchange_mod
id_rough_moist, id_rough_heat, id_rough_mom, &
id_land_mask, id_ice_mask, &
id_u_star, id_b_star, id_q_star, id_u_flux, id_v_flux, &
id_t_surf, id_t_flux, id_r_flux, id_q_flux, id_slp, &
id_t_surf, id_t_ocean, id_t_flux, id_r_flux, id_q_flux, id_slp, &
id_t_atm, id_u_atm, id_v_atm, id_wind, &
id_thv_atm, id_thv_surf, &
id_t_ref, id_rh_ref, id_u_ref, id_v_ref, id_wind_ref, &
Expand Down Expand Up @@ -546,6 +546,7 @@ subroutine atm_land_ice_flux_exchange_init(Time, Atm, Land, Ice, atmos_ice_bound
!allocate land_ice_atmos_boundary
call fms_mpp_domains_get_compute_domain( Atm%domain, is, ie, js, je )
allocate( land_ice_atmos_boundary%t(is:ie,js:je) )
allocate( land_ice_atmos_boundary%t_ocean(is:ie,js:je) )! Joseph: surf ocean temp
allocate( land_ice_atmos_boundary%u_ref(is:ie,js:je) ) ! bqx
allocate( land_ice_atmos_boundary%v_ref(is:ie,js:je) ) ! bqx
allocate( land_ice_atmos_boundary%t_ref(is:ie,js:je) ) ! cjg: PBL depth mods
Expand Down Expand Up @@ -573,9 +574,11 @@ subroutine atm_land_ice_flux_exchange_init(Time, Atm, Land, Ice, atmos_ice_bound
allocate( land_ice_atmos_boundary%thv_atm(is:ie,js:je) )
allocate( land_ice_atmos_boundary%thv_surf(is:ie,js:je) )
allocate( land_ice_atmos_boundary%rough_mom(is:ie,js:je) )
allocate( land_ice_atmos_boundary%rough_heat(is:ie,js:je) ) ! Kun
allocate( land_ice_atmos_boundary%frac_open_sea(is:ie,js:je) )
! initialize boundary values for override experiments (mjh)
land_ice_atmos_boundary%t=273.0
land_ice_atmos_boundary%t_ocean=200.0
land_ice_atmos_boundary%u_ref=0.0 ! bqx
land_ice_atmos_boundary%v_ref=0.0 ! bqx
land_ice_atmos_boundary%t_ref=273.0 ! cjg: PBL depth mods
Expand Down Expand Up @@ -603,6 +606,7 @@ subroutine atm_land_ice_flux_exchange_init(Time, Atm, Land, Ice, atmos_ice_bound
land_ice_atmos_boundary%thv_atm=0.0
land_ice_atmos_boundary%thv_surf=0.0
land_ice_atmos_boundary%rough_mom=0.01
land_ice_atmos_boundary%rough_heat=0.01
land_ice_atmos_boundary%frac_open_sea=0.0

! allocate fields for extra tracers
Expand Down Expand Up @@ -1424,6 +1428,8 @@ subroutine sfc_boundary_layer ( dt, Time, Atm, Land, Ice, Land_Ice_Atmos_Boundar
enddo

! [6.2] put relevant quantities onto atmospheric boundary
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%t_ocean, 'ATM', &
ex_t_surf , xmap_sfc, complete=.false.) !joseph
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%t, 'ATM', ex_t_surf4 , xmap_sfc, complete=.false.)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%frac_open_sea,'ATM',ex_frac_open_sea, xmap_sfc)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%albedo, 'ATM', ex_albedo , xmap_sfc, complete=.false.)
Expand All @@ -1436,6 +1442,8 @@ subroutine sfc_boundary_layer ( dt, Time, Atm, Land, Ice, Land_Ice_Atmos_Boundar
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%albedo_nir_dif, 'ATM', &
ex_albedo_nir_dif , xmap_sfc, complete=.false.)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%rough_mom, 'ATM', ex_rough_mom, xmap_sfc, complete=.false.)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%rough_heat,'ATM', &
ex_rough_heat, xmap_sfc, complete=.false.) !kgao
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%land_frac, 'ATM', ex_land_frac, xmap_sfc, complete=.false.)

call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%u_flux, 'ATM', ex_flux_u, xmap_sfc, complete=.false.)
Expand All @@ -1449,9 +1457,10 @@ subroutine sfc_boundary_layer ( dt, Time, Atm, Land, Ice, Land_Ice_Atmos_Boundar
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%u_ref, 'ATM', ex_ref_u , xmap_sfc, complete=.false.) !bqx
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%v_ref, 'ATM', ex_ref_v , xmap_sfc, complete=.true.) !bqx

! kgao: for shield+mom6 coupling; used by shield pbl schemes (am5 with tke-edmf should do the same)
#ifndef use_AM3_physics
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%shflx, 'ATM', ex_flux_t , xmap_sfc)
JosephMouallem marked this conversation as resolved.
Show resolved Hide resolved
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%lhflx, 'ATM', ex_flux_tr(:,isphum), xmap_sfc)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%shflx,'ATM', ex_flux_t , xmap_sfc, complete=.false.)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%lhflx,'ATM', ex_flux_tr(:,isphum), xmap_sfc, complete=.true.)
#endif
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%wind, 'ATM', ex_wind , xmap_sfc)
call fms_xgrid_get_from_xgrid (Land_Ice_Atmos_Boundary%thv_atm, 'ATM', ex_thv_atm , xmap_sfc)
Expand Down Expand Up @@ -1702,6 +1711,9 @@ subroutine sfc_boundary_layer ( dt, Time, Atm, Land, Ice, Land_Ice_Atmos_Boundar
if ( id_psl > 0 ) used = fms_diag_send_data ( id_psl, diag_atm, Time )
endif

if ( id_t_ocean > 0 ) then
used = fms_diag_send_data ( id_t_ocean, Land_Ice_Atmos_Boundary%t_ocean, Time )
endif
!-----------------------------------------------------------------------
!--------- diagnostics for fields at reference level ---------
!cjg
Expand Down Expand Up @@ -3499,6 +3511,11 @@ subroutine diag_field_init ( Time, atmos_axes, land_axes, land_pe )
fms_diag_register_diag_field ( mod_name, 'tau_y', atmos_axes, Time, &
'meridional wind stress', 'pa' )

id_t_ocean = &
fms_diag_register_diag_field ( mod_name, 't_ocean', atmos_axes, Time, &
'surface temperature from ocean output', 'deg_k', &
range=trange )

id_t_surf = &
fms_diag_register_diag_field ( mod_name, 't_surf', atmos_axes, Time, &
'surface temperature', 'deg_k', &
Expand Down
Loading