/[lmdze]/trunk/Sources/phylmd/Mobidic/regr_pr_int.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Mobidic/regr_pr_int.f

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

revision 167 by guez, Fri Jul 24 16:54:30 2015 UTC revision 168 by guez, Wed Sep 9 10:41:47 2015 UTC
# Line 1  Line 1 
1  module regr_pr_int_m  module regr_pr_int_m
2    
3      ! Author: Lionel GUEZ
4    
5    implicit none    implicit none
6    
7  contains  contains
# Line 25  contains Line 27  contains
27    
28      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
29      use dimphy, only: klon      use dimphy, only: klon
30      use grid_change, only: dyn_phy      use grid_change, only: gr_dyn_phy
31      use netcdf95, only: nf95_inq_varid, nf95_get_var      use netcdf95, only: nf95_inq_varid, nf95_get_var
32      use nr_util, only: assert      use nr_util, only: assert
33      use numer_rec_95, only: regr1_lint      use numer_rec_95, only: regr1_lint
# Line 57  contains Line 59  contains
59      ! Field on the "physics" horizontal grid.  "v2(i, k >= 1)" is at      ! Field on the "physics" horizontal grid.  "v2(i, k >= 1)" is at
60      ! longitude "xlon(i)", latitude "xlat(i)" and pressure "plev(k)".)      ! longitude "xlon(i)", latitude "xlat(i)" and pressure "plev(k)".)
61    
62      integer i, k      integer i
63    
64      !--------------------------------------------      !--------------------------------------------
65    
# Line 74  contains Line 76  contains
76      ! Complete "v1" with the value at 0 pressure:      ! Complete "v1" with the value at 0 pressure:
77      v1(:, 0) = top_value      v1(:, 0) = top_value
78    
79      forall (k = 0:size(plev)) v2(:, k) = pack(spread(v1(:, k), dim = 1, &      v2 = gr_dyn_phy(spread(v1, dim = 1, ncopies = iim + 1))
          ncopies = iim + 1), dyn_phy)  
80    
81      ! Regrid in pressure at each horizontal position:      ! Regrid in pressure at each horizontal position:
82      do i = 1, klon      do i = 1, klon
Legend:
Removed from v.167  
changed lines
  Added in v.168
  ViewVC Help
Powered by ViewVC 1.1.21