Skip to content

swap signs on imported fields #4

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
BinLiu-NOAA merged 3 commits into from
Aug 11, 2023
Merged

swap signs on imported fields #4

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
c772b80
reverse sign for taux,tauy to match atm export
DeniseWorthen Jun 21, 2023
55b61c2
change sign for sensible heat to match ATM export
DeniseWorthen Jun 22, 2023
81edd3c
swap sign of latent heat to match ATM export
DeniseWorthen Jun 26, 2023
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
15 changes: 4 additions & 11 deletions NUOPC/hycom_couple.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -538,10 +538,7 @@ subroutine import_to_hycom_deb(tlb, tub, impData, fill_value, fieldName, &
do i=1, ii
! imp_taux(i,j,1)=mgrid(i,j)
if (ishlf(i,j).eq.1) then
! The momentum flux received from CMEPS/FV3ATM are from sea to air (friction to air).
! However, the momentum flux expected and used in HYCOM should be air to sea (wind stress).
! Thus, we need to reverse the sign when retrieving taux/y.
imp_taux(i,j,1)=-impData(i+i0,j+j0)
imp_taux(i,j,1)=impData(i+i0,j+j0)
else
imp_taux(i,j,1)=0.0
endif
Expand All @@ -557,10 +554,7 @@ subroutine import_to_hycom_deb(tlb, tub, impData, fill_value, fieldName, &
do j=1, jja
do i=1, ii
if (ishlf(i,j).eq.1) then
! The momentum flux received from CMEPS/FV3ATM are from sea to air (friction to air).
! However, the momentum flux expected and used in HYCOM should be air to sea (wind stress).
! Thus, we need to reverse the sign when retrieving taux/y.
imp_tauy(i,j,1)=-impData(i+i0,j+j0)
imp_tauy(i,j,1)=impData(i+i0,j+j0)
if (impData(i+i0,j+j0).ne.fill_value) then
! rotate taux and tauy to (x,y)ward
! assumes rotation only needed for impData
Expand Down Expand Up @@ -882,7 +876,7 @@ subroutine import_to_hycom_deb(tlb, tub, impData, fill_value, fieldName, &
do j=1, jja
do i=1, ii
if (ishlf(i,j).eq.1) then
imp_latflx(i,j,1)=impData(i+i0,j+j0)
imp_latflx(i,j,1)=-impData(i+i0,j+j0)
else
imp_latflx(i,j,1)=0.0
endif
Expand All @@ -898,7 +892,7 @@ subroutine import_to_hycom_deb(tlb, tub, impData, fill_value, fieldName, &
do j=1, jja
do i=1, ii
if (ishlf(i,j).eq.1) then
imp_sensflx(i,j,1)=impData(i+i0,j+j0)
imp_sensflx(i,j,1)=-impData(i+i0,j+j0)
else
imp_sensflx(i,j,1)=0.0
endif
Expand Down Expand Up @@ -1289,4 +1283,3 @@ end subroutine hycom_couple_final
!===============================================================================
end module hycom_couple
!===============================================================================