/[lmdze]/trunk/libf/bibio/regr_pr.f90
ViewVC logotype

Diff of /trunk/libf/bibio/regr_pr.f90

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

revision 13 by guez, Fri Jul 25 19:59:34 2008 UTC revision 18 by guez, Thu Aug 7 12:29:13 2008 UTC
# Line 23  contains Line 23  contains
23      use numer_rec, only: assert      use numer_rec, only: assert
24      use regr1_step_av_m, only: regr1_step_av      use regr1_step_av_m, only: regr1_step_av
25      use pressure_var, only: p3d      use pressure_var, only: p3d
26        use grid_change, only: dyn_phy
27    
28      real, intent(in):: v(:, :)      real, intent(in):: v(:, :)
29      ! ("v(j, l)" is at latitude "rlatu(j)" and for pressure interval      ! ("v(j, l)" is at latitude "rlatu(j)" and for pressure interval
# Line 46  contains Line 47  contains
47      call assert(size(press_in) == size(v, 2) + 1, "regr_pr_av 2")      call assert(size(press_in) == size(v, 2) + 1, "regr_pr_av 2")
48    
49      ! Regrid in pressure by averaging a step function of pressure.      ! Regrid in pressure by averaging a step function of pressure.
50        do j = 1, jjm + 1
     ! Poles ("p3d" does not depend on longitude):  
     do j = 1, jjm + 1, jjm  
        ! Average on pressure, only for first longitude:  
        regr_pr_av(1, j, llm:1:-1) &  
             = regr1_step_av(v(j, :), press_in, 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:  
51         do i = 1, iim         do i = 1, iim
52            regr_pr_av(i, j, llm:1:-1) &            if (dyn_phy(i, j)) then
53                 = regr1_step_av(v(j, :), press_in, p3d(i, j, llm+1:1:-1))               regr_pr_av(i, j, llm:1:-1) &
54            ! (invert order of indices because "p3d" is decreasing)                    = regr1_step_av(v(j, :), press_in, p3d(i, j, llm+1:1:-1))
55                 ! (invert order of indices because "p3d" is decreasing)
56              end if
57         end do         end do
58      end do      end do
59    
# Line 96  contains Line 88  contains
88      use numer_rec, only: assert      use numer_rec, only: assert
89      use regr1_lint_m, only: regr1_lint      use regr1_lint_m, only: regr1_lint
90      use pressure_var, only: pls      use pressure_var, only: pls
91        use grid_change, only: dyn_phy
92    
93      real, intent(in):: v(:, :)      real, intent(in):: v(:, :)
94      ! ("v(j, l)" is at latitude "rlatu(j)" and pressure level "press_in(l)".)      ! ("v(j, l)" is at latitude "rlatu(j)" and pressure level "press_in(l)".)
# Line 117  contains Line 110  contains
110      call assert(size(press_in) == size(v, 2), "regr_pr_int 2")      call assert(size(press_in) == size(v, 2), "regr_pr_int 2")
111    
112      ! Regrid in pressure by linear interpolation      ! Regrid in pressure by linear interpolation
113        do j = 1, jjm + 1
     ! Poles ("pls" does not depend on longitude):  
     do j = 1, jjm + 1, jjm  
        ! Interpolate in pressure, only for first longitude:  
        regr_pr_int(1, j, llm:1:-1) &  
             = regr1_lint(v(j, :), press_in, 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:  
114         do i = 1, iim         do i = 1, iim
115            regr_pr_int(i, j, llm:1:-1) &            if (dyn_phy(i, j)) then
116                 = regr1_lint(v(j, :), press_in, pls(i, j, llm:1:-1))               regr_pr_int(i, j, llm:1:-1) &
117            ! (invert order of indices because "pls" is decreasing)                    = regr1_lint(v(j, :), press_in, pls(i, j, llm:1:-1))
118                 ! (invert order of indices because "pls" is decreasing)
119              end if
120         end do         end do
121      end do      end do
122    
Legend:
Removed from v.13  
changed lines
  Added in v.18
  ViewVC Help
Powered by ViewVC 1.1.21