source: lmdz_wrf/trunk/WRFV3/lmdz/regr_pr_int_m.F90 @ 1554

Last change on this file since 1554 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 3.6 KB
Line 
1! $Id$
2module regr_pr_int_m
3
4  ! Author: Lionel GUEZ
5
6  implicit none
7
8contains
9
10  subroutine regr_pr_int(ncid, name, julien, plev, pplay, top_value, v3)
11
12    ! "regr_pr_int" stands for "regrid pressure interpolation".
13    ! In this procedure:
14    ! -- the root process reads a 2D latitude-pressure field from a
15    !    NetCDF file, at a given day.
16    ! -- the field is packed to the LMDZ horizontal "physics"
17    !    grid and scattered to all threads of all processes;
18    ! -- in all the threads of all the processes, the field is regridded in
19    !    pressure to the LMDZ vertical grid.
20    ! We assume that, in the input file, the field has 3 dimensions:
21    ! latitude, pressure, julian day.
22    ! We assume that latitudes are in ascending order in the input file.
23    ! The target vertical LMDZ grid is the grid of mid-layers.
24    ! Regridding is by linear interpolation.
25
26    use dimphy, only: klon
27    use netcdf95, only: nf95_inq_varid, handle_err
28    use netcdf, only: nf90_get_var
29    use assert_m, only: assert
30    use regr1_lint_m, only: regr1_lint
31    use mod_phys_lmdz_mpi_data, only: is_mpi_root
32
33    use mod_phys_lmdz_transfert_para, only: scatter2d
34    ! (pack to the LMDZ horizontal "physics" grid and scatter)
35
36    integer, intent(in):: ncid ! NetCDF ID of the file
37    character(len=*), intent(in):: name ! of the NetCDF variable
38    integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
39
40    real, intent(in):: plev(:)
41    ! (pressure level of input data, in Pa, in strictly ascending order)
42
43    real, intent(in):: pplay(:, :) ! (klon, llm)
44    ! (pression pour le mileu de chaque couche, en Pa)
45
46    real, intent(in):: top_value
47    ! (extra value of field at 0 pressure)
48
49    real, intent(out):: v3(:, :) ! (klon, llm)
50    ! (regridded field on the partial "physics" grid)
51    ! ("v3(i, k)" is at longitude "xlon(i)", latitude
52    ! "xlat(i)", middle of layer "k".)
53
54    ! Variables local to the procedure:
55
56    include "dimensions.h"
57    integer varid, ncerr ! for NetCDF
58
59    real  v1(iim, jjm + 1, 0:size(plev))
60    ! (input field at day "julien", on the global "dynamics" horizontal grid)
61    ! (First dimension is for longitude.
62    ! The value is the same for all longitudes.
63    ! "v1(:, j, k >=1)" is at latitude "rlatu(j)" and pressure "plev(k)".)
64
65    real v2(klon, 0:size(plev))
66    ! (field scattered to the partial "physics" horizontal grid)
67    ! "v2(i, k >= 1)" is at longitude "xlon(i)", latitude "xlat(i)"
68    ! and pressure "plev(k)".)
69
70    integer i
71
72    !--------------------------------------------
73
74    call assert(shape(v3) == (/klon, llm/), "regr_pr_int v3")
75    call assert(shape(pplay) == (/klon, llm/), "regr_pr_int pplay")
76
77    !$omp master
78    if (is_mpi_root) then
79       call nf95_inq_varid(ncid, name, varid)
80
81       ! Get data at the right day from the input file:
82       ncerr = nf90_get_var(ncid, varid, v1(1, :, 1:), start=(/1, 1, julien/))
83       call handle_err("regr_pr_int nf90_get_var " // name, ncerr, ncid)
84       ! Latitudes are in ascending order in the input file while
85       ! "rlatu" is in descending order so we need to invert order:
86       v1(1, :, 1:) = v1(1, jjm+1:1:-1, 1:)
87
88       ! Complete "v1" with the value at 0 pressure:
89       v1(1, :, 0) = top_value
90
91       ! Duplicate on all longitudes:
92       v1(2:, :, :) = spread(v1(1, :, :), dim=1, ncopies=iim-1)
93    end if
94    !$omp end master
95
96    call scatter2d(v1, v2)
97
98    ! Regrid in pressure at each horizontal position:
99    do i = 1, klon
100       v3(i, llm:1:-1) = regr1_lint(v2(i, :), (/0., plev/), pplay(i, llm:1:-1))
101       ! (invert order of indices because "pplay" is in descending order)
102    end do
103
104  end subroutine regr_pr_int
105
106end module regr_pr_int_m
Note: See TracBrowser for help on using the repository browser.