Changeset 2440

Timestamp:

Feb 9, 2016, 3:45:31 PM (9 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.

Location:

LMDZ5/trunk/libf

Files:

• misc/regr1_conserv_m.F90 (added)
• misc/slopes_m.F90 (added)
• obsolete/LIST.txt (modified) (1 diff)
• obsolete/regr1_step_av_m.F90 (moved) (moved from LMDZ5/trunk/libf/misc/regr1_step_av_m.F90)
• phylmd/regr_lat_time_climoz_m.F90 (modified) (9 diffs)
• phylmd/regr_lat_time_coefoz_m.F90 (modified) (3 diffs)
• phylmd/regr_pr_av_m.F90 (modified) (4 diffs)
• phylmd/regr_pr_o3_m.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/obsolete/LIST.txt

@@ -6 +6 @@
 phylmd/mkstat.F90             2320
 phylmd/inistats.F90           2320
+misc/regr1_step_av_m.F90      2439

LMDZ5/trunk/libf/phylmd/regr_lat_time_climoz_m.F90

@@ -67 +67 @@
 
     use mod_grid_phy_lmdz, ONLY : nbp_lat
-    use regr1_step_av_m, only: regr1_step_av
+    use regr1_conserv_m, only: regr1_conserv
     use regr3_lint_m, only: regr3_lint
     use netcdf95, only: handle_err, nf95_close, nf95_get_att, nf95_gw_var, &
@@ -76 +76 @@
     use regular_lonlat_mod, only : boundslat_reg, south
     use nrtype, only: pi
+    use slopes_m, only: slopes
 
     integer, intent(in):: read_climoz ! read ozone climatology
@@ -92 +93 @@
     ! (of input data, converted to rad, sorted in strictly ascending order)
 
-    real, allocatable:: lat_in_edg(:)
-    ! (edges of latitude intervals for input data, in rad, in strictly
+    real, allocatable:: sin_lat_in_edg(:)
+    ! (sine of edges of latitude intervals for input data, in rad, in strictly
     ! ascending order)
 
@@ -115 +116 @@
 
     real, allocatable:: o3_regr_lat(:, :, :, :)
-    ! (jjm + 1, n_plev, 0:13, read_climoz)
+    ! (nbp_lat, n_plev, 0:13, read_climoz)
     ! mean of "o3_in" over a latitude interval of LMDZ
     ! First dimension is latitude interval.
     ! The latitude interval for "o3_regr_lat(j,:, :, :)" contains "rlatu(j)".
-    ! If "j" is between 2 and "jjm" then the interval is:
+    ! If "j" is between 2 and "nbp_lat - 1" then the interval is:
     ! [rlatv(j), rlatv(j-1)]
-    ! If "j" is 1 or "jjm + 1" then the interval is:
+    ! If "j" is 1 or "nbp_lat" then the interval is:
     ! [rlatv(1), pi / 2]
     ! or:
-    ! [- pi / 2, rlatv(jjm)]
+    ! [- pi / 2, rlatv(nbp_lat - 1)]
     ! respectively.
     ! "o3_regr_lat(:, k, :, :)" is for pressure level "plev(k)".
@@ -132 +133 @@
 
     real, allocatable:: o3_out(:, :, :, :)
-    ! (jjm + 1, n_plev, 360, read_climoz)
+    ! (nbp_lat, n_plev, 360, read_climoz)
     ! regridded ozone climatology
     ! "o3_out(j, k, l, :)" is at latitude "rlatu(j)", pressure
@@ -175 +176 @@
     latitude = latitude / 180. * pi
     n_lat = size(latitude)
-    ! We need to supply the latitudes to "regr1_step_av" in
+    ! We need to supply the latitudes to "regr1_conserv" in
     ! ascending order, so invert order if necessary:
     desc_lat = latitude(1) > latitude(n_lat)
@@ -181 +182 @@
 
     ! Compute edges of latitude intervals:
-    allocate(lat_in_edg(n_lat + 1))
-    lat_in_edg(1) = - pi / 2
-    forall (j = 2:n_lat) lat_in_edg(j) = (latitude(j - 1) + latitude(j)) / 2
-    lat_in_edg(n_lat + 1) = pi / 2
+    allocate(sin_lat_in_edg(n_lat + 1))
+    sin_lat_in_edg(1) = - 1.
+    forall (j = 2:n_lat) sin_lat_in_edg(j) = sin((latitude(j - 1) &
+         + latitude(j)) / 2.)
+    sin_lat_in_edg(n_lat + 1) = 1.
     deallocate(latitude) ! pointer
 
@@ -292 +294 @@
        print *, &
             "Found 12 months in ozone climatologies, assuming periodicity..."
-       o3_regr_lat(nbp_lat:1:-1, :, 1:12, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, boundslat_reg(nbp_lat-1:1:-1,south), pi / 2/)))
+       call regr1_conserv(o3_in, xs = sin_lat_in_edg, &
+            xt = (/- 1., sin(boundslat_reg(nbp_lat - 1:1:- 1, south)), 1./), &
+            vt = o3_regr_lat(nbp_lat:1:- 1, :, 1:12, :), &
+            slope = slopes(o3_in, sin_lat_in_edg))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)
@@ -303 +307 @@
     else
        print *, "Using 14 months in ozone climatologies..."
