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 1026 del file volgrid6d_class.F90. #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 SUBROUTINE import_from_gridinfovv(this, gridinfov, dup_mode, clone, decode, & time_definition, anavar, categoryappend) TYPE(volgrid6d),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=*),INTENT(IN),OPTIONAL :: anavar(:) CHARACTER(len=*),INTENT(in),OPTIONAL :: categoryappend 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) & .AND. .NOT.isanavar(:)) IF (j <= 0) THEN CALL l4f_category_log(category, l4f_fatal, 'grid n.'//t2c(i)// & ' has only constant data, this is not allowed') CALL l4f_category_log(category, l4f_fatal, 'please check anavar argument') CALL raise_fatal_error() ENDIF ENDIF IF (ltime_definition == 1 .OR. ltime_definition == 2) THEN ! verification time ntime = count_distinct(correctedtime, &