Changeset 2346 for LMDZ5/trunk/libf/phylmd/regr_pr_int_m.F90

Timestamp:

Aug 21, 2015, 5:13:46 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Physics/dynamics separation:

• remove all references to dimensions.h from physics. nbp_lon (==iim) , nbp_lat (==jjm+1) and nbp_lev (==llm) from mod_grid_phy_lmdz should be used instead.
• added module regular_lonlat_mod in phy_common to store information about the global (lon-lat) grid cell boundaries and centers.

EM

File:

• LMDZ5/trunk/libf/phylmd/regr_pr_int_m.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/regr_pr_int_m.F90

@@ -30 +30 @@
     use regr1_lint_m, only: regr1_lint
     use mod_phys_lmdz_mpi_data, only: is_mpi_root
-
+    use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat, nbp_lev
     use mod_phys_lmdz_transfert_para, only: scatter2d
     ! (pack to the LMDZ horizontal "physics" grid and scatter)
@@ -41 +41 @@
     ! (pressure level of input data, in Pa, in strictly ascending order)
 
-    real, intent(in):: pplay(:, :) ! (klon, llm)
+    real, intent(in):: pplay(:, :) ! (klon, nbp_lev)
     ! (pression pour le mileu de chaque couche, en Pa)
 
@@ -47 +47 @@
     ! (extra value of field at 0 pressure)
 
-    real, intent(out):: v3(:, :) ! (klon, llm)
+    real, intent(out):: v3(:, :) ! (klon, nbp_lev)
     ! (regridded field on the partial "physics" grid)
     ! ("v3(i, k)" is at longitude "xlon(i)", latitude
@@ -54 +54 @@
     ! Variables local to the procedure:
 
-    include "dimensions.h"
     integer varid, ncerr ! for NetCDF
 
-    real  v1(iim, jjm + 1, 0:size(plev))
+    real  v1(nbp_lon, nbp_lat, 0:size(plev))
     ! (input field at day "julien", on the global "dynamics" horizontal grid)
     ! (First dimension is for longitude.
@@ -72 +71 @@
     !--------------------------------------------
 
-    call assert(shape(v3) == (/klon, llm/), "regr_pr_int v3")
-    call assert(shape(pplay) == (/klon, llm/), "regr_pr_int pplay")
+    call assert(shape(v3) == (/klon, nbp_lev/), "regr_pr_int v3")
+    call assert(shape(pplay) == (/klon, nbp_lev/), "regr_pr_int pplay")
 
     !$omp master
@@ -84 +83 @@
        ! Latitudes are in ascending order in the input file while
        ! "rlatu" is in descending order so we need to invert order:
-       v1(1, :, 1:) = v1(1, jjm+1:1:-1, 1:)
+       v1(1, :, 1:) = v1(1, nbp_lat:1:-1, 1:)
 
        ! Complete "v1" with the value at 0 pressure:
@@ -90 +89 @@
 
        ! Duplicate on all longitudes:
-       v1(2:, :, :) = spread(v1(1, :, :), dim=1, ncopies=iim-1)
+       v1(2:, :, :) = spread(v1(1, :, :), dim=1, ncopies=nbp_lon-1)
     end if
     !$omp end master
@@ -98 +97 @@
     ! Regrid in pressure at each horizontal position:
     do i = 1, klon
-       v3(i, llm:1:-1) = regr1_lint(v2(i, :), (/0., plev/), pplay(i, llm:1:-1))
+       v3(i, nbp_lev:1:-1) = regr1_lint(v2(i, :), (/0., plev/), pplay(i, nbp_lev:1:-1))
        ! (invert order of indices because "pplay" is in descending order)
     end do