/[lmdze]/trunk/libf/phylmd/Mobidic/regr_lat_time_coefoz.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/Mobidic/regr_lat_time_coefoz.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/regr_coefoz.f90 revision 4 by guez, Thu Feb 28 18:05:06 2008 UTC trunk/libf/dyn3d/rrltwc.f90 revision 7 by guez, Mon Mar 31 12:24:17 2008 UTC
# Line 1  Line 1 
1  module regr_coefoz_m  module regr_lat_time_coefoz_m
2    
3    ! This module is clean: no C preprocessor directive, no include line.    ! This module is clean: no C preprocessor directive, no include line.
4    ! Author: Lionel GUEZ    ! Author: Lionel GUEZ
# Line 6  module regr_coefoz_m Line 6  module regr_coefoz_m
6    implicit none    implicit none
7    
8    private    private
9    public regr_coefoz    public read_regr_lat_time_write_coefoz
10    
11  contains  contains
12    
13    subroutine regr_coefoz    subroutine regr_lat_time_coefoz
14    
15      ! "regr_coefoz" stands for "regrid coefficients ozone".      ! "regr_lat_time_coefoz" stands for "regrid latitude time
16        ! coefficients ozone".
17    
18        ! 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    
     ! This procedure reads from a file parameters for ozone  
     ! chemistry and regrids them for LMDZ.  
22      ! The input fields depend on time, pressure level and      ! The input fields depend on time, pressure level and
23      ! latitude.      ! latitude.
24      ! We assume that the input fields are step functions      ! We assume that the input fields are step functions
25      ! of latitude.      ! of latitude.
26      ! Horizontal regridding is made by averaging on latitude, with a      ! Regridding in latitude is made by averaging, with a cosine of
27      ! cosine of latitude factor.      ! latitude factor.
28      ! The target horizontal LMDZ grid is the "scalar" grid: "rlonv", "rlatu".      ! The target LMDZ latitude grid is the "scalar" grid: "rlatu".
29      ! The values of "rlatu" are taken to be the centers of intervals.      ! The values of "rlatu" are taken to be the centers of intervals.
30      ! The target vertical LMDZ grid is the grid of mid-layers.      ! Regridding in time is by linear interpolation.
31      ! The input latitude values are different from the latitude values      ! Monthly values are processed to get daily values, on the basis
32      ! of the LMDZ "scalar" grid.      ! of a 360-day calendar.
33      ! The input data does not depend on longitude, but the pressure  
34      ! at LMDZ mid-layers does.      ! We assume that in the input file:
35      ! Therefore, the values on the LMDZ grid do depend on longitude.      ! -- the latitude is in degrees and strictly monotonic (as all
36      ! The regridded fields are written to a file.      ! NetCDF coordinate variables should be);
37        ! -- time increases from January to December (even though we do
38      ! We assume that in the input file the latitude is in degrees      ! not use values of the input time coordinate).
     ! and the pressure level is in hPa, and that both are strictly  
     ! monotonic (as all NetCDF coordinate variables should be).  
39    
40      use dimens_m, only: iim, jjm, llm      use dimens_m, only: jjm
41      use comgeom, only: rlatv      use comgeom, only: rlatv
42      use comconst, only: pi      use comconst, only: pi
     use pressure_m, only: p3d, pls  
43      use regr1_step_av_m, only: regr1_step_av      use regr1_step_av_m, only: regr1_step_av
44      use regr1_lint_m, only: regr1_lint      use regr3_lint_m, only: regr3_lint
45      use netcdf95, only: nf95_open, nf90_nowrite, coordin, nf95_close, &      use netcdf95, only: nf95_open, nf95_get_coord, nf95_close, &
46           nf95_inq_varid, handle_err, nf90_put_var, nf90_get_var           nf95_inq_varid, handle_err, nf95_put_var
47        use netcdf, only: nf90_nowrite, nf90_get_var
48    
49      ! Variables local to the procedure:      ! Variables local to the procedure:
50    
51      integer ncid_in, ncid_out ! NetCDF IDs for input and output files      integer ncid_in, ncid_out ! NetCDF IDs for input and output files
   
