|
libsim Versione 8.0.0
|
◆ grid_id_display()
Display on standard output a description of the grid_id object provided. Also the grib key names and values are printed; the set of keys returned can be controlled with the input variable namespace. Available namespaces are
Definizione alla linea 840 del file grid_id_class.F90. 841 RETURN
842 ENDIF
843
844else if (editionnumber /= 1) then
845
846 CALL l4f_log(l4f_error, &
847 "grib_api GribEditionNumber not supported: "//t2c(editionnumber))
848 call raise_error()
849 field(:,:) = rmiss
850 RETURN
851
852end if
853
854CALL grib_get(gaid,'iScansNegatively',iscansnegatively,ierr)
855IF (ierr /= grib_success) iscansnegatively=0
856CALL grib_get(gaid,'jScansPositively',jscanspositively,ierr)
857IF (ierr /= grib_success) jscanspositively=1
858CALL grib_get(gaid,'jPointsAreConsecutive',jpointsareconsecutive,ierr)
859IF (ierr /= grib_success) jpointsareconsecutive=0
860
861call grib_get(gaid,'numberOfPoints',numberofpoints)
862
863IF (numberofpoints /= SIZE(field)) THEN
864 CALL l4f_log(l4f_error, 'grid_id_decode_data_gribapi numberOfPoints and grid size different')
865 CALL l4f_log(l4f_error, 'grid_id_decode_data_gribapi numberOfPoints: ' &
866 //t2c(numberofpoints)//', nx,ny: '&
867 //t2c(SIZE(field,1))//','//t2c(SIZE(field,2)))
868 CALL raise_error()
869 field(:,:) = rmiss
870 RETURN
871ENDIF
872
873! get data values
874#ifdef DEBUG
875call grib_get(gaid,'numberOfValues',numberofvalues)
876call l4f_log(l4f_info,'grib_api number of values: '//to_char(numberofvalues))
877call l4f_log(l4f_info,'grib_api number of points: '//to_char(numberofpoints))
878#endif
879
880CALL grib_set(gaid,'missingValue',rmiss)
881
882IF (iscansnegatively == 0 .AND. jscanspositively == 0 .AND. jpointsareconsecutive == 0) THEN
883
884 fieldp(1:SIZE(field)) => field
|