/[lmdze]/trunk/phylmd/Mobidic/regr_lat_time_coefoz.f
ViewVC logotype

Diff of /trunk/phylmd/Mobidic/regr_lat_time_coefoz.f

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/phylmd/Mobidic/regr_lat_time_coefoz.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC
# Line 1  Line 1 
1  module regr_coefoz_m  module regr_lat_time_coefoz_m
2    
   ! This module is clean: no C preprocessor directive, no include line.  
3    ! Author: Lionel GUEZ    ! Author: Lionel GUEZ
4    
5    implicit none    implicit none
6    
7    private    private
8    public regr_coefoz    public regr_lat_time_coefoz
9    
10  contains  contains
11    
12    subroutine regr_coefoz    subroutine regr_lat_time_coefoz
13    
14      ! "regr_coefoz" stands for "regrid coefficients ozone".      ! "regr_lat_time_coefoz" stands for "regrid latitude time
15        ! coefficients ozone".
16    
17        ! This procedure reads from a NetCDF file parameters for ozone
18        ! chemistry, regrids them in latitude and time, and writes the
19        ! regridded fields to a new NetCDF file.
20    
     ! This procedure reads from a file parameters for ozone  
     ! chemistry and regrids them for LMDZ.  
21      ! The input fields depend on time, pressure level and      ! The input fields depend on time, pressure level and
22      ! latitude.      ! latitude.
23      ! We assume that the input fields are step functions      ! We assume that the input fields are step functions
24      ! of latitude.      ! of latitude.
25      ! Horizontal regridding is made by averaging on latitude, with a      ! Regridding in latitude is made by averaging, with a cosine of
26      ! cosine of latitude factor.      ! latitude factor.
27      ! The target horizontal LMDZ grid is the "scalar" grid: "rlonv", "rlatu".      ! The target LMDZ latitude grid is the "scalar" grid: "rlatu".
28      ! The values of "rlatu" are taken to be the centers of intervals.      ! The values of "rlatu" are taken to be the centers of intervals.
29      ! The target vertical LMDZ grid is the grid of mid-layers.      ! Regridding in time is by linear interpolation.
30      ! The input latitude values are different from the latitude values      ! Monthly values are processed to get daily values, on the basis
31      ! of the LMDZ "scalar" grid.      ! of a 360-day calendar.
32      ! The input data does not depend on longitude, but the pressure  
33      ! at LMDZ mid-layers does.      ! We assume that in the input file:
34      ! Therefore, the values on the LMDZ grid do depend on longitude.      ! -- the latitude is in degrees and strictly monotonic (as all
35      ! The regridded fields are written to a file.      ! NetCDF coordinate variables should be);
36        ! -- time increases from January to December (even though we do
37      ! 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).  
38    
39      use dimens_m, only: iim, jjm, llm      use dimens_m, only: jjm
40      use comgeom, only: rlatv      use comgeom, only: rlatv
41      use comconst, only: pi      use nr_util, only: pi
42      use pressure_m, only: p3d, pls      use numer_rec_95, only: regr1_step_av, regr3_lint
43      use regr1_step_av_m, only: regr1_step_av      use netcdf95, only: nf95_open, nf95_gw_var, nf95_close, &
44      use regr1_lint_m, only: regr1_lint           nf95_inq_varid, handle_err, nf95_put_var
45      use netcdf95, only: nf95_open, nf90_nowrite, coordin, nf95_close, &      use netcdf, only: nf90_nowrite, nf90_get_var
          nf95_inq_varid, handle_err, nf90_put_var, nf90_get_var  
46    
47      ! Variables local to the procedure:      ! Variables local to the procedure:
48    
49      integer ncid_in, ncid_out ! NetCDF IDs for input and output files      integer ncid_in, ncid_out ! NetCDF IDs for input and output files
   
50      integer n_plev ! number of pressure levels in the input data      integer n_plev ! number of pressure levels in the input data
51      integer n_lat! number of latitudes in the input data      integer n_lat! number of latitudes in the input data
52    
# Line 59  contains Line 57  contains
57      ! (edges of latitude intervals for input data, in rad, in strictly      ! (edges of latitude intervals for input data, in rad, in strictly
58      ! increasing order)      ! increasing order)
59    
60      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)  
   