52      integer n_plev ! number of pressure levels in the input data      integer n_plev ! number of pressure levels in the input data
53      integer n_lat! number of latitudes in the input data      integer n_lat! number of latitudes in the input data
54    
# Line 59  contains Line 59  contains
59      ! (edges of latitude intervals for input data, in rad, in strictly      ! (edges of latitude intervals for input data, in rad, in strictly
60      ! increasing order)      ! increasing order)
61    
62      real, pointer:: plev(:)      real, pointer:: plev(:) ! pressure level of input data
     ! (pressure level of input data, converted to Pa, sorted  
     ! in strictly increasing order)  
   
     real, allocatable:: press_in_edg(:)  
     ! (edges of pressure intervals for input data, in Pa, in strictly  
     ! increasing order)  
   
63      logical decr_lat ! decreasing latitude in the input file      logical decr_lat ! decreasing latitude in the input file
     logical decr_plev ! decreasing pressure level in the input file  
64    
65      real, allocatable:: o3_par_in(:, :, :)      real, allocatable:: o3_par_in(:, :, :) ! (n_lat, n_plev, 12)
66      ! (ozone parameter from the input file)      ! (ozone parameter from the input file)
67      ! ("o3_par_in(j, l, month)" is at latitude "latitude(j)" and pressure      ! ("o3_par_in(j, l, month)" is at latitude "latitude(j)" and pressure
68      ! level "plev(l)". "month" is between 1 and 12.)      ! level "plev(l)". "month" is between 1 and 12.)
69    
70      real o3_par_out(iim + 1, jjm + 1, llm, 12)      real, allocatable:: v_regr_lat(:, :, :) ! (jjm + 1, n_plev, 0:13)
     ! (ozone parameter adapted to the LMDZ grid)  
     ! (Last dimension is month number.  
     ! "o3_par_out(i, j, l, month)" is at longitude "rlonv(i)", latitude  
     ! "rlatu(j)" and pressure level "pls(i, j, l)")  
   
     integer i, 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 variable in the input file)  
   
     character(len=9) name_out(n_o3_param)  
     ! (name of NetCDF variable in the output file)  
   
     logical:: stepav_choice(n_o3_param) = .true.  
     ! (vertical regridding by step average, otherwise linear interpolation)  
   
     real top_value(n_o3_param)  
     ! (value at 0 pressure, only used for linear interpolation)  
   
     integer varid_in(n_o3_param), varid_out(n_o3_param), ncerr ! for NetCDF  
   
     real, allocatable:: v_regr_lat(:, :, :) ! (jjm + 1, 0:n_plev, 12)  
71      ! (mean of a variable "v" over a latitude interval)      ! (mean of a variable "v" over a latitude interval)
72      ! First dimension is latitude interval.      ! (First dimension is latitude interval.
73      ! The latitude interval for "v_regr_lat(j,:, :)" contains "rlatu(j)".      ! The latitude interval for "v_regr_lat(j,:, :)" contains "rlatu(j)".
74      ! If "j" is between 2 and "jjm" then the interval is:      ! If "j" is between 2 and "jjm" then the interval is:
75      ! [rlatv(j), rlatv(j-1)]      ! [rlatv(j), rlatv(j-1)]
# Line 110  contains Line 78  contains
78      ! or:      ! or:
79      ! [- pi / 2, rlatv(jjm)]      ! [- pi / 2, rlatv(jjm)]
80      ! respectively.      ! respectively.
81      ! "v_regr_lat(:, l, :)" is for pressure interval      ! "v_regr_lat(:, l, :)" is for pressure level "plev(l)".
     ! "[press_in_edg(l), press_in_edg(l+1)]" or pressure level "plev(l)",  
     ! depending on the type of vertical regridding, step average or linear  
     ! interpolation.  
82      ! Last dimension is month number.)      ! Last dimension is month number.)
83    
84        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        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        ! (name of NetCDF primary variable in the input file)
96    
97        character(len=9) name_out(n_o3_param)
98        ! (name of NetCDF primary variable in the output file)
99    
100        integer varid_in(n_o3_param), varid_out(n_o3_param), varid_plev, varid_time
101        integer ncerr
102        ! (for NetCDF)
103    
104        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    
110        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    
114      !---------------------------------      !---------------------------------
115    
116      print *, "Call sequence information: regr_coefoz"      print *, "Call sequence information: regr_coefoz"
117    
118      ! Details for each ozone parameter:      ! Names of ozone parameters:
119      i_v = 0      i_v = 0
120    
121      i_v = i_v + 1      i_v = i_v + 1
# Line 150  contains Line 145  contains
145      i_v = i_v + 1      i_v = i_v + 1
146      name_in(i_v) = "Sigma"      name_in(i_v) = "Sigma"
147      name_out(i_v) = "Sigma_Mob"      name_out(i_v) = "Sigma_Mob"
     stepav_choice(i_v) = .false.  
     top_value(i_v) = 0.  
