Skip to content

Enable exchange of runoff DOC from Land to Ice #178

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

Open
nikizadehgfdl wants to merge 1 commit into
base: ESM45branch
Choose a base branch
from
+12 −1
Open

Enable exchange of runoff DOC from Land to Ice #178

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
5 changes: 5 additions & 0 deletions full/ice_ocean_flux_exchange.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -136,6 +136,7 @@ subroutine ice_ocean_flux_exchange_init(Time, Ice, Ocean, Ocean_state, ice_ocean
allocate( ice_ocean_boundary%lprec (is:ie,js:je) ) ; ice_ocean_boundary%lprec = 0.0
allocate( ice_ocean_boundary%fprec (is:ie,js:je) ) ; ice_ocean_boundary%fprec = 0.0
allocate( ice_ocean_boundary%runoff (is:ie,js:je) ) ; ice_ocean_boundary%runoff = 0.0
allocate( ice_ocean_boundary%runoff_carbon (is:ie,js:je) ) ; ice_ocean_boundary%runoff_carbon = 0.0
allocate( ice_ocean_boundary%calving (is:ie,js:je) ) ; ice_ocean_boundary%calving = 0.0
allocate( ice_ocean_boundary%runoff_hflx (is:ie,js:je) ) ; ice_ocean_boundary%runoff_hflx = 0.0
allocate( ice_ocean_boundary%calving_hflx (is:ie,js:je) ) ; ice_ocean_boundary%calving_hflx = 0.0
Expand Down Expand Up @@ -292,6 +293,9 @@ subroutine flux_ice_to_ocean ( Ice, Ocean, Ice_Ocean_Boundary )
if(ASSOCIATED(Ice_Ocean_Boundary%runoff) ) call flux_ice_to_ocean_redistribute( Ice, Ocean, &
Ice%runoff, Ice_Ocean_Boundary%runoff, Ice_Ocean_Boundary%xtype, do_area_weighted_flux )

if(ASSOCIATED(Ice_Ocean_Boundary%runoff_carbon) ) call flux_ice_to_ocean_redistribute( Ice, Ocean, &
Ice%runoff_carbon, Ice_Ocean_Boundary%runoff_carbon, Ice_Ocean_Boundary%xtype, do_area_weighted_flux )

if(ASSOCIATED(Ice_Ocean_Boundary%calving) ) call flux_ice_to_ocean_redistribute( Ice, Ocean, &
Ice%calving, Ice_Ocean_Boundary%calving, Ice_Ocean_Boundary%xtype, do_area_weighted_flux )

Expand Down Expand Up @@ -340,6 +344,7 @@ subroutine flux_ice_to_ocean_finish ( Time, Ice_Ocean_Boundary )
call fms_data_override('OCN', 'lprec', Ice_Ocean_Boundary%lprec , Time )
call fms_data_override('OCN', 'fprec', Ice_Ocean_Boundary%fprec , Time )
call fms_data_override('OCN', 'runoff', Ice_Ocean_Boundary%runoff , Time )
call fms_data_override('OCN', 'runoff_carbon', Ice_Ocean_Boundary%runoff_carbon , Time )
call fms_data_override('OCN', 'calving', Ice_Ocean_Boundary%calving , Time )
call fms_data_override('OCN', 'runoff_hflx', Ice_Ocean_Boundary%runoff_hflx , Time )
call fms_data_override('OCN', 'calving_hflx', Ice_Ocean_Boundary%calving_hflx , Time )
Expand Down
8 changes: 7 additions & 1 deletion full/land_ice_flux_exchange.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -77,11 +77,13 @@ subroutine land_ice_flux_exchange_init(Land, Ice, land_ice_boundary, Dt_cpl_in,

!allocate land_ice_boundary
allocate( land_ice_boundary%runoff(is:ie,js:je) )
allocate( land_ice_boundary%runoff_carbon(is:ie,js:je) )
allocate( land_ice_boundary%calving(is:ie,js:je) )
allocate( land_ice_boundary%runoff_hflx(is:ie,js:je) )
allocate( land_ice_boundary%calving_hflx(is:ie,js:je) )
! initialize values for override experiments (mjh)
land_ice_boundary%runoff=0.0
land_ice_boundary%runoff_carbon=0.0
land_ice_boundary%calving=0.0
land_ice_boundary%runoff_hflx=0.0
land_ice_boundary%calving_hflx=0.0
Expand Down Expand Up @@ -109,7 +111,7 @@ subroutine flux_land_to_ice( Time, Land, Ice, Land_Ice_Boundary )
!! fluxes passed from land to ice

integer :: ier
real, dimension(n_xgrid_runoff) :: ex_runoff, ex_calving, ex_runoff_hflx, ex_calving_hflx
real, dimension(n_xgrid_runoff) :: ex_runoff, ex_calving, ex_runoff_hflx, ex_calving_hflx, ex_runoff_DOC
real, dimension(size(Land_Ice_Boundary%runoff,1),size(Land_Ice_Boundary%runoff,2),1) :: ice_buf

!Balaji
Expand All @@ -121,11 +123,14 @@ subroutine flux_land_to_ice( Time, Land, Ice, Land_Ice_Boundary )

if (do_runoff) then
call fms_xgrid_put_to_xgrid ( Land%discharge, 'LND', ex_runoff, xmap_runoff)
call fms_xgrid_put_to_xgrid ( Land%discharge_DOC, 'LND', ex_runoff_DOC, xmap_runoff)
call fms_xgrid_put_to_xgrid ( Land%discharge_snow, 'LND', ex_calving, xmap_runoff)
call fms_xgrid_put_to_xgrid ( Land%discharge_heat, 'LND', ex_runoff_hflx, xmap_runoff)
call fms_xgrid_put_to_xgrid ( Land%discharge_snow_heat, 'LND', ex_calving_hflx, xmap_runoff)
call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_runoff, xmap_runoff)
Land_Ice_Boundary%runoff = ice_buf(:,:,1);
call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_runoff_DOC, xmap_runoff)
Land_Ice_Boundary%runoff_carbon = ice_buf(:,:,1);
call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_calving, xmap_runoff)
Land_Ice_Boundary%calving = ice_buf(:,:,1);
call fms_xgrid_get_from_xgrid (ice_buf, 'OCN', ex_runoff_hflx, xmap_runoff)
Expand All @@ -134,6 +139,7 @@ subroutine flux_land_to_ice( Time, Land, Ice, Land_Ice_Boundary )
Land_Ice_Boundary%calving_hflx = ice_buf(:,:,1);
!Balaji
call fms_data_override('ICE', 'runoff' , Land_Ice_Boundary%runoff , Time)
call fms_data_override('ICE', 'runoff_DOC' , Land_Ice_Boundary%runoff_carbon , Time)
call fms_data_override('ICE', 'calving', Land_Ice_Boundary%calving, Time)
call fms_data_override('ICE', 'runoff_hflx' , Land_Ice_Boundary%runoff_hflx , Time)
call fms_data_override('ICE', 'calving_hflx', Land_Ice_Boundary%calving_hflx, Time)
Expand Down
Loading