WriteDataToNETCDF Subroutine

public subroutine WriteDataToNETCDF(dim1, vecdim, nVal, ndims, DimNames, VarNames, Coord, Values, FileString, coord1, coord2, coord3, CoordNames, attr_values, attr_coords)

Uses

  • proc~~writedatatonetcdf~~UsesGraph proc~writedatatonetcdf WriteDataToNETCDF module~modgvec_globals MODgvec_Globals proc~writedatatonetcdf->module~modgvec_globals module~modgvec_io_netcdf MODgvec_IO_NETCDF proc~writedatatonetcdf->module~modgvec_io_netcdf iso_fortran_env iso_fortran_env module~modgvec_globals->iso_fortran_env module~modgvec_io_netcdf->module~modgvec_globals netcdf netcdf module~modgvec_io_netcdf->netcdf

Subroutine to write multidimensional data to netCDF format

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: dim1

dimension of the data (either 1D,2D or 3D)
integer, intent(in) :: vecdim

dimension of coordinates
integer, intent(in) :: nVal

Number of nodal output variables
integer, intent(in) :: ndims(1:dim1)

size of the data in each dimension
character(len=*), intent(in) :: DimNames(1:dim1)

Names of dimensions of multi-dimensional array
character(len=*), intent(in) :: VarNames(nVal)

Names of all variables that will be written out
real(kind=wp), intent(in) :: Coord(vecdim,1:PRODUCT(ndims))
real(kind=wp), intent(in) :: Values(nVal,1:PRODUCT(ndims))

Statevector
character(len=*), intent(in) :: FileString

Output file name (without .nc ending)
real(kind=wp), intent(in), optional :: coord1(:)

Netcdf coordinate values e.g. rho, theta and zeta
real(kind=wp), intent(in), optional :: coord2(:)

Netcdf coordinate values e.g. rho, theta and zeta
real(kind=wp), intent(in), optional :: coord3(:)

Netcdf coordinate values e.g. rho, theta and zeta
character(len=*), intent(in), optional :: CoordNames(1:dim1)

Names of the dimensions
character(len=255), intent(in), optional :: attr_values(nVal,2)
character(len=255), intent(in), optional :: attr_coords(dim1,2)

Calls

proc~~writedatatonetcdf~~CallsGraph proc~writedatatonetcdf WriteDataToNETCDF proc~ncfile_def_dim t_ncfile%ncfile_def_dim proc~writedatatonetcdf->proc~ncfile_def_dim proc~ncfile_end_def_mode t_ncfile%ncfile_end_def_mode proc~writedatatonetcdf->proc~ncfile_end_def_mode proc~ncfile_init ncfile_init proc~writedatatonetcdf->proc~ncfile_init proc~ncfile_put_array t_ncfile%ncfile_put_array proc~writedatatonetcdf->proc~ncfile_put_array proc~ncfile_put_attributes_char t_ncfile%ncfile_put_attributes_char proc~writedatatonetcdf->proc~ncfile_put_attributes_char proc~ncfile_put_char t_ncfile%ncfile_put_char proc~writedatatonetcdf->proc~ncfile_put_char nf90_def_dim nf90_def_dim proc~ncfile_def_dim->nf90_def_dim proc~mpi_check_single_access mpi_check_single_access proc~ncfile_def_dim->proc~mpi_check_single_access proc~ncfile_handle_error t_ncfile%ncfile_handle_error proc~ncfile_def_dim->proc~ncfile_handle_error nf90_enddef nf90_enddef proc~ncfile_end_def_mode->nf90_enddef proc~ncfile_end_def_mode->proc~mpi_check_single_access proc~ncfile_end_def_mode->proc~ncfile_handle_error proc~ncfile_init->proc~mpi_check_single_access proc~ncfile_openfile t_ncfile%ncfile_openfile proc~ncfile_init->proc~ncfile_openfile nf90_def_var nf90_def_var proc~ncfile_put_array->nf90_def_var nf90_inq_varid nf90_inq_varid proc~ncfile_put_array->nf90_inq_varid nf90_put_var nf90_put_var proc~ncfile_put_array->nf90_put_var proc~ncfile_put_array->proc~mpi_check_single_access proc~ncfile_put_array->proc~ncfile_handle_error proc~ncfile_put_array->proc~ncfile_openfile proc~ncfile_put_attributes_char->nf90_inq_varid nf90_put_att nf90_put_att proc~ncfile_put_attributes_char->nf90_put_att proc~ncfile_put_attributes_char->proc~mpi_check_single_access proc~ncfile_put_attributes_char->proc~ncfile_handle_error proc~ncfile_put_char->nf90_def_var proc~ncfile_put_char->nf90_inq_varid proc~ncfile_put_char->nf90_put_var proc~ncfile_put_char->proc~mpi_check_single_access proc~ncfile_put_char->proc~ncfile_handle_error proc~ncfile_handle_error->proc~mpi_check_single_access nf90_strerror nf90_strerror proc~ncfile_handle_error->nf90_strerror proc~ncfile_openfile->proc~mpi_check_single_access proc~ncfile_openfile->proc~ncfile_handle_error nf90_create nf90_create proc~ncfile_openfile->nf90_create nf90_open nf90_open proc~ncfile_openfile->nf90_open