148    
149      i_v = i_v + 1      i_v = i_v + 1
150      name_in(i_v) = "R_Het"      name_in(i_v) = "R_Het"
# Line 161  contains Line 154  contains
154    
155      ! Get coordinates from the input file:      ! Get coordinates from the input file:
156    
157      latitude => coordin(ncid_in, "latitude")      call nf95_get_coord(ncid_in, "latitude", latitude)
158      ! Convert from degrees to rad, because "rlatv" is in rad:      ! Convert from degrees to rad, because "rlatv" is in rad:
159      latitude = latitude / 180. * pi      latitude = latitude / 180. * pi
160      n_lat = size(latitude)      n_lat = size(latitude)
161      ! We need to supply the latitudes to "stepav" in      ! We need to supply the latitudes to "regr1_step_av" in
162      ! increasing order, so invert order if necessary:      ! increasing order, so invert order if necessary:
163      decr_lat = latitude(1) > latitude(n_lat)      decr_lat = latitude(1) > latitude(n_lat)
164      if (decr_lat) latitude = latitude(n_lat:1:-1)      if (decr_lat) latitude = latitude(n_lat:1:-1)
# Line 177  contains Line 170  contains
170      lat_in_edg(n_lat + 1) = pi / 2      lat_in_edg(n_lat + 1) = pi / 2
171      deallocate(latitude) ! pointer      deallocate(latitude) ! pointer
172    
173      plev => coordin(ncid_in, "plev")      call nf95_get_coord(ncid_in, "plev", plev)
     ! Convert from hPa to Pa because "p3d" and "pls" are in Pa:  
     plev = plev * 100.  
174      n_plev = size(plev)      n_plev = size(plev)
175      ! We need to supply the pressure levels to "stepav" in      ! (We only need the pressure coordinate to copy it to the output file.)
     ! increasing order, so invert order if necessary:  
     decr_plev = plev(1) > plev(n_plev)  
     if (decr_plev) plev = plev(n_plev:1:-1)  
   
     ! Compute edges of pressure intervals:  
     allocate(press_in_edg(n_plev + 1))  
     press_in_edg(1) = 0.  
     ! We choose edges halfway in logarithm:  
     forall (j = 2:n_plev) press_in_edg(j) = sqrt(plev(j - 1) * plev(j))  
     press_in_edg(n_plev + 1) = huge(0.)  
     ! (infinity, but any value guaranteed to be greater than the  
     ! surface pressure would do)  
