Sources/misc/createnewfield.f

module createnewfield_m

  implicit none

  integer, parameter:: MaxWriteField = 100
  character(len=255), save:: FieldName(MaxWriteField)
  integer:: NbField = 0
  integer, save:: Ncid(MaxWriteField)

contains

  subroutine CreateNewField(name, my_shape)

    USE netcdf, ONLY: nf90_clobber, nf90_double, nf90_unlimited
    USE netcdf95, ONLY: nf95_create, nf95_def_dim, nf95_def_var, nf95_enddef

    character(len = *), intent(in):: name
    integer, intent(in):: my_shape(:)

    ! Local:
    integer Dimid(size(my_shape) + 1), n_dim, varid

    !--------------------------------------------------------------------

    NbField = NbField + 1
    FieldName(NbField) = ADJUSTL(name)
    n_dim = size(my_shape)

    call NF95_CREATE(TRIM(FieldName(NbField)) // '.nc', NF90_CLOBBER, &
         Ncid(NbField))
    call NF95_DEF_DIM(Ncid(NbField), 'X', my_shape(1), Dimid(1))
    if (n_dim >= 2) then
       call NF95_DEF_DIM(Ncid(NbField), 'Y', my_shape(2), Dimid(2))
       if (n_dim >= 3) then
          call NF95_DEF_DIM(Ncid(NbField), 'Z', my_shape(3), Dimid(3))
       end if
    end if
    call NF95_DEF_DIM(Ncid(NbField), 'iter', NF90_UNLIMITED, Dimid(n_dim + 1))
    call NF95_DEF_VAR(Ncid(NbField), FieldName(NbField), NF90_DOUBLE, Dimid, &
         Varid)
    call NF95_ENDDEF(Ncid(NbField))

  end subroutine CreateNewField

end module createnewfield_m
1	module createnewfield_m
2
3	implicit none
4
5	integer, parameter:: MaxWriteField = 100
6	character(len=255), save:: FieldName(MaxWriteField)
7	integer:: NbField = 0
8	integer, save:: Ncid(MaxWriteField)
9
10	contains
11
12	subroutine CreateNewField(name, my_shape)
13
14	USE netcdf, ONLY: nf90_clobber, nf90_double, nf90_unlimited
15	USE netcdf95, ONLY: nf95_create, nf95_def_dim, nf95_def_var, nf95_enddef
16
17	character(len = *), intent(in):: name
18	integer, intent(in):: my_shape(:)
19
20	! Local:
21	integer Dimid(size(my_shape) + 1), n_dim, varid
22
23	!--------------------------------------------------------------------
24
25	NbField = NbField + 1
26	FieldName(NbField) = ADJUSTL(name)
27	n_dim = size(my_shape)
28
29	call NF95_CREATE(TRIM(FieldName(NbField)) // '.nc', NF90_CLOBBER, &
30	Ncid(NbField))
31	call NF95_DEF_DIM(Ncid(NbField), 'X', my_shape(1), Dimid(1))
32	if (n_dim >= 2) then
33	call NF95_DEF_DIM(Ncid(NbField), 'Y', my_shape(2), Dimid(2))
34	if (n_dim >= 3) then
35	call NF95_DEF_DIM(Ncid(NbField), 'Z', my_shape(3), Dimid(3))
36	end if
37	end if
38	call NF95_DEF_DIM(Ncid(NbField), 'iter', NF90_UNLIMITED, Dimid(n_dim + 1))
39	call NF95_DEF_VAR(Ncid(NbField), FieldName(NbField), NF90_DOUBLE, Dimid, &
40	Varid)
41	call NF95_ENDDEF(Ncid(NbField))
42
43	end subroutine CreateNewField
44
45	end module createnewfield_m