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

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

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

revision 18 by guez, Tue Aug 5 13:31:32 2008 UTC revision 19 by guez, Thu Aug 7 15:46:20 2008 UTC
# Line 1  Line 1 
1  module regr_pr_coefoz  module regr_pr_coefoz
2    
3      ! Both procedures of this module read a single Mobidic ozone
4      ! coefficient from "coefoz_LMDZ.nc", at the current day, regrid this
5      ! coefficient in pressure to the LMDZ vertical grid and pack it to the LMDZ
6      ! horizontal "physics" grid.
7      ! The input data is a 2D latitude -- pressure field.
8      ! The target horizontal LMDZ grid is the "scalar" grid: "rlonv", "rlatu".
9      ! We assume that the input data is already on the LMDZ "rlatu"
10      ! latitude grid.
11    
12    implicit none    implicit none
13    
14  contains  contains
# Line 8  contains Line 17  contains
17    
18      ! "regr_pr_av_coefoz" stands for "regrid pressure averaging      ! "regr_pr_av_coefoz" stands for "regrid pressure averaging
19      ! coefficient ozone".      ! coefficient ozone".
20      ! This procedure reads a single Mobidic ozone coefficient from      ! The target vertical LMDZ grid is the grid of layer boundaries.
21      ! "coefoz_LMDZ.nc", at the current day, regrids this parameter in      ! The input data does not depend on longitude, but the pressure
22      ! pressure to the LMDZ vertical grid and packs it to the LMDZ      ! at LMDZ layers does.
23      ! horizontal "physics" grid.      ! Therefore, the values on the LMDZ grid do depend on longitude.
24      ! Regridding in pressure is done by averaging a step function.      ! Regridding in pressure is done by averaging a step function.
25    
26      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
# Line 19  contains Line 28  contains
28      use netcdf95, only: nf95_inq_varid, handle_err      use netcdf95, only: nf95_inq_varid, handle_err
29      use netcdf, only: nf90_get_var      use netcdf, only: nf90_get_var
30      use grid_change, only: dyn_phy      use grid_change, only: dyn_phy
     use regr_pr, only: regr_pr_av  
31      use numer_rec, only: assert      use numer_rec, only: assert
32      use press_coefoz_m, only: press_in_edg      use press_coefoz_m, only: press_in_edg
33        use regr1_step_av_m, only: regr1_step_av
34        use pressure_var, only: p3d
35    
36      integer, intent(in):: ncid ! NetCDF ID of the file      integer, intent(in):: ncid ! NetCDF ID of the file
37      character(len=*), intent(in):: name ! of the NetCDF variable      character(len=*), intent(in):: name ! of the NetCDF variable
# Line 34  contains Line 44  contains
44    
45      ! Variables local to the procedure:      ! Variables local to the procedure:
46      integer varid, ncerr      integer varid, ncerr
47      integer k      integer i, j, k
48    
49      real  v1(jjm + 1, size(press_in_edg) - 1)      real  v1(jjm + 1, size(press_in_edg) - 1)
50      ! (ozone coefficient from "coefoz_LMDZ.nc" at day "julien")      ! (ozone coefficient from "coefoz_LMDZ.nc" at day "julien")
# Line 42  contains Line 52  contains
52      ! pressure interval "[press_in_edg(k), press_in_edg(k+1)]".)      ! pressure interval "[press_in_edg(k), press_in_edg(k+1)]".)
53    
54      real v2(iim + 1, jjm + 1, llm)      real v2(iim + 1, jjm + 1, llm)
55      ! (ozone parameter from Mobidic on the "dynamics" grid)      ! (ozone coefficient from Mobidic on the "dynamics" grid)
56      ! "v2(i, j, k)" is at longitude "rlonv(i)", latitude      ! "v2(i, j, k)" is at longitude "rlonv(i)", latitude
57      ! "rlatu(j)", middle of layer "k".)      ! "rlatu(j)", middle of layer "k".)
58    
# Line 60  contains Line 70  contains
70      v1 = v1(jjm+1:1:-1, :)      v1 = v1(jjm+1:1:-1, :)
71    
72      ! Regrid in pressure at each horizontal position:      ! Regrid in pressure at each horizontal position:
73      v2 = regr_pr_av(v1, press_in_edg)      do j = 1, jjm + 1
74           do i = 1, iim
75              if (dyn_phy(i, j)) then
76                 v2(i, j, llm:1:-1) &
77                      = regr1_step_av(v1(j, :), press_in_edg, &
78                      p3d(i, j, llm+1:1:-1))
79                 ! (invert order of indices because "p3d" is decreasing)
80              end if
81           end do
82        end do
83    
84      forall (k = 1:llm) v3(:, k) = pack(v2(:, :, k), dyn_phy)      forall (k = 1:llm) v3(:, k) = pack(v2(:, :, k), dyn_phy)
85    
# Line 70  contains Line 89  contains
89    
90    subroutine regr_pr_int_coefoz(ncid, name, julien, top_value, v3)    subroutine regr_pr_int_coefoz(ncid, name, julien, top_value, v3)
91    
92      ! This procedure reads a single Mobidic ozone coefficient from      ! "regr_pr_int_coefoz" stands for "regrid pressure interpolation
93      !"coefoz_LMDZ.nc", at the current day, regrids this parameter in      ! coefficient ozone".
94      ! pressure to the LMDZ vertical grid and packs it to the LMDZ      ! The target vertical LMDZ grid is the grid of mid-layers.
95      ! horizontal "physics" grid.      ! The input data does not depend on longitude, but the pressure
96        ! at LMDZ mid-layers does.
97        ! Therefore, the values on the LMDZ grid do depend on longitude.
98      ! Regridding is by linear interpolation.      ! Regridding is by linear interpolation.
99    
100      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
# Line 81  contains Line 102  contains
102      use netcdf95, only: nf95_inq_varid, handle_err      use netcdf95, only: nf95_inq_varid, handle_err
103      use netcdf, only: nf90_get_var      use netcdf, only: nf90_get_var
104      use grid_change, only: dyn_phy      use grid_change, only: dyn_phy
     use regr_pr, only: regr_pr_int  