176    
177      ! Get the IDs of ozone parameters in the input file:      ! Get the IDs of ozone parameters in the input file:
178      do i_v = 1, n_o3_param      do i_v = 1, n_o3_param
179         call nf95_inq_varid(ncid_in, trim(name_in(i_v)), varid_in(i_v))         call nf95_inq_varid(ncid_in, trim(name_in(i_v)), varid_in(i_v))
180      end do      end do
181    
182      call prepare_out(ncid_in, varid_in, name_out, ncid_out, varid_out)      ! Create the output file and get the variable IDs:
183      allocate(o3_par_in(n_lat, n_plev, 12), v_regr_lat(jjm + 1, 0:n_plev, 12))      call prepare_out(ncid_in, varid_in, n_plev, name_out, ncid_out, &
184             varid_out, varid_plev, varid_time)
185    
186        ! Write remaining coordinate variables:
187        call nf95_put_var(ncid_out, varid_time, tmidday)
188        call nf95_put_var(ncid_out, varid_plev, plev)
189    
190        deallocate(plev) ! pointer
191    
192        allocate(o3_par_in(n_lat, n_plev, 12))
193        allocate(v_regr_lat(jjm + 1, n_plev, 0:13))
194        allocate(o3_par_out(jjm + 1, n_plev, 360))
195    
196      do i_v = 1, n_o3_param      do i_v = 1, n_o3_param
197         ! Process ozone parameter "name_in(i_v)"         ! Process ozone parameter "name_in(i_v)"
# Line 210  contains Line 200  contains
200         call handle_err("nf90_get_var", ncerr, ncid_in)         call handle_err("nf90_get_var", ncerr, ncid_in)
201    
202         if (decr_lat) o3_par_in = o3_par_in(n_lat:1:-1, :, :)         if (decr_lat) o3_par_in = o3_par_in(n_lat:1:-1, :, :)
        if (decr_plev) o3_par_in = o3_par_in(:, n_plev:1:-1, :)  
203    
204         ! Regrid in latitude:         ! Regrid in latitude:
205         ! We average with respect to sine of latitude, which is         ! We average with respect to sine of latitude, which is
206         ! equivalent to weighting by cosine of latitude:         ! equivalent to weighting by cosine of latitude:
207         v_regr_lat(jjm+1:1:-1, 1:, :) = regr1_step_av(o3_par_in, &         v_regr_lat(jjm+1:1:-1, :, 1:12) = regr1_step_av(o3_par_in, &
208              xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))              xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
209         ! (invert order of indices in "v_regr_lat" because "rlatu" is         ! (invert order of indices in "v_regr_lat" because "rlatu" is
210         ! decreasing)         ! decreasing)
211    
212         ! Regrid in pressure at each horizontal position:         ! Duplicate January and December values, in preparation of time
213           ! interpolation:
214         if (stepav_choice(i_v)) then         v_regr_lat(:, :, 0) = v_regr_lat(:, :, 12)
215            ! Regrid in pressure by averaging a step function of pressure         v_regr_lat(:, :, 13) = v_regr_lat(:, :, 1)
216    
217            ! Poles ("p3d" does not depend on longitude):         ! Regrid in time by linear interpolation:
218            do j = 1, jjm + 1, jjm         o3_par_out = regr3_lint(v_regr_lat, tmidmonth, tmidday)
              ! Average on pressure, only for first longitude:  
              o3_par_out(1, j, llm:1:-1, :) &  
                   = regr1_step_av(v_regr_lat(j, 1:, :), press_in_edg, &  
                   p3d(1, j, llm+1:1:-1))  
              ! (invert order of indices because "p3d" is decreasing)  
           end do  
   
           ! Latitudes other than poles ("p3d" depends on longitude):  
           do j = 2, jjm  
              ! Average on pressure at each longitude:  
              do i = 1, iim  
                 o3_par_out(i, j, llm:1:-1, :) &  
                      = regr1_step_av(v_regr_lat(j, 1:, :), press_in_edg, &  
                      p3d(i, j, llm+1:1:-1))  
                 ! (invert order of indices because "p3d" is decreasing)  
              end do  
           end do  
        else  
           ! Regrid in pressure by linear interpolation  
   
           ! Complete "v_regr_lat" with the value at 0 pressure:  
           v_regr_lat(:, 0, :) = top_value(i_v)  
   
           ! Poles ("pls" does not depend on longitude):  
           do j = 1, jjm + 1, jjm  
              ! Interpolate in pressure, only for first longitude:  
              o3_par_out(1, j, llm:1:-1, :) = regr1_lint(v_regr_lat(j, :, :), &  
                   xs=(/0., plev/), xt=pls(1, j, llm:1:-1))  
              ! (invert order of indices because "pls" is decreasing)  
           end do  
   
           ! Latitudes other than poles ("pls" depends on longitude):  
           do j = 2, jjm  
              ! Average on pressure at each longitude:  
              do i = 1, iim  
                 o3_par_out(i, j, llm:1:-1, :) &  
                      = regr1_lint(v_regr_lat(j, :, :), xs=(/0., plev/), &  
                      xt=pls(i, j, llm:1:-1))  
                 ! (invert order of indices because "pls" is decreasing)  
              end do  
           end do  
        end if  
   
        ! Duplicate pole values on all longitudes:  
        o3_par_out(2:, 1, :, :) &  
             = spread(o3_par_out(1, 1, :, :), dim=1, ncopies=iim)  
        o3_par_out(2:, jjm + 1, :, :) &  
             = spread(o3_par_out(1, jjm + 1, :, :), dim=1, ncopies=iim)  
   
        ! Duplicate first longitude to last longitude:  
        o3_par_out(iim + 1, 2:jjm, :, :) = o3_par_out(1, 2:jjm, :, :)  
