bibio/Writefield/writefield.f

module WriteField_m

  use CreateNewField_m, only: CreateNewField, ncid, nbfield, MaxWriteField
  use GetFieldIndex_m, only: GetFieldIndex
  USE netcdf95, ONLY: nf95_put_var

  implicit none

  integer, save:: Record(MaxWriteField)

  interface WriteField
     module procedure WriteField3d,WriteField2d,WriteField1d
  end interface WriteField

  private
  public WriteField

contains

  subroutine WriteField1d(name,Field)

    character(len=*), intent(in):: name
    real, intent(in):: Field(:)

    ! Local:
    integer index

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

    Index=GetFieldIndex(name)

    if (Index==-1) then
       call CreateNewField(name, shape(field))
       Index=nbfield
       Record(Index) = 1
    else
       Record(Index)=Record(Index)+1
    endif

    call NF95_PUT_VAR(Ncid(Index), Varid = 1, values = Field, &
         start = (/1, Record(Index)/))

  end subroutine WriteField1d

  !****************************************************************

  subroutine WriteField2d(name,Field)

    character(len=*), intent(in):: name
    real, intent(in):: Field(:, :)

    ! Local:
    integer index

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

    Index=GetFieldIndex(name)

    if (Index==-1) then
       call CreateNewField(name, shape(field))
       Index=nbfield
       Record(Index) = 1
    else
       Record(Index)=Record(Index)+1
    endif

    call NF95_PUT_VAR(Ncid(Index), Varid = 1, values = Field, &
         start = (/1, 1, Record(Index)/))

  end subroutine WriteField2d

  !****************************************************************

  subroutine WriteField3d(name,Field)

    character(len=*), intent(in):: name
    real, intent(in):: Field(:, :, :)

    ! Local:
    integer index

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

    Index=GetFieldIndex(name)

    if (Index==-1) then
       call CreateNewField(name, shape(field))
       Index=nbfield
       Record(Index) = 1
    else
       Record(Index)=Record(Index)+1
    endif

    call NF95_PUT_VAR(Ncid(Index), Varid = 1, values = Field, &
         start = (/1, 1, 1, Record(Index)/))

  end subroutine WriteField3d

end module WriteField_m
1	module WriteField_m
2
3	use CreateNewField_m, only: CreateNewField, ncid, nbfield, MaxWriteField
4	use GetFieldIndex_m, only: GetFieldIndex
5	USE netcdf95, ONLY: nf95_put_var
6
7	implicit none
8
9	integer, save:: Record(MaxWriteField)
10
11	interface WriteField
12	module procedure WriteField3d,WriteField2d,WriteField1d
13	end interface WriteField
14
15	private
16	public WriteField
17
18	contains
19
20	subroutine WriteField1d(name,Field)
21
22	character(len=*), intent(in):: name
23	real, intent(in):: Field(:)
24
25	! Local:
26	integer index
27
28	!-------------------------------------------
29
30	Index=GetFieldIndex(name)
31
32	if (Index==-1) then
33	call CreateNewField(name, shape(field))
34	Index=nbfield
35	Record(Index) = 1
36	else
37	Record(Index)=Record(Index)+1
38	endif
39
40	call NF95_PUT_VAR(Ncid(Index), Varid = 1, values = Field, &
41	start = (/1, Record(Index)/))
42
43	end subroutine WriteField1d
44
45	!****************************************************************
46
47	subroutine WriteField2d(name,Field)
48
49	character(len=*), intent(in):: name
50	real, intent(in):: Field(:, :)
51
52	! Local:
53	integer index
54
55	!-------------------------------------------
56
57	Index=GetFieldIndex(name)
58
59	if (Index==-1) then
60	call CreateNewField(name, shape(field))
61	Index=nbfield
62	Record(Index) = 1
63	else
64	Record(Index)=Record(Index)+1
65	endif
66
67	call NF95_PUT_VAR(Ncid(Index), Varid = 1, values = Field, &
68	start = (/1, 1, Record(Index)/))
69
70	end subroutine WriteField2d
71
72	!****************************************************************
73
74	subroutine WriteField3d(name,Field)
75
76	character(len=*), intent(in):: name
77	real, intent(in):: Field(:, :, :)
78
79	! Local:
80	integer index
81
82	!-------------------------------------------
83
84	Index=GetFieldIndex(name)
85
86	if (Index==-1) then
87	call CreateNewField(name, shape(field))
88	Index=nbfield
89	Record(Index) = 1
90	else
91	Record(Index)=Record(Index)+1
92	endif
93
94	call NF95_PUT_VAR(Ncid(Index), Varid = 1, values = Field, &
95	start = (/1, 1, 1, Record(Index)/))
96
97	end subroutine WriteField3d
98
99	end module WriteField_m