core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Testing: Fortran key setting

This commit is contained in:
Shahram Najm 2024-01-10 21:05:38 +00:00
parent 45f47d1b90
commit d3e7e2ecb2
1 changed files with 10 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
examples/F90/grib_set_gvc.f90
View File

@ -15,8 +15,12 @@
program set program set
use eccodes use eccodes
implicit none implicit none
integer :: infile, outfile integer :: infile, outfile
integer :: igrib integer :: igrib
real(4) :: nlev_real4 = 12.21
real(8) :: nlev_real8 = 12.21
integer(4) :: PDTN_int = 11
integer(8) :: PDTN_long = 11
call codes_open_file(infile, '../../data/sample.grib2', 'r') call codes_open_file(infile, '../../data/sample.grib2', 'r')
@ -25,13 +29,15 @@ program set
call codes_grib_new_from_file(infile, igrib) call codes_grib_new_from_file(infile, igrib)
! individual ensemble forecast ! individual ensemble forecast
call codes_set(igrib, 'productDefinitionTemplateNumber', 11) call codes_set(igrib, 'productDefinitionTemplateNumber', PDTN_int)
call codes_set(igrib, 'productDefinitionTemplateNumber', PDTN_long)
! select level type as Generalized Vertical Height Coordinate ! select level type as Generalized Vertical Height Coordinate
call codes_set(igrib, 'typeOfLevel', 'generalVertical') call codes_set(igrib, 'typeOfLevel', 'generalVertical')
! now set keys specific to this level type ! now set keys specific to this level type
call codes_set(igrib, 'nlev', 12.21) call codes_set(igrib, 'nlev', nlev_real4)
call codes_set(igrib, 'nlev', nlev_real8)
call codes_set(igrib, 'numberOfVGridUsed', 13.55) call codes_set(igrib, 'numberOfVGridUsed', 13.55)
! check integrity of GRIB message ! check integrity of GRIB message