display_volgrid6d display_volgrid6dv

display_volgrid6d() subroutine display_volgrid6d ( type(volgrid6d), intent(in) this ) Display object on screen. Show brief content on screen. Parametri [in] this object to display Definizione alla linea 3497 del file volgrid6d_class.F90. type(volgrid6d),intent(out) :: vg6dout !< output volume type(vol7d_level),intent(in) :: roundlevel(:) !< new level list to use for rounding type(vol7d_timerange),intent(in) :: roundtimerange(:) !< new timerange list to use for rounding logical,intent(in),optional :: merge !< if there are data on equal levels or timeranges will be merged POINT BY POINT with priority for the first data found ordered by icreasing var index (require to decode all the data) integer :: nlevel,ntime,ntimerange,nvar,ilevel,itimerange,ivar,indl,indt,itime,nx,ny real,allocatable :: vol2d(:,:) nx=vg6din%griddim%dim%nx ny=vg6din%griddim%dim%ny nlevel=count_distinct(roundlevel,back=.true.) ntime=size(vg6din%time) ntimerange=count_distinct(roundtimerange,back=.true.) nvar=size(vg6din%var) call init(vg6dout, vg6din%griddim, vg6din%time_definition, categoryappend="generated by vg6d_reduce") call volgrid6d_alloc(vg6dout, vg6din%griddim%dim, ntime, nlevel, ntimerange, nvar) if ( ASSOCIATED(vg6din%voldati) .or. optio_log(merge)) then call volgrid6d_alloc_vol(vg6dout,inivol=.true.,decode=.true.) allocate(vol2d(nx,ny)) else call volgrid6d_alloc_vol(vg6dout,inivol=.true.,decode=.false.) end if vg6dout%time=vg6din%time vg6dout%var=vg6din%var vg6dout%timerange=pack_distinct(roundtimerange,ntimerange,back=.true.) vg6dout%level=pack_distinct(roundlevel,nlevel,back=.true.) ! sort modified dimensions CALL sort(vg6dout%timerange) CALL sort(vg6dout%level) do ilevel=1,size(vg6din%level) indl=index(vg6dout%level,roundlevel(ilevel)) do itimerange=1,size(vg6din%timerange) indt=index(vg6dout%timerange,roundtimerange(itimerange)) do ivar=1, nvar do itime=1,ntime if ( ASSOCIATED(vg6din%voldati)) then vol2d=vg6din%voldati(:,:,ilevel,itime,itimerange,ivar) end if if (optio_log(merge)) then if ( .not. ASSOCIATED(vg6din%voldati)) then CALL grid_id_decode_data(vg6din%gaid(ilevel,itime,itimerange,ivar), vol2d) end if !! merge present data point by point where (.not. c_e(vg6dout%voldati(:,:,indl,itime,indt,ivar)))