source: LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90 @ 408

Last change on this file since 408 was 395, checked in by lmdzadmin, 22 years ago

Changement dans les sorties histoires (rain devient precip pour eviter toute
confusion) + unites, commentaires des fichiers histoire
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.5 KB
Line 
1!
2! $Header$
3!
4
5  subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
6 &                     fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &
7 &                     ratqsbas,ratqshaut,if_ebil)
8
9   use IOIPSL
10   implicit none
11
12#include "conema3.h"
13#include "fisrtilp.inc"
14#include "nuage.h"
15!
16! Configuration de la "physique" de LMDZ a l'aide de la fonction
17! GETIN de IOIPSL
18!
19! LF 05/2001
20!
21
22!
23! ocean:      type d'ocean (force, slab, couple)
24! ok_veget:   type de modele de vegetation
25! ok_journe:  sorties journalieres
26! ok_mensuel: sorties mensuelles
27! ok_instan:  sorties instantanees
28
29
30
31! Sortie:
32  character (len = 6)  :: ocean
33  logical              :: ok_veget, ok_newmicro
34  logical              :: ok_journe, ok_mensuel, ok_instan
35  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
36  integer              :: iflag_cldcon, if_ebil
37! Local
38  integer              :: numout = 6
39
40!
41!
42!
43
44
45!Config Key  = OCEAN
46!Config Desc = Type d'ocean
47!Config Def  = force
48!Config Help = Type d'ocean utilise: force, slab,couple
49!
50  ocean = 'force '
51  call getin('OCEAN', ocean)
52!
53!Config Key  = VEGET
54!Config Desc = Type de modele de vegetation
55!Config Def  = .false.
56!Config Help = Type de modele de vegetation utilise
57!
58  ok_veget = .false.
59  call getin('VEGET', ok_veget)
60!
61!Config Key  = OK_journe
62!Config Desc = Pour des sorties journalieres
63!Config Def  = .false.
64!Config Help = Pour creer le fichier histday contenant les sorties
65!              journalieres
66!
67  ok_journe = .false.
68  call getin('OK_journe', ok_journe)
69!
70!Config Key  = OK_mensuel
71!Config Desc = Pour des sorties mensuelles
72!Config Def  = .true.
73!Config Help = Pour creer le fichier histmth contenant les sorties
74!              mensuelles
75!
76  ok_mensuel = .true.
77  call getin('OK_mensuel', ok_mensuel)
78!
79!Config Key  = OK_instan
80!Config Desc = Pour des sorties instantanees
81!Config Def  = .false.
82!Config Help = Pour creer le fichier histins contenant les sorties
83!              instantanees
84!
85  ok_instan = .false.
86  call getin('OK_instan', ok_instan)
87!
88!Config Key  = if_ebil
89!Config Desc = Niveau de sortie pour les diags bilan d'energie
90!Config Def  = 0
91!Config Help =
92!               
93!
94  if_ebil = 0
95  call getin('if_ebil', if_ebil)
96!!
97!! KE
98!
99!Config Key  = epmax
100!Config Desc = Efficacite precip
101!Config Def  = 0.993
102!Config Help =
103!
104  epmax = .993
105  call getin('epmax', epmax)
106!
107!Config Key  = ok_adj_ema
108!Config Desc = 
109!Config Def  = false
110!Config Help =
111!
112  ok_adj_ema = .false.
113  call getin('ok_adj_ema',ok_adj_ema)
114!
115!Config Key  = iflag_clw
116!Config Desc = 
117!Config Def  = 0
118!Config Help =
119!
120  iflag_clw = 0
121  call getin('iflag_clw',iflag_clw)
122!
123!Config Key  = cld_lc_lsc
124!Config Desc = 
125!Config Def  = 2.6e-4
126!Config Help =
127!
128  cld_lc_lsc = 2.6e-4
129  call getin('cld_lc_lsc',cld_lc_lsc)
130!
131!Config Key  = cld_lc_con
132!Config Desc = 
133!Config Def  = 2.6e-4
134!Config Help =
135!
136  cld_lc_con = 2.6e-4
137  call getin('cld_lc_con',cld_lc_con)
138!
139!Config Key  = cld_tau_lsc
140!Config Desc = 
141!Config Def  = 3600.
142!Config Help =
143!
144  cld_tau_lsc = 3600.
145  call getin('cld_tau_lsc',cld_tau_lsc)
146!
147!Config Key  = cld_tau_con
148!Config Desc = 
149!Config Def  = 3600.
150!Config Help =
151!
152  cld_tau_con = 3600.
153  call getin('cld_tau_con',cld_tau_con)
154!
155!Config Key  = ffallv_lsc
156!Config Desc = 
157!Config Def  = 1.
158!Config Help =
159!
160  ffallv_lsc = 1.
161  call getin('ffallv_lsc',ffallv_lsc)
162!
163!Config Key  = ffallv_con
164!Config Desc = 
165!Config Def  = 1.
166!Config Help =
167!
168  ffallv_con = 1.
169  call getin('ffallv_con',ffallv_con)
170!
171!Config Key  = coef_eva
172!Config Desc = 
173!Config Def  = 2.e-5
174!Config Help =
175!
176  coef_eva = 2.e-5
177  call getin('coef_eva',coef_eva)
178!
179!Config Key  = reevap_ice
180!Config Desc = 
181!Config Def  = .false.
182!Config Help =
183!
184  reevap_ice = .false.
185  call getin('reevap_ice',reevap_ice)
186!
187!Config Key  = iflag_cldcon
188!Config Desc = 
189!Config Def  = 1
190!Config Help =
191!
192  iflag_cldcon = 1
193  call getin('iflag_cldcon',iflag_cldcon)
194
195!
196!Config Key  = iflag_pdf
197!Config Desc = 
198!Config Def  = 0
199!Config Help =
200!
201  iflag_pdf = 0
202  call getin('iflag_pdf',iflag_pdf)
203!
204!Config Key  = fact_cldcon
205!Config Desc = 
206!Config Def  = 0.375
207!Config Help =
208!
209  fact_cldcon = 0.375
210  call getin('fact_cldcon',fact_cldcon)
211
212!
213!Config Key  = facttemps
214!Config Desc = 
215!Config Def  = 1.e-4
216!Config Help =
217!
218  facttemps = 1.e-4
219  call getin('facttemps',facttemps)
220
221!
222!Config Key  = ok_newmicro
223!Config Desc = 
224!Config Def  = .true.
225!Config Help =
226!
227  ok_newmicro = .true.
228  call getin('ok_newmicro',ok_newmicro)
229!
230!Config Key  = ratqsbas
231!Config Desc = 
232!Config Def  = 0.01
233!Config Help =
234!
235  ratqsbas = 0.01
236  call getin('ratqsbas',ratqsbas)
237!
238!Config Key  = ratqshaut
239!Config Desc = 
240!Config Def  = 0.3
241!Config Help =
242!
243  ratqshaut = 0.3
244  call getin('ratqshaut',ratqshaut)
245
246!
247!Config Key  = rad_froid
248!Config Desc = 
249!Config Def  = 35.0
250!Config Help =
251!
252  rad_froid = 35.0
253  call getin('rad_froid',rad_froid)
254
255!
256!Config Key  = rad_chau1
257!Config Desc = 
258!Config Def  = 13.0
259!Config Help =
260!
261  rad_chau1 = 13.0
262  call getin('rad_chau1',rad_chau1)
263
264!
265!Config Key  = rad_chau2
266!Config Desc = 
267!Config Def  = 9.0
268!Config Help =
269!
270  rad_chau2 = 9.0
271  call getin('rad_chau2',rad_chau2)
272
273!
274!Config Key  =
275!Config Desc = 
276!Config Def  =
277!Config Help =
278!
279!   =
280!  call getin('',)
281!
282!
283!
284!
285
286  write(numout,*)' ##############################################'
287  write(numout,*)' Configuration des parametres de la physique: '
288  write(numout,*)' Config ocean = ', ocean
289  write(numout,*)' Config veget = ', ok_veget
290  write(numout,*)' Sortie journaliere = ', ok_journe
291  write(numout,*)' Sortie mensuelle = ', ok_mensuel
292  write(numout,*)' Sortie instantanee = ', ok_instan
293  write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil
294  write(numout,*)' epmax = ', epmax
295  write(numout,*)' ok_adj_ema = ', ok_adj_ema
296  write(numout,*)' iflag_clw = ', iflag_clw
297  write(numout,*)' cld_lc_lsc = ', cld_lc_lsc
298  write(numout,*)' cld_lc_con = ', cld_lc_con
299  write(numout,*)' cld_tau_lsc = ', cld_tau_lsc
300  write(numout,*)' cld_tau_con = ', cld_tau_con
301  write(numout,*)' ffallv_lsc = ', ffallv_lsc
302  write(numout,*)' ffallv_con = ', ffallv_con
303  write(numout,*)' coef_eva = ', coef_eva
304  write(numout,*)' reevap_ice = ', reevap_ice
305  write(numout,*)' iflag_pdf = ', iflag_pdf
306  write(numout,*)' iflag_cldcon = ', iflag_cldcon
307  write(numout,*)' fact_cldcon = ', fact_cldcon
308  write(numout,*)' facttemps = ', facttemps
309  write(numout,*)' ok_newmicro = ',ok_newmicro
310  write(numout,*)' ratqsbas = ',ratqsbas
311  write(numout,*)' ratqshaut = ',ratqshaut
312
313  return
314
315  end subroutine conf_phys
316
Note: See TracBrowser for help on using the repository browser.