import_from_gridinfo import_from_gridinfovv volgrid6d_import_from_file volgrid6d_read_from_file

import_from_gridinfovv() subroutine import_from_gridinfovv ( type(volgrid6d), dimension(:), pointer this, type(arrayof_gridinfo), intent(in) gridinfov, integer, intent(in), optional dup_mode, logical, intent(in), optional clone, logical, intent(in), optional decode, integer, intent(in), optional time_definition, character(len=*), dimension(:), intent(in), optional anavar, character(len=*), intent(in), optional categoryappend ) Import an array of gridinfo objects into an array of volgrid6d objects. This method imports an array of gridded fields from an arrayof_gridinfo object into a suitable number of volgrid6d objects. The number of volgrid6d allocated is determined by the number of different grids encountered in arrayof_gridinfo. Unlike the import for a single gridinfo, here the volgrid6d object is a non-associated pointer to a 1-d array of uninitialized objects, and all the dimension descriptors in each of the objects are allocated and assigned within the method according to the data contained in gridinfov. If the anavar array argument is provided, all the input messages whose variable maps to one of the B-table variables contained in anavar are treated as time-independent (AKA anagraphic data, station data, etc.), thus their time and timerange are ignored and they are replicated for every time and timerange present in the corresponding data volume. Parametri this object in which to import [in] gridinfov array of gridinfo objects to be imported [in] dup_mode determines the behavior in case of duplicate metadata: if dup_mode is not provided or 0, a duplicate field overwrites, if dup_mode is 1, duplicate fields are merged with priority to the last [in] clone if provided and .TRUE. , clone the gaid's from gridinfo to this [in] decode if provided and .FALSE. the data volume in the elements of this is not allocated and successive work will be performed on grid_id's [in] time_definition 0=time is reference time; 1=time is validity time [in] anavar list of variables (B-table code equivalent) to be treated as time-independent data [in] categoryappend append this suffix to log4fortran namespace category Definizione alla linea 1264 del file volgrid6d_class.F90. !> Export a \a volgrid6d object to an \a arrayof_gridinfo object. !! The multidimensional \a volgrid6d structure is serialized into a !! one-dimensional array of gridinfo_def objects, which is allocated !! to the proper size if not already allocated, or it is extended !! keeping the old data if any. SUBROUTINE export_to_gridinfov(this, gridinfov, gaid_template, clone) TYPE(volgrid6d),INTENT(inout) :: this !< volume to be exported TYPE(arrayof_gridinfo),INTENT(inout) :: gridinfov !< output array of gridinfo_def objects TYPE(grid_id),INTENT(in),OPTIONAL :: gaid_template !< \a grid_id template to be used for output data replacing the one contained in \a this LOGICAL,INTENT(in),OPTIONAL :: clone !< if provided and \c .TRUE., clone the grid_id included in \a this rather than making a shallow copy INTEGER :: i ,itime, itimerange, ilevel, ivar INTEGER :: ntime, ntimerange, nlevel, nvar TYPE(gridinfo_def) :: gridinfol #ifdef DEBUG CALL l4f_category_log(this%category,l4f_debug,"start export_to_gridinfov") #endif ! this is necessary in order not to repeatedly and uselessly copy the ! same array of coordinates for each 2d grid during export, the ! side-effect is that the computed projection in this is lost CALL dealloc(this%griddim%dim) i=0 ntime=size(this%time) ntimerange=size(this%timerange) nlevel=size(this%level) nvar=size(this%var) DO itime=1,ntime DO itimerange=1,ntimerange DO ilevel=1,nlevel DO ivar=1,nvar CALL init(gridinfol) CALL export(this, gridinfol, itime, itimerange, ilevel, ivar, & gaid_template=gaid_template, clone=clone) IF (c_e(gridinfol%gaid)) THEN CALL insert(gridinfov, gridinfol) ELSE CALL delete(gridinfol) ENDIF ENDDO ENDDO ENDDO ENDDO END SUBROUTINE export_to_gridinfov !> Export an array of \a volgrid6d objects to an \a arrayof_gridinfo object. !! The multidimensional \a volgrid6d structures are serialized into a !! one-dimensional array of gridinfo_def objects, which is allocated !! to the proper size if not already allocated, or it is extended !! keeping the old data if any. SUBROUTINE export_to_gridinfovv(this, gridinfov, gaid_template, clone) !, & ! categoryappend) TYPE(volgrid6d),INTENT(inout) :: this(:) !< volume array to be exported TYPE(arrayof_gridinfo),INTENT(inout) :: gridinfov !< output array of gridinfo_def objects TYPE(grid_id),INTENT(in),OPTIONAL :: gaid_template !< \a grid_id template to be used for output data replacing the one contained in \a this LOGICAL,INTENT(in),OPTIONAL :: clone !< if provided and \c .TRUE., clone the grid_id included in \a this rather than making a shallow copy INTEGER :: i DO i = 1, SIZE(this) #ifdef DEBUG CALL l4f_category_log(this(i)%category,l4f_debug, & "export_to_gridinfovv grid index: "//t2c(i)) #endif CALL export(this(i), gridinfov, gaid_template=gaid_template, clone=clone) ENDDO END SUBROUTINE export_to_gridinfovv !> Import the content of a supported file (like grib or gdal-supported !! format) into an array of \a volgrid6d objects. This method imports !! a set of gridded fields from a file object into a suitable number !! of \a volgrid6d objects. The data are imported by creating a !! temporary \a gridinfo object, importing it into the \a volgrid6d !! object cloning the gaid's and then destroying the gridinfo, so it !! works similarly to volgrid6d_class::import_from_gridinfovv() method. !! For a detailed explanation of the \a anavar argument, see the !! documentation of volgrid6d_class::import_from_gridinfovv() method. SUBROUTINE volgrid6d_import_from_file(this, filename, dup_mode, decode, & time_definition, anavar, categoryappend) TYPE(volgrid6d),POINTER :: this(:) !< object in which to import CHARACTER(len=*),INTENT(in) :: filename !< name of file from which to import INTEGER,INTENT(in),OPTIONAL :: dup_mode !< determines the behavior in case of duplicate metadata: if \a dup_mode is not provided or 0, a duplicate field overwrites, if \a dup_mode is 1, duplicate fields are merged with priority to the last LOGICAL,INTENT(in),OPTIONAL :: decode !< if provided and \a .FALSE. the data volume in the elements of \a this is not allocated and successive work will be performed on grid_id's INTEGER,INTENT(IN),OPTIONAL :: time_definition !< 0=time is reference time; 1=time is validity time CHARACTER(len=*),INTENT(IN),OPTIONAL :: anavar(:) !< list of variables (B-table code equivalent) to be treated as time-independent data character(len=*),INTENT(in),OPTIONAL :: categoryappend !< append this suffix to log4fortran namespace category TYPE(arrayof_gridinfo) :: gridinfo TYPE(l4f_handle) :: category CHARACTER(len=512) :: a_name NULLIFY(this) IF (PRESENT(categoryappend))THEN CALL l4f_launcher(a_name,a_name_append= & trim(subcategory)//"."//trim(categoryappend)) ELSE CALL l4f_launcher(a_name,a_name_append=trim(subcategory)) ENDIF category=l4f_category_get_handle(a_name) CALL import(gridinfo, filename=filename, categoryappend=categoryappend) IF (gridinfo%arraysize > 0) THEN CALL import(this, gridinfo, dup_mode=dup_mode, clone=.true., decode=decode, & time_definition=time_definition, anavar=anavar, & categoryappend=categoryappend) CALL l4f_category_log(category,l4f_info,"deleting gridinfo") CALL delete(gridinfo) ELSE CALL l4f_category_log(category,l4f_info,"file does not contain gridded data") ENDIF ! close logger CALL l4f_category_delete(category) END SUBROUTINE volgrid6d_import_from_file !> High level method for exporting a volume array to file. !! All the information contained into an array of \a volgrid6d !! objects, i.e. dimension descriptors and data, is exported to a file !! using the proper output driver (typically grib_api for grib !! format). If a template is provided, it will determine the !! characteristic of the output file, otherwise the \a grid_id !! descriptors contained in the volgrid6d object will be used SUBROUTINE volgrid6d_export_to_file(this, filename, gaid_template, categoryappend) TYPE(volgrid6d) :: this(:) !< volume(s) to be exported CHARACTER(len=*),INTENT(in) :: filename !< output file name TYPE(grid_id),INTENT(in),OPTIONAL :: gaid_template !< template for the output file, if provided the grid_id information stored in the volgrid6d objects will be ignored character(len=*),INTENT(in),OPTIONAL :: categoryappend !< append this suffix to log4fortran namespace category TYPE(arrayof_gridinfo) :: gridinfo TYPE(l4f_handle) :: category CHARACTER(len=512) :: a_name IF (PRESENT(categoryappend)) THEN CALL l4f_launcher(a_name,a_name_append=trim(subcategory)//"."//trim(categoryappend)) ELSE CALL l4f_launcher(a_name,a_name_append=trim(subcategory)) ENDIF category=l4f_category_get_handle(a_name) #ifdef DEBUG CALL l4f_category_log(category,l4f_debug,"start export to file") #endif CALL l4f_category_log(category,l4f_info,"writing volgrid6d to grib file: "//trim(filename)) !IF (ASSOCIATED(this)) THEN CALL export(this, gridinfo, gaid_template=gaid_template, clone=.true.) IF (gridinfo%arraysize > 0) THEN CALL export(gridinfo, filename) CALL delete(gridinfo) ENDIF !ELSE ! CALL l4f_category_log(category,L4F_INFO,"volume volgrid6d is not associated") !ENDIF ! close logger CALL l4f_category_delete(category) END SUBROUTINE volgrid6d_export_to_file !> Array destructor for \a volgrid6d class. !! Delete an array of \a volgrid6d objects and deallocate the array !! itself. SUBROUTINE volgrid6dv_delete(this) TYPE(volgrid6d),POINTER :: this(:) !< vector of volgrid6d object INTEGER :: i