Changeset 2440 for LMDZ5/trunk/libf/phylmd/regr_pr_av_m.F90

Timestamp:

Feb 9, 2016, 3:45:31 PM (8 years ago)

Author:

lguez

Message:

For read_climoz = 1 or 2, replaced first order conservative regridding
of ozone by second order conservative regridding, with Van Leer
slope-limiting. The replacement is done for both latitude and pressure
regridding. The replacement is beneficial if the resolution of the
input data is coarser than the resolution of LMDZ. If the resolution
of the input data is finer, then the replacement is neutral, it does
not change much.

File:

• LMDZ5/trunk/libf/phylmd/regr_pr_av_m.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/regr_pr_av_m.F90

@@ -26 +26 @@
 
     ! The target vertical LMDZ grid is the grid of layer boundaries.
-    ! Regridding in pressure is done by averaging a step function of pressure.
+    ! Regridding in pressure is conservative, second order.
 
     ! All the fields are regridded as a single multi-dimensional array
@@ -38 +38 @@
     use assert_m, only: assert
     use assert_eq_m, only: assert_eq
-    use regr1_step_av_m, only: regr1_step_av
+    use regr1_conserv_m, only: regr1_conserv
+    use slopes_m, only: slopes
     use mod_phys_lmdz_mpi_data, only: is_mpi_root
     use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, nbp_lev
@@ -83 +84 @@
     !--------------------------------------------
 
-    call assert(size(v3, 1) == klon, size(v3, 2) == nbp_lev, "regr_pr_av v3 klon")
+    call assert(size(v3, 1) == klon, size(v3, 2) == nbp_lev, &
+         "regr_pr_av v3 klon")
     n_var = assert_eq(size(name), size(v3, 3), "regr_pr_av v3 n_var")
     call assert(shape(paprs) == (/klon, nbp_lev+1/), "regr_pr_av paprs")
@@ -112 +114 @@
     ! Regrid in pressure at each horizontal position:
     do i = 1, klon
-       v3(i, nbp_lev:1:-1, :) = regr1_step_av(v2(i, :, :), press_in_edg, &
-            paprs(i, nbp_lev+1:1:-1))
+       call regr1_conserv(v2(i, :, :), press_in_edg, &
+            paprs(i, nbp_lev + 1:1:-1), v3(i, nbp_lev:1:-1, :), &
+            slopes(v2(i, :, :), press_in_edg))
        ! (invert order of indices because "paprs" is in descending order)
     end do