phylmd/Mobidic/regr_lat_time_coefoz.f90

module regr_lat_time_coefoz_m

  ! This module is clean: no C preprocessor directive, no include line.
  ! Author: Lionel GUEZ

  implicit none

  private
  public regr_lat_time_coefoz

contains

  subroutine regr_lat_time_coefoz

    ! "regr_lat_time_coefoz" stands for "regrid latitude time
    ! coefficients ozone".

    ! This procedure reads from a NetCDF file parameters for ozone
    ! chemistry, regrids them in latitude and time, and writes the
    ! regridded fields to a new NetCDF file.

    ! The input fields depend on time, pressure level and
    ! latitude.
    ! We assume that the input fields are step functions
    ! of latitude.
    ! Regridding in latitude is made by averaging, with a cosine of
    ! latitude factor.
    ! The target LMDZ latitude grid is the "scalar" grid: "rlatu".
    ! The values of "rlatu" are taken to be the centers of intervals.
    ! Regridding in time is by linear interpolation.
    ! Monthly values are processed to get daily values, on the basis
    ! of a 360-day calendar.

    ! We assume that in the input file:
    ! -- the latitude is in degrees and strictly monotonic (as all
    ! NetCDF coordinate variables should be);
    ! -- time increases from January to December (even though we do
    ! not use values of the input time coordinate).

    use dimens_m, only: jjm
    use comgeom, only: rlatv
    use comconst, only: pi
    use regr1_step_av_m, only: regr1_step_av
    use regr3_lint_m, only: regr3_lint
    use netcdf95, only: nf95_open, nf95_gw_var, nf95_close, &
         nf95_inq_varid, handle_err, nf95_put_var
    use netcdf, only: nf90_nowrite, nf90_get_var

    ! Variables local to the procedure:

    integer ncid_in, ncid_out ! NetCDF IDs for input and output files
    integer n_plev ! number of pressure levels in the input data
    integer n_lat! number of latitudes in the input data

    real, pointer:: latitude(:)
    ! (of input data, converted to rad, sorted in strictly increasing order)

    real, allocatable:: lat_in_edg(:)
    ! (edges of latitude intervals for input data, in rad, in strictly
    ! increasing order)

    real, pointer:: plev(:) ! pressure level of input data
    logical decr_lat ! decreasing latitude in the input file

    real, allocatable:: o3_par_in(:, :, :) ! (n_lat, n_plev, 12)
    ! (ozone parameter from the input file)
    ! ("o3_par_in(j, l, month)" is at latitude "latitude(j)" and pressure
    ! level "plev(l)". "month" is between 1 and 12.)

    real, allocatable:: v_regr_lat(:, :, :) ! (jjm + 1, n_plev, 0:13)
    ! (mean of a variable "v" over a latitude interval)
    ! (First dimension is latitude interval.
    ! The latitude interval for "v_regr_lat(j,:, :)" contains "rlatu(j)".
    ! If "j" is between 2 and "jjm" then the interval is:
    ! [rlatv(j), rlatv(j-1)]
    ! If "j" is 1 or "jjm + 1" then the interval is:
    ! [rlatv(1), pi / 2]
    ! or:
    ! [- pi / 2, rlatv(jjm)]
    ! respectively.
    ! "v_regr_lat(:, l, :)" is for pressure level "plev(l)".
    ! Last dimension is month number.)

    real, allocatable:: o3_par_out(:, :, :) ! (jjm + 1, n_plev, 360)
    ! (regridded ozone parameter)
    ! ("o3_par_out(j, l, day)" is at latitude "rlatu(j)", pressure
    ! level "plev(l)" and date "January 1st 0h" + "tmidday(day)", in a
    ! 360-day calendar.)

    integer j
    integer i_v ! index of ozone parameter
    integer, parameter:: n_o3_param = 8 ! number of ozone parameters

    character(len=11) name_in(n_o3_param)
    ! (name of NetCDF primary variable in the input file)

    character(len=9) name_out(n_o3_param)
    ! (name of NetCDF primary variable in the output file)

    integer varid_in(n_o3_param), varid_out(n_o3_param), varid_plev, varid_time
    integer ncerr, varid
    ! (for NetCDF)

    real, parameter:: tmidmonth(0:13) = (/(-15. + 30. * j, j = 0, 13)/)
    ! (time to middle of month, in days since January 1st 0h, in a
    ! 360-day calendar)
    ! (We add values -15 and 375 so that, for example, day 3 of the year is
    ! interpolated between the December and the January value.)

    real, parameter:: tmidday(360) = (/(j + 0.5, j = 0, 359)/)
    ! (time to middle of day, in days since January 1st 0h, in a
    ! 360-day calendar)

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

    print *, "Call sequence information: regr_lat_time_coefoz"

    ! Names of ozone parameters:
    i_v = 0

    i_v = i_v + 1
    name_in(i_v) = "P_net"
    name_out(i_v) = "P_net_Mob"

    i_v = i_v + 1
    name_in(i_v) = "a2"
    name_out(i_v) = "a2"

    i_v = i_v + 1
    name_in(i_v) = "tro3"
    name_out(i_v) = "r_Mob"

    i_v = i_v + 1
    name_in(i_v) = "a4"
    name_out(i_v) = "a4"

    i_v = i_v + 1
    name_in(i_v) = "temperature"
    name_out(i_v) = "temp_Mob"

    i_v = i_v + 1
    name_in(i_v) = "a6"
    name_out(i_v) = "a6"

    i_v = i_v + 1
    name_in(i_v) = "Sigma"
    name_out(i_v) = "Sigma_Mob"

    i_v = i_v + 1
    name_in(i_v) = "R_Het"
    name_out(i_v) = "R_Het"

    call nf95_open("coefoz.nc", nf90_nowrite, ncid_in)

    ! Get coordinates from the input file:

    call nf95_inq_varid(ncid_in, "latitude", varid)
    call nf95_gw_var(ncid_in, varid, latitude)
    ! Convert from degrees to rad, because "rlatv" is in rad:
    latitude = latitude / 180. * pi
    n_lat = size(latitude)
    ! We need to supply the latitudes to "regr1_step_av" in
    ! increasing order, so invert order if necessary:
    decr_lat = latitude(1) > latitude(n_lat)
    if (decr_lat) latitude = latitude(n_lat:1:-1)

    ! Compute edges of latitude intervals:
    allocate(lat_in_edg(n_lat + 1))
    lat_in_edg(1) = - pi / 2
    forall (j = 2:n_lat) lat_in_edg(j) = (latitude(j - 1) + latitude(j)) / 2
    lat_in_edg(n_lat + 1) = pi / 2
    deallocate(latitude) ! pointer

    call nf95_inq_varid(ncid_in, "plev", varid)
    call nf95_gw_var(ncid_in, varid, plev)
    n_plev = size(plev)
    ! (We only need the pressure coordinate to copy it to the output file.)

    ! Get the IDs of ozone parameters in the input file:
    do i_v = 1, n_o3_param
       call nf95_inq_varid(ncid_in, trim(name_in(i_v)), varid_in(i_v))
    end do

    ! Create the output file and get the variable IDs:
    call prepare_out(ncid_in, varid_in, n_plev, name_out, ncid_out, &
         varid_out, varid_plev, varid_time)

    ! Write remaining coordinate variables:
    call nf95_put_var(ncid_out, varid_time, tmidday)
    call nf95_put_var(ncid_out, varid_plev, plev)

    deallocate(plev) ! pointer

    allocate(o3_par_in(n_lat, n_plev, 12))
    allocate(v_regr_lat(jjm + 1, n_plev, 0:13))
    allocate(o3_par_out(jjm + 1, n_plev, 360))

    do i_v = 1, n_o3_param
       ! Process ozone parameter "name_in(i_v)"

       ncerr = nf90_get_var(ncid_in, varid_in(i_v), o3_par_in)
       call handle_err("nf90_get_var", ncerr, ncid_in)

       if (decr_lat) o3_par_in = o3_par_in(n_lat:1:-1, :, :)

       ! Regrid in latitude:
       ! We average with respect to sine of latitude, which is
       ! equivalent to weighting by cosine of latitude:
       v_regr_lat(jjm+1:1:-1, :, 1:12) = regr1_step_av(o3_par_in, &
            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
       ! (invert order of indices in "v_regr_lat" because "rlatu" is
       ! decreasing)

       ! Duplicate January and December values, in preparation of time
       ! interpolation:
       v_regr_lat(:, :, 0) = v_regr_lat(:, :, 12)
       v_regr_lat(:, :, 13) = v_regr_lat(:, :, 1)

       ! Regrid in time by linear interpolation:
       o3_par_out = regr3_lint(v_regr_lat, tmidmonth, tmidday)

       ! Write to file:
       call nf95_put_var(ncid_out, varid_out(i_v), &
            o3_par_out(jjm+1:1:-1, :, :))
       ! (The order of "rlatu" is inverted in the output file)
    end do

    call nf95_close(ncid_out)
    call nf95_close(ncid_in)

  end subroutine regr_lat_time_coefoz

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

  subroutine prepare_out(ncid_in, varid_in, n_plev, name_out, ncid_out, &
       varid_out, varid_plev, varid_time)

    ! This subroutine creates the NetCDF output file, defines
    ! dimensions and variables, and writes one of the coordinate variables.

    use dimens_m, only: jjm
    use comgeom, only: rlatu
    use comconst, only: pi
    use numer_rec, only: assert_eq

    use netcdf95, only: nf95_create, nf95_def_dim, nf95_def_var, &
         nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
    use netcdf, only: nf90_clobber, nf90_float, nf90_copy_att, nf90_global

    integer, intent(in):: ncid_in, varid_in(:), n_plev
    character(len=*), intent(in):: name_out(:) ! of NetCDF variables
    integer, intent(out):: ncid_out, varid_out(:), varid_plev, varid_time

    ! Variables local to the procedure:
    integer ncerr
    integer dimid_rlatu, dimid_plev, dimid_time
    integer varid_rlatu
    integer i, n_o3_param

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

    print *, "Call sequence information: prepare_out"
    n_o3_param = assert_eq(size(varid_in), size(varid_out), &
         size(name_out), "prepare_out") 

    call nf95_create("coefoz_LMDZ.nc", nf90_clobber, ncid_out)

    ! Dimensions:
    call nf95_def_dim(ncid_out, "time", 360, dimid_time)
    call nf95_def_dim(ncid_out, "plev", n_plev, dimid_plev)
    call nf95_def_dim(ncid_out, "rlatu", jjm + 1, dimid_rlatu)

    ! Define coordinate variables:

    call nf95_def_var(ncid_out, "time", nf90_float, dimid_time, varid_time)
    call nf95_put_att(ncid_out, varid_time, "units", "days since 2000-1-1")
    call nf95_put_att(ncid_out, varid_time, "calendar", "360_day")
    call nf95_put_att(ncid_out, varid_time, "standard_name", "time")

    call nf95_def_var(ncid_out, "plev", nf90_float, dimid_plev, varid_plev)
    call nf95_put_att(ncid_out, varid_plev, "units", "millibar")
    call nf95_put_att(ncid_out, varid_plev, "standard_name", "air_pressure")
    call nf95_put_att(ncid_out, varid_plev, "long_name", "air pressure")

    call nf95_def_var(ncid_out, "rlatu", nf90_float, dimid_rlatu, varid_rlatu)
    call nf95_put_att(ncid_out, varid_rlatu, "units", "degrees_north")
    call nf95_put_att(ncid_out, varid_rlatu, "standard_name", "latitude")

    ! Define primary variables:

    do i = 1, n_o3_param
       call nf95_def_var(ncid_out, name_out(i), nf90_float, &
            (/dimid_rlatu, dimid_plev, dimid_time/), varid_out(i))

       ! The following commands may fail. That is OK. It should just
       ! mean that the attribute is not defined in the input file.

       ncerr = nf90_copy_att(ncid_in, varid_in(i), "long_name",&
            & ncid_out, varid_out(i))
       call handle_err_copy_att("long_name")

       ncerr = nf90_copy_att(ncid_in, varid_in(i), "units", ncid_out,&
            & varid_out(i))
       call handle_err_copy_att("units")

       ncerr = nf90_copy_att(ncid_in, varid_in(i), "standard_name", ncid_out,&
            & varid_out(i))
       call handle_err_copy_att("standard_name")
    end do

    ! Global attributes:
    call nf95_copy_att(ncid_in, nf90_global, "Conventions", ncid_out, &
         nf90_global)
    call nf95_copy_att(ncid_in, nf90_global, "title", ncid_out, nf90_global)
    call nf95_copy_att(ncid_in, nf90_global, "source", ncid_out, nf90_global)
    call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")

    call nf95_enddef(ncid_out)

    ! Write one of the coordinate variables:
    call nf95_put_var(ncid_out, varid_rlatu, rlatu(jjm+1:1:-1) / pi * 180.)
    ! (convert from rad to degrees and sort in increasing order)

  contains

    subroutine handle_err_copy_att(att_name)

      use netcdf, only: nf90_noerr, nf90_strerror

      character(len=*), intent(in):: att_name

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

      if (ncerr /= nf90_noerr) then
         print *, "prepare_out " // trim(name_out(i)) &
              // " nf90_copy_att " // att_name // " -- " &
              // trim(nf90_strerror(ncerr))
      end if

    end subroutine handle_err_copy_att

  end subroutine prepare_out

end module regr_lat_time_coefoz_m
1	guez	7	module regr_lat_time_coefoz_m
2	guez	3
3			! This module is clean: no C preprocessor directive, no include line.
4	guez	4	! Author: Lionel GUEZ
5	guez	3
6			implicit none
7
8			private
9	guez	9	public regr_lat_time_coefoz
10	guez	3
11			contains
12
13	guez	7	subroutine regr_lat_time_coefoz
14	guez	3
15	guez	7	! "regr_lat_time_coefoz" stands for "regrid latitude time
16			! coefficients ozone".
17	guez	3
18	guez	7	! This procedure reads from a NetCDF file parameters for ozone
19			! chemistry, regrids them in latitude and time, and writes the
20			! regridded fields to a new NetCDF file.
21
22	guez	3	! The input fields depend on time, pressure level and
23			! latitude.
24			! We assume that the input fields are step functions
25			! of latitude.
26	guez	7	! Regridding in latitude is made by averaging, with a cosine of
27			! latitude factor.
28			! The target LMDZ latitude grid is the "scalar" grid: "rlatu".
29	guez	3	! The values of "rlatu" are taken to be the centers of intervals.
30	guez	7	! Regridding in time is by linear interpolation.
31			! Monthly values are processed to get daily values, on the basis
32			! of a 360-day calendar.
33	guez	3
34	guez	7	! We assume that in the input file:
35			! -- the latitude is in degrees and strictly monotonic (as all
36			! NetCDF coordinate variables should be);
37			! -- time increases from January to December (even though we do
38			! not use values of the input time coordinate).
39	guez	3
40	guez	7	use dimens_m, only: jjm
41	guez	3	use comgeom, only: rlatv
42			use comconst, only: pi
43			use regr1_step_av_m, only: regr1_step_av
44	guez	7	use regr3_lint_m, only: regr3_lint
45	guez	22	use netcdf95, only: nf95_open, nf95_gw_var, nf95_close, &
46	guez	7	nf95_inq_varid, handle_err, nf95_put_var
47			use netcdf, only: nf90_nowrite, nf90_get_var
48	guez	3
49			! Variables local to the procedure:
50
51			integer ncid_in, ncid_out ! NetCDF IDs for input and output files
52			integer n_plev ! number of pressure levels in the input data
53			integer n_lat! number of latitudes in the input data
54
55			real, pointer:: latitude(:)
56			! (of input data, converted to rad, sorted in strictly increasing order)
57
58			real, allocatable:: lat_in_edg(:)
59			! (edges of latitude intervals for input data, in rad, in strictly
60			! increasing order)
61
62	guez	7	real, pointer:: plev(:) ! pressure level of input data
63	guez	3	logical decr_lat ! decreasing latitude in the input file
64
65	guez	7	real, allocatable:: o3_par_in(:, :, :) ! (n_lat, n_plev, 12)
66	guez	3	! (ozone parameter from the input file)
67			! ("o3_par_in(j, l, month)" is at latitude "latitude(j)" and pressure
68			! level "plev(l)". "month" is between 1 and 12.)
69
70	guez	7	real, allocatable:: v_regr_lat(:, :, :) ! (jjm + 1, n_plev, 0:13)
71			! (mean of a variable "v" over a latitude interval)
72			! (First dimension is latitude interval.
73			! The latitude interval for "v_regr_lat(j,:, :)" contains "rlatu(j)".
74			! If "j" is between 2 and "jjm" then the interval is:
75			! [rlatv(j), rlatv(j-1)]
76			! If "j" is 1 or "jjm + 1" then the interval is:
77			! [rlatv(1), pi / 2]
78			! or:
79			! [- pi / 2, rlatv(jjm)]
80			! respectively.
81			! "v_regr_lat(:, l, :)" is for pressure level "plev(l)".
82			! Last dimension is month number.)
83	guez	3
84	guez	7	real, allocatable:: o3_par_out(:, :, :) ! (jjm + 1, n_plev, 360)
85			! (regridded ozone parameter)
86			! ("o3_par_out(j, l, day)" is at latitude "rlatu(j)", pressure
87			! level "plev(l)" and date "January 1st 0h" + "tmidday(day)", in a
88			! 360-day calendar.)
89
90			integer j
91	guez	3	integer i_v ! index of ozone parameter
92			integer, parameter:: n_o3_param = 8 ! number of ozone parameters
93
94			character(len=11) name_in(n_o3_param)
95	guez	7	! (name of NetCDF primary variable in the input file)
96	guez	4
97	guez	3	character(len=9) name_out(n_o3_param)
98	guez	7	! (name of NetCDF primary variable in the output file)
99	guez	3
100	guez	7	integer varid_in(n_o3_param), varid_out(n_o3_param), varid_plev, varid_time
101	guez	22	integer ncerr, varid
102	guez	7	! (for NetCDF)
103	guez	3
104	guez	7	real, parameter:: tmidmonth(0:13) = (/(-15. + 30. * j, j = 0, 13)/)
105			! (time to middle of month, in days since January 1st 0h, in a
106			! 360-day calendar)
107			! (We add values -15 and 375 so that, for example, day 3 of the year is
108			! interpolated between the December and the January value.)
109	guez	3
110	guez	7	real, parameter:: tmidday(360) = (/(j + 0.5, j = 0, 359)/)
111			! (time to middle of day, in days since January 1st 0h, in a
112			! 360-day calendar)
113	guez	3
114			!---------------------------------
115
116	guez	10	print *, "Call sequence information: regr_lat_time_coefoz"
117	guez	3
118	guez	7	! Names of ozone parameters:
119	guez	3	i_v = 0
120
121			i_v = i_v + 1
122			name_in(i_v) = "P_net"
123			name_out(i_v) = "P_net_Mob"
124
125			i_v = i_v + 1
126			name_in(i_v) = "a2"
127			name_out(i_v) = "a2"
128
129			i_v = i_v + 1
130	guez	23	name_in(i_v) = "tro3"
131	guez	3	name_out(i_v) = "r_Mob"
132
133			i_v = i_v + 1
134			name_in(i_v) = "a4"
135			name_out(i_v) = "a4"
136
137			i_v = i_v + 1
138			name_in(i_v) = "temperature"
139			name_out(i_v) = "temp_Mob"
140
141			i_v = i_v + 1
142			name_in(i_v) = "a6"
143			name_out(i_v) = "a6"
144
145			i_v = i_v + 1
146			name_in(i_v) = "Sigma"
147			name_out(i_v) = "Sigma_Mob"
148
149			i_v = i_v + 1
150			name_in(i_v) = "R_Het"
151			name_out(i_v) = "R_Het"
152
153	guez	22	call nf95_open("coefoz.nc", nf90_nowrite, ncid_in)
154	guez	3
155			! Get coordinates from the input file:
156
157	guez	22	call nf95_inq_varid(ncid_in, "latitude", varid)
158			call nf95_gw_var(ncid_in, varid, latitude)
159	guez	3	! Convert from degrees to rad, because "rlatv" is in rad:
160	guez	4	latitude = latitude / 180. * pi
161	guez	3	n_lat = size(latitude)
162	guez	7	! We need to supply the latitudes to "regr1_step_av" in
163	guez	3	! increasing order, so invert order if necessary:
164			decr_lat = latitude(1) > latitude(n_lat)
165	guez	4	if (decr_lat) latitude = latitude(n_lat:1:-1)
166	guez	3
167			! Compute edges of latitude intervals:
168			allocate(lat_in_edg(n_lat + 1))
169			lat_in_edg(1) = - pi / 2
170			forall (j = 2:n_lat) lat_in_edg(j) = (latitude(j - 1) + latitude(j)) / 2
171			lat_in_edg(n_lat + 1) = pi / 2
172			deallocate(latitude) ! pointer
173
174	guez	22	call nf95_inq_varid(ncid_in, "plev", varid)
175			call nf95_gw_var(ncid_in, varid, plev)
176	guez	3	n_plev = size(plev)
177	guez	7	! (We only need the pressure coordinate to copy it to the output file.)
178	guez	3
179			! Get the IDs of ozone parameters in the input file:
180			do i_v = 1, n_o3_param
181			call nf95_inq_varid(ncid_in, trim(name_in(i_v)), varid_in(i_v))
182			end do
183
184	guez	7	! Create the output file and get the variable IDs:
185			call prepare_out(ncid_in, varid_in, n_plev, name_out, ncid_out, &
186			varid_out, varid_plev, varid_time)
187	guez	3
188	guez	7	! Write remaining coordinate variables:
189			call nf95_put_var(ncid_out, varid_time, tmidday)
190			call nf95_put_var(ncid_out, varid_plev, plev)
191
192			deallocate(plev) ! pointer
193
194			allocate(o3_par_in(n_lat, n_plev, 12))
195			allocate(v_regr_lat(jjm + 1, n_plev, 0:13))
196			allocate(o3_par_out(jjm + 1, n_plev, 360))
197
198	guez	3	do i_v = 1, n_o3_param
199			! Process ozone parameter "name_in(i_v)"
200
201			ncerr = nf90_get_var(ncid_in, varid_in(i_v), o3_par_in)
202			call handle_err("nf90_get_var", ncerr, ncid_in)
203
204			if (decr_lat) o3_par_in = o3_par_in(n_lat:1:-1, :, :)
205
206			! Regrid in latitude:
207			! We average with respect to sine of latitude, which is
208			! equivalent to weighting by cosine of latitude:
209	guez	7	v_regr_lat(jjm+1:1:-1, :, 1:12) = regr1_step_av(o3_par_in, &
210	guez	3	xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
211	guez	4	! (invert order of indices in "v_regr_lat" because "rlatu" is
212			! decreasing)
213	guez	3
214	guez	7	! Duplicate January and December values, in preparation of time
215			! interpolation:
216			v_regr_lat(:, :, 0) = v_regr_lat(:, :, 12)
217			v_regr_lat(:, :, 13) = v_regr_lat(:, :, 1)
218	guez	3
219	guez	7	! Regrid in time by linear interpolation:
220			o3_par_out = regr3_lint(v_regr_lat, tmidmonth, tmidday)
221	guez	3
222			! Write to file:
223	guez	7	call nf95_put_var(ncid_out, varid_out(i_v), &
224			o3_par_out(jjm+1:1:-1, :, :))
225	guez	3	! (The order of "rlatu" is inverted in the output file)
226			end do
227
228			call nf95_close(ncid_out)
229			call nf95_close(ncid_in)
230
231	guez	7	end subroutine regr_lat_time_coefoz
232	guez	3
233			!********************************************
234
235	guez	7	subroutine prepare_out(ncid_in, varid_in, n_plev, name_out, ncid_out, &
236			varid_out, varid_plev, varid_time)
237	guez	3
238			! This subroutine creates the NetCDF output file, defines
239	guez	7	! dimensions and variables, and writes one of the coordinate variables.
240	guez	3
241	guez	7	use dimens_m, only: jjm
242			use comgeom, only: rlatu
243	guez	3	use comconst, only: pi
244	guez	13	use numer_rec, only: assert_eq
245	guez	3
246	guez	7	use netcdf95, only: nf95_create, nf95_def_dim, nf95_def_var, &
247			nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
248			use netcdf, only: nf90_clobber, nf90_float, nf90_copy_att, nf90_global
249	guez	3
250	guez	7	integer, intent(in):: ncid_in, varid_in(:), n_plev
251	guez	3	character(len=*), intent(in):: name_out(:) ! of NetCDF variables
252	guez	7	integer, intent(out):: ncid_out, varid_out(:), varid_plev, varid_time
253	guez	3
254			! Variables local to the procedure:
255			integer ncerr
256	guez	7	integer dimid_rlatu, dimid_plev, dimid_time
257			integer varid_rlatu
258	guez	3	integer i, n_o3_param
259
260			!---------------------------
261
262			print *, "Call sequence information: prepare_out"
263			n_o3_param = assert_eq(size(varid_in), size(varid_out), &
264			size(name_out), "prepare_out")
265
266			call nf95_create("coefoz_LMDZ.nc", nf90_clobber, ncid_out)
267
268			! Dimensions:
269	guez	7	call nf95_def_dim(ncid_out, "time", 360, dimid_time)
270			call nf95_def_dim(ncid_out, "plev", n_plev, dimid_plev)
271	guez	3	call nf95_def_dim(ncid_out, "rlatu", jjm + 1, dimid_rlatu)
272
273	guez	7	! Define coordinate variables:
274	guez	3
275	guez	7	call nf95_def_var(ncid_out, "time", nf90_float, dimid_time, varid_time)
276			call nf95_put_att(ncid_out, varid_time, "units", "days since 2000-1-1")
277			call nf95_put_att(ncid_out, varid_time, "calendar", "360_day")
278			call nf95_put_att(ncid_out, varid_time, "standard_name", "time")
279	guez	3
280	guez	7	call nf95_def_var(ncid_out, "plev", nf90_float, dimid_plev, varid_plev)
281			call nf95_put_att(ncid_out, varid_plev, "units", "millibar")
282			call nf95_put_att(ncid_out, varid_plev, "standard_name", "air_pressure")
283			call nf95_put_att(ncid_out, varid_plev, "long_name", "air pressure")
284	guez	3
285			call nf95_def_var(ncid_out, "rlatu", nf90_float, dimid_rlatu, varid_rlatu)
286			call nf95_put_att(ncid_out, varid_rlatu, "units", "degrees_north")
287	guez	7	call nf95_put_att(ncid_out, varid_rlatu, "standard_name", "latitude")
288	guez	3
289	guez	7	! Define primary variables:
290	guez	3
291			do i = 1, n_o3_param
292			call nf95_def_var(ncid_out, name_out(i), nf90_float, &
293	guez	7	(/dimid_rlatu, dimid_plev, dimid_time/), varid_out(i))
294
295	guez	3	! The following commands may fail. That is OK. It should just
296			! mean that the attribute is not defined in the input file.
297	guez	7
298	guez	3	ncerr = nf90_copy_att(ncid_in, varid_in(i), "long_name",&
299			& ncid_out, varid_out(i))
300			call handle_err_copy_att("long_name")
301	guez	7
302	guez	3	ncerr = nf90_copy_att(ncid_in, varid_in(i), "units", ncid_out,&
303			& varid_out(i))
304			call handle_err_copy_att("units")
305	guez	7
306	guez	3	ncerr = nf90_copy_att(ncid_in, varid_in(i), "standard_name", ncid_out,&
307			& varid_out(i))
308			call handle_err_copy_att("standard_name")
309			end do
310
311			! Global attributes:
312	guez	7	call nf95_copy_att(ncid_in, nf90_global, "Conventions", ncid_out, &
313			nf90_global)
314	guez	3	call nf95_copy_att(ncid_in, nf90_global, "title", ncid_out, nf90_global)
315			call nf95_copy_att(ncid_in, nf90_global, "source", ncid_out, nf90_global)
316			call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")
317
318			call nf95_enddef(ncid_out)
319
320	guez	7	! Write one of the coordinate variables:
321			call nf95_put_var(ncid_out, varid_rlatu, rlatu(jjm+1:1:-1) / pi * 180.)
322	guez	3	! (convert from rad to degrees and sort in increasing order)
323
324			contains
325
326			subroutine handle_err_copy_att(att_name)
327
328	guez	7	use netcdf, only: nf90_noerr, nf90_strerror
329	guez	3
330			character(len=*), intent(in):: att_name
331
332			!----------------------------------------
333
334			if (ncerr /= nf90_noerr) then
335			print *, "prepare_out " // trim(name_out(i)) &
336			// " nf90_copy_att " // att_name // " -- " &
337			// trim(nf90_strerror(ncerr))
338			end if
339
340			end subroutine handle_err_copy_att
341
342			end subroutine prepare_out
343
344	guez	7	end module regr_lat_time_coefoz_m