-       o3_regr_lat(nbp_lat:1:-1, :, :, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, boundslat_reg(nbp_lat-1:1:-1,south), pi / 2/)))
+       call regr1_conserv(o3_in, xs = sin_lat_in_edg, &
+            xt = (/- 1., sin(boundslat_reg(nbp_lat - 1:1:- 1, south)), 1./), &
+            vt = o3_regr_lat(nbp_lat:1:- 1, :, :, :), &
+            slope = slopes(o3_in, sin_lat_in_edg))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)

LMDZ5/trunk/libf/phylmd/regr_lat_time_coefoz_m.F90

@@ -41 +41 @@
 
     use mod_grid_phy_lmdz, ONLY : nbp_lat
-    use regr1_step_av_m, only: regr1_step_av
+    use regr1_conserv_m, only: regr1_conserv
     use regr3_lint_m, only: regr3_lint
     use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, handle_err, &
@@ -162 +162 @@
     latitude = latitude / 180. * pi
     n_lat = size(latitude)
-    ! We need to supply the latitudes to "regr1_step_av" in
+    ! We need to supply the latitudes to "regr1_conserv" in
     ! ascending order, so invert order if necessary:
     desc_lat = latitude(1) > latitude(n_lat)
@@ -209 +209 @@
        ! We average with respect to sine of latitude, which is
        ! equivalent to weighting by cosine of latitude:
-       v_regr_lat(nbp_lat:1:-1, :, 1:12) = regr1_step_av(o3_par_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, boundslat_reg(nbp_lat-1:1:-1,south), pi / 2/)))
+       call regr1_conserv(o3_par_in, xs = sin(lat_in_edg), &
+            xt = (/-1., sin((/boundslat_reg(nbp_lat-1:1:-1,south)/)), 1./), &
+            vt = v_regr_lat(nbp_lat:1:-1, :, 1:12))
        ! (invert order of indices in "v_regr_lat" because "rlatu" is
        ! in descending order)

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

LMDZ5/trunk/libf/phylmd/regr_pr_o3_m.F90

@@ -28 +28 @@
     use netcdf, only:  nf90_nowrite, nf90_get_var
     use assert_m, only: assert
-    use regr1_step_av_m, only: regr1_step_av
+    use regr1_conserv_m, only: regr1_conserv
     use press_coefoz_m, only: press_in_edg
     use time_phylmdz_mod, only: day_ref
@@ -75 +75 @@
     ! Poles:
     do j = 1, nbp_lat, nbp_lat-1
-       o3_mob_regr(1, j, nbp_lev:1:-1) &
-            = regr1_step_av(r_mob(j, :), press_in_edg, p3d(1, j, nbp_lev+1:1:-1))
+       call regr1_conserv(r_mob(j, :), press_in_edg, &
+            p3d(1, j, nbp_lev + 1:1:-1), o3_mob_regr(1, j, nbp_lev:1:-1))
        ! (invert order of indices because "p3d" is in descending order)
     end do
@@ -83 +83 @@
     do j = 2, nbp_lat-1
        do i = 1, nbp_lon
-          o3_mob_regr(i, j, nbp_lev:1:-1) &
-               = regr1_step_av(r_mob(j, :), press_in_edg, &
-               p3d(i, j, nbp_lev+1:1:-1))
-             ! (invert order of indices because "p3d" is in descending order)
+          call regr1_conserv(r_mob(j, :), press_in_edg, &
+               p3d(i, j, nbp_lev + 1:1:-1), o3_mob_regr(i, j, nbp_lev:1:-1))
+          ! (invert order of indices because "p3d" is in descending order)
        end do
     end do