61      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  
62    
63      real, allocatable:: o3_par_in(:, :, :)      real, allocatable:: o3_par_in(:, :, :) ! (n_lat, n_plev, 12)
64      ! (ozone parameter from the input file)      ! (ozone parameter from the input file)
65      ! ("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
66      ! level "plev(l)". "month" is between 1 and 12.)      ! level "plev(l)". "month" is between 1 and 12.)
67    
68      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)  
69      ! (mean of a variable "v" over a latitude interval)      ! (mean of a variable "v" over a latitude interval)
70      ! First dimension is latitude interval.      ! (First dimension is latitude interval.
71      ! The latitude interval for "v_regr_lat(j,:, :)" contains "rlatu(j)".      ! The latitude interval for "v_regr_lat(j,:, :)" contains "rlatu(j)".
72      ! If "j" is between 2 and "jjm" then the interval is:      ! If "j" is between 2 and "jjm" then the interval is:
73      ! [rlatv(j), rlatv(j-1)]      ! [rlatv(j), rlatv(j-1)]
# Line 110  contains Line 76  contains
76      ! or:      ! or:
77      ! [- pi / 2, rlatv(jjm)]      ! [- pi / 2, rlatv(jjm)]
78      ! respectively.      ! respectively.
79      ! "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.  
80      ! Last dimension is month number.)      ! Last dimension is month number.)
81    
82        real, allocatable:: o3_par_out(:, :, :) ! (jjm + 1, n_plev, 360)
83        ! (regridded ozone parameter)
84        ! ("o3_par_out(j, l, day)" is at latitude "rlatu(j)", pressure
85        ! level "plev(l)" and date "January 1st 0h" + "tmidday(day)", in a
86        ! 360-day calendar.)
87    
88        integer j
89        integer i_v ! index of ozone parameter
90        integer, parameter:: n_o3_param = 8 ! number of ozone parameters
91    
92        character(len=11) name_in(n_o3_param)
93        ! (name of NetCDF primary variable in the input file)
94    
95        character(len=9) name_out(n_o3_param)
96        ! (name of NetCDF primary variable in the output file)
97    
98        integer varid_in(n_o3_param), varid_out(n_o3_param), varid_plev, varid_time
99        integer ncerr, varid
100        ! (for NetCDF)
101    
102        real, parameter:: tmidmonth(0:13) = (/(-15. + 30. * j, j = 0, 13)/)
103        ! (time to middle of month, in days since January 1st 0h, in a
104        ! 360-day calendar)
105        ! (We add values -15 and 375 so that, for example, day 3 of the year is
106        ! interpolated between the December and the January value.)
107    
108        real, parameter:: tmidday(360) = (/(j + 0.5, j = 0, 359)/)
109        ! (time to middle of day, in days since January 1st 0h, in a
110        ! 360-day calendar)
111    
112      !---------------------------------      !---------------------------------
113    
114      print *, "Call sequence information: regr_coefoz"      print *, "Call sequence information: regr_lat_time_coefoz"
115    
116      ! Details for each ozone parameter:      ! Names of ozone parameters:
117      i_v = 0      i_v = 0
118    
119      i_v = i_v + 1      i_v = i_v + 1
# Line 132  contains Line 125  contains
125      name_out(i_v) = "a2"      name_out(i_v) = "a2"
126    
127      i_v = i_v + 1      i_v = i_v + 1
128      name_in(i_v) = "r"      name_in(i_v) = "tro3"
129      name_out(i_v) = "r_Mob"      name_out(i_v) = "r_Mob"
130    
131      i_v = i_v + 1      i_v = i_v + 1
# Line 150  contains Line 143  contains
143      i_v = i_v + 1      i_v = i_v + 1
144      name_in(i_v) = "Sigma"      name_in(i_v) = "Sigma"
145      name_out(i_v) = "Sigma_Mob"      name_out(i_v) = "Sigma_Mob"
     stepav_choice(i_v) = .false.  
     top_value(i_v) = 0.  
