Changeset 3819 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_comb_coefoz_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_comb_coefoz_m.F90 (modified) (7 diffs)

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

@@ -39 +39 @@
     ! It allocates module variables.
 
-    use dimphy, only: klon
+    use dimphy, only: klon, klev
 
     ! Variables local to the procedure:
-    include "dimensions.h"
 
     !---------------------------------------
@@ -49 +48 @@
     print *, "Call sequence information: alloc_coefoz"
     !$omp end master
-    allocate(c_Mob(klon, llm), a2(klon, llm), a4_mass(klon, llm))
-    allocate(a6_mass(klon, llm), r_het_interm(klon, llm))
+    allocate(c_Mob(klon, klev), a2(klon, klev), a4_mass(klon, klev))
+    allocate(a6_mass(klon, klev), r_het_interm(klon, klev))
 
   end subroutine alloc_coefoz
@@ -74 +73 @@
     use netcdf, only: nf90_nowrite
     use assert_m, only: assert
-    use dimphy, only: klon
+    use dimphy, only: klon, klev
     use mod_phys_lmdz_mpi_data, only: is_mpi_root
     use regr_pr_av_m, only: regr_pr_av
@@ -85 +84 @@
     ! (latitude on the partial "physics" grid, in degrees)
 
-    real, intent(in):: paprs(:, :) ! (klon, llm + 1)
+    real, intent(in):: paprs(:, :) ! (klon, klev + 1)
     ! (pression pour chaque inter-couche, en Pa)
 
-    real, intent(in):: pplay(:, :) ! (klon, llm)
+    real, intent(in):: pplay(:, :) ! (klon, klev)
     ! (pression pour le mileu de chaque couche, en Pa)
 
     ! Variables local to the procedure:
 
-    include "dimensions.h"
     integer ncid ! for NetCDF
 
-    real coefoz(klon, llm, 7)
+    real coefoz(klon, klev, 7)
     ! (temporary storage for 7 ozone coefficients)
     ! (On the partial "physics" grid.
@@ -102 +100 @@
     ! middle of layer "k".)
 
-    real a6(klon, llm)
+    real a6(klon, klev)
     ! (derivative of "P_net_Mob" with respect to column-density of ozone
     ! above, in cm2 s-1)
@@ -121 +119 @@
     call assert((/size(rlat), size(paprs, 1), size(pplay, 1)/) == klon, &
          "regr_pr_comb_coefoz klon")
-    call assert((/size(paprs, 2) - 1, size(pplay, 2)/) == llm, &
-         "regr_pr_comb_coefoz llm")
+    call assert((/size(paprs, 2) - 1, size(pplay, 2)/) == klev, &
+         "regr_pr_comb_coefoz klev")
 
     !$omp master
@@ -150 +148 @@
     r_het_interm = coefoz(:, :, 7)
     ! Heterogeneous chemistry is only at high latitudes:
-    forall (k = 1: llm)
+    forall (k = 1: klev)
        where (abs(rlat) <= 45.) r_het_interm(:, k) = 0.
     end forall