display_volgrid6d display_volgrid6dv export_to_gridinfo export_to_gridinfov export_to_gridinfovv import_from_gridinfo import_from_gridinfovv v7d_v7d_transform v7d_volgrid6d_transform vg6d_c2a vg6d_reduce vg6d_rounding vg6d_wind_rot vg6d_wind_unrot vg6dv_rounding volgrid6d_alloc volgrid6d_alloc_vol volgrid6d_delete volgrid6d_export_to_file volgrid6d_import_from_file volgrid6d_init volgrid6d_read_from_file volgrid6d_transform volgrid6d_v7d_transform volgrid6d_v7d_transform_compute volgrid6d_write_on_file volgrid6dv_delete volgrid6dv_transform volgrid6dv_v7d_transform volgrid_get_vol_2d volgrid_get_vol_3d volgrid_set_vol_2d volgrid_set_vol_3d

volgrid6d_write_on_file() subroutine volgrid6d_class::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  ) private 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)