Changeset 3819 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_av_m.F90

Timestamp:

Apr 22, 2015, 6:28:50 PM (10 years ago)

Author:

ymipsl

Message:

Removed all iim et jjm depedency. Replaced by nbp_lon and nbp_lat.
Supress gr_fi_ecrit, replaced by grid1dTo2d_glo

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_av_m.F90 (modified) (6 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_av_m.F90

@@ -33 +33 @@
     ! NetCDF variable.
 
-    use dimphy, only: klon
+    use dimphy, only: klon, klev
+    use mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
     use netcdf95, only: nf95_inq_varid, handle_err
     use netcdf, only: nf90_get_var
@@ -52 +53 @@
     ! ascending order
 
-    real, intent(in):: paprs(:, :) ! (klon, llm + 1)
+    real, intent(in):: paprs(:, :) ! (klon, klev + 1)
     ! (pression pour chaque inter-couche, en Pa)
 
-    real, intent(out):: v3(:, :, :) ! (klon, llm, size(name))
+    real, intent(out):: v3(:, :, :) ! (klon, klev, size(name))
     ! regridded fields on the partial "physics" grid
     ! "v3(i, k, l)" is at longitude "xlon(i)", latitude
@@ -63 +64 @@
     ! Variables local to the procedure:
 
-    include "dimensions.h"
     integer varid, ncerr ! for NetCDF
 
-    real  v1(iim, jjm + 1, size(press_in_edg) - 1, size(name))
+    real  v1(nbp_lon, nbp_lat, size(press_in_edg) - 1, size(name))
     ! input fields at day "julien", on the global "dynamics" horizontal grid
     ! First dimension is for longitude.
@@ -84 +84 @@
     !--------------------------------------------
 
-    call assert(size(v3, 1) == klon, size(v3, 2) == llm, "regr_pr_av v3 klon")
+    call assert(size(v3, 1) == klon, size(v3, 2) == klev, "regr_pr_av v3 klon")
     n_var = assert_eq(size(name), size(v3, 3), "regr_pr_av v3 n_var")
-    call assert(shape(paprs) == (/klon, llm+1/), "regr_pr_av paprs")
+    call assert(shape(paprs) == (/klon, klev+1/), "regr_pr_av paprs")
 
     !$omp master
@@ -102 +102 @@
        ! Latitudes are in ascending order in the input file while
        ! "rlatu" is in descending order so we need to invert order:
-       v1(1, :, :, :) = v1(1, jjm+1:1:-1, :, :)
+       v1(1, :, :, :) = v1(1, nbp_lat:1:-1, :, :)
 
        ! 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
@@ -113 +113 @@
     ! Regrid in pressure at each horizontal position:
     do i = 1, klon
-       v3(i, llm:1:-1, :) = regr1_step_av(v2(i, :, :), press_in_edg, &
-            paprs(i, llm+1:1:-1))
+       v3(i, klev:1:-1, :) = regr1_step_av(v2(i, :, :), press_in_edg, &
+            paprs(i, klev+1:1:-1))
        ! (invert order of indices because "paprs" is in descending order)
     end do