source: LMDZ5/branches/LMDZ5_AR5/libf/phylmd/ini_histday_seri.h @ 2072

Last change on this file since 2072 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: 4.5 KB
Line 
1c
2c $Id: ini_histday_seri.h 1403 2010-07-01 09:02:53Z lguez $
3c
4cym Ne fonctionnera pas en mode parallele
5      IF (is_sequential) THEN
6     
7      IF (type_run.EQ."AMIP") THEN
8c
9       zstophy = dtime
10       zout = ecrit_day
11c
12         idayref = day_ref
13         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14c
15         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
16         DO i = 1, iim
17            zx_lon(i,1) = rlon(i+1)
18            zx_lon(i,jjmp1) = rlon(i+1)
19         ENDDO
20         DO ll=1,klev
21            znivsig(ll)=REAL(ll)
22         ENDDO
23         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
24c
25         imin_debut=1 
26         nbpti=1
27         jmin_debut=1 
28         nbptj=1
29c
30         CALL histbeg("histday_seri.nc", 
31     .                 iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
32     .                 imin_debut,nbpti,jmin_debut,nbptj,
33     .                 itau_phy, zjulian, dtime,
34     .                 nhori, nid_day_seri)
35c
36         CALL histvert(nid_day_seri, "presnivs", 
37     .                "Vertical levels","mb",
38     .                 klev, presnivs/100., nvert)
39c
40         CALL histdef(nid_day_seri, "bilTOA", 
41     .                "Net radiation at model top", "W/m2",
42     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
43     .                "ave(X)", zstophy,zout)
44c
45         CALL histdef(nid_day_seri, "bils", 
46     .                "Net downward energy flux at surface","W/m2",
47     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
48     .                "ave(X)", zstophy,zout)
49c
50         CALL histdef(nid_day_seri, "ecin", 
51     .                "Total kinetic energy (per unit area)","J/m2",
52     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
53     .                "ave(X)", zstophy,zout)
54c
55cIM 151004 BEG
56         IF(1.EQ.0) THEN
57c
58         CALL histdef(nid_day_seri, "momang", 
59     .               "Total relative angular momentum (per unit area)",
60     .               "kg/s",
61     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
62     .                "ave(X)", zstophy,zout)
63c
64         CALL histdef(nid_day_seri, "frictor", 
65     .               "Friction torque (per unit area)", "N/m",
66     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
67     .                "ave(X)", zstophy,zout)
68c
69         CALL histdef(nid_day_seri, "mountor", 
70     .               "Mountain torque (per unit area)", "N/m",
71     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
72     .                "ave(X)", zstophy,zout)
73c
74         ENDIF !(1.EQ.0) THEN
75c
76         CALL histdef(nid_day_seri, "momang", 
77     .               "Axial angular momentum (per unit area)",
78     .               "kg/s",
79     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
80     .                "ave(X)", zstophy,zout)
81c
82         CALL histdef(nid_day_seri, "torsfc", 
83     .        "Total surface torque (including mountain torque)", "N/m",
84     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
85     .                "ave(X)", zstophy,zout)
86c
87cIM 151004 END       
88c
89         CALL histdef(nid_day_seri, "tamv", 
90     .                "Temperature (mass-weighted vert. ave)", "K",
91     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
92     .                "ave(X)", zstophy,zout)
93c
94         CALL histdef(nid_day_seri, "psol", 
95     .                "Surface pressure", "Pa",
96     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
97     .                "ave(X)", zstophy,zout)
98c
99         CALL histdef(nid_day_seri, "evap", 
100     .                "Evaporation and sublimation (per unit area)", 
101     .                "kg/(m2*s)",
102     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
103     .                "ave(X)", zstophy,zout)
104c
105c          call histdef(nid_day_seri, 
106c    .         "SnowFrac", 
107c    .         "Snow-covered area ", "%", 
108c    .         iim,jjmp1,nhori, 1,1,1, -99, 32,
109c    .         "ave(X)", zstophy,zout)
110c
111c        CALL histdef(nid_day_seri, "snow_depth", 
112cIM 080904  .                "Snow Depth (water equivalent)", "m",
113cIM 191104  .                "Snow Depth (water equivalent)", "kg/m2",
114c    .                "Snow Mass", "kg/m2",
115c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
116c    .               "ave(X)", zstophy,zout)
117c
118           call histdef(nid_day_seri, 
119     .         "tsol_"//clnsurf(is_oce),
120     .         "SST over open (ice-free) ocean ", "K", 
121     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
122     .         "ave(X)", zstophy,zout)
123c
124c=================================================================
125c
126         CALL histend(nid_day_seri)
127c
128c=================================================================
129      ENDIF ! fin de test sur type_run.EQ.AMIP
130     
131      ENDIF ! is_sequential
Note: See TracBrowser for help on using the repository browser.