Changeset 1379 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/regr_pr_av_m.F90

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/phylmd/regr_pr_av_m.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/regr_pr_av_m.F90

@@ -21 +21 @@
     ! latitude, pressure, julian day.
     ! We assume that the input fields are already on the "rlatu"
-    ! latitudes, excepth that latitudes are in ascending order in the input
+    ! latitudes, except that latitudes are in ascending order in the input
     ! file.
-    ! We assume that the inputs fields have the same pressure coordinate.
+    ! We assume that all the inputs fields have the same coordinates.
 
     ! The target vertical LMDZ grid is the grid of layer boundaries.
@@ -91 +91 @@
     if (is_mpi_root) then
        do i = 1, n_var
-          call nf95_inq_varid(ncid, name(i), varid)
+          call nf95_inq_varid(ncid, trim(name(i)), varid)
           
           ! Get data at the right day from the input file:
           ncerr = nf90_get_var(ncid, varid, v1(1, :, :, i), &
                start=(/1, 1, julien/))
-          call handle_err("regr_pr_av nf90_get_var " // name(i), ncerr, ncid)
+          call handle_err("regr_pr_av nf90_get_var " // trim(name(i)), ncerr, &
+               ncid)
        end do