export_to_gridinfo export_to_gridinfov export_to_gridinfovv volgrid6d_export_to_file volgrid6d_write_on_file

volgrid6d_write_on_file() subroutine volgrid6d_write_on_file ( type(volgrid6d), intent(in)  this, integer, intent(inout), optional  unit, character(len=*), intent(in), optional  description, character(len=*), intent(in), optional  filename, character(len=*), intent(out), optional  filename_auto  ) Scrittura su file di un volume Volgrid6d. Scrittura su file unformatted di un intero volume Volgrid6d. Il volume viene serializzato e scritto su file. Il file puo' essere aperto opzionalmente dall'utente. Si possono passare opzionalmente unità e nome del file altrimenti assegnati internamente a dei default; se assegnati internamente tali parametri saranno in output. Se non viene fornito il nome file viene utilizzato un file di default con nome pari al nome del programma in esecuzione con postfisso ".vg6d". Come parametro opzionale c'è la description che insieme alla data corrente viene inserita nell'header del file. Parametri [in] this volume volgrid6d da scrivere [in,out] unit unità su cui scrivere; se passata =0 ritorna il valore rielaborato (default =rielaborato internamente con getlun ) [in] filename nome del file su cui scrivere [out] filename_auto nome del file generato se "filename" è omesso [in] description descrizione del volume Definizione alla linea 822 del file volgrid6d_class.F90. !! descriptor of the grid in the \a volgrid object is assigned to the !! descriptor contained in \a gridinfo if it is missing in \a volgrid, !! otherwise it is checked and the object is rejected if grids do not !! match. SUBROUTINE import_from_gridinfo(this, gridinfo, force, dup_mode, clone, & isanavar) 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 CHARACTER(len=255) :: type INTEGER :: itime0, itimerange0, itime1, itimerange1, itime, itimerange, & ilevel, ivar, ldup_mode LOGICAL :: dup TYPE(datetime) :: correctedtime TYPE(vol7d_timerange) :: correctedtimerange REAL,ALLOCATABLE :: tmpgrid(:,:) IF (PRESENT(dup_mode)) THEN ldup_mode = dup_mode ELSE ldup_mode = 0 ENDIF call get_val(this%griddim,proj_type=type) #ifdef DEBUG call l4f_category_log(this%category,l4f_debug,"import_from_gridinfo: "//trim(type)) #endif if (.not. c_e(type))then call copy(gridinfo%griddim, this%griddim) ! ho gia` fatto init, altrimenti non potrei fare la get_val(this%griddim) ! per cui meglio non ripetere ! call init(this,gridinfo%griddim,categoryappend) CALL volgrid6d_alloc_vol(this, ini=.true.) ! decode? else if (.not. (this%griddim == gridinfo%griddim ))then CALL l4f_category_log(this%category,l4f_error, & "volgrid and gridinfo grid type or size are different, gridinfo rejected") CALL raise_error() RETURN end if ! Cerco gli indici del campo da inserire, se non trovo metto nel primo missing ilevel = index(this%level, gridinfo%level) IF (ilevel == 0 .AND. optio_log(force)) THEN ilevel = index(this%level, vol7d_level_miss) IF (ilevel /= 0) this%level(ilevel) = gridinfo%level ENDIF IF (ilevel == 0) THEN CALL l4f_category_log(this%category,l4f_error, & "volgrid6d: level not valid for volume, gridinfo rejected") CALL raise_error() RETURN ENDIF IF (optio_log(isanavar)) THEN ! assign to all times and timeranges itime0 = 1 itime1 = SIZE(this%time) itimerange0 = 1 itimerange1 = SIZE(this%timerange) ELSE ! usual case correctedtime = gridinfo%time IF (this%time_definition == 1 .OR. this%time_definition == 2) correctedtime = correctedtime + & timedelta_new(sec=gridinfo%timerange%p1) itime0 = index(this%time, correctedtime) IF (itime0 == 0 .AND. optio_log(force)) THEN itime0 = index(this%time, datetime_miss) IF (itime0 /= 0) this%time(itime0) = correctedtime ENDIF IF (itime0 == 0) THEN CALL l4f_category_log(this%category,l4f_error, & "volgrid6d: time not valid for volume, gridinfo rejected") CALL raise_error() RETURN ENDIF itime1 = itime0 correctedtimerange = gridinfo%timerange IF (this%time_definition == 2) correctedtimerange%p1 = 0 itimerange0 = index(this%timerange, correctedtimerange)