Changeset 1379 for LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/etat0_netcdf.F

Timestamp:

May 6, 2010, 2:19:18 PM (14 years ago)

Author:

lguez

Message:

Added optional ozone tracer with chemistry parameterized by Daniel
Cariolle. This tracer is passive: it has no influence on the rest of
the simulation.

Added variable "zmasse" in file "histrac.nc". Corrected long name of
variable "pplay" in "histrac.nc". Changed name of variable "t" to "T"
in "histrac.nc", corrected long name and unit.

In "phytrac", moved definition of "zmasse" toward the beginning of the
procedure, so that "zmasse" can be given as argument to
"traclmdz". Also added arguments "julien", "gmtime" and "xlon" to
"traclmdz". The four additional arguments are required for the ozone
tracer.

In module "traclmdz_mod", factorized declaration "implicit none" that
was in each procedure. There is now an equivalent single declaration
at the module level.

In procedure "traclmdz", removed variable "delp". Use "zmasse * rg"
instead since we now have "zmasse" as an argument.

Tests. Compilations on Brodie only, with optimization options "-debug"
and "-dev", parallelization options "none", "mpi", "omp" and
"mpi_omp", this revision and revision 1373. Run cases with and without
ozone tracer, 1 and 2 MPI processes, 1 and 2 OpenMP
threads. Comparisons of all cases are ok, except for strange
variations in variables "d_tr_cl_RN" and "d_tr_cl_PB" of file
"histrac.nc", variables "RN" and "PB" of "restart.nc", variables
"trs_RN" and "trs_PB" of "restartphy.nc". Relative variations of these
variables between cases are of order 1e-7 or less, after one day of
simulation.

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/etat0_netcdf.F (modified) (3 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/etat0_netcdf.F

@@ -17 +17 @@
       use regr_lat_time_climoz_m, only: regr_lat_time_climoz
       use conf_phys_m, only: conf_phys
+!     For parameterization of ozone chemistry:
+      use regr_lat_time_coefoz_m, only: regr_lat_time_coefoz
+      use press_coefoz_m, only: press_coefoz
+      use regr_pr_o3_m, only: regr_pr_o3
 #endif
 !#endif of #ifdef CPP_EARTH
@@ -145 +149 @@
       real :: seuil_inversion
       real :: alp_offset
+      logical found
 
       integer  read_climoz ! read ozone climatology
@@ -426 +431 @@
       q3d(:,:,:,1) = qd(:,:,:)
       !
+!     Parameterization of ozone chemistry:
+C     Look for ozone tracer:
+      i = 1
+      DO
+         found = tname(i)=="O3" .OR. tname(i)=="o3"
+         if (found .or. i == nqtot) exit
+         i = i + 1
+      end do
+      if (found) then
+         call regr_lat_time_coefoz
+         call press_coefoz
+         call regr_pr_o3(p3d, q3d(:, :, :, i))
+C     Convert from mole fraction to mass fraction:
+         q3d(:, :, :, i) = q3d(:, :, :, i)  * 48. / 29.
+      end if
 
 !     Ozone climatology: