source: LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/printflag.F @ 3722

Last change on this file since 3722 was 1311, checked in by Laurent Fairhead, 15 years ago

Modifications to thermals for TKE transport


Modifications aux thermiques pour le transport de la TKE

pbl_surface_mode.F90 : ok_flux_surf=.false. seulement pour klon>1
physiq.F : option iflag_pbl=10 pour transporter la TKE avec les thermiques.
calltherm.F90 : passage de iflag_thermals_ed en argument pour thermcell_plume
thermcell_main.F90 : Appel a plusieurs version de thermcell_plume en option
thermcell_plume.F90 : plusieurs versions dans le meme fichier (temporaire)
thermcell_height.F90 : verrue pour les cas ou les thermiques montent tout

en haut

yamada4 : inclusion de la diffusion verticale en option iflag_pbl=9

+ variables anciennement common, puis save/allocatable, remises en local

  • 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.