/[lmdze]/trunk/phylmd/Mobidic/regr_pr_o3.f

Diff of /trunk/phylmd/Mobidic/regr_pr_o3.f

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/phylmd/Mobidic/regr_pr_o3.f90
revision 8 by guez,
Mon Mar 31 12:51:21 2008 UTC
+trunk/phylmd/Mobidic/regr_pr_o3.f
revision 265 by guez,
Tue Mar 20 09:35:59 2018 UTC
 Line 4 
 module regr_pr_o3_m
  contains
-   subroutine regr_pr_o3(o3_mob_regr)
+   subroutine regr_pr_o3(p3d, o3_mob_regr)
      ! "regr_pr_o3" stands for "regrid pressure ozone".
      ! This procedure reads Mobidic ozone mole fraction from
      ! "coefoz_LMDZ.nc" at the initial day and regrids it in pressure.
+     ! Ozone mole fraction from "coefoz_LMDZ.nc" is a 2D latitude --
+     ! pressure variable.
+     ! The target horizontal LMDZ grid is the "scalar" grid: "rlonv", "rlatu".
+     ! The target vertical LMDZ grid is the grid of layer boundaries.
+     ! We assume that the input variable is already on the LMDZ "rlatu"
+     ! latitude grid.
+     ! The input variable does not depend on longitude, but the
+     ! pressure at LMDZ layers does.
+     ! Therefore, the values on the LMDZ grid do depend on longitude.
      ! Regridding is by averaging, assuming a step function.
-     ! We assume that, in the input file, the pressure levels are in hPa
+     ! We assume that, in the input file, the pressure levels are in
-     ! and strictly increasing.
+     ! hPa and strictly increasing.
-     use conf_gcm_m, only: dayref
+     use dimensions, only: iim, jjm, llm
-     use dimens_m, only: iim, jjm, llm
+     use dynetat0_m, only: day_ref
-     use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, handle_err, &
+     use grid_change, only: dyn_phy
-          nf95_get_coord
      use netcdf, only:  nf90_nowrite, nf90_get_var
-     use regr_pr, only: regr_pr_av
+     use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, handle_err, &
-     use nrutil, only: assert
+          nf95_gw_var
+     use nr_util, only: assert
+     use numer_rec_95, only: regr1_step_av
+     REAL, intent(in):: p3d(:, :, :) ! (iim + 1, jjm + 1, llm+1)
+     ! pressure at layer interfaces, in Pa
+     ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",
+     ! for interface "l")
      real, intent(out):: o3_mob_regr(:, :, :) ! (iim + 1, jjm + 1, llm)
      ! (ozone mole fraction from Mobidic adapted to the LMDZ grid)
-Line 28 
 contains
+Line 43 
 contains
      ! Variables local to the procedure:
-     real, pointer:: plev(:)
+     real, allocatable:: plev(:)
      ! (pressure levels of Mobidic data, in Pa, in strictly increasing order)
      real, allocatable:: press_in_edg(:)
-Line 37 
 contains
+Line 52 
 contains
      integer ncid, varid, ncerr ! for NetCDF
      integer n_plev ! number of pressure levels in Mobidic data
-     integer j
+     integer i, j
      real, allocatable:: r_mob(:, :)! (jjm + 1, n_plev)
-     ! (ozone mole fraction from Mobidic at day "dayref")
+     ! (ozone mole fraction from Mobidic at day "day_ref")
      ! (r_mob(j, k) is at latitude "rlatu(j)" and pressure level "plev(k)".)
      !------------------------------------------------------------
-     call assert(shape(o3_mob_regr) == (/iim + 1, jjm + 1, llm/), "regr_pr_o3")
+     print *, "Call sequence information: regr_pr_o3"
+     call assert(shape(o3_mob_regr) == (/iim + 1, jjm + 1, llm/), &
+          "regr_pr_o3 o3_mob_regr")
+     call assert(shape(p3d) == (/iim + 1, jjm + 1, llm + 1/), &
+          "regr_pr_o3 p3d")
      call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)
-     call nf95_get_coord(ncid, "plev", plev)
+     call nf95_inq_varid(ncid, "plev", varid)
+     call nf95_gw_var(ncid, varid, plev)
      ! Convert from hPa to Pa because "regr_pr_av" requires so:
      plev = plev * 100.
      n_plev = size(plev)
-Line 63 
 contains
+Line 84 
 contains
      ! (infinity, but any value guaranteed to be greater than the
      ! surface pressure would do)
-     deallocate(plev) ! pointer
      call nf95_inq_varid(ncid, "r_Mob", varid)
      allocate(r_mob(jjm + 1, n_plev))
      ! Get data at the right day from the input file:
-     ncerr = nf90_get_var(ncid, varid, r_mob, start=(/1, 1, dayref/))
+     ncerr = nf90_get_var(ncid, varid, r_mob, start=(/1, 1, day_ref/))
      call handle_err("nf90_get_var r_Mob", ncerr)
      ! Latitudes are in increasing order in the input file while
      ! "rlatu" is in decreasing order so we need to invert order:
-Line 77 
 contains
+Line 96 
 contains
      call nf95_close(ncid)
-     o3_mob_regr = regr_pr_av(r_mob, press_in_edg)
+     ! Regrid in pressure by averaging a step function of pressure:
+     do j = 1, jjm + 1
+        do i = 1, iim
+           if (dyn_phy(i, j)) then
+              o3_mob_regr(i, j, llm:1:-1) &
+                   = regr1_step_av(r_mob(j, :), press_in_edg, &
+                   p3d(i, j, llm+1:1:-1))
+              ! (invert order of indices because "p3d" is decreasing)
+           end if
+        end do
+     end do
+     ! Duplicate pole values on all longitudes:
+     o3_mob_regr(2:, 1, :) = spread(o3_mob_regr(1, 1, :), dim=1, ncopies=iim)
+     o3_mob_regr(2:, jjm + 1, :) &
+          = spread(o3_mob_regr(1, jjm + 1, :), dim=1, ncopies=iim)
+     ! Duplicate first longitude to last longitude:
+     o3_mob_regr(iim + 1, 2:jjm, :) = o3_mob_regr(1, 2:jjm, :)
    end subroutine regr_pr_o3

 Legend:



Removed from v.8
 


changed lines


 
Added in v.265
 Legend:



Removed from v.8
 


changed lines


 
Added in v.265
-Removed from v.8
+Added in v.265

	ViewVC Help
Powered by ViewVC 1.1.21