source: LMDZ5/branches/LF-private/libf/phylmd/regr_pr_o3_m.F90

Last change on this file was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

  • IPSLCM5A configuration (comparison between trunk and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent
  • MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
    • numerical convergence on dynamical fields over 3 days
    • start files are equivalent (except for RN and PB fields)
    • daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

  • IPSLCM5A configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent
  • MH07 configuration:
    • start files are equivalent (except for RN and PB fields)
    • monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files


Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

  • pour une config IPSLCM5A (comparaison tronc/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux (à part sur RN et Pb)
    • fichiers histoire égaux
  • pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
    • convergence numérique sur les champs dynamiques après 3 jours
    • restart et restartphy égaux
    • fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

  • pour une config IPSLCM5A:
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc
  • pour une config MH07
    • restart et restartphy égaux (à part sur RN et PB)
    • pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90

File size: 3.7 KB
Line 
1! $Id$
2module regr_pr_o3_m
3
4  implicit none
5
6contains
7
8  subroutine regr_pr_o3(p3d, o3_mob_regr)
9
10    ! "regr_pr_o3" stands for "regrid pressure ozone".
11    ! This procedure reads Mobidic ozone mole fraction from
12    ! "coefoz_LMDZ.nc" at the initial day of the run and regrids it in
13    ! pressure.
14    ! Ozone mole fraction from "coefoz_LMDZ.nc" at the initial day is
15    ! a 2D latitude -- pressure variable.
16    ! The target horizontal LMDZ grid is the "scalar" grid: "rlonv", "rlatu".
17    ! The target vertical LMDZ grid is the grid of layer boundaries.
18    ! We assume that the input variable is already on the LMDZ "rlatu"
19    ! latitude grid.
20    ! The input variable does not depend on longitude, but the
21    ! pressure at LMDZ layers does.
22    ! Therefore, the values on the LMDZ grid do depend on longitude.
23    ! Regridding is by averaging, assuming a step function.
24    ! We assume that, in the input file, the pressure levels are in
25    ! hPa and strictly increasing.
26
27    use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, handle_err
28    use netcdf, only:  nf90_nowrite, nf90_get_var
29    use assert_m, only: assert
30    use regr1_step_av_m, only: regr1_step_av
31    use press_coefoz_m, only: press_in_edg
32    use control_mod, only: dayref
33
34    REAL, intent(in):: p3d(:, :, :) ! pressure at layer interfaces, in Pa
35    ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",
36    ! for interface "l")
37
38    real, intent(out):: o3_mob_regr(:, :, :) ! (iim + 1, jjm + 1, llm)
39    ! (ozone mole fraction from Mobidic adapted to the LMDZ grid)
40    ! ("o3_mob_regr(i, j, l)" is at longitude "rlonv(i)", latitude
41    ! "rlatu(j)" and pressure level "pls(i, j, l)")
42
43    ! Variables local to the procedure:
44
45    include "dimensions.h"
46
47    integer ncid, varid, ncerr ! for NetCDF
48    integer i, j
49
50    real r_mob(jjm + 1, size(press_in_edg) - 1)
51    ! (ozone mole fraction from Mobidic at day "dayref")
52    ! (r_mob(j, k) is at latitude "rlatu(j)", in pressure interval
53    ! "[press_in_edg(k), press_in_edg(k+1)]".)
54
55    !------------------------------------------------------------
56
57    print *, "Call sequence information: regr_pr_o3"
58    call assert(shape(o3_mob_regr) == (/iim + 1, jjm + 1, llm/), &
59         "regr_pr_o3 o3_mob_regr")
60    call assert(shape(p3d) == (/iim + 1, jjm + 1, llm + 1/), "regr_pr_o3 p3d")
61
62    call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)
63
64    call nf95_inq_varid(ncid, "r_Mob", varid)
65    ! Get data at the right day from the input file:
66    ncerr = nf90_get_var(ncid, varid, r_mob, start=(/1, 1, dayref/))
67    call handle_err("nf90_get_var r_Mob", ncerr)
68    ! Latitudes are in ascending order in the input file while
69    ! "rlatu" is in descending order so we need to invert order:
70    r_mob = r_mob(jjm+1:1:-1, :)
71
72    call nf95_close(ncid)
73
74    ! Regrid in pressure by averaging a step function of pressure:
75
76    ! Poles:
77    do j = 1, jjm + 1, jjm
78       o3_mob_regr(1, j, llm:1:-1) &
79            = regr1_step_av(r_mob(j, :), press_in_edg, p3d(1, j, llm+1:1:-1))
80       ! (invert order of indices because "p3d" is in descending order)
81    end do
82
83    ! Other latitudes:
84    do j = 2, jjm
85       do i = 1, iim
86          o3_mob_regr(i, j, llm:1:-1) &
87               = regr1_step_av(r_mob(j, :), press_in_edg, &
88               p3d(i, j, llm+1:1:-1))
89             ! (invert order of indices because "p3d" is in descending order)
90       end do
91    end do
92
93    ! Duplicate pole values on all longitudes:
94    o3_mob_regr(2:, 1, :) = spread(o3_mob_regr(1, 1, :), dim=1, ncopies=iim)
95    o3_mob_regr(2:, jjm + 1, :) &
96         = spread(o3_mob_regr(1, jjm + 1, :), dim=1, ncopies=iim)
97
98    ! Duplicate first longitude to last longitude:
99    o3_mob_regr(iim + 1, 2:jjm, :) = o3_mob_regr(1, 2:jjm, :)
100
101  end subroutine regr_pr_o3
102
103end module regr_pr_o3_m
Note: See TracBrowser for help on using the repository browser.