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

Add valid time to the 6 tile cubed sphere grid history files #550

Merged
merged 5 commits into from
Jun 24, 2022
Merged

Add valid time to the 6 tile cubed sphere grid history files #550

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
e6895cc
Add 'time_iso' variable to 6-tile history files
DusanJovic-NOAA Jun 8, 2022
78224b5
add '_Encoding' attribute to time_iso
jswhit2 Mar 29, 2022
a80a13e
add '_Encoding' attribute to time_iso
DusanJovic-NOAA Jun 8, 2022
a9fe9c1
Merge remote-tracking branch 'origin/develop' into time_iso
DusanJovic-NOAA Jun 9, 2022
520801c
Merge remote-tracking branch 'origin/develop' into time_iso
DusanJovic-NOAA Jun 21, 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_netcdf.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -236,6 +236,7 @@ subroutine write_netcdf(wrtfb, filename, &
ncerr = nf90_def_var(ncid, "time_iso", NF90_CHAR, [ch_dimid,time_dimid], timeiso_varid); NC_ERR_STOP(ncerr)
ncerr = nf90_put_att(ncid, timeiso_varid, "long_name", "valid time"); NC_ERR_STOP(ncerr)
ncerr = nf90_put_att(ncid, timeiso_varid, "description", "ISO 8601 datetime string"); NC_ERR_STOP(ncerr)
ncerr = nf90_put_att(ncid, timeiso_varid, "_Encoding", "UTF-8"); NC_ERR_STOP(ncerr)

! coordinate variable attributes based on output_grid type
if (trim(output_grid(grid_id)) == 'gaussian_grid' .or. &
Expand Down
21 changes: 21 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 @@ -3130,6 +3130,8 @@ subroutine ioCompRun(comp, importState, exportState, clock, rc)
real(ESMF_KIND_R4) :: valueR4
real(ESMF_KIND_R8) :: valueR8
logical :: thereAreVerticals
integer :: ch_dimid, timeiso_varid
character(len=ESMF_MAXSTR) :: time_iso

rc = ESMF_SUCCESS

Expand Down Expand Up @@ -3271,6 +3273,10 @@ subroutine ioCompRun(comp, importState, exportState, clock, rc)
name="time", value=time, rc=rc)
if (ESMF_LogFoundError(rcToCheck=rc, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__)) return

call ESMF_AttributeGet(grid, convention="NetCDF", purpose="FV3", &
name="time_iso", value=time_iso, rc=rc)
if (ESMF_LogFoundError(rcToCheck=rc, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__)) return

ncerr = nf90_redef(ncid=ncid)
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return

Expand All @@ -3285,6 +3291,17 @@ subroutine ioCompRun(comp, importState, exportState, clock, rc)
dimids=(/dimid/), varid=varid)
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return

ncerr = nf90_def_dim(ncid, "nchars", 20, ch_dimid)
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return
ncerr = nf90_def_var(ncid, "time_iso", NF90_CHAR, [ch_dimid,dimid], timeiso_varid)
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return
ncerr = nf90_put_att(ncid, timeiso_varid, "long_name", "valid time")
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return
ncerr = nf90_put_att(ncid, timeiso_varid, "description", "ISO 8601 datetime string")
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return
ncerr = nf90_put_att(ncid, timeiso_varid, "_Encoding", "UTF-8")
if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return

ncerr = nf90_enddef(ncid=ncid)

if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return
Expand All @@ -3293,6 +3310,10 @@ subroutine ioCompRun(comp, importState, exportState, clock, rc)

if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return

ncerr = nf90_put_var(ncid, timeiso_varid, values=[trim(time_iso)])

if (ESMF_LogFoundNetCDFError(ncerr, msg=ESMF_LOGERR_PASSTHRU, line=__LINE__, file=__FILE__, rcToReturn=rc)) return

! loop over all the grid attributes that start with "time:", and
! put them on the "time" variable in the NetCDF file

Expand Down