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

Jump to bottom

Upp inline post for RRFS and HAFS #567

Merged
merged 5 commits into from
Aug 18, 2022
Merged

Upp inline post for RRFS and HAFS #567

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
f4032d0
1)update upp revision to 8ed92cb; 2)read instant model top longwave;
WenMeng-NOAA Jul 30, 2022
37b3287
Correct a tybo
WenMeng-NOAA Jul 30, 2022
2ec1826
Merge remote-tracking branch 'upstream/develop' into upp_update_aug
WenMeng-NOAA Aug 10, 2022
0c95ee9
Merge remote-tracking branch 'upstream/develop' into upp_update_aug
WenMeng-NOAA Aug 15, 2022
68cc00a
Merge remote-tracking branch 'upstream/develop' into upp_update_aug
WenMeng-NOAA Aug 16, 2022
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
1 change: 1 addition & 0 deletions io/module_write_internal_state.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,6 +49,7 @@ module write_internal_state
integer :: lat_start, lon_start
integer :: lat_end, lon_end
real :: latstart, latlast, lonstart, lonlast
real :: latse, latnw, lonse, lonnw
integer,dimension(:),allocatable :: lat_start_wrtgrp, lon_start_wrtgrp
integer,dimension(:),allocatable :: lat_end_wrtgrp, lon_end_wrtgrp
real,dimension(:,:),allocatable :: lonPtr, latPtr
Expand Down
24 changes: 24 additions & 0 deletions io/module_wrt_grid_comp.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -756,6 +756,30 @@ subroutine wrt_initialize_p1(wrt_comp, imp_state_write, exp_state_write, clock,
enddo
enddo
wrt_int_state%post_maptype = 207
rot_lon = lon1(n)
rot_lat = lat1(n)
call rtll(rot_lon, rot_lat, geo_lon, geo_lat, dble(cen_lon(n)), dble(cen_lat(n)))
if (geo_lon < 0.0) geo_lon = geo_lon + 360.0
wrt_int_state%lonstart = geo_lon
wrt_int_state%latstart = geo_lat
rot_lon = lon2(n)
rot_lat = lat1(n)
call rtll(rot_lon, rot_lat, geo_lon, geo_lat, dble(cen_lon(n)), dble(cen_lat(n)))
if (geo_lon < 0.0) geo_lon = geo_lon + 360.0
wrt_int_state%lonse = geo_lon
wrt_int_state%latse = geo_lat
rot_lon = lon1(n)
rot_lat = lat2(n)
call rtll(rot_lon, rot_lat, geo_lon, geo_lat, dble(cen_lon(n)), dble(cen_lat(n)))
if (geo_lon < 0.0) geo_lon = geo_lon + 360.0
wrt_int_state%lonnw = geo_lon
wrt_int_state%latnw = geo_lat
rot_lon = lon2(n)
rot_lat = lat2(n)
call rtll(rot_lon, rot_lat, geo_lon, geo_lat, dble(cen_lon(n)), dble(cen_lat(n)))
if (geo_lon < 0.0) geo_lon = geo_lon + 360.0
wrt_int_state%lonlast = geo_lon
wrt_int_state%latlast = geo_lat
else if ( trim(output_grid(n)) == 'rotated_latlon_moving' ) then
! Do not compute lonPtr, latPtr here. Will be done in the run phase
wrt_int_state%post_maptype = 207
Expand Down
45 changes: 39 additions & 6 deletions io/post_fv3.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,6 +37,9 @@ subroutine post_run_fv3(wrt_int_state,mypei,mpicomp,lead_write, &
! 6)read 3D cloud fraction from cld_amt for GFDL MP,
! and from cldfra for other MPs.
! Jun 2022 J. Meng 2D decomposition
! Jul 2022 W. Meng 1)output lat/lon of four corner point for rotated
! lat-lon grid.
! 2)read instant model top logwave
!
!-----------------------------------------------------------------------
!*** run post on write grid comp
Expand Down Expand Up @@ -249,7 +252,8 @@ subroutine post_getattr_fv3(wrt_int_state,grid_id)
lonstartv, lonlastv, cenlonv, &
latstartv, latlastv, cenlatv, &
latstart_r,latlast_r,lonstart_r, &
lonlast_r, STANDLON, maptype, gridtype
lonlast_r, STANDLON, maptype, gridtype, &
latse,lonse,latnw,lonnw
!
implicit none
!
Expand Down Expand Up @@ -356,6 +360,14 @@ subroutine post_getattr_fv3(wrt_int_state,grid_id)
lonstart_r = lonstart
latlast_r = latlast
lonlast_r = lonlast
lonstart = nint(wrt_int_state%lonstart*gdsdegr)
latstart = nint(wrt_int_state%latstart*gdsdegr)
lonse = nint(wrt_int_state%lonse*gdsdegr)
latse = nint(wrt_int_state%latse*gdsdegr)
lonnw = nint(wrt_int_state%lonnw*gdsdegr)
latnw = nint(wrt_int_state%latnw*gdsdegr)
lonlast = nint(wrt_int_state%lonlast*gdsdegr)
latlast = nint(wrt_int_state%latlast*gdsdegr)