Called by

proc~~writedatatonetcdf~~CalledByGraph proc~writedatatonetcdf WriteDataToNETCDF proc~visu_axisnb Visu_axisNB proc~visu_axisnb->proc~writedatatonetcdf proc~visufrenet VisuFrenet proc~visufrenet->proc~writedatatonetcdf proc~hmap_axisnb_init_params hmap_axisNB_init_params proc~hmap_axisnb_init_params->proc~visu_axisnb proc~hmap_frenet_init_params hmap_frenet_init_params proc~hmap_frenet_init_params->proc~visufrenet interface~t_hmap_axisnb t_hmap_axisNB interface~t_hmap_axisnb->proc~hmap_axisnb_init_params proc~hmap_axisnb_init hmap_axisNB_init interface~t_hmap_axisnb->proc~hmap_axisnb_init interface~t_hmap_frenet t_hmap_frenet interface~t_hmap_frenet->proc~hmap_frenet_init_params proc~hmap_frenet_init hmap_frenet_init interface~t_hmap_frenet->proc~hmap_frenet_init proc~hmap_axisnb_init->proc~hmap_axisnb_init_params proc~hmap_frenet_init->proc~hmap_frenet_init_params

Source Code

SUBROUTINE WriteDataToNETCDF(dim1,vecdim,nVal,ndims,Dimnames,VarNames,Coord,Values,FileString,coord1,coord2,coord3, CoordNames, attr_values, attr_coords)
! MODULES
USE MODgvec_Globals, ONLY:wp,abort,UNIT_stdOut
USE MODgvec_io_netcdf, ONLY:t_ncfile,ncfile_init
! IMPLICIT VARIABLE HANDLING
IMPLICIT NONE
!-----------------------------------------------------------------------------------------------------------------------------------
! INPUT VARIABLES
INTEGER,INTENT(IN)            :: dim1                    !! dimension of the data (either 1D,2D or 3D)
INTEGER,INTENT(IN)            :: vecdim                  !! dimension of coordinates
INTEGER,INTENT(IN)            :: nVal                    !! Number of nodal output variables
INTEGER,INTENT(IN)            :: ndims(1:dim1)           !! size of the data in each dimension
CHARACTER(LEN=*),INTENT(IN)   :: DimNames(1:dim1)        !! Names of dimensions of multi-dimensional array
REAL(wp),INTENT(IN)           :: Coord(vecdim,1:PRODUCT(ndims))      ! CoordinatesVector
CHARACTER(LEN=*),INTENT(IN)   :: VarNames(nVal)          !! Names of all variables that will be written out
REAL(wp),INTENT(IN)           :: Values(nVal,1:PRODUCT(ndims))   !! Statevector
CHARACTER(LEN=*),INTENT(IN)   :: FileString              !! Output file name (without .nc ending)
REAL(wp),INTENT(IN),OPTIONAL  :: coord1(:),coord2(:),coord3(:) !! Netcdf coordinate values e.g. rho, theta and zeta
CHARACTER(LEN=*),INTENT(IN),OPTIONAL :: CoordNames(1:dim1) !! Names of the dimensions
CHARACTER(LEN=255),INTENT(IN),OPTIONAL :: attr_values(nVal,2)
CHARACTER(LEN=255),INTENT(IN),OPTIONAL :: attr_coords(dim1,2)
!-----------------------------------------------------------------------------------------------------------------------------------
! OUTPUT VARIABLES
!-----------------------------------------------------------------------------------------------------------------------------------
! LOCAL VARIABLES
CLASS(t_ncfile),ALLOCATABLE  :: nc  !! container for netcdf-file
CHARACTER(LEN=255) :: tmpVarName,tmpVarNameY,tmpVarNameZ, coord_attr
CHARACTER(LEN=255) :: tmpCoordNames(1:dim1)
INTEGER :: def_put_mode,i,iVal,StrLen,dimids(1:dim1),vecdimid
LOGICAL            :: isVector,maybeVector
!===================================================================================================================================
#if NETCDF
WRITE(UNIT_stdOut,'(A)',ADVANCE='NO')'   WRITE DATA TO NETCDF FILE "'//TRIM(FileString)//'.nc" ...'
#else
WRITE(UNIT_stdOut,'(A)')'   OMIT WRITING DATA TO NETCDF FILE "'//TRIM(FileString)//'.nc" (not compiled with netcdf)'
RETURN
#endif
CALL ncfile_init(nc,TRIM(fileString)//".nc","o")
CALL nc%def_dim("xyz",vecdim,vecdimid)
DO i=1,dim1
  CALL nc%def_dim(TRIM(DimNames(i)),ndims(i),dimids(i))
END DO
IF (PRESENT(CoordNames)) THEN
  tmpCoordNames = CoordNames
ELSE
  tmpCoordNames = DimNames
END IF
DO def_put_mode=1,2
  IF(def_put_mode.EQ.2) CALL nc%end_def_mode()
  CALL nc%put_array("pos",1+dim1,(/vecdim,ndims/),(/vecdimid,dimids/),def_put_mode,real_in=Coord)
  IF(def_put_mode.EQ.1) THEN
    coord_attr = "xyz"
    DO i=1,dim1
      coord_attr = TRIM(coord_attr)//" "//TRIM(tmpCoordNames(i))
    END DO
    CALL nc%put_attr_char("pos",1,(/"coordinates"/),(/coord_attr/))
    CALL nc%put_attr_char("pos",2,(/"long_name","symbol   "/),(/"position vector", "\mathbf{x}     "/))
  END IF
  !accout for vectors:
  ! if Variable Name ends with an X and the following have the same name with Y and Z
  ! then it forms a vector variable (X is omitted for the name)

  iVal=0 !scalars
  IF(def_put_mode.EQ.1) THEN
    coord_attr = ""
    DO i=1,dim1
      coord_attr = TRIM(coord_attr)//" "//TRIM(tmpCoordNames(i))
    END DO
  END IF
  DO WHILE(iVal.LT.nVal)
    iVal=iVal+1
    tmpVarName=TRIM(VarNames(iVal))
    StrLen=LEN(TRIM(tmpVarName))
    maybeVector=(iVal+vecdim-1.LE.nVal)
    isVector=.FALSE.
    IF(maybeVector)THEN
      SELECT CASE(vecdim)
      CASE(2)
        tmpVarNameY=TRIM(VarNames(iVal+1))
        isVector=((iVal+2.LE.nVal).AND.(INDEX(tmpVarName( StrLen:StrLen),"X").NE.0) &
                                  .AND.(INDEX(tmpVarNameY(:StrLen),TRIM(tmpVarName(:StrLen-1))//"Y").NE.0))
      CASE(3)
        tmpVarNameY=TRIM(VarNames(iVal+1))
        tmpVarNameZ=TRIM(VarNames(iVal+2))
        isVector=((iVal+2.LE.nVal).AND.(INDEX(tmpVarName( StrLen:StrLen),"X").NE.0) &
                                  .AND.(INDEX(tmpVarNameY(:StrLen),TRIM(tmpVarName(:StrLen-1))//"Y").NE.0) &
                                  .AND.(INDEX(tmpVarNameZ(:StrLen),TRIM(tmpVarName(:StrLen-1))//"Z").NE.0))

      END SELECT
    END IF !maybevector

    IF(isvector)THEN !variable is a vector!
      tmpVarName=tmpVarName(:StrLen-1)
      CALL nc%put_array(TRIM(tmpVarName),1+dim1,(/vecdim,ndims/),(/vecdimid,dimids/),def_put_mode, &
                        real_in=Values(iVal:ival-1+vecdim,:))
      iVal=iVal+vecdim-1 !skip the Y (& Z) components
    ELSE
      CALL nc%put_array(TRIM(tmpVarName),dim1,ndims,dimids,def_put_mode,real_in=Values(iVal,:))
    END IF !isvector
    IF (def_put_mode.EQ.1) THEN
      CALL nc%put_attr_char(TRIM(tmpVarName),1,(/"coordinates"/),(/coord_attr/))
      IF (PRESENT(attr_values)) THEN
        CALL nc%put_attr_char(TRIM(tmpVarName),2,(/"long_name","symbol   "/),attr_values(iVal,:))
      END IF
    END IF
  END DO !iVal <=nVal
  IF (PRESENT(coord1)) THEN
    CALL nc%put_array(TRIM(tmpCoordNames(1)),1,(/ndims(1)/),(/dimids(1)/),def_put_mode,real_in=coord1)
  END IF
  IF (PRESENT(coord2)) THEN
    CALL nc%put_array(TRIM(tmpCoordNames(2)),1,(/ndims(2)/),(/dimids(2)/),def_put_mode,real_in=coord2)
  END IF
  IF (PRESENT(coord3)) THEN
    CALL nc%put_array(TRIM(tmpCoordNames(3)),1,(/ndims(3)/),(/dimids(3)/),def_put_mode,real_in=coord3)
  END IF
  IF (def_put_mode.EQ.1) THEN
    IF (PRESENT(attr_coords)) THEN
      DO i=1,dim1
        CALL nc%put_attr_char(TRIM(tmpCoordNames(i)),2,(/"long_name","symbol   "/),attr_coords(i,:))
      END DO
    END IF
  END IF
  CALL nc%put_char("xyz",vecdimid,def_put_mode,char_in="xyz")
  IF (def_put_mode.EQ.1) THEN
    CALL nc%put_attr_char("xyz",2,(/"long_name","symbol   "/),(/"cartesian vector components","\mathbf{x}                 "/))
  END IF
END DO !mode
CALL nc%free()
WRITE(UNIT_stdOut,'(A)',ADVANCE='YES')"   DONE"
END SUBROUTINE WriteDataToNETCDF