source: LMDZ5/branches/LMDZ5-dev032011/libf/phylmd/printflag.F @ 4225

Last change on this file since 4225 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: 5.4 KB
Line 
1!
2! $Header$
3!
4       SUBROUTINE  printflag( tabcntr0, radpas,
5     ,                        ok_journe,ok_instan,ok_region        )
6c
7
8c
9c      Auteur :  P. Le Van
10
11       IMPLICIT NONE
12
13       REAL tabcntr0( 100 )
14       LOGICAL cycle_diurn0,soil_model0,new_oliq0,ok_orodr0
15       LOGICAL ok_orolf0,ok_limitvr0
16       LOGICAL ok_journe,ok_instan,ok_region
17       INTEGER radpas , radpas0
18c
19#include "clesphys.h"
20c
21c
22       PRINT 100
23       PRINT *,' *******************************************************
24     ,************'
25       PRINT *,' ********   Choix  des principales  cles de la physique
26     ,   *********'
27       PRINT *,' *******************************************************
28     ,************'
29       PRINT 100
30       PRINT 10, cycle_diurne,  soil_model 
31       PRINT 100
32
33       IF   (    iflag_con.EQ. 1 )   THEN
34           PRINT *,' *****           Shema  convection   LMD           
35     ,          ******'
36       ELSE IF ( iflag_con.EQ. 2 )   THEN
37           PRINT *,' *****           Shema  convection  Tiedtke 
38     ,          ******'
39       ELSE IF ( iflag_con.GE. 3 )   THEN
40           PRINT *,' *****           Shema  convection    Emanuel     
41     ,          ******'
42       ENDIF
43       PRINT 100
44
45       PRINT 11, new_oliq, ok_orodr, ok_orolf   
46       PRINT 100
47
48       PRINT 7,  ok_limitvrai   
49       PRINT 100
50
51       PRINT 12, nbapp_rad
52       PRINT 100
53
54       PRINT 8, radpas
55       PRINT 100
56
57       PRINT 4,ok_journe,ok_instan,ok_region
58       PRINT 100
59       PRINT 100
60c
61c
62        cycle_diurn0  = .FALSE.
63        soil_model0   = .FALSE.
64        new_oliq0     = .FALSE.
65        ok_orodr0     = .FALSE.
66        ok_orolf0     = .FALSE.
67        ok_limitvr0   = .FALSE.
68
69        IF( tabcntr0( 7 ).EQ. 1. )   cycle_diurn0 = .TRUE.
70        IF( tabcntr0( 8 ).EQ. 1. )    soil_model0 = .TRUE.
71        IF( tabcntr0( 9 ).EQ. 1. )      new_oliq0 = .TRUE.
72        IF( tabcntr0(10 ).EQ. 1. )      ok_orodr0 = .TRUE.
73        IF( tabcntr0(11 ).EQ. 1. )      ok_orolf0 = .TRUE.
74        IF( tabcntr0(12 ).EQ. 1. )    ok_limitvr0 = .TRUE.
75
76        PRINT *,' $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
77     ,$$$$$$$$$$$$$'
78        PRINT 100
79c
80       IF( INT( tabcntr0( 5 ) ) .NE. iflag_con  )   THEN
81        PRINT 20, INT(tabcntr0(5)), iflag_con
82        PRINT 100
83       ENDIF
84
85       IF( INT( tabcntr0( 6 ) ) .NE. nbapp_rad  )   THEN
86        PRINT 21,  INT(tabcntr0(6)), nbapp_rad
87!        radpas0  = NINT( 86400./tabcntr0(1)/INT( tabcntr0(6) ) )
88        PRINT 100
89        PRINT 22, radpas0, radpas
90        PRINT 100
91       ENDIF
92
93       IF( cycle_diurn0.AND..NOT.cycle_diurne.OR..NOT.cycle_diurn0.AND.
94     ,        cycle_diurne )     THEN
95        PRINT 13, cycle_diurn0, cycle_diurne
96        PRINT 100
97       ENDIF
98
99       IF( soil_model0.AND..NOT.soil_model.OR..NOT.soil_model0.AND.
100     ,        soil_model )     THEN
101        PRINT 14, soil_model0, soil_model
102        PRINT 100
103       ENDIF
104
105       IF( new_oliq0.AND..NOT.new_oliq.OR..NOT.new_oliq0.AND.
106     ,        new_oliq )     THEN
107        PRINT 16, new_oliq0, new_oliq
108        PRINT 100
109       ENDIF
110
111       IF( ok_orodr0.AND..NOT.ok_orodr.OR..NOT.ok_orodr0.AND.
112     ,        ok_orodr )     THEN
113        PRINT 15, ok_orodr0, ok_orodr
114        PRINT 100
115       ENDIF
116
117       IF( ok_orolf0.AND..NOT.ok_orolf.OR..NOT.ok_orolf0.AND.
118     ,        ok_orolf )     THEN
119        PRINT 17, ok_orolf0, ok_orolf
120        PRINT 100
121       ENDIF
122
123       IF( ok_limitvr0.AND..NOT.ok_limitvrai.OR..NOT.ok_limitvr0.
124     ,     AND.ok_limitvrai )     THEN
125        PRINT 18, ok_limitvr0, ok_limitvrai
126        PRINT 100
127       ENDIF
128
129       PRINT 100
130       PRINT *,' *******************************************************
131     ,************'
132       PRINT 100
133
134 4    FORMAT(2x,5("*"),'  ok_journe= ',l3,3x,',ok_instan = ',
135     , l3,3x,',ok_region = ',l3,3x,5("*") )
136
137 7     FORMAT(2x,5("*"),15x,'      ok_limitvrai   = ',l3,16x,5("*") )
138
139 8     FORMAT(2x,'*****             radpas    =                      ' ,
140     , i4,6x,' *****')
141
142 10    FORMAT(2x,5("*"),'    Cycle_diurne = ',l3,4x,', Soil_model = ',
143     , l3,12x,6("*") )
144
145
146 11    FORMAT(2x,5("*"),'  new_oliq = ',l3,3x,', Ok_orodr = ',
147     , l3,3x,', Ok_orolf = ',l3,3x,5("*") )
148
149
150 12    FORMAT(2x,'*****  Nb d appels /jour des routines de rayonn. = ' ,
151     , i4,6x,' *****')
152
153 13    FORMAT(2x,'$$$$$$$$   Attention !!  cycle_diurne  different  sur',
154     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
155
156 14    FORMAT(2x,'$$$$$$$$   Attention !!    soil_model  different  sur',
157     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
158
159 15    FORMAT(2x,'$$$$$$$$   Attention !!      ok_orodr  different  sur',
160     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
161
162 16    FORMAT(2x,'$$$$$$$$   Attention !!      new_oliq  different  sur',
163     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
164
165 17    FORMAT(2x,'$$$$$$$$   Attention !!      ok_orolf  different  sur',
166     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
167
168 18    FORMAT(2x,'$$$$$$$$   Attention !!  ok_limitvrai  different  sur',
169     , /1x,10x,' startphy = ',l3,2x,' et  run.def = ',l3)
170
171 20    FORMAT(/2x,'$$$$$$$$   Attention !!    iflag_con  different  sur',
172     , /1x,10x,' startphy = ',i3,2x,' et  run.def = ',i3 )
173
174 21    FORMAT(2x,'$$$$$$$$   Attention !!     nbapp_rad  different  sur',
175     , /1x,10x,' startphy = ',i3,2x,' et  run.def = ',i3 )
176
177 22    FORMAT(2x,'$$$$$$$$   Attention !!        radpas  different  sur',
178     , /1x,10x,' startphy = ',i3,2x,' et  run.def = ',i3 )
179
180 100   FORMAT(/)
181
182       RETURN
183       END
Note: See TracBrowser for help on using the repository browser.