import_from_gridinfo import_from_gridinfovv volgrid6d_import_from_file volgrid6d_read_from_file

import_from_gridinfo() subroutine import_from_gridinfo ( type(volgrid6d), intent(inout) this, type(gridinfo_def), intent(in) gridinfo, logical, intent(in), optional force, integer, intent(in), optional dup_mode, logical, intent(in), optional clone, logical, intent(in), optional isanavar ) Import a single gridinfo object into a volgrid6d object. This methods imports a single gridded field from a gridinfo object into a volgrid6d object, inserting it into the multidimensional structure of volgrid6d. The volgrid6d object must have been already initialized and the dimensions specified with volgrid6d_alloc(). If the force argument is missing or .FALSE. , the volgrid6d object dimension descriptors (time, timerange, vertical level, physical variable) must already have space for the corresponding values coming from gridinfo, otherwise the object will be rejected; this means that all the volgrid6d dimension descriptors should be correctly assigned. If force is .TRUE. , the gridinfo dimension descriptors that do not fit into available descriptors in the volgrid6d structure, will be accomodated in a empty (i.e. equal to missing value) descriptor, if available, otherwise the gridinfo will be rejected. The descriptor of the grid in the volgrid object is assigned to the descriptor contained in gridinfo if it is missing in volgrid, otherwise it is checked and the object is rejected if grids do not match. Parametri [in,out] this object in which to import [in] gridinfo gridinfo object to be imported [in] force if provided and .TRUE., the gridinfo is forced into an empty element of this, if required and possible [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] isanavar if provides and .TRUE., the gridinfo object is treated as time-independent and replicated for every time and timerange Definizione alla linea 1014 del file volgrid6d_class.F90. CALL copy(gaid_template, gaid) #ifdef DEBUG CALL l4f_category_log(this%category,l4f_debug,"template cloned to a new gaid") #endif usetemplate = c_e(gaid) ENDIF IF (.NOT.usetemplate) THEN IF (optio_log(clone)) THEN CALL copy(this%gaid(ilevel,itime,itimerange,ivar), gaid) #ifdef DEBUG CALL l4f_category_log(this%category,l4f_debug,"original gaid cloned to a new one") #endif ELSE gaid = this%gaid(ilevel,itime,itimerange,ivar) ENDIF ENDIF IF (this%time_definition == 1 .OR. this%time_definition == 2) THEN correctedtime = this%time(itime) - & timedelta_new(sec=this%timerange(itimerange)%p1) ELSE correctedtime = this%time(itime) ENDIF CALL init(gridinfo,gaid, this%griddim, correctedtime, this%timerange(itimerange), & this%level(ilevel), this%var(ivar)) ! reset the gridinfo, bad but necessary at this point for encoding the field CALL export(gridinfo%griddim, gridinfo%gaid) ! encode the field IF (ASSOCIATED(this%voldati)) THEN CALL encode_gridinfo(gridinfo, this%voldati(:,:,ilevel,itime,itimerange,ivar)) ELSE IF (usetemplate) THEN ! field must be forced into template in this case NULLIFY(voldati) CALL volgrid_get_vol_2d(this, ilevel, itime, itimerange, ivar, voldati) CALL encode_gridinfo(gridinfo, voldati) DEALLOCATE(voldati) ENDIF END SUBROUTINE export_to_gridinfo !> Import an array of \a gridinfo objects into an array of \a !! volgrid6d objects. !! This method imports an array of gridded fields from an \a !! arrayof_gridinfo object into a suitable number of \a volgrid6d !! objects. The number of \a volgrid6d allocated is determined by the !! number of different grids encountered in \a arrayof_gridinfo. !! Unlike the import for a single \a gridinfo, here the \a 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 \a gridinfov. !! If the \a anavar array argument is provided, all the input messages !! whose variable maps to one of the B-table variables contained in \a !! 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. SUBROUTINE import_from_gridinfovv(this, gridinfov, dup_mode, clone, decode, & time_definition, anavar, categoryappend) TYPE(volgrid6d),POINTER :: this(:) !< object in which to import TYPE(arrayof_gridinfo),INTENT(in) :: gridinfov !< array of gridinfo objects to be imported 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 :: clone !< if provided and \c .TRUE. , clone the gaid's from \a gridinfo to \a this 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 INTEGER :: i, j, stallo INTEGER :: ngrid, ntime, ntimerange, nlevel, nvar, ltime_definition TYPE(l4f_handle) :: category CHARACTER(len=512) :: a_name TYPE(datetime),ALLOCATABLE :: correctedtime(:) LOGICAL,ALLOCATABLE :: isanavar(:) TYPE(vol7d_var) :: lvar TYPE(vol7d_timerange),ALLOCATABLE :: correctedtimerange(:) ! category temporanea (altrimenti non possiamo loggare) 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 import_from_gridinfovv") #endif IF (PRESENT(time_definition)) THEN ltime_definition = max(min(time_definition, 2), 0) ELSE ltime_definition = 0 ENDIF ngrid=count_distinct(gridinfov%array(1:gridinfov%arraysize)%griddim,back=.true.) CALL l4f_category_log(category,l4f_info, t2c(ngrid)// & ' different grid definition(s) found in input data') ALLOCATE(this(ngrid),stat=stallo) IF (stallo /= 0)THEN CALL l4f_category_log(category,l4f_fatal,"allocating memory") CALL raise_fatal_error() ENDIF DO i = 1, ngrid IF (PRESENT(categoryappend))THEN CALL init(this(i), time_definition=ltime_definition, categoryappend=trim(categoryappend)//"-vol"//t2c(i)) ELSE CALL init(this(i), time_definition=ltime_definition, categoryappend="vol"//t2c(i)) ENDIF ENDDO this(:)%griddim=pack_distinct(gridinfov%array(1:gridinfov%arraysize)%griddim, & ngrid, back=.true.) ! mark elements as ana variables (time-independent) ALLOCATE(isanavar(gridinfov%arraysize)) isanavar(:) = .false. IF (PRESENT(anavar)) THEN DO i = 1, gridinfov%arraysize DO j = 1, SIZE(anavar) lvar = convert(gridinfov%array(i)%var) IF (lvar%btable == anavar(j)) THEN isanavar(i) = .true. EXIT ENDIF ENDDO ENDDO CALL l4f_category_log(category,l4f_info,t2c(count(isanavar))//'/'// & t2c(gridinfov%arraysize)//' constant-data messages found in input data') ENDIF IF (ltime_definition == 1 .OR. ltime_definition == 2) THEN ! verification time ALLOCATE(correctedtime(gridinfov%arraysize)) correctedtime(:) = gridinfov%array(1:gridinfov%arraysize)%time DO i = 1, gridinfov%arraysize correctedtime(i) = correctedtime(i) + & timedelta_new(sec=gridinfov%array(i)%timerange%p1) ENDDO ENDIF IF (ltime_definition == 2) THEN ! set all to analysis ALLOCATE(correctedtimerange(gridinfov%arraysize)) correctedtimerange(:) = gridinfov%array(1:gridinfov%arraysize)%timerange correctedtimerange(:)%p1 = 0 ENDIF DO i = 1, ngrid IF (PRESENT(anavar)) THEN j = count((this(i)%griddim == gridinfov%array(1:gridinfov%arraysize)%griddim) &