/[lmdze]/trunk/bibio/Writefield/createnewfield.f
ViewVC logotype

Contents of /trunk/bibio/Writefield/createnewfield.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 109 - (show annotations)
Wed Sep 17 10:08:00 2014 UTC (9 years, 7 months ago) by guez
File size: 1297 byte(s) 
Moved a call to writefield from guide to tau2alpha. (dxdys does not
change with itau.) So dxdys does not need to be a module variable any
longer. Other variables of modules tau2alpha_m downgraded to local
variables of tau2alpha, since they were not used elsewhere.

Procedures write_field[13]d and formcoord were never called. Could
then remove int2str.

Inline writefield_gen into writefield.

CreateNewField takes an integer array argument instead of 3 scalar
integers. CreateNewField now creates a number of dimensions adapted to
the rank of the output field, instead of always 4 dimensions.

Changed names of variables of module write_field: fieldid to
ncid, fieldindex to record, fieldvarid to varid.

In writefield_gen, if index == -1, no need to call GetFieldIndex
again, we know that the result is nbfield.

In guide, moved calls to writefield for some variables inside if
first_call: those variables do not change with time. Removed ztau:
computed only to be output, does not seem meaningful. Removed
writefield for aire: does not change with time and is already in
"grilles_gcm.nc".
1 module createnewfield_m
2
3 implicit none
4
5 contains
6
7 subroutine CreateNewField(name, my_shape)
8
9 USE netcdf, ONLY: nf90_clobber, nf90_double, nf90_unlimited
10 USE netcdf95, ONLY: nf95_create, nf95_def_dim, nf95_def_var, nf95_enddef
11 USE write_field, ONLY: ncid, record, fieldname, varid, nbfield
12
13 character(len = *), intent(in):: name
14 integer, intent(in):: my_shape(:)
15
16 ! Local:
17 integer Dimid(size(my_shape) + 1), n_dim
18
19 !--------------------------------------------------------------------
20
21 NbField = NbField + 1
22 FieldName(NbField) = ADJUSTL(name)
23 Record(NbField) = 1
24 n_dim = size(my_shape)
25
26 call NF95_CREATE(TRIM(FieldName(NbField)) // '.nc', NF90_CLOBBER, &
27 Ncid(NbField))
28 call NF95_DEF_DIM(Ncid(NbField), 'X', my_shape(1), Dimid(1))
29 if (n_dim >= 2) then
30 call NF95_DEF_DIM(Ncid(NbField), 'Y', my_shape(2), Dimid(2))
31 if (n_dim >= 3) then
32 call NF95_DEF_DIM(Ncid(NbField), 'Z', my_shape(3), Dimid(3))
33 end if
34 end if
35 call NF95_DEF_DIM(Ncid(NbField), 'iter', NF90_UNLIMITED, Dimid(n_dim + 1))
36 call NF95_DEF_VAR(Ncid(NbField), FieldName(NbField), NF90_DOUBLE, &
37 Dimid, Varid(NbField))
38 call NF95_ENDDEF(Ncid(NbField))
39
40 end subroutine CreateNewField
41
42 end module createnewfield_m
  ViewVC Help
Powered by ViewVC 1.1.21