Changeset 2440 for LMDZ5/trunk/libf/phylmd/regr_lat_time_climoz_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_lat_time_climoz_m.F90 (modified) (9 diffs)

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)