/[lmdze]/trunk/libf/phylmd/Mobidic/regr_pr_comb_coefoz.f90

Diff of /trunk/libf/phylmd/Mobidic/regr_pr_comb_coefoz.f90

Parent Directory | Revision Log | View Patch Patch

-trunk/libf/phylmd/read_coefoz_m.f90
revision 3 by guez,
Wed Feb 27 13:16:39 2008 UTC
+trunk/libf/phylmd/Mobidic/regr_pr_comb_coefoz.f90
revision 17 by guez,
Tue Aug  5 13:31:32 2008 UTC
 Line 1
- module read_coefoz_m
+ module regr_pr_comb_coefoz_m
    ! This module is clean: no C preprocessor directive, no include line.
 Line 7 
 module read_coefoz_m
    implicit none
-   real, save:: c_Mob(klon, llm, 12)
+   ! The five module variables declared here are on the "physics" grid.
+   ! The value of each variable for index "(i, k)" is at longitude
+   ! "rlon(i)", latitude "rlat(i)" and middle of layer "k".
+   real, save:: c_Mob(klon, llm)
    ! (sum of Mobidic terms in the net mass production rate of ozone
    ! by chemistry, per unit mass of air, in s-1)
-   ! (On the "physics" grid.
-   ! Third dimension is the number of the month in the year.
-   ! "c_Mob(i, k, month)" is at longitude "rlon(i)", latitude "rlat(i)",
-   ! middle of layer "k".)
-   real, save:: a2(klon, llm, 12)
+   real, save:: a2(klon, llm)
    ! (derivative of mass production rate of ozone per unit mass of
    ! air with respect to ozone mass fraction, in s-1)
-   ! (On the "physics" grid.
-   ! Third dimension is the number of the month in the year.
-   ! "a2(i, k, month)" is at longitude "rlon(i)", latitude "rlat(i)",
-   ! middle of layer "k".)
-   real, save:: a4_mass(klon, llm, 12)
+   real, save:: a4_mass(klon, llm)
    ! (derivative of mass production rate of ozone per unit mass of
    ! air with respect to temperature, in s-1 K-1)
-   ! (On the "physics" grid.
-   ! Third dimension is the number of the month in the year.
-   ! "a4_mass(i, k, month)" is at longitude "rlon(i)", latitude "rlat(i)",
-   ! middle of layer "k".)
-   real, save:: a6_mass(klon, llm, 12)
+   real, save:: a6_mass(klon, llm)
    ! (derivative of mass production rate of ozone per unit mass of
    ! air with respect to mass column-density of ozone above, in m2 s-1 kg-1)
-   ! (On the "physics" grid.
-   ! Third dimension is the number of the month in the year.
-   ! "a6_mass(i, k, month)" is at longitude "rlon(i)", latitude "rlat(i)",
-   ! middle of layer "k".)
-   real, save:: r_het_interm(klon, llm, 12)
+   real, save:: r_het_interm(klon, llm)
    ! (net mass production rate by heterogeneous chemistry, per unit
    ! mass of ozone, corrected for chlorine content and latitude, but
    ! not for temperature and sun direction, in s-1)
-   ! (On the "physics" grid.
-   ! Third dimension is the number of the month in the year.
-   ! "r_het_interm(i, k, month)" is at longitude "rlon(i)", latitude "rlat(i)",
-   ! middle of layer "k".)
    private klon, llm
  contains
-   subroutine read_coefoz
+   subroutine regr_pr_comb_coefoz(julien)
-     ! This subroutine reads from a file all eight parameters for ozone
+     ! "regr_pr_comb_coefoz" stands for "regrid pressure combine
-     ! chemistry, at all months.
+     ! coefficients ozone".
-     ! The parameters are packed to the "physics" grid.
-     ! Finally, the eight parameters are combined to define the five
-     ! module variables.
-     use netcdf95, only: nf95_open, nf95_close, nf90_nowrite
+     ! This subroutine :
-     use o3_Mob_ph_m, only: o3_Mob_ph
+     ! -- reads from a file all eight coefficients for ozone chemistry,
+     !    at the current day ;
+     ! -- regrids the coefficients in pressure to the LMDZ vertical grid ;
+     ! -- packs the coefficients to the "physics" horizontal grid ;
+     ! -- combines the eight coefficients to define the five module variables.
+     ! We assume that, in "coefoz_LMDZ.nc", the pressure levels are in hPa
+     ! and strictly increasing.
+     use netcdf95, only: nf95_open, nf95_close, nf95_get_coord
+     use netcdf, only: nf90_nowrite
+     use regr_pr_coefoz, only: regr_pr_av_coefoz, regr_pr_int_coefoz
      use phyetat0_m, only: rlat
