Changeset 2408 for LMDZ5/branches/testing/libf/phylmd/regr_pr_o3_m.F90

Timestamp:

Dec 14, 2015, 11:43:09 AM (9 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2298:2396 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/regr_pr_o3_m.F90 (modified) (6 diffs)

LMDZ5/branches/testing/libf/phylmd/regr_pr_o3_m.F90

@@ -30 +30 @@
     use regr1_step_av_m, only: regr1_step_av
     use press_coefoz_m, only: press_in_edg
-    use control_mod, only: dayref
+    use time_phylmdz_mod, only: day_ref
+    use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, nbp_lev
 
     REAL, intent(in):: p3d(:, :, :) ! pressure at layer interfaces, in Pa
@@ -43 +44 @@
     ! Variables local to the procedure:
 
-    include "dimensions.h"
-
     integer ncid, varid, ncerr ! for NetCDF
     integer i, j
 
-    real r_mob(jjm + 1, size(press_in_edg) - 1)
-    ! (ozone mole fraction from Mobidic at day "dayref")
+    real r_mob(nbp_lat, size(press_in_edg) - 1)
+    ! (ozone mole fraction from Mobidic at day "day_ref")
     ! (r_mob(j, k) is at latitude "rlatu(j)", in pressure interval
     ! "[press_in_edg(k), press_in_edg(k+1)]".)
@@ -56 +55 @@
 
     print *, "Call sequence information: regr_pr_o3"
-    call assert(shape(o3_mob_regr) == (/iim + 1, jjm + 1, llm/), &
+    call assert(shape(o3_mob_regr) == (/nbp_lon + 1, nbp_lat, nbp_lev/), &
          "regr_pr_o3 o3_mob_regr")
-    call assert(shape(p3d) == (/iim + 1, jjm + 1, llm + 1/), "regr_pr_o3 p3d")
+    call assert(shape(p3d) == (/nbp_lon + 1, nbp_lat, nbp_lev + 1/), "regr_pr_o3 p3d")
 
     call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)
@@ -64 +63 @@
     call nf95_inq_varid(ncid, "r_Mob", varid)
     ! 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 ascending order in the input file while
     ! "rlatu" is in descending order so we need to invert order:
-    r_mob = r_mob(jjm+1:1:-1, :)
+    r_mob = r_mob(nbp_lat:1:-1, :)
 
     call nf95_close(ncid)
@@ -75 +74 @@
 
     ! Poles:
-    do j = 1, jjm + 1, jjm
-       o3_mob_regr(1, j, llm:1:-1) &
-            = regr1_step_av(r_mob(j, :), press_in_edg, p3d(1, j, llm+1:1:-1))
+    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))
        ! (invert order of indices because "p3d" is in descending order)
     end do
 
     ! Other latitudes:
-    do j = 2, jjm
-       do i = 1, iim
-          o3_mob_regr(i, j, llm:1:-1) &
+    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, llm+1:1:-1))
+               p3d(i, j, nbp_lev+1:1:-1))
              ! (invert order of indices because "p3d" is in descending order)
        end do
@@ -92 +91 @@
 
     ! 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)
+    o3_mob_regr(2:, 1, :) = spread(o3_mob_regr(1, 1, :), dim=1, ncopies=nbp_lon)
+    o3_mob_regr(2:, nbp_lat, :) &
+         = spread(o3_mob_regr(1, nbp_lat, :), dim=1, ncopies=nbp_lon)
 
     ! Duplicate first longitude to last longitude:
-    o3_mob_regr(iim + 1, 2:jjm, :) = o3_mob_regr(1, 2:jjm, :)
+    o3_mob_regr(nbp_lon + 1, 2:nbp_lat-1, :) = o3_mob_regr(1, 2:nbp_lat-1, :)
 
   end subroutine regr_pr_o3