Skip to content

Include surface ocean current in the computation of air-sea fluxes #311

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

Merged
bensonr merged 6 commits into from
Feb 1, 2024
Merged

Include surface ocean current in the computation of air-sea fluxes #311

Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
1ea4e44
Send sea surface current field to nested grids.
binli2337 Dec 27, 2023
5906ba3
Update driver/fvGFS/atmosphere.F90.
binli2337 Dec 30, 2023
cc50611
Change the variable name for zonal ocean current from ssu to usfco.
binli2337 Jan 5, 2024
64ba0b0
Merge branch 'dev/emc' into feature/hafsv2_baseline_with_ssc
binli2337 Jan 18, 2024
005a1b0
Revise fv_arrays.F90.
binli2337 Jan 23, 2024
15e9fe7
Merge branch 'dev/emc' into feature/hafsv2_baseline_with_ssc
binli2337 Jan 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
66 changes: 66 additions & 0 deletions driver/fvGFS/atmosphere.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2390,6 +2390,72 @@ subroutine atmosphere_fill_nest_cpl(Atm_block, IPD_control, IPD_data)
endif
endif

! Deal with usfco and vsfco
if (IPD_control%cplocn2atm .and. IPD_control%icplocn2atm==1) then
! Extract the coupling field
do nb = 1,Atm_block%nblks
blen = Atm_block%blksz(nb)
do ix = 1, blen
i = Atm_block%index(nb)%ii(ix)
j = Atm_block%index(nb)%jj(ix)
Atm(mygrid)%parent2nest_2d(i,j) = IPD_Data(nb)%Sfcprop%usfco(ix)
enddo
enddo
! Loop through and fill all nested grids
do n=2,ngrids
if (n==mygrid .or. mygrid==Atm(n)%parent_grid%grid_number) then
call fill_nested_grid_cpl(n, n==mygrid)
endif
enddo
! Update the nested grids
if (Atm(mygrid)%neststruct%nested) then
do nb = 1,Atm_block%nblks
blen = Atm_block%blksz(nb)
do ix = 1, blen
i = Atm_block%index(nb)%ii(ix)
j = Atm_block%index(nb)%jj(ix)
if (IPD_data(nb)%Sfcprop%oceanfrac(ix) > 0.) then
IPD_data(nb)%Sfcprop%usfco(ix) = Atm(mygrid)%parent2nest_2d(i,j)
else
IPD_data(nb)%Sfcprop%usfco(ix) = 0.0_kind_phys
endif
enddo
enddo
endif

! Extract the coupling field
do nb = 1,Atm_block%nblks
blen = Atm_block%blksz(nb)
do ix = 1, blen
i = Atm_block%index(nb)%ii(ix)
j = Atm_block%index(nb)%jj(ix)
Atm(mygrid)%parent2nest_2d(i,j) = IPD_Data(nb)%Sfcprop%vsfco(ix)
enddo
enddo
! Loop through and fill all nested grids
do n=2,ngrids
if (n==mygrid .or. mygrid==Atm(n)%parent_grid%grid_number) then
call fill_nested_grid_cpl(n, n==mygrid)
endif
enddo
! Update the nested grids
if (Atm(mygrid)%neststruct%nested) then
do nb = 1,Atm_block%nblks
blen = Atm_block%blksz(nb)
do ix = 1, blen
i = Atm_block%index(nb)%ii(ix)
j = Atm_block%index(nb)%jj(ix)
if (IPD_data(nb)%Sfcprop%oceanfrac(ix) > 0.) then
IPD_data(nb)%Sfcprop%vsfco(ix) = Atm(mygrid)%parent2nest_2d(i,j)
else
IPD_data(nb)%Sfcprop%vsfco(ix) = 0.0_kind_phys
endif
enddo
enddo
endif

endif

! Deal with zorlwav (sea surface roughness length)
if (IPD_control%cplwav2atm) then
! Extract the coupling field
Expand Down