source: LMDZ4/trunk/libf/phylmd/o3cm.F @ 2852

Last change on this file since 2852 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.5 KB
RevLine 
[524]1!
[1403]2! $Id: o3cm.F 1403 2010-07-01 09:02:53Z jyg $
[524]3!
4      SUBROUTINE o3cm (amb, bmb, sortie, ntab)
5      IMPLICIT none
6c======================================================================
7c Auteur(s): Z.X. Li (LMD/CNRS) date: 19930818
8c Objet: Ce programme calcule le contenu en ozone "sortie"
9c        (unite: cm.atm) entre deux niveaux "amb" et "bmb" (unite: mb)
10c        "ntab" est le nombre d'intervalles pour l'integration, sa
11c        valeur depend bien sur de l'epaisseur de la couche et de
12c        la precision qu'on souhaite a obtenir
13c======================================================================
14      REAL amb, bmb, sortie
15      INTEGER ntab
16c======================================================================
17      INTEGER n
18      REAL xtab(500), xa, xb, ya, yb, xincr
19c======================================================================
20      external mbtozm
[1403]21      CHARACTER (LEN=20) :: modname=''
22      CHARACTER (LEN=80) :: abort_message
[524]23c======================================================================
24c la fonction en ligne w(x) donne le profil de l'ozone en fonction
25c de l'altitude (unite: cm.atm / km)
26c (Green 1964, Appl. Opt. 3: 203-208)
27      REAL wp, xp, h, x, w, con
28      PARAMETER (wp=0.218, xp=23.25, h=4.63, con=1.0)
29      w(x) = wp/h * EXP((x-xp)/h)/ (con+EXP((x-xp)/h))**2
30c======================================================================
[1403]31      IF (ntab .GT. 499) THEN
32        abort_message = 'BIG ntab'
33        CALL abort_gcm (modname,abort_message,1)
34      ENDIF
35      xincr = (bmb-amb) / REAL(ntab)
[524]36      xtab(1) = amb
37      DO n = 2, ntab
38         xtab(n) = xtab(n-1) + xincr
39      ENDDO
40      xtab(ntab+1) = bmb
41      sortie = 0.0
42      DO n = 1, ntab
43         CALL mbtozm(xtab(n), xa)
44         CALL mbtozm(xtab(n+1), xb)
45         xa = xa / 1000.
46         xb = xb / 1000.
47         ya = w(xa)
48         yb = w(xb)
49         sortie = sortie + (ya+yb)/2.0 * ABS(xb-xa)
50      ENDDO
51      RETURN
52      END
53      SUBROUTINE mbtozm(rmb,zm)
54      IMPLICIT none
55c======================================================================
56c Auteur(s): Z.X. Li (LMD/CNRS)
57c Objet: transformer une hauteur de mb (rmb) en metre (zm)
58c======================================================================
59      REAL rmb, zm
60c======================================================================
61      REAL gama, tzero, pzero, g, r
62      PARAMETER (gama=6.5e-3, tzero=288., pzero=1013.25)
63      PARAMETER (g=9.81, r=287.0)
64      zm = tzero/gama * ( 1.-(rmb/pzero)**(r*gama/g) )
65      RETURN
66      END
Note: See TracBrowser for help on using the repository browser.