+     integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
      ! Variables local to the procedure:
      integer ncid ! for NetCDF
-     real a6(klon, llm, 12)
+     real coefoz(klon, llm)
-     ! (derivative of P_net_Mob with respect to column-density of ozone
+     ! (temporary storage for an ozone coefficient)
+     ! (On the "physics" grid.
+     ! "coefoz(i, k)" is at longitude "rlon(i)", latitude "rlat(i)",
+     ! middle of layer "k".)
+     real a6(klon, llm)
+     ! (derivative of "P_net_Mob" with respect to column-density of ozone
      ! above, in cm2 s-1)
      ! (On the "physics" grid.
-     ! Third dimension is the number of the month in the year.
+     ! "a6(i, k)" is at longitude "rlon(i)", latitude "rlat(i)",
-     ! "a6(i, k, month)" is at longitude "rlon(i)", latitude "rlat(i)",
      ! middle of layer "k".)
      real, parameter:: amu = 1.6605402e-27 ! atomic mass unit, in kg
-Line 80 
 contains
+Line 79 
 contains
      real, parameter:: Clx = 3.8e-9
      ! (total chlorine content in the upper stratosphere)
-     integer k, month
+     integer k
      !------------------------------------
-     print *, "Call sequence information: read_coefoz"
+     print *, "Call sequence information: regr_pr_comb_coefoz"
      call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)
-     a2 = o3_Mob_ph(ncid, "a2")
+     call regr_pr_av_coefoz(ncid, "a2", julien, a2)
-     a4_mass = o3_Mob_ph(ncid, "a4") * 48. / 29.
-     a6 = o3_Mob_ph(ncid, "a6")
+     call regr_pr_av_coefoz(ncid, "a4", julien, a4_mass)
+     a4_mass = a4_mass * 48. / 29.
+     call regr_pr_av_coefoz(ncid, "a6", julien, a6)
      ! Compute "a6_mass" avoiding underflow, do not divide by 1e4
      ! before dividing by molecular mass:
      a6_mass = a6 / (1e4 * 29. * amu)
      ! (factor 1e4: conversion from cm2 to m2)
-     c_Mob = 48. / 29. * (o3_Mob_ph(ncid, "P_net_Mob") &
+     ! Combine coefficients to get "c_Mob":
-          - a2 * o3_Mob_ph(ncid, "r_Mob") - a6 * o3_Mob_ph(ncid, "Sigma_Mob")) &
+     ! (We use as few local variables as possible, in order to spare
-          - a4_mass * o3_Mob_ph(ncid, "temp_Mob")
+     ! main memory.)
+     call regr_pr_av_coefoz(ncid, "P_net_Mob", julien, c_Mob)
+     call regr_pr_av_coefoz(ncid, "r_Mob", julien, coefoz)
+     c_mob = c_mob - a2 * coefoz
+     call regr_pr_int_coefoz(ncid, "Sigma_Mob", julien, top_value=0., v3=coefoz)
+     c_mob = (c_mob - a6 * coefoz) * 48. / 29.
+     call regr_pr_av_coefoz(ncid, "temp_Mob", julien, coefoz)
+     c_mob = c_mob - a4_mass * coefoz
-     r_het_interm = o3_Mob_ph(ncid, "R_Het") * (Clx / 3.8e-9)**2
+     call regr_pr_av_coefoz(ncid, "R_Het", julien, r_het_interm)
      ! Heterogeneous chemistry is only at high latitudes:
-     forall (k = 1: llm, month = 1: 12)
+     forall (k = 1: llm)
-        where (abs(rlat) <= 45.) r_het_interm(:, k, month) = 0.
+        where (abs(rlat) <= 45.) r_het_interm(:, k) = 0.
      end forall
+     r_het_interm = r_het_interm * (Clx / 3.8e-9)**2
      call nf95_close(ncid)
-   end subroutine read_coefoz
+   end subroutine regr_pr_comb_coefoz
- end module read_coefoz_m
+ end module regr_pr_comb_coefoz_m

 Legend:



Removed from v.3
 


changed lines


 
Added in v.17
 Legend:



Removed from v.3
 


changed lines


 
Added in v.17
-Removed from v.3
+Added in v.17

	ViewVC Help
Powered by ViewVC 1.1.21