gridinfo_init

gridinfo_init() subroutine gridinfo_init ( type(gridinfo_def), intent(out)  this, type(grid_id), intent(in), optional  gaid, type(griddim_def), intent(in), optional  griddim, type(datetime), intent(in), optional  time, type(vol7d_timerange), intent(in), optional  timerange, type(vol7d_level), intent(in), optional  level, type(volgrid6d_var), intent(in), optional  var, logical, intent(in), optional  clone, character(len=*), intent(in), optional  categoryappend  ) Constructor, it creates a new instance of the object. All the additional parameters are optional and they will be initialised to the corresponding missing value if not provided. Parametri [out] this object to be initialized [in] gaid identificator of the grid to be described [in] griddim grid descriptor [in] time time dimension descriptor [in] timerange timerange (forecast, analysis, statistically processed) dimension descriptor [in] level vertical level dimension descriptor [in] var physical variable dimension descriptor [in] clone if provided and .TRUE., the gaid will be cloned and not simply copied into the gridinfo object [in] categoryappend append this suffix to log4fortran namespace category Definizione alla linea 702 del file gridinfo_class.F90. CALL grib_get(gaid,'GRIBEditionNumber',EditionNumber) IF (editionnumber == 1) THEN CALL code_referencetime(this) ELSE IF (editionnumber == 2 )THEN IF (timerange%p1 >= timerange%p2) THEN ! forecast-like CALL code_referencetime(this) ELSE IF (timerange%p1 == 0) THEN ! analysis-like ! ready for coding with general convention CALL getenv('LIBSIM_G2COSMO_BEHAVIOR', env_var) g2cosmo_behavior = len_trim(env_var) > 0 CALL grib_get(gaid,'centre',centre) IF (g2cosmo_behavior .AND. centre == 78) THEN ! DWD analysis exception CALL code_referencetime(this) ELSE ! cosmo or old simc convention CALL code_referencetime(this-timedelta_new(sec=timerange%p2)) ENDIF ELSE ! bad timerange CALL l4f_log( l4f_error, 'Timerange with 0>p1>p2 cannot be exported in grib2') CALL raise_error() ENDIF ELSE CALL l4f_log(l4f_error,'GribEditionNumber '//t2c(editionnumber)//' not supported') CALL raise_error() ENDIF CONTAINS SUBROUTINE code_referencetime(reftime) TYPE(datetime),INTENT(in) :: reftime INTEGER :: date,time CHARACTER(len=17) :: date_time ! datetime is AAAAMMGGhhmmssmsc CALL getval(reftime, simpledate=date_time) READ(date_time(:8),'(I8)')date READ(date_time(9:12),'(I4)')time CALL grib_set(gaid,'dataDate',date) CALL grib_set(gaid,'dataTime',time) END SUBROUTINE code_referencetime END SUBROUTINE time_export_gribapi SUBROUTINE level_import_gribapi(this, gaid) TYPE(vol7d_level),INTENT(out) :: this ! vol7d_level object INTEGER,INTENT(in) :: gaid ! grib_api id of the grib loaded in memory to import INTEGER :: EditionNumber,level1,l1,level2,l2 INTEGER :: ltype,ltype1,scalef1,scalev1,ltype2,scalef2,scalev2 call grib_get(gaid,'GRIBEditionNumber',EditionNumber) if (editionnumber == 1)then call grib_get(gaid,'indicatorOfTypeOfLevel',ltype) call grib_get(gaid,'topLevel',l1) call grib_get(gaid,'bottomLevel',l2)