146    
147      i_v = i_v + 1      i_v = i_v + 1
148      name_in(i_v) = "R_Het"      name_in(i_v) = "R_Het"
149      name_out(i_v) = "R_Het"      name_out(i_v) = "R_Het"
150    
151      call nf95_open("coefoz_v2_3.nc", nf90_nowrite, ncid_in)      call nf95_open("coefoz.nc", nf90_nowrite, ncid_in)
152    
153      ! Get coordinates from the input file:      ! Get coordinates from the input file:
154    
155      latitude => coordin(ncid_in, "latitude")      call nf95_inq_varid(ncid_in, "latitude", varid)
156        call nf95_gw_var(ncid_in, varid, latitude)
157      ! Convert from degrees to rad, because "rlatv" is in rad:      ! Convert from degrees to rad, because "rlatv" is in rad:
158      latitude = latitude / 180. * pi      latitude = latitude / 180. * pi
159      n_lat = size(latitude)      n_lat = size(latitude)
160      ! We need to supply the latitudes to "stepav" in      ! We need to supply the latitudes to "regr1_step_av" in
161      ! increasing order, so invert order if necessary:      ! increasing order, so invert order if necessary:
162      decr_lat = latitude(1) > latitude(n_lat)      decr_lat = latitude(1) > latitude(n_lat)
163      if (decr_lat) latitude = latitude(n_lat:1:-1)      if (decr_lat) latitude = latitude(n_lat:1:-1)
# Line 177  contains Line 169  contains
169      lat_in_edg(n_lat + 1) = pi / 2      lat_in_edg(n_lat + 1) = pi / 2
170      deallocate(latitude) ! pointer      deallocate(latitude) ! pointer
171    
172      plev => coordin(ncid_in, "plev")      call nf95_inq_varid(ncid_in, "plev", varid)
173      ! Convert from hPa to Pa because "p3d" and "pls" are in Pa:      call nf95_gw_var(ncid_in, varid, plev)
     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           v_regr_lat(:, :, 0) = v_regr_lat(:, :, 12)
215           v_regr_lat(:, :, 13) = v_regr_lat(:, :, 1)
216    
217         if (stepav_choice(i_v)) then         ! Regrid in time by linear interpolation:
218            ! Regrid in pressure by averaging a step function of pressure         o3_par_out = regr3_lint(v_regr_lat, tmidmonth, tmidday)
   
           ! Poles ("p3d" does not depend on longitude):  
           do j = 1, jjm + 1, jjm  
              ! 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: jjm
240        use comgeom, only: rlatu
241        use nr_util, only: assert_eq, pi
242    
243      use dimens_m, only: iim, jjm, llm      use netcdf95, only: nf95_create, nf95_def_dim, nf95_def_var, &
244      use comgeom, only: rlatu, rlonv           nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
245      use pressure_m, only: pls      use netcdf, only: nf90_clobber, nf90_float, nf90_copy_att, nf90_global
     use comconst, only: pi  
     use nrutil, only: assert_eq  
   
     use netcdf95, only: nf95_create, nf90_clobber, nf95_def_dim, &  
          nf95_def_var, nf90_float, nf90_int, nf95_put_att, nf95_enddef, &  
          nf90_put_var, handle_err, nf90_copy_att, nf95_copy_att, nf90_global  
246    
247      integer, intent(in):: ncid_in, varid_in(:)      integer, intent(in):: ncid_in, varid_in(:), n_plev
248      character(len=*), intent(in):: name_out(:) ! of NetCDF variables      character(len=*), intent(in):: name_out(:) ! of NetCDF variables
249      integer, intent(out):: ncid_out, varid_out(:)      integer, intent(out):: ncid_out, varid_out(:), varid_plev, varid_time
250    
251      ! Variables local to the procedure:      ! Variables local to the procedure:
252      integer ncerr      integer ncerr
253      integer dimid_rlonv, dimid_rlatu, dimid_sigs, dimid_month      integer dimid_rlatu, dimid_plev, dimid_time
254      integer varid_rlonv, varid_rlatu, varid_sigs, varid_month      integer varid_rlatu
     integer varid_pls  
255      integer i, n_o3_param      integer i, n_o3_param
256    
257      !---------------------------      !---------------------------
# Line 330  contains Line 263  contains
263      call nf95_create("coefoz_LMDZ.nc", nf90_clobber, ncid_out)      call nf95_create("coefoz_LMDZ.nc", nf90_clobber, ncid_out)
264    
265      ! Dimensions:      ! Dimensions:
266      call nf95_def_dim(ncid_out, "month", 12, dimid_month)      call nf95_def_dim(ncid_out, "time", 360, dimid_time)
267      call nf95_def_dim(ncid_out, "sigs", llm, dimid_sigs)      call nf95_def_dim(ncid_out, "plev", n_plev, dimid_plev)
268      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)  
   
     ! Coordinate variables:  