if(dlon(grid_id)<spval) then
dxval = dlon(grid_id)*gdsdegr
Expand Down Expand Up @@ -530,7 +542,7 @@ subroutine set_postvars_fv3(wrt_int_state,mpicomp,setvar_atmfile, &
dxval, dyval, truelat2, truelat1, psmapf, cenlat, &
lonstartv, lonlastv, cenlonv, latstartv, latlastv, &
cenlatv,latstart_r,latlast_r,lonstart_r,lonlast_r, &
maptype, gridtype, STANDLON
maptype, gridtype, STANDLON,latse,lonse,latnw,lonnw
use lookup_mod, only: thl, plq, ptbl, ttbl, rdq, rdth, rdp, rdthe, pl, &
qs0, sqs, sthe, ttblq, rdpq, rdtheq, stheq, the0q, the0
use physcons, only: grav => con_g, fv => con_fvirt, rgas => con_rd, &
Expand Down Expand Up @@ -571,7 +583,7 @@ subroutine set_postvars_fv3(wrt_int_state,mpicomp,setvar_atmfile, &
real,external::FPVSNEW
real,dimension(:,:),allocatable :: dummy, p2d, t2d, q2d, qs2d, &
cw2d, cfr2d
character(len=80) :: fieldname, wrtFBName
character(len=80) :: fieldname, wrtFBName, flatlon
type(ESMF_Grid) :: wrtGrid
type(ESMF_Field) :: theField
type(ESMF_Field), allocatable :: fcstField(:)
Expand Down Expand Up @@ -714,22 +726,20 @@ subroutine set_postvars_fv3(wrt_int_state,mpicomp,setvar_atmfile, &
! time averaged cloud fraction, set acfrst to spval, ncfrst to 1
! UNDERGROUND RUNOFF, bgroff
! inst incoming sfc longwave
! inst model top outgoing longwave,rlwtoa
! inst incoming sfc shortwave, rswin
! inst incoming clear sky sfc shortwave, rswinc
! inst outgoing sfc shortwave, rswout
! snow phase change heat flux, snopcx
! GFS does not use total momentum flux,sfcuvx
!$omp parallel do default(none),private(i,j),shared(jsta,jend,im,spval,ista,iend), &
!$omp& shared(acfrcv,ncfrcv,acfrst,ncfrst,bgroff,rlwtoa,rswin,rswinc,rswout,snopcx,sfcuvx)
!$omp& shared(acfrcv,ncfrcv,acfrst,ncfrst,bgroff,rswin,rswinc,rswout,snopcx,sfcuvx)
do j=jsta,jend
do i=ista,iend
acfrcv(i,j) = spval
ncfrcv(i,j) = 1.0
acfrst(i,j) = spval
ncfrst(i,j) = 1.0
bgroff(i,j) = spval
rlwtoa(i,j) = spval
rswinc(i,j) = spval
enddo
enddo
Expand Down Expand Up @@ -1678,6 +1688,17 @@ subroutine set_postvars_fv3(wrt_int_state,mpicomp,setvar_atmfile, &
enddo
endif

! outgoing model top logwave
if(trim(fieldname)=='ulwrf_toa') then
!$omp parallel do default(none) private(i,j) shared(jsta,jend,ista,iend,rlwtoa,arrayr42d,fillValue,spval)
do j=jsta,jend
do i=ista, iend
rlwtoa(i,j) = arrayr42d(i,j)
if( abs(arrayr42d(i,j)-fillValue)< small) rlwtoa(i,j) = spval
enddo
enddo
endif

! time averaged incoming sfc shortwave
if(trim(fieldname)=='dswrf_ave') then
!$omp parallel do default(none) private(i,j) shared(jsta,jend,ista,iend,aswin,arrayr42d,fillValue,spval)
Expand Down Expand Up @@ -3265,6 +3286,18 @@ subroutine set_postvars_fv3(wrt_int_state,mpicomp,setvar_atmfile, &
!
!more fields need to be computed
!

! write lat/lon of the four corner point for rotated lat-lon grid
if(mype == 0 .and. maptype == 207)then
write(flatlon,1001)ifhr
open(112,file=trim(flatlon),form='formatted',status='unknown')
jkbk2004 marked this conversation as resolved.
Show resolved Hide resolved
write(112,1002)latstart/1000,lonstart/1000,&
latse/1000,lonse/1000,latnw/1000,lonnw/1000, &
latlast/1000,lonlast/1000
1001 format('latlons_corners.txt.f',I3.3)
1002 format(4(I6,I7,X))
close(112)
endif
end subroutine set_postvars_fv3


Expand Down