219    
220         ! Write to file:         ! Write to file:
221           call nf95_put_var(ncid_out, varid_out(i_v), &
222         ncerr = nf90_put_var(ncid_out, varid_out(i_v), &              o3_par_out(jjm+1:1:-1, :, :))
             o3_par_out(:,jjm+1:1:-1, :, :))  
223         ! (The order of "rlatu" is inverted in the output file)         ! (The order of "rlatu" is inverted in the output file)
        call handle_err("nf90_put_var", ncerr, ncid_out)  
224      end do      end do
225    
     deallocate(plev) ! pointer  
226      call nf95_close(ncid_out)      call nf95_close(ncid_out)
227      call nf95_close(ncid_in)      call nf95_close(ncid_in)
228    
229    end subroutine regr_coefoz    end subroutine regr_lat_time_coefoz
230    
231    !********************************************    !********************************************
232    
233    subroutine prepare_out(ncid_in, varid_in, name_out, ncid_out, varid_out)    subroutine prepare_out(ncid_in, varid_in, n_plev, name_out, ncid_out, &
234           varid_out, varid_plev, varid_time)
235    
236      ! This subroutine creates the NetCDF output file, defines      ! This subroutine creates the NetCDF output file, defines
237      ! dimensions and variables, and writes coordinate variables and "pls".      ! dimensions and variables, and writes one of the coordinate variables.
238    
239      use dimens_m, only: iim, jjm, llm      use dimens_m, only: jjm
240      use comgeom, only: rlatu, rlonv      use comgeom, only: rlatu
     use pressure_m, only: pls  
241      use comconst, only: pi      use comconst, only: pi
242      use nrutil, only: assert_eq      use nrutil, only: assert_eq
243    
244      use netcdf95, only: nf95_create, nf90_clobber, nf95_def_dim, &      use netcdf95, only: nf95_create, nf95_def_dim, nf95_def_var, &
245           nf95_def_var, nf90_float, nf90_int, nf95_put_att, nf95_enddef, &           nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
246           nf90_put_var, handle_err, nf90_copy_att, nf95_copy_att, nf90_global      use netcdf, only: nf90_clobber, nf90_float, nf90_copy_att, nf90_global
247    
248      integer, intent(in):: ncid_in, varid_in(:)      integer, intent(in):: ncid_in, varid_in(:), n_plev
249      character(len=*), intent(in):: name_out(:) ! of NetCDF variables      character(len=*), intent(in):: name_out(:) ! of NetCDF variables
250      integer, intent(out):: ncid_out, varid_out(:)      integer, intent(out):: ncid_out, varid_out(:), varid_plev, varid_time
251    
252      ! Variables local to the procedure:      ! Variables local to the procedure:
253      integer ncerr      integer ncerr
254      integer dimid_rlonv, dimid_rlatu, dimid_sigs, dimid_month      integer dimid_rlatu, dimid_plev, dimid_time
255      integer varid_rlonv, varid_rlatu, varid_sigs, varid_month      integer varid_rlatu
     integer varid_pls  