269    
270      call nf95_def_var(ncid_out, "month", nf90_float, dimid_month, varid_month)      ! Define coordinate variables:
     call nf95_put_att(ncid_out, varid_month, "units", &  
          "calendar_month as %m.%f")  
     call nf95_put_att(ncid_out, varid_month, "long_name", "seasonal phase")  
271    
272      call nf95_def_var(ncid_out, "sigs", nf90_int, dimid_sigs, varid_sigs)      call nf95_def_var(ncid_out, "time", nf90_float, dimid_time, varid_time)
273      call nf95_put_att(ncid_out, varid_sigs, "long_name", "s-layer index")      call nf95_put_att(ncid_out, varid_time, "units", "days since 2000-1-1")
274        call nf95_put_att(ncid_out, varid_time, "calendar", "360_day")
275        call nf95_put_att(ncid_out, varid_time, "standard_name", "time")
276    
277        call nf95_def_var(ncid_out, "plev", nf90_float, dimid_plev, varid_plev)
278        call nf95_put_att(ncid_out, varid_plev, "units", "millibar")
279        call nf95_put_att(ncid_out, varid_plev, "standard_name", "air_pressure")
280        call nf95_put_att(ncid_out, varid_plev, "long_name", "air pressure")
281    
282      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)
283      call nf95_put_att(ncid_out, varid_rlatu, "units", "degrees_north")      call nf95_put_att(ncid_out, varid_rlatu, "units", "degrees_north")
284      call nf95_put_att(ncid_out, varid_rlatu, "long_name", "latitude")      call nf95_put_att(ncid_out, varid_rlatu, "standard_name", "latitude")
285    
286      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")  
287    
288      do i = 1, n_o3_param      do i = 1, n_o3_param
289         call nf95_def_var(ncid_out, name_out(i), nf90_float, &         call nf95_def_var(ncid_out, name_out(i), nf90_float, &
290              (/dimid_rlonv, dimid_rlatu, dimid_sigs, dimid_month/),&              (/dimid_rlatu, dimid_plev, dimid_time/), varid_out(i))
291              & varid_out(i))  
292         ! The following commands may fail. That is OK. It should just         ! The following commands may fail. That is OK. It should just
293         ! mean that the attribute is not defined in the input file.         ! mean that the attribute is not defined in the input file.
294    
295         ncerr = nf90_copy_att(ncid_in, varid_in(i), "long_name",&         ncerr = nf90_copy_att(ncid_in, varid_in(i), "long_name",&
296              & ncid_out, varid_out(i))              & ncid_out, varid_out(i))
297         call handle_err_copy_att("long_name")         call handle_err_copy_att("long_name")
298    
299         ncerr = nf90_copy_att(ncid_in, varid_in(i), "units", ncid_out,&         ncerr = nf90_copy_att(ncid_in, varid_in(i), "units", ncid_out,&
300              & varid_out(i))              & varid_out(i))
301         call handle_err_copy_att("units")         call handle_err_copy_att("units")
302    
303         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,&
304              & varid_out(i))              & varid_out(i))
305         call handle_err_copy_att("standard_name")         call handle_err_copy_att("standard_name")
306      end do      end do
307    
308      ! Global attributes:      ! Global attributes:
309        call nf95_copy_att(ncid_in, nf90_global, "Conventions", ncid_out, &
310             nf90_global)
311      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)
312      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)
313      call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")      call nf95_put_att(ncid_out, nf90_global, "comment", "Regridded for LMDZ")
314    
315      call nf95_enddef(ncid_out)      call nf95_enddef(ncid_out)
316    
317      ! Coordinate variables:      ! Write one of the coordinate variables:
318        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.)  
319      ! (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)  
320    
321    contains    contains
322    
323      subroutine handle_err_copy_att(att_name)      subroutine handle_err_copy_att(att_name)
324    
325        use netcdf95, only: nf90_noerr, nf90_strerror        use netcdf, only: nf90_noerr, nf90_strerror
326    
327        character(len=*), intent(in):: att_name        character(len=*), intent(in):: att_name
328    
# Line 427  contains Line 338  contains
338    
339    end subroutine prepare_out    end subroutine prepare_out
340    
341  end module regr_coefoz_m  end module regr_lat_time_coefoz_m
Legend:
Removed from v.4  
changed lines
  Added in v.76
  ViewVC Help
Powered by ViewVC 1.1.21