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

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.

File:

• LMDZ5/trunk/libf/phylmd/regr_pr_o3_m.F90 (modified) (3 diffs)

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