256      integer i, n_o3_param      integer i, n_o3_param
257    
258      !---------------------------      !---------------------------
# Line 330  contains Line 264  contains
264      call nf95_create("coefoz_LMDZ.nc", nf90_clobber, ncid_out)      call nf95_create("coefoz_LMDZ.nc", nf90_clobber, ncid_out)
265    
266      ! Dimensions:      ! Dimensions:
267      call nf95_def_dim(ncid_out, "month", 12, dimid_month)      call nf95_def_dim(ncid_out, "time", 360, dimid_time)
268      call nf95_def_dim(ncid_out, "sigs", llm, dimid_sigs)      call nf95_def_dim(ncid_out, "plev", n_plev, dimid_plev)
269      call nf95_def_dim(ncid_out, "rlatu", jjm + 1, dimid_rlatu)      call nf95_def_dim(ncid_out, "rlatu", jjm + 1, dimid_rlatu)
     call nf95_def_dim(ncid_out, "rlonv", iim + 1, dimid_rlonv)  
270    
271      ! Coordinate variables:      ! Define coordinate variables:
272    
273      call nf95_def_var(ncid_out, "month", nf90_float, dimid_month, varid_month)      call nf95_def_var(ncid_out, "time", nf90_float, dimid_time, varid_time)
274      call nf95_put_att(ncid_out, varid_month, "units", &      call nf95_put_att(ncid_out, varid_time, "units", "days since 2000-1-1")
275           "calendar_month as %m.%f")      call nf95_put_att(ncid_out, varid_time, "calendar", "360_day")
276      call nf95_put_att(ncid_out, varid_month, "long_name", "seasonal phase")      call nf95_put_att(ncid_out, varid_time, "standard_name", "time")
277    
278      call nf95_def_var(ncid_out, "sigs", nf90_int, dimid_sigs, varid_sigs)      call nf95_def_var(ncid_out, "plev", nf90_float, dimid_plev, varid_plev)
279      call nf95_put_att(ncid_out, varid_sigs, "long_name", "s-layer index")      call nf95_put_att(ncid_out, varid_plev, "units", "millibar")
280        call nf95_put_att(ncid_out, varid_plev, "standard_name", "air_pressure")
281        call nf95_put_att(ncid_out, varid_plev, "long_name", "air pressure")
282    
283      call nf95_def_var(ncid_out, "rlatu", nf90_float, dimid_rlatu, varid_rlatu)      call nf95_def_var(ncid_out, "rlatu", nf90_float, dimid_rlatu, varid_rlatu)
284      call nf95_put_att(ncid_out, varid_rlatu, "units", "degrees_north")      call nf95_put_att(ncid_out, varid_rlatu, "units", "degrees_north")
285      call nf95_put_att(ncid_out, varid_rlatu, "long_name", "latitude")      call nf95_put_att(ncid_out, varid_rlatu, "standard_name", "latitude")
286    
287      call nf95_def_var(ncid_out, "rlonv", nf90_float, dimid_rlonv, varid_rlonv)      ! Define primary variables:
     call nf95_put_att(ncid_out, varid_rlonv, "units", "degrees_east")  
     call nf95_put_att(ncid_out, varid_rlonv, "long_name", "longitude")  
   
     ! Primary variables:  
   
     call nf95_def_var(ncid_out, "pls", nf90_float, &  
          (/dimid_rlonv, dimid_rlatu, dimid_sigs/), varid_pls)  
     call nf95_put_att(ncid_out, varid_pls, "units", "millibar")  
     call nf95_put_att(ncid_out, varid_pls, "long_name", &  
          "pressure at LMDZ mid-layers")  
