Skip to content

Commit

Permalink
Merge pull request #143 from JosephMouallem/shiemom_fluxes
Browse files Browse the repository at this point in the history 
Add ocean temp and rough heat in for shield coupling
  • Loading branch information
@thomas-robinson
thomas-robinson authored Sep 27, 2024
2 parents f306ac1 + eadc455 commit 008399d
Showing 1 changed file with 20 additions and 3 deletions.
23 changes: 20 additions & 3 deletions full/atm_land_ice_flux_exchange.F90
View file
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)
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

0 comments on commit 008399d

Please sign in to comment.