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_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 811 del file volgrid6d_class.F90. !! volgrid6d_alloc(). If the \a force argument is missing or \a !! .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 \a force is !! \a .TRUE. , the gridinfo dimension descriptors that do not fit into !! available descriptors in the \a 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 \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 !< object in which to import TYPE(gridinfo_def),INTENT(in) :: gridinfo !< gridinfo object to be imported LOGICAL,INTENT(in),OPTIONAL :: force !< if provided and \c .TRUE., the gridinfo is forced into an empty element of \a this, if required and possible 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 :: isanavar !< if provides and \a .TRUE., the gridinfo object is treated as time-independent and replicated for every time and timerange 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