105      use numer_rec, only: assert      use numer_rec, only: assert
106      use press_coefoz_m, only: plev      use press_coefoz_m, only: plev
107        use regr1_lint_m, only: regr1_lint
108        use pressure_var, only: pls
109    
110      integer, intent(in):: ncid ! NetCDF ID of the file      integer, intent(in):: ncid ! NetCDF ID of the file
111      character(len=*), intent(in):: name ! of the NetCDF variable      character(len=*), intent(in):: name ! of the NetCDF variable
# Line 93  contains Line 115  contains
115      ! (extra value of ozone coefficient at 0 pressure)      ! (extra value of ozone coefficient at 0 pressure)
116    
117      real, intent(out):: v3(:, :) ! (klon, llm)      real, intent(out):: v3(:, :) ! (klon, llm)
118      ! (ozone parameter from Mobidic on the "physics" grid)      ! (ozone coefficient from Mobidic on the "physics" grid)
119      ! ("v3(i, k)" is at longitude "xlon(i)", latitude      ! ("v3(i, k)" is at longitude "xlon(i)", latitude
120      ! "xlat(i)", middle of layer "k".)      ! "xlat(i)", middle of layer "k".)
121    
122      ! Variables local to the procedure:      ! Variables local to the procedure:
123      integer varid, ncerr      integer varid, ncerr
124      integer k      integer i, j, k
125    
126      real  v1(jjm + 1, 0:size(plev))      real  v1(jjm + 1, 0:size(plev))
127      ! (ozone coefficient from "coefoz_LMDZ.nc" at day "julien")      ! (ozone coefficient from "coefoz_LMDZ.nc" at day "julien")
128      ! ("v1(j, k >=1)" is at latitude "rlatu(j)" and pressure "plev(k)".)      ! ("v1(j, k >=1)" is at latitude "rlatu(j)" and pressure "plev(k)".)
129    
130      real v2(iim + 1, jjm + 1, llm)      real v2(iim + 1, jjm + 1, llm)
131      ! (ozone parameter from Mobidic on the "dynamics" grid)      ! (ozone coefficient from Mobidic on the "dynamics" grid)
132      ! "v2(i, j, k)" is at longitude "rlonv(i)", latitude      ! "v2(i, j, k)" is at longitude "rlonv(i)", latitude
133      ! "rlatu(j)", middle of layer "k".)      ! "rlatu(j)", middle of layer "k".)
134    
# Line 127  contains Line 149  contains
149      v1(:, 0) = top_value      v1(:, 0) = top_value
150    
151      ! Regrid in pressure at each horizontal position:      ! Regrid in pressure at each horizontal position:
152      v2 = regr_pr_int(v1, (/0., plev/))      do j = 1, jjm + 1
153           do i = 1, iim
154              if (dyn_phy(i, j)) then
155                 v2(i, j, llm:1:-1) &
156                      = regr1_lint(v1(j, :), (/0., plev/), pls(i, j, llm:1:-1))
157                 ! (invert order of indices because "pls" is decreasing)
158              end if
159           end do
160        end do
161    
162      forall (k = 1:llm) v3(:, k) = pack(v2(:, :, k), dyn_phy)      forall (k = 1:llm) v3(:, k) = pack(v2(:, :, k), dyn_phy)
163    
Legend:
Removed from v.18  
changed lines
  Added in v.19
  ViewVC Help
Powered by ViewVC 1.1.21