288    
289      do i = 1, n_o3_param      do i = 1, n_o3_param
290         call nf95_def_var(ncid_out, name_out(i), nf90_float, &         call nf95_def_var(ncid_out, name_out(i), nf90_float, &
291              (/dimid_rlonv, dimid_rlatu, dimid_sigs, dimid_month/),&              (/dimid_rlatu, dimid_plev, dimid_time/), varid_out(i))
292              & varid_out(i))  
293         ! The following commands may fail. That is OK. It should just         ! The following commands may fail. That is OK. It should just
294         ! mean that the attribute is not defined in the input file.         ! mean that the attribute is not defined in the input file.
295    
296         ncerr = nf90_copy_att(ncid_in, varid_in(i), "long_name",&         ncerr = nf90_copy_att(ncid_in, varid_in(i), "long_name",&
297              & ncid_out, varid_out(i))              & ncid_out, varid_out(i))
298         call handle_err_copy_att("long_name")         call handle_err_copy_att("long_name")
299    
300         ncerr = nf90_copy_att(ncid_in, varid_in(i), "units", ncid_out,&         ncerr = nf90_copy_att(ncid_in, varid_in(i), "units", ncid_out,&
301              & varid_out(i))              & varid_out(i))
302         call handle_err_copy_att("units")         call handle_err_copy_att("units")
303    
304         ncerr = nf90_copy_att(ncid_in, varid_in(i), "standard_name", ncid_out,&         ncerr = nf90_copy_att(ncid_in, varid_in(i), "standard_name", ncid_out,&
305              & varid_out(i))              & varid_out(i))
306         call handle_err_copy_att("standard_name")         call handle_err_copy_att("standard_name")
307      end do      end do
308    
309      ! Global attributes:      ! Global attributes:
310        call nf95_copy_att(ncid_in, nf90_global, "Conventions", ncid_out, &
311             nf90_global)
312      call nf95_copy_att(ncid_in, nf90_global, "title", ncid_out, nf90_global)      call nf95_copy_att(ncid_in, nf90_global, "title", ncid_out, nf90_global)
313      call nf95_copy_att(ncid_in, nf90_global, "source", ncid_out, nf90_global)      call nf95_copy_att(ncid_in, nf90_global, "source", ncid_out, nf90_global)
314      call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")      call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")
315    
316      call nf95_enddef(ncid_out)      call nf95_enddef(ncid_out)
317    
318      ! Coordinate variables:      ! Write one of the coordinate variables:
319        call nf95_put_var(ncid_out, varid_rlatu, rlatu(jjm+1:1:-1) / pi * 180.)
     ncerr = nf90_put_var(ncid_out, varid_month, (/(i + 0.5, i = 1, 12)/))  
     call handle_err("nf90_put_var", ncerr, ncid_out)  
   
     ncerr = nf90_put_var(ncid_out, varid_sigs, (/(i, i = 1, llm)/))  
     call handle_err("nf90_put_var", ncerr, ncid_out)  
   
     ncerr = nf90_put_var(ncid_out, varid_rlatu, rlatu(jjm+1:1:-1) / pi * 180.)  
320      ! (convert from rad to degrees and sort in increasing order)      ! (convert from rad to degrees and sort in increasing order)
     call handle_err("nf90_put_var", ncerr, ncid_out)  
   
     ncerr = nf90_put_var(ncid_out, varid_rlonv, rlonv / pi * 180.)  
     ! (convert from rad to degrees)  
     call handle_err("nf90_put_var", ncerr, ncid_out)  
   
     ! Primary variable:  
   
     ncerr = nf90_put_var(ncid_out, varid_pls, pls(:, jjm+1:1:-1, :) / 100.)  
     ! (convert from Pa to hPa)  
     call handle_err("nf90_put_var", ncerr, ncid_out)  
321    
322    contains    contains
323    
324      subroutine handle_err_copy_att(att_name)      subroutine handle_err_copy_att(att_name)
325    
326        use netcdf95, only: nf90_noerr, nf90_strerror        use netcdf, only: nf90_noerr, nf90_strerror
327    
328        character(len=*), intent(in):: att_name        character(len=*), intent(in):: att_name
329    
# Line 427  contains Line 339  contains
339    
340    end subroutine prepare_out    end subroutine prepare_out
341    
342  end module regr_coefoz_m  end module regr_lat_time_coefoz_m
Legend:
Removed from v.4  
changed lines
  Added in v.7
  ViewVC Help
Powered by ViewVC 1.1.21