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

Feature/enkf q2 #568

Merged
merged 13 commits into from
Jun 27, 2023
Merged

Feature/enkf q2 #568

Changes from 6 commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
0a6d8c4
trivial changes in ush/build.sh to pull from upstream
TingLei-daprediction Dec 21, 2022
1642e55
Merge branch 'NOAA-EMC:develop' into develop
TingLei-daprediction Feb 14, 2023
3f0d680
Merge branch 'NOAA-EMC:develop' into develop
TingLei-daprediction Mar 24, 2023
dee603d
Merge branch 'NOAA-EMC:develop' into develop
TingLei-daprediction May 3, 2023
243d01c
GitHub Issue NOAA-EMC/GSI#566 A quick enhancement for clipping negati…
TingLei-daprediction May 4, 2023
7385d79
GitHub Issue NOAA-EMC/GSI#566 A quick enhancement for clipping negati…
TingLei-daprediction May 4, 2023
2a82e40
GitHub Issue NOAA-EMC/GSI#566 A quick enhancement for clipping negati…
TingLei-daprediction May 4, 2023
027c06d
Merge branch 'feature/enkf_q2' of https://github.com/TingLei-dapredic…
TingLei-daprediction May 5, 2023
620bf3a
Merge branch 'NOAA-EMC:develop' into develop
TingLei-daprediction May 15, 2023
a67dad1
Merge branch 'develop' into feature/enkf_q2
TingLei-daprediction May 15, 2023
5328776
Merge branch 'NOAA-EMC:develop' into develop
TingLei-daprediction Jun 2, 2023
ccf2c90
Merge branch 'develop' into feature/enkf_q2
TingLei-daprediction Jun 2, 2023
40a6fd0
Merge branch 'NOAA-EMC:develop' into feature/enkf_q2
TingLei-daprediction Jun 8, 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
16 changes: 16 additions & 0 deletions src/enkf/gridio_fv3reg.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -853,6 +853,10 @@ subroutine writegriddata(nanal1,nanal2,vars3d,vars2d,n3d,n2d,levels,ndim,vargrid
varstrname = 'sphum'
call fv3lamfile%get_idfn(varstrname,file_id,fv3filename)
call read_fv3_restart_data3d(varstrname,fv3filename,file_id,qworkvar3d)
if ( cliptracers ) then ! set cliptracers to remove negative hydrometers
clip = tiny(qworkvar3d(1,1,1))
where (qworkvar3d < clip) qworkvar3d = clip
end if
tvworkvar3d=tsenworkvar3d*(one+fv*qworkvar3d)
tvworkvar3d=tvworkvar3d+workinc3d
if(q_ind > 0) then
Expand All @@ -866,6 +870,10 @@ subroutine writegriddata(nanal1,nanal2,vars3d,vars2d,n3d,n2d,levels,ndim,vargrid
enddo
enddo
qworkvar3d=qworkvar3d+workinc3d
if ( cliptracers ) then ! set cliptracers to remove negative hydrometers
clip = tiny(qworkvar3d(1,1,1))
where (qworkvar3d < clip) qworkvar3d = clip
end if
endif
tsenworkvar3d=tvworkvar3d/(one+fv*qworkvar3d)
varstrname = 'T'
Expand Down Expand Up @@ -2173,6 +2181,10 @@ subroutine writegriddata_pnc(vars3d,vars2d,n3d,n2d,levels,ndim,vargrid,no_inflat
enddo

if(iope ==0 ) then
if ( cliptracers ) then ! set cliptracers to remove negative hydrometers
clip = tiny(qworkvar3d(1,1,1))
where (qworkvar3d < clip) qworkvar3d = clip
end if
tvworkvar3d=tsenworkvar3d*(one+fv*qworkvar3d)
tvworkvar3d=tvworkvar3d+workinc3d
if(q_ind > 0) then
Expand All @@ -2186,6 +2198,10 @@ subroutine writegriddata_pnc(vars3d,vars2d,n3d,n2d,levels,ndim,vargrid,no_inflat
enddo
enddo
qworkvar3d=qworkvar3d+workinc3d
if ( cliptracers ) then ! set cliptracers to remove negative hydrometers
clip = tiny(qworkvar3d(1,1,1))
where (qworkvar3d < clip) qworkvar3d = clip
end if
endif
tsenworkvar3d=tvworkvar3d/(one+fv*qworkvar3d)
endif
Expand Down