Changeset 2408 for LMDZ5/branches/testing/libf/phylmd/o3_chem_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/o3_chem_m.F90 (modified) (11 diffs)

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

@@ -14 +14 @@
 
     ! All the 2-dimensional arrays are on the partial "physics" grid.
-    ! Their shape is "(/klon, llm/)".
+    ! Their shape is "(/klon, nbp_lev/)".
     ! Index "(i, :)" is for longitude "rlon(i)", latitude "rlat(i)".
 
@@ -20 +20 @@
     use dimphy, only: klon
     use regr_pr_comb_coefoz_m, only: c_Mob, a4_mass, a2, r_het_interm
+    use mod_grid_phy_lmdz, only: nbp_lev
+    use nrtype, only: pi
 
     integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
     real, intent(in):: gmtime ! heure de la journée en fraction de jour
-    real, intent(in):: t_seri(:, :) ! (klon, llm) temperature, in K
+    real, intent(in):: t_seri(:, :) ! (klon, nbp_lev) temperature, in K
 
-    real, intent(in):: zmasse(:, :) ! (klon, llm)
+    real, intent(in):: zmasse(:, :) ! (klon, nbp_lev)
     ! (column-density of mass of air in a cell, in kg m-2)
     ! "zmasse(:, k)" is for layer "k".)
@@ -34 +36 @@
     ! (longitude and latitude of each horizontal position, in degrees)
 
-    real, intent(inout):: q(:, :) ! (klon, llm) mass fraction of ozone
+    real, intent(inout):: q(:, :) ! (klon, nbp_lev) mass fraction of ozone
     ! "q(:, k)" is at middle of layer "k".)
 
     ! Variables local to the procedure:
-    include "dimensions.h"
-    include "comconst.h"
     ! (for "pi")
     integer k
 
-    real c(klon, llm)
+    real c(klon, nbp_lev)
     ! (constant term during a time step in the net mass production
     ! rate of ozone by chemistry, per unit mass of air, in s-1)
     ! "c(:, k)" is at middle of layer "k".)
 
-    real b(klon, llm)
+    real b(klon, nbp_lev)
     ! (coefficient of "q" in the net mass production
     ! rate of ozone by chemistry, per unit mass of air, in s-1)
     ! "b(:, k)" is at middle of layer "k".)
 
-    real dq_o3_chem(klon, llm)
+    real dq_o3_chem(klon, nbp_lev)
     ! (variation of ozone mass fraction due to chemistry during a time step)
     ! "dq_o3_chem(:, k)" is at middle of layer "k".)
@@ -69 +69 @@
     call assert(klon == (/size(q, 1), size(t_seri, 1), size(zmasse, 1), &
          size(rlat), size(rlon)/), "o3_chem klon")
-    call assert(llm == (/size(q, 2), size(t_seri, 2), size(zmasse, 2)/), &
-         "o3_chem llm")
+    call assert(nbp_lev == (/size(q, 2), size(t_seri, 2), size(zmasse, 2)/), &
+         "o3_chem nbp_lev")
 
     c = c_Mob + a4_mass * t_seri
@@ -86 +86 @@
     call orbite(real(julien), earth_long, trash1)
     call zenang(earth_long, gmtime, pdtphys, rlat, rlon, pmu0, trash2)
-    forall (k = 1: llm)
+    forall (k = 1: nbp_lev)
        where (pmu0 <= cos(87. / 180. * pi)) b(:, k) = 0.
     end forall
@@ -112 +112 @@
 
     ! All the 2-dimensional arrays are on the partial "physics" grid.
-    ! Their shape is "(/klon, llm/)".
+    ! Their shape is "(/klon, nbp_lev/)".
     ! Index "(i, :)" is for longitude "rlon(i)", latitude "rlat(i)".
 
@@ -118 +118 @@
     use assert_m, only: assert
     use dimphy, only: klon
+    use mod_grid_phy_lmdz, only: nbp_lev
 
     real, intent(in):: q(:, :) ! mass fraction of ozone
@@ -136 +137 @@
     ! ("b(:, k)" is at middle of layer "k".)
 
-    include "dimensions.h"
-
-    real o3_prod(klon, llm)
+    real o3_prod(klon, nbp_lev)
     ! (net mass production rate of ozone by chemistry, per unit mass
     ! of air, in s-1)
@@ -145 +144 @@
     ! Variables local to the procedure:
 
-    real sigma_mass(klon, llm)
+    real sigma_mass(klon, nbp_lev)
     ! (mass column-density of ozone above point, in kg m-2)
     ! ("sigma_mass(:, k)" is at middle of layer "k".)
@@ -155 +154 @@
     call assert(klon == (/size(q, 1), size(zmasse, 1), size(c, 1), &
          size(b, 1)/), "o3_prod 1")
-    call assert(llm == (/size(q, 2), size(zmasse, 2), size(c, 2), &
+    call assert(nbp_lev == (/size(q, 2), size(zmasse, 2), size(c, 2), &
          size(b, 2)/), "o3_prod 2")
 
@@ -161 +160 @@
     ! "k", and, as a first approximation, take it as column-density
     ! above the middle of layer "k":
-    sigma_mass(:, llm) = zmasse(:, llm) * q(:, llm) ! top layer
-    do k =  llm - 1, 1, -1
+    sigma_mass(:, nbp_lev) = zmasse(:, nbp_lev) * q(:, nbp_lev) ! top layer
+    do k =  nbp_lev - 1, 1, -1
        sigma_mass(:, k) = sigma_mass(:, k+1) + zmasse(:, k) * q(:, k)
     end do