source: LMDZ5/branches/testing/libf/phylmd/conf_phys_m.F90 @ 1856

Last change on this file since 1856 was 1795, checked in by Ehouarn Millour, 11 years ago

Version testing basee sur la r1794


Testing release based on r1794

File size: 52.8 KB
Line 
1!
2! $Id: conf_phys.F90 1668 2012-10-12 10:47:37Z idelkadi $
3!
4!
5!
6module conf_phys_m
7
8   implicit none
9
10contains
11
12  subroutine conf_phys(ok_journe, ok_mensuel, ok_instan, ok_hf, &
13                       ok_LES,&
14                       callstats,&
15                       solarlong0,seuil_inversion, &
16                       fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
17                       iflag_cldcon, &
18                       iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
19                       ok_ade, ok_aie, ok_cdnc, aerosol_couple, &
20                       flag_aerosol, flag_aerosol_strat, new_aod, &
21                       bl95_b0, bl95_b1,&
22                       read_climoz, &
23                       alp_offset)
24
25   use IOIPSL
26   USE surface_data
27   USE phys_cal_mod
28   USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl
29   use control_mod
30
31 include "conema3.h"
32 include "fisrtilp.h"
33 include "nuage.h"
34 include "YOMCST.h"
35 include "YOMCST2.h"
36
37 include "thermcell.h"
38 include "iniprint.h"
39
40!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
41include "clesphys.h"
42include "compbl.h"
43include "comsoil.h"
44!
45! Configuration de la "physique" de LMDZ a l'aide de la fonction
46! GETIN de IOIPSL
47!
48! LF 05/2001
49!
50
51!
52! type_ocean:      type d'ocean (force, slab, couple)
53! version_ocean:   version d'ocean (opa8/nemo pour type_ocean=couple ou
54!                                   sicOBS pour type_ocean=slab)
55! ok_veget:   type de modele de vegetation
56! ok_journe:  sorties journalieres
57! ok_hf:  sorties haute frequence
58! ok_mensuel: sorties mensuelles
59! ok_instan:  sorties instantanees
60! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
61! ok_cdnc, ok cloud droplet number concentration
62! flag_aerosol_strat : flag pour les aerosols stratos
63! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
64!
65
66
67! Sortie:
68  logical              :: ok_newmicro
69  integer              :: iflag_radia
70  logical              :: ok_journe, ok_mensuel, ok_instan, ok_hf
71  logical              :: ok_LES
72  LOGICAL              :: callstats
73  LOGICAL              :: ok_ade, ok_aie, ok_cdnc, aerosol_couple
74  INTEGER              :: flag_aerosol
75  LOGICAL              :: flag_aerosol_strat
76  LOGICAL              :: new_aod
77  REAL                 :: bl95_b0, bl95_b1
78  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
79  integer              :: iflag_cldcon
80  integer              :: iflag_ratqs
81
82  character (len = 6),SAVE  :: type_ocean_omp, version_ocean_omp, ocean_omp
83  character (len = 10),SAVE  :: type_veget_omp
84  CHARACTER(len = 8),SAVE   :: aer_type_omp
85  logical,SAVE              :: ok_newmicro_omp
86  logical,SAVE        :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp       
87  logical,SAVE        :: ok_LES_omp   
88  LOGICAL,SAVE        :: callstats_omp
89  LOGICAL,SAVE        :: ok_ade_omp, ok_aie_omp, ok_cdnc_omp, aerosol_couple_omp
90  INTEGER, SAVE       :: flag_aerosol_omp
91  LOGICAL, SAVE       :: flag_aerosol_strat_omp
92  LOGICAL, SAVE       :: new_aod_omp
93  REAL,SAVE           :: bl95_b0_omp, bl95_b1_omp
94  REAL,SAVE           :: freq_ISCCP_omp, ecrit_ISCCP_omp
95  REAL,SAVE           :: freq_COSP_omp
96  real,SAVE           :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
97  real,SAVE           :: ratqshaut_omp
98  real,SAVE           :: tau_ratqs_omp
99  integer,SAVE        :: iflag_radia_omp
100  integer,SAVE        :: iflag_rrtm_omp
101  integer,SAVE        :: iflag_cldcon_omp, ip_ebil_phy_omp
102  integer,SAVE        :: iflag_ratqs_omp
103
104  Real,SAVE           :: f_cdrag_ter_omp,f_cdrag_oce_omp
105  Real,SAVE           :: f_rugoro_omp   
106
107! Local
108  real                 :: zzz
109
110  real :: seuil_inversion
111  real,save :: seuil_inversion_omp
112
113  integer,SAVE :: iflag_thermals_ed_omp,iflag_thermals_optflux_omp
114  integer,SAVE :: iflag_thermals_omp,nsplit_thermals_omp
115  real,save :: tau_thermals_omp,alp_bl_k_omp
116!!! nrlmd le 10/04/2012
117  integer,SAVE :: iflag_trig_bl_omp,iflag_clos_bl_omp
118  integer,SAVE :: tau_trig_shallow_omp,tau_trig_deep_omp
119  real,SAVE    :: s_trig_omp
120!!! fin nrlmd le 10/04/2012
121  real :: alp_offset
122  REAL, SAVE :: alp_offset_omp
123  integer,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp
124  integer,SAVE :: iflag_cvl_sigd_omp
125  REAL, SAVE :: supcrit1_omp, supcrit2_omp
126  INTEGER, SAVE :: iflag_mix_omp
127  real, save :: scut_omp, qqa1_omp, qqa2_omp, gammas_omp, Fmax_omp, alphas_omp
128
129  REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp
130  REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp
131  REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp
132  REAL,SAVE :: N2O_ppb_omp, RN2O_omp, N2O_ppb_per_omp, RN2O_per_omp
133  REAL,SAVE :: CFC11_ppt_omp,RCFC11_omp,CFC11_ppt_per_omp,RCFC11_per_omp
134  REAL,SAVE :: CFC12_ppt_omp,RCFC12_omp,CFC12_ppt_per_omp,RCFC12_per_omp
135  LOGICAL,SAVE :: ok_4xCO2atm_omp
136  REAL,SAVE :: epmax_omp
137  LOGICAL,SAVE :: ok_adj_ema_omp
138  INTEGER,SAVE :: iflag_clw_omp
139  REAL,SAVE :: cld_lc_lsc_omp,cld_lc_con_omp,cld_tau_lsc_omp,cld_tau_con_omp
140  REAL,SAVE :: ffallv_lsc_omp, ffallv_con_omp,coef_eva_omp
141  LOGICAL,SAVE :: reevap_ice_omp
142  INTEGER,SAVE :: iflag_pdf_omp
143  REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
144  REAL,SAVE :: t_glace_min_omp, t_glace_max_omp
145  REAL,SAVE :: rei_min_omp, rei_max_omp
146  REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
147  REAL,SAVE :: qsol0_omp
148  REAL      :: solarlong0
149  REAL,SAVE :: solarlong0_omp
150  INTEGER,SAVE :: top_height_omp,overlap_omp
151  REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp
152  LOGICAL,SAVE :: ok_kzmin_omp
153  REAL, SAVE ::  fmagic_omp, pmagic_omp
154  INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
155  Integer, save :: lev_histins_omp, lev_histLES_omp
156  INTEGER, SAVE :: lev_histdayNMC_omp
157  LOGICAL, SAVE :: ok_histNMC_omp(3)
158  REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3)
159  CHARACTER*4, SAVE :: type_run_omp
160  LOGICAL,SAVE :: ok_isccp_omp
161  LOGICAL,SAVE :: ok_cosp_omp
162  LOGICAL,SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
163  REAL,SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
164  REAL,SAVE :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp
165  REAL,SAVE :: ecrit_ins_omp
166  REAL,SAVE :: ecrit_LES_omp
167  REAL,SAVE :: ecrit_tra_omp
168  REAL,SAVE :: cvl_corr_omp
169  LOGICAL,SAVE :: ok_lic_melt_omp
170!
171  LOGICAL,SAVE :: cycle_diurne_omp,soil_model_omp,new_oliq_omp
172  LOGICAL,SAVE :: ok_orodr_omp, ok_orolf_omp, ok_limitvrai_omp
173  INTEGER, SAVE :: nbapp_rad_omp, iflag_con_omp
174  INTEGER, SAVE :: iflag_ener_conserv_omp
175  LOGICAL,SAVE :: ok_strato_omp
176  LOGICAL,SAVE :: ok_hines_omp
177  LOGICAL,SAVE      :: carbon_cycle_tr_omp
178  LOGICAL,SAVE      :: carbon_cycle_cpl_omp
179
180  integer, intent(out):: read_climoz ! read ozone climatology, OpenMP shared
181  ! Allowed values are 0, 1 and 2
182  ! 0: do not read an ozone climatology
183  ! 1: read a single ozone climatology that will be used day and night
184  ! 2: read two ozone climatologies, the average day and night
185  ! climatology and the daylight climatology
186
187!$OMP MASTER
188!Config Key  = type_ocean
189!Config Desc = Type d'ocean
190!Config Def  = force
191!Config Help = Type d'ocean utilise: force, slab,couple
192!
193  type_ocean_omp = 'force '
194  call getin('type_ocean', type_ocean_omp)
195!
196!Config Key  = version_ocean
197!Config Desc = Version d'ocean
198!Config Def  = xxxxxx
199!Config Help = Version d'ocean utilise: opa8/nemo/sicOBS/xxxxxx
200!
201  version_ocean_omp = 'xxxxxx'
202  call getin('version_ocean', version_ocean_omp)
203
204!Config Key  = OCEAN
205!Config Desc = Old parameter name for type_ocean
206!Config Def  = yyyyyy
207!Config Help = This is only for testing purpose
208!
209  ocean_omp = 'yyyyyy'
210  call getin('OCEAN', ocean_omp)
211  IF (ocean_omp /= 'yyyyyy') THEN
212     WRITE(lunout,*)'ERROR!! Old variable name OCEAN used in parmeter file.'
213     WRITE(lunout,*)'Variable OCEAN has been replaced by the variable type_ocean.'
214     WRITE(lunout,*)'You have to update your parameter file physiq.def to succed running'
215     CALL abort_gcm('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1)
216  END IF
217
218!
219!Config Key  = VEGET
220!Config Desc = Type de modele de vegetation
221!Config Def  = .false.
222!Config Help = Type de modele de vegetation utilise
223!
224  type_veget_omp ='orchidee'
225  call getin('VEGET', type_veget_omp)
226!
227!Config Key  = OK_journe
228!Config Desc = Pour des sorties journalieres
229!Config Def  = .false.
230!Config Help = Pour creer le fichier histday contenant les sorties
231!              journalieres
232!
233  ok_journe_omp = .false.
234  call getin('OK_journe', ok_journe_omp)
235!
236!Config Key  = ok_hf
237!Config Desc = Pour des sorties haute frequence
238!Config Def  = .false.
239!Config Help = Pour creer le fichier histhf contenant les sorties
240!              haute frequence ( 3h ou 6h)
241!
242  ok_hf_omp = .false.
243  call getin('ok_hf', ok_hf_omp)
244!
245!Config Key  = OK_mensuel
246!Config Desc = Pour des sorties mensuelles
247!Config Def  = .true.
248!Config Help = Pour creer le fichier histmth contenant les sorties
249!              mensuelles
250!
251  ok_mensuel_omp = .true.
252  call getin('OK_mensuel', ok_mensuel_omp)
253!
254!Config Key  = OK_instan
255!Config Desc = Pour des sorties instantanees
256!Config Def  = .false.
257!Config Help = Pour creer le fichier histins contenant les sorties
258!              instantanees
259!
260  ok_instan_omp = .false.
261  call getin('OK_instan', ok_instan_omp)
262!
263!Config Key  = ok_ade
264!Config Desc = Aerosol direct effect or not?
265!Config Def  = .false.
266!Config Help = Used in radlwsw.F
267!
268  ok_ade_omp = .false.
269  call getin('ok_ade', ok_ade_omp)
270
271!
272!Config Key  = ok_aie
273!Config Desc = Aerosol indirect effect or not?
274!Config Def  = .false.
275!Config Help = Used in nuage.F and radlwsw.F
276!
277  ok_aie_omp = .false.
278  call getin('ok_aie', ok_aie_omp)
279
280!
281!Config Key  = ok_cdnc
282!Config Desc = ok cloud droplet number concentration
283!Config Def  = .false.
284!Config Help = Used in newmicro.F
285!
286  ok_cdnc_omp = .false.
287  call getin('ok_cdnc', ok_cdnc_omp)
288!
289!Config Key  = aerosol_couple
290!Config Desc = read aerosol in file or calcul by inca
291!Config Def  = .false.
292!Config Help = Used in physiq.F
293!
294  aerosol_couple_omp = .false.
295  CALL getin('aerosol_couple',aerosol_couple_omp)
296!
297!Config Key  = flag_aerosol
298!Config Desc = which aerosol is use for coupled model
299!Config Def  = 1
300!Config Help = Used in physiq.F
301!
302! - flag_aerosol=0 => no aerosol
303! - flag_aerosol=1 => so4 only (defaut)
304! - flag_aerosol=2 => bc  only
305! - flag_aerosol=3 => pom only
306! - flag_aerosol=4 => seasalt only
307! - flag_aerosol=5 => dust only
308! - flag_aerosol=6 => all aerosol
309
310  flag_aerosol_omp = 0
311  CALL getin('flag_aerosol',flag_aerosol_omp)
312!
313!Config Key  = flag_aerosol_strat
314!Config Desc = use stratospheric aerosols T/F
315!Config Def  = false
316!Config Help = Used in physiq.F
317!
318!
319  flag_aerosol_strat_omp = .false.
320  CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
321
322! Temporary variable for testing purpose!!
323!Config Key  = new_aod
324!Config Desc = which calcul of aeropt
325!Config Def  = false
326!Config Help = Used in physiq.F
327!
328  new_aod_omp = .true.
329  CALL getin('new_aod',new_aod_omp)
330
331!
332!Config Key  = aer_type
333!Config Desc = Use a constant field for the aerosols
334!Config Def  = scenario
335!Config Help = Used in readaerosol.F90
336!
337  aer_type_omp = 'scenario'
338  call getin('aer_type', aer_type_omp)
339
340!
341!Config Key  = bl95_b0
342!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
343!Config Def  = .false.
344!Config Help = Used in nuage.F
345!
346  bl95_b0_omp = 2.
347  call getin('bl95_b0', bl95_b0_omp)
348
349!Config Key  = bl95_b1
350!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
351!Config Def  = .false.
352!Config Help = Used in nuage.F
353!
354  bl95_b1_omp = 0.2
355  call getin('bl95_b1', bl95_b1_omp)
356
357!Config Key  = freq_ISCCP
358!Config Desc = Frequence d'appel du simulateur ISCCP en secondes;
359!              par defaut 10800, i.e. 3 heures
360!Config Def  = 10800.
361!Config Help = Used in ini_histISCCP.h
362!
363  freq_ISCCP_omp = 10800.
364  call getin('freq_ISCCP', freq_ISCCP_omp)
365!
366!Config Key  = ecrit_ISCCP
367!Config Desc = Frequence d'ecriture des resultats du simulateur ISCCP en nombre de jours;
368!              par defaut 1., i.e. 1 jour
369!Config Def  = 1.
370!Config Help = Used in ini_histISCCP.h
371!
372!
373  ecrit_ISCCP_omp = 1.
374  call getin('ecrit_ISCCP', ecrit_ISCCP_omp)
375
376!Config Key  = freq_COSP
377!Config Desc = Frequence d'appel du simulateur COSP en secondes;
378!              par defaut 10800, i.e. 3 heures
379!Config Def  = 10800.
380!Config Help = Used in ini_histdayCOSP.h
381!
382  freq_COSP_omp = 10800.
383  call getin('freq_COSP', freq_COSP_omp)
384
385!
386!Config Key  = ip_ebil_phy
387!Config Desc = Niveau de sortie pour les diags bilan d'energie
388!Config Def  = 0
389!Config Help =
390!               
391  ip_ebil_phy_omp = 0
392  call getin('ip_ebil_phy', ip_ebil_phy_omp)
393!
394!Config Key  = seuil_inversion
395!Config Desc = Seuil ur dTh pour le choix entre les schemas de CL
396!Config Def  = -0.1
397!Config Help =
398!               
399  seuil_inversion_omp = -0.1
400  call getin('seuil_inversion', seuil_inversion_omp)
401
402!!
403!! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
404!!
405!Config Key  = R_ecc
406!Config Desc = Excentricite
407!Config Def  = 0.016715
408!Config Help =
409!               
410!valeur AMIP II
411  R_ecc_omp = 0.016715
412  call getin('R_ecc', R_ecc_omp)
413!!
414!Config Key  = R_peri
415!Config Desc = Equinoxe
416!Config Def  =
417!Config Help =
418!               
419!
420!valeur AMIP II
421  R_peri_omp = 102.7
422  call getin('R_peri', R_peri_omp)
423!!
424!Config Key  = R_incl
425!Config Desc = Inclinaison
426!Config Def  =
427!Config Help =
428!               
429!
430!valeur AMIP II
431  R_incl_omp = 23.441
432  call getin('R_incl', R_incl_omp)
433!!
434!Config Key  = solaire
435!Config Desc = Constante solaire en W/m2
436!Config Def  = 1365.
437!Config Help =
438!               
439!
440!valeur AMIP II
441  solaire_omp = 1365.
442  call getin('solaire', solaire_omp)
443!!
444!Config Key  = co2_ppm
445!Config Desc = concentration du gaz carbonique en ppmv
446!Config Def  = 348.
447!Config Help =
448!               
449!
450!valeur AMIP II
451  co2_ppm_omp = 348.
452  call getin('co2_ppm', co2_ppm_omp)
453!!
454!Config Key  = RCO2
455!Config Desc = Concentration du CO2
456!Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
457!Config Def  = 348. * 1.0e-06  * 44.011/28.97
458!Config Help =
459!               
460! RCO2 = 5.286789092164308E-04
461!ancienne valeur
462  RCO2_omp = co2_ppm_omp * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.
463
464!!  call getin('RCO2', RCO2)
465!!
466!Config Key  = RCH4
467!Config Desc = Concentration du CH4
468!Config Def  = 1.65E-06* 16.043/28.97
469!Config Help =
470!               
471!
472!valeur AMIP II
473!OK  RCH4 = 1.65E-06* 16.043/28.97
474! RCH4 = 9.137366240938903E-07
475!
476!ancienne valeur
477! RCH4 = 1.72E-06* 16.043/28.97
478!OK call getin('RCH4', RCH4)
479  zzz = 1650.
480  call getin('CH4_ppb', zzz)
481  CH4_ppb_omp = zzz
482  RCH4_omp = CH4_ppb_omp * 1.0E-09 * 16.043/28.97
483!!
484!Config Key  = RN2O
485!Config Desc = Concentration du N2O
486!Config Def  = 306.E-09* 44.013/28.97
487!Config Help =
488!               
489!
490!valeur AMIP II
491!OK  RN2O = 306.E-09* 44.013/28.97
492! RN2O = 4.648939592682085E-07
493!
494!ancienne valeur
495! RN2O = 310.E-09* 44.013/28.97
496!OK  call getin('RN2O', RN2O)
497  zzz=306.
498  call getin('N2O_ppb', zzz)
499  N2O_ppb_omp = zzz
500  RN2O_omp = N2O_ppb_omp * 1.0E-09 * 44.013/28.97
501!!
502!Config Key  = RCFC11
503!Config Desc = Concentration du CFC11
504!Config Def  = 280.E-12* 137.3686/28.97
505!Config Help =
506!               
507!
508!OK RCFC11 = 280.E-12* 137.3686/28.97
509  zzz = 280.
510  call getin('CFC11_ppt',zzz)
511  CFC11_ppt_omp = zzz
512  RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * 137.3686/28.97
513! RCFC11 = 1.327690990680013E-09
514!OK call getin('RCFC11', RCFC11)
515!!
516!Config Key  = RCFC12
517!Config Desc = Concentration du CFC12
518!Config Def  = 484.E-12* 120.9140/28.97
519!Config Help =
520!               
521!
522!OK RCFC12 = 484.E-12* 120.9140/28.97
523  zzz = 484.
524  call getin('CFC12_ppt',zzz)
525  CFC12_ppt_omp = zzz
526  RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * 120.9140/28.97
527! RCFC12 = 2.020102726958923E-09
528!OK call getin('RCFC12', RCFC12)
529
530!ajout CFMIP begin
531!!
532!Config Key  = co2_ppm_per
533!Config Desc = concentration du co2_ppm_per
534!Config Def  = 348.
535!Config Help =
536!               
537  co2_ppm_per_omp = co2_ppm_omp
538  call getin('co2_ppm_per', co2_ppm_per_omp)
539!!
540!Config Key  = RCO2_per
541!Config Desc = Concentration du CO2_per
542!Config Def  = co2_ppm_per * 1.0e-06  * 44.011/28.97
543!Config Def  = 348. * 1.0e-06  * 44.011/28.97
544!Config Help =
545!               
546  RCO2_per_omp = co2_ppm_per_omp * 1.0e-06  * 44.011/28.97
547
548!Config Key  = ok_4xCO2atm
549!Config Desc = Calcul ou non effet radiatif 4xco2
550!Config Def  = .false.
551!Config Help =
552
553  ok_4xCO2atm_omp = .false.
554  call getin('ok_4xCO2atm',ok_4xCO2atm_omp)
555
556!Config Key  = RCH4_per
557!Config Desc = Concentration du CH4_per
558!Config Def  = 1.65E-06* 16.043/28.97
559!Config Help =
560!               
561  zzz = CH4_ppb_omp
562  call getin('CH4_ppb_per', zzz)
563  CH4_ppb_per_omp = zzz
564  RCH4_per_omp = CH4_ppb_per_omp * 1.0E-09 * 16.043/28.97
565!!
566!Config Key  = RN2O_per
567!Config Desc = Concentration du N2O_per
568!Config Def  = 306.E-09* 44.013/28.97
569!Config Help =
570!               
571  zzz = N2O_ppb_omp
572  call getin('N2O_ppb_per', zzz)
573  N2O_ppb_per_omp = zzz
574  RN2O_per_omp = N2O_ppb_per_omp * 1.0E-09 * 44.013/28.97
575!!
576!Config Key  = RCFC11_per
577!Config Desc = Concentration du CFC11_per
578!Config Def  = 280.E-12* 137.3686/28.97
579!Config Help =
580!               
581  zzz = CFC11_ppt_omp
582  call getin('CFC11_ppt_per',zzz)
583  CFC11_ppt_per_omp = zzz
584  RCFC11_per_omp=CFC11_ppt_per_omp* 1.0E-12 * 137.3686/28.97
585!!
586!Config Key  = RCFC12_per
587!Config Desc = Concentration du CFC12_per
588!Config Def  = 484.E-12* 120.9140/28.97
589!Config Help =
590!               
591  zzz = CFC12_ppt_omp
592  call getin('CFC12_ppt_per',zzz)
593  CFC12_ppt_per_omp = zzz
594  RCFC12_per_omp = CFC12_ppt_per_omp * 1.0E-12 * 120.9140/28.97
595!ajout CFMIP end
596
597!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
598! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
599! Constantes precedemment dans dyn3d/conf_gcm
600
601!Config  Key  = cycle_diurne
602!Config  Desc = Cycle ddiurne
603!Config  Def  = y
604!Config  Help = Cette option permet d'eteidre le cycle diurne.
605!Config         Peut etre util pour accelerer le code !
606       cycle_diurne_omp = .TRUE.
607       CALL getin('cycle_diurne',cycle_diurne_omp)
608
609!Config  Key  = soil_model
610!Config  Desc = Modele de sol
611!Config  Def  = y
612!Config  Help = Choix du modele de sol (Thermique ?)
613!Config         Option qui pourait un string afin de pouvoir
614!Config         plus de choix ! Ou meme une liste d'options !
615       soil_model_omp = .TRUE.
616       CALL getin('soil_model',soil_model_omp)
617
618!Config  Key  = new_oliq
619!Config  Desc = Nouvelle eau liquide
620!Config  Def  = y
621!Config  Help = Permet de mettre en route la
622!Config         nouvelle parametrisation de l'eau liquide !
623       new_oliq_omp = .TRUE.
624       CALL getin('new_oliq',new_oliq_omp)
625
626!Config  Key  = ok_orodr
627!Config  Desc = Orodr ???
628!Config  Def  = y
629!Config  Help = Y en a pas comprendre !
630!Config         
631       ok_orodr_omp = .TRUE.
632       CALL getin('ok_orodr',ok_orodr_omp)
633
634!Config  Key  =  ok_orolf
635!Config  Desc = Orolf ??
636!Config  Def  = y
637!Config  Help = Connais pas !
638       ok_orolf_omp = .TRUE.
639       CALL getin('ok_orolf', ok_orolf_omp)
640
641!Config  Key  = ok_limitvrai
642!Config  Desc = Force la lecture de la bonne annee
643!Config  Def  = n
644!Config  Help = On peut forcer le modele a lire le
645!Config         fichier SST de la bonne annee. C'est une tres bonne
646!Config         idee, pourquoi ne pas mettre toujours a y ???
647       ok_limitvrai_omp = .FALSE.
648       CALL getin('ok_limitvrai',ok_limitvrai_omp)
649
650!Config  Key  = nbapp_rad
651!Config  Desc = Frequence d'appel au rayonnement
652!Config  Def  = 12
653!Config  Help = Nombre  d'appels des routines de rayonnements
654!Config         par jour.
655       nbapp_rad_omp = 12
656       CALL getin('nbapp_rad',nbapp_rad_omp)
657
658!Config  Key  = iflag_con
659!Config  Desc = Flag de convection
660!Config  Def  = 2
661!Config  Help = Flag  pour la convection les options suivantes existent :
662!Config         1 pour LMD,
663!Config         2 pour Tiedtke,
664!Config         3 pour CCM(NCAR) 
665       iflag_con_omp = 2
666       CALL getin('iflag_con',iflag_con_omp)
667
668!Config  Key  = iflag_ener_conserv
669!Config  Desc = Flag de convection
670!Config  Def  = 1
671!Config  Help = Flag  pour la convection les options suivantes existent :
672!Config         -1 pour Kinetic energy correction
673!Config         1  conservation kinetic and enthalpy
674       iflag_ener_conserv_omp = -1
675       CALL getin('iflag_ener_conserv',iflag_ener_conserv_omp)
676
677!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
678!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
679!!
680!! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
681!!
682!! KE
683!
684
685!Config key  = cvl_corr
686!Config Desc = Facteur multiplication des precip convectives dans KE
687!Config Def  = 1.00
688!Config Help = 1.02 pour un moderne ou un pre-ind. A ajuster pour un glaciaire
689  cvl_corr_omp = 1.00
690  CALL getin('cvl_corr', cvl_corr_omp)
691
692
693!Config Key  = epmax
694!Config Desc = Efficacite precip
695!Config Def  = 0.993
696!Config Help =
697!
698  epmax_omp = .993
699  call getin('epmax', epmax_omp)
700!
701!Config Key  = ok_adj_ema
702!Config Desc = 
703!Config Def  = false
704!Config Help =
705!
706  ok_adj_ema_omp = .false.
707  call getin('ok_adj_ema',ok_adj_ema_omp)
708!
709!Config Key  = iflag_clw
710!Config Desc = 
711!Config Def  = 0
712!Config Help =
713!
714  iflag_clw_omp = 0
715  call getin('iflag_clw',iflag_clw_omp)
716!
717!Config Key  = cld_lc_lsc
718!Config Desc = 
719!Config Def  = 2.6e-4
720!Config Help =
721!
722  cld_lc_lsc_omp = 2.6e-4
723  call getin('cld_lc_lsc',cld_lc_lsc_omp)
724!
725!Config Key  = cld_lc_con
726!Config Desc = 
727!Config Def  = 2.6e-4
728!Config Help =
729!
730  cld_lc_con_omp = 2.6e-4
731  call getin('cld_lc_con',cld_lc_con_omp)
732!
733!Config Key  = cld_tau_lsc
734!Config Desc = 
735!Config Def  = 3600.
736!Config Help =
737!
738  cld_tau_lsc_omp = 3600.
739  call getin('cld_tau_lsc',cld_tau_lsc_omp)
740!
741!Config Key  = cld_tau_con
742!Config Desc = 
743!Config Def  = 3600.
744!Config Help =
745!
746  cld_tau_con_omp = 3600.
747  call getin('cld_tau_con',cld_tau_con_omp)
748!
749!Config Key  = ffallv_lsc
750!Config Desc = 
751!Config Def  = 1.
752!Config Help =
753!
754  ffallv_lsc_omp = 1.
755  call getin('ffallv_lsc',ffallv_lsc_omp)
756!
757!Config Key  = ffallv_con
758!Config Desc = 
759!Config Def  = 1.
760!Config Help =
761!
762  ffallv_con_omp = 1.
763  call getin('ffallv_con',ffallv_con_omp)
764!
765!Config Key  = coef_eva
766!Config Desc = 
767!Config Def  = 2.e-5
768!Config Help =
769!
770  coef_eva_omp = 2.e-5
771  call getin('coef_eva',coef_eva_omp)
772!
773!Config Key  = reevap_ice
774!Config Desc = 
775!Config Def  = .false.
776!Config Help =
777!
778  reevap_ice_omp = .false.
779  call getin('reevap_ice',reevap_ice_omp)
780
781!Config Key  = iflag_ratqs
782!Config Desc =
783!Config Def  = 1
784!Config Help =
785!
786  iflag_ratqs_omp = 1
787  call getin('iflag_ratqs',iflag_ratqs_omp)
788
789!
790!Config Key  = iflag_radia
791!Config Desc = 
792!Config Def  = 1
793!Config Help =
794!
795  iflag_radia_omp = 1
796  call getin('iflag_radia',iflag_radia_omp)
797
798!
799!Config Key  = iflag_rrtm
800!Config Desc = 
801!Config Def  = 0
802!Config Help =
803!
804  iflag_rrtm_omp = 0
805  call getin('iflag_rrtm',iflag_rrtm_omp)
806
807!
808!Config Key  = iflag_cldcon
809!Config Desc = 
810!Config Def  = 1
811!Config Help =
812!
813  iflag_cldcon_omp = 1
814  call getin('iflag_cldcon',iflag_cldcon_omp)
815
816!
817!Config Key  = iflag_pdf
818!Config Desc = 
819!Config Def  = 0
820!Config Help =
821!
822  iflag_pdf_omp = 0
823  call getin('iflag_pdf',iflag_pdf_omp)
824!
825!Config Key  = fact_cldcon
826!Config Desc = 
827!Config Def  = 0.375
828!Config Help =
829!
830  fact_cldcon_omp = 0.375
831  call getin('fact_cldcon',fact_cldcon_omp)
832
833!
834!Config Key  = facttemps
835!Config Desc = 
836!Config Def  = 1.e-4
837!Config Help =
838!
839  facttemps_omp = 1.e-4
840  call getin('facttemps',facttemps_omp)
841
842!
843!Config Key  = ok_newmicro
844!Config Desc = 
845!Config Def  = .true.
846!Config Help =
847!
848  ok_newmicro_omp = .true.
849  call getin('ok_newmicro',ok_newmicro_omp)
850!
851!Config Key  = ratqsbas
852!Config Desc = 
853!Config Def  = 0.01
854!Config Help =
855!
856  ratqsbas_omp = 0.01
857  call getin('ratqsbas',ratqsbas_omp)
858!
859!Config Key  = ratqshaut
860!Config Desc = 
861!Config Def  = 0.3
862!Config Help =
863!
864  ratqshaut_omp = 0.3
865  call getin('ratqshaut',ratqshaut_omp)
866
867!Config Key  = tau_ratqs
868!Config Desc = 
869!Config Def  = 1800.
870!Config Help =
871!
872  tau_ratqs_omp = 1800.
873  call getin('tau_ratqs',tau_ratqs_omp)
874
875!
876!-----------------------------------------------------------------------
877! Longitude solaire pour le calcul de l'ensoleillement en degre
878! si on veut imposer la saison. Sinon, solarlong0=-999.999
879!Config Key  = solarlong0
880!Config Desc = 
881!Config Def  = -999.999
882!Config Help =
883!
884  solarlong0_omp = -999.999
885  call getin('solarlong0',solarlong0_omp)
886!
887!-----------------------------------------------------------------------
888!  Valeur imposee de l'humidite du sol pour le modele bucket.
889!Config Key  = qsol0
890!Config Desc = 
891!Config Def  = -1.
892!Config Help =
893!
894  qsol0_omp = -1.
895  call getin('qsol0',qsol0_omp)
896!
897!-----------------------------------------------------------------------
898!
899!Config Key  = inertie_ice
900!Config Desc = 
901!Config Def  = 2000.
902!Config Help =
903!
904  inertie_ice_omp = 2000.
905  call getin('inertie_ice',inertie_ice_omp)
906!
907!Config Key  = inertie_sno
908!Config Desc = 
909!Config Def  = 2000.
910!Config Help =
911!
912  inertie_sno_omp = 2000.
913  call getin('inertie_sno',inertie_sno_omp)
914!
915!Config Key  = inertie_sol
916!Config Desc = 
917!Config Def  = 2000.
918!Config Help =
919!
920  inertie_sol_omp = 2000.
921  call getin('inertie_sol',inertie_sol_omp)
922
923!
924!Config Key  = rad_froid
925!Config Desc = 
926!Config Def  = 35.0
927!Config Help =
928!
929  rad_froid_omp = 35.0
930  call getin('rad_froid',rad_froid_omp)
931
932!
933!Config Key  = rad_chau1
934!Config Desc = 
935!Config Def  = 13.0
936!Config Help =
937!
938  rad_chau1_omp = 13.0
939  call getin('rad_chau1',rad_chau1_omp)
940
941!
942!Config Key  = rad_chau2
943!Config Desc = 
944!Config Def  = 9.0
945!Config Help =
946!
947  rad_chau2_omp = 9.0
948  call getin('rad_chau2',rad_chau2_omp)
949
950!
951!Config Key  = t_glace_min
952!Config Desc = 
953!Config Def  = 258.
954!Config Help =
955!
956  t_glace_min_omp = 258.
957  call getin('t_glace_min',t_glace_min_omp)
958
959!
960!Config Key  = t_glace_max
961!Config Desc = 
962!Config Def  = 273.13
963!Config Help =
964!
965  t_glace_max_omp = 273.13
966  call getin('t_glace_max',t_glace_max_omp)
967
968!Config Key  = rei_min
969!Config Desc = 
970!Config Def  = 3.5
971!Config Help =
972!
973  rei_min_omp = 3.5
974  call getin('rei_min',rei_min_omp)
975
976!
977!Config Key  = rei_max
978!Config Desc = 
979!Config Def  = 61.29
980!Config Help =
981!
982  rei_max_omp = 61.29
983  call getin('rei_max',rei_max_omp)
984
985!
986!Config Key  = top_height
987!Config Desc =
988!Config Def  = 3
989!Config Help =
990!
991  top_height_omp = 3
992  call getin('top_height',top_height_omp)
993
994!
995!Config Key  = overlap
996!Config Desc =
997!Config Def  = 3
998!Config Help =
999!
1000  overlap_omp = 3
1001  call getin('overlap',overlap_omp)
1002
1003
1004!
1005!
1006!Config Key  = cdmmax
1007!Config Desc =
1008!Config Def  = 1.3E-3
1009!Config Help =
1010!
1011  cdmmax_omp = 1.3E-3
1012  call getin('cdmmax',cdmmax_omp)
1013
1014!
1015!Config Key  = cdhmax
1016!Config Desc =
1017!Config Def  = 1.1E-3
1018!Config Help =
1019!
1020  cdhmax_omp = 1.1E-3
1021  call getin('cdhmax',cdhmax_omp)
1022
1023!261103
1024!
1025!Config Key  = ksta
1026!Config Desc =
1027!Config Def  = 1.0e-10
1028!Config Help =
1029!
1030  ksta_omp = 1.0e-10
1031  call getin('ksta',ksta_omp)
1032
1033!
1034!Config Key  = ksta_ter
1035!Config Desc =
1036!Config Def  = 1.0e-10
1037!Config Help =
1038!
1039  ksta_ter_omp = 1.0e-10
1040  call getin('ksta_ter',ksta_ter_omp)
1041
1042!
1043!Config Key  = ok_kzmin
1044!Config Desc =
1045!Config Def  = .true.
1046!Config Help =
1047!
1048  ok_kzmin_omp = .true.
1049  call getin('ok_kzmin',ok_kzmin_omp)
1050
1051!
1052!Config Key  = fmagic
1053!Config Desc = additionnal multiplicator factor used for albedo
1054!Config Def  = 1.
1055!Config Help = additionnal multiplicator factor used in albedo.F
1056!
1057  fmagic_omp = 1.
1058  call getin('fmagic',fmagic_omp)
1059
1060!
1061!Config Key  = pmagic
1062!Config Desc = additional factor used for albedo
1063!Config Def  = 0.
1064!Config Help = additional factor used in albedo.F
1065!
1066  pmagic_omp = 0.
1067  call getin('pmagic',pmagic_omp)
1068
1069
1070!Config Key = ok_lic_melt
1071!Config Desc = Prise en compte de la fonte de la calotte dans le bilan d'eau
1072!Config Def  = .false.
1073!Config Help = mettre a .false. pour assurer la conservation en eau
1074  ok_lic_melt_omp = .false.
1075  call getin('ok_lic_melt', ok_lic_melt_omp)
1076
1077!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1078! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
1079!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1080
1081!Config Key  = iflag_pbl
1082!Config Desc =
1083!Config Def  = 1
1084!Config Help =
1085!
1086  iflag_pbl_omp = 1
1087  call getin('iflag_pbl',iflag_pbl_omp)
1088!
1089!Config Key  = iflag_thermals
1090!Config Desc =
1091!Config Def  = 0
1092!Config Help =
1093!
1094  iflag_thermals_omp = 0
1095  call getin('iflag_thermals',iflag_thermals_omp)
1096!
1097!
1098!Config Key  = iflag_thermals_ed
1099!Config Desc =
1100!Config Def  = 0
1101!Config Help =
1102!
1103  iflag_thermals_ed_omp = 0
1104  call getin('iflag_thermals_ed',iflag_thermals_ed_omp)
1105!
1106!
1107!Config Key  = iflag_thermals_optflux
1108!Config Desc =
1109!Config Def  = 0
1110!Config Help =
1111!
1112  iflag_thermals_optflux_omp = 0
1113  call getin('iflag_thermals_optflux',iflag_thermals_optflux_omp)
1114!
1115!
1116!Config Key  = nsplit_thermals
1117!Config Desc =
1118!Config Def  = 1
1119!Config Help =
1120!
1121  nsplit_thermals_omp = 1
1122  call getin('nsplit_thermals',nsplit_thermals_omp)
1123
1124!Config Key  = alp_bl_k
1125!Config Desc =
1126!Config Def  = 0.
1127!Config Help =
1128!
1129  alp_bl_k_omp = 1.
1130  call getin('alp_bl_k',alp_bl_k_omp)
1131
1132!!! nrlmd le 10/04/2012
1133
1134!Config Key  = iflag_trig_bl
1135!Config Desc = 
1136!Config Def  = 0
1137!Config Help =
1138!
1139  iflag_trig_bl_omp = 0
1140  call getin('iflag_trig_bl',iflag_trig_bl_omp)
1141
1142!Config Key  = s_trig_bl
1143!Config Desc = 
1144!Config Def  = 0
1145!Config Help =
1146!
1147  s_trig_omp = 2e7
1148  call getin('s_trig',s_trig_omp)
1149
1150!Config Key  = tau_trig_shallow
1151!Config Desc = 
1152!Config Def  = 0
1153!Config Help =
1154!
1155  tau_trig_shallow_omp = 600
1156  call getin('tau_trig_shallow',tau_trig_shallow_omp)
1157
1158!Config Key  = tau_trig_deep
1159!Config Desc = 
1160!Config Def  = 0
1161!Config Help =
1162!
1163  tau_trig_deep_omp = 1800
1164  call getin('tau_trig_deep',tau_trig_deep_omp)
1165
1166!Config Key  = iflag_clos_bl
1167!Config Desc = 
1168!Config Def  = 0
1169!Config Help =
1170!
1171  iflag_clos_bl_omp = 0
1172  call getin('iflag_clos_bl',iflag_clos_bl_omp)
1173
1174!!! fin nrlmd le 10/04/2012
1175
1176!
1177!Config Key  = tau_thermals
1178!Config Desc =
1179!Config Def  = 0.
1180!Config Help =
1181!
1182  tau_thermals_omp = 0.
1183  call getin('tau_thermals',tau_thermals_omp)
1184
1185!
1186!Config Key  = iflag_coupl
1187!Config Desc =
1188!Config Def  = 0
1189!Config Help =
1190!
1191  iflag_coupl_omp = 0
1192  call getin('iflag_coupl',iflag_coupl_omp)
1193
1194!
1195!Config Key  = iflag_clos
1196!Config Desc = 
1197!Config Def  = 0
1198!Config Help =
1199!
1200  iflag_clos_omp = 1
1201  call getin('iflag_clos',iflag_clos_omp)
1202!
1203!Config Key  = iflag_cvl_sigd
1204!Config Desc = 
1205!Config Def  = 0
1206!Config Help =
1207!
1208  iflag_cvl_sigd_omp = 0
1209  call getin('iflag_cvl_sigd',iflag_cvl_sigd_omp)
1210
1211!Config Key  = iflag_wake
1212!Config Desc = 
1213!Config Def  = 0
1214!Config Help =
1215!
1216  iflag_wake_omp = 0
1217  call getin('iflag_wake',iflag_wake_omp)
1218
1219!Config Key  = alp_offset
1220!Config Desc = 
1221!Config Def  = 0
1222!Config Help =
1223!
1224  alp_offset_omp = 0.
1225  call getin('alp_offset',alp_offset_omp)
1226
1227!
1228!Config Key  = lev_histhf
1229!Config Desc =
1230!Config Def  = 1
1231!Config Help =
1232!
1233  lev_histhf_omp = 1
1234  call getin('lev_histhf',lev_histhf_omp)
1235
1236!
1237!Config Key  = lev_histday
1238!Config Desc =
1239!Config Def  = 1
1240!Config Help =
1241!
1242  lev_histday_omp = 1
1243  call getin('lev_histday',lev_histday_omp)
1244
1245!
1246!Config Key  = lev_histmth
1247!Config Desc =
1248!Config Def  = 2
1249!Config Help =
1250!
1251  lev_histmth_omp = 2
1252  call getin('lev_histmth',lev_histmth_omp)
1253!
1254!Config Key  = lev_histins
1255!Config Desc =
1256!Config Def  = 1
1257!Config Help =
1258!
1259  lev_histins_omp = 1
1260  call getin('lev_histins',lev_histins_omp)
1261  !
1262!Config Key  = lev_histLES
1263!Config Desc =
1264!Config Def  = 1
1265!Config Help =
1266!
1267  lev_histLES_omp = 1
1268  call getin('lev_histLES',lev_histLES_omp)
1269!
1270!Config Key  = lev_histdayNMC
1271!Config Desc =
1272!Config Def  = 8
1273!Config Help =
1274!
1275  lev_histdayNMC_omp = 8
1276  call getin('lev_histdayNMC',lev_histdayNMC_omp)
1277!
1278!histNMC BEG
1279!Config Key  = ok_histNMC
1280!Config Desc = ok_histNMC(1) = frequence de sortie fichiers histmthNMC
1281!Config Desc = ok_histNMC(2) = frequence de sortie fichiers histdayNMC
1282!Config Desc = ok_histNMC(3) = frequence de sortie fichiers histhfNMC
1283!Config Def  = n, n, n
1284!Config Help =
1285!
1286  ok_histNMC_omp(1) = .false.
1287  ok_histNMC_omp(2) = .false.
1288  ok_histNMC_omp(3) = .false.
1289  call getin('ok_histNMC',ok_histNMC_omp)
1290!
1291!Config Key  = freq_outNMC
1292!Config Desc = freq_outNMC(1) = frequence de sortie fichiers histmthNMC
1293!Config Desc = freq_outNMC(2) = frequence de sortie fichiers histdayNMC
1294!Config Desc = freq_outNMC(3) = frequence de sortie fichiers histhfNMC
1295!Config Def  = 2592000., 86400., 21600.
1296!Config Help =
1297!
1298! freq_outNMC_omp(1) = 2592000.
1299  freq_outNMC_omp(1) = mth_len*86400.
1300  freq_outNMC_omp(2) = 86400.
1301  freq_outNMC_omp(3) = 21600.
1302  call getin('freq_outNMC',freq_outNMC_omp)
1303!
1304!Config Key  = freq_calNMC
1305!Config Desc = freq_calNMC(1) = frequence de calcul fichiers histmthNMC
1306!Config Desc = freq_calNMC(2) = frequence de calcul fichiers histdayNMC
1307!Config Desc = freq_calNMC(3) = frequence de calcul fichiers histhfNMC
1308!Config Def  = pasphys
1309!Config Help =
1310!
1311  freq_calNMC_omp(1) = pasphys
1312  freq_calNMC_omp(2) = pasphys
1313  freq_calNMC_omp(3) = pasphys
1314  call getin('freq_calNMC',freq_calNMC_omp)
1315!
1316!Config Key  = type_run
1317!Config Desc =
1318!Config Def  = 'AMIP'/'CFMIP'  ou 'CLIM'/'ENSP'
1319!Config Help =
1320!
1321  type_run_omp = 'AMIP'
1322  call getin('type_run',type_run_omp)
1323
1324!
1325!Config Key  = ok_isccp
1326!Config Desc =
1327!Config Def  = .true.
1328!Config Help =
1329!
1330! ok_isccp = .true.
1331  ok_isccp_omp = .false.
1332  call getin('ok_isccp',ok_isccp_omp)
1333
1334!
1335!Config Key  = ok_cosp
1336!Config Desc =
1337!Config Def  = .false.
1338!Config Help =
1339!
1340  ok_cosp_omp = .false.
1341  call getin('ok_cosp',ok_cosp_omp)
1342
1343!
1344!Config Key  = ok_mensuelCOSP
1345!Config Desc =
1346!Config Def  = .true.
1347!Config Help =
1348!
1349  ok_mensuelCOSP_omp = .true.
1350  call getin('ok_mensuelCOSP',ok_mensuelCOSP_omp)
1351
1352!
1353!Config Key  = ok_journeCOSP
1354!Config Desc =
1355!Config Def  = .true.
1356!Config Help =
1357!
1358  ok_journeCOSP_omp = .true.
1359  call getin('ok_journeCOSP',ok_journeCOSP_omp)
1360
1361!
1362!Config Key  = ok_hfCOSP
1363!Config Desc =
1364!Config Def  = .false.
1365!Config Help =
1366!
1367  ok_hfCOSP_omp = .false.
1368  call getin('ok_hfCOSP',ok_hfCOSP_omp)
1369
1370!
1371! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone
1372! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
1373!
1374!Config Key  = lonmin_ins
1375!Config Desc = 100. 
1376!Config Def  = longitude minimale sorties "bilKP_ins"
1377!Config Help =
1378!
1379   lonmin_ins_omp = 100.
1380   call getin('lonmin_ins',lonmin_ins_omp)
1381!
1382!Config Key  = lonmax_ins
1383!Config Desc = 130.
1384!Config Def  = longitude maximale sorties "bilKP_ins"
1385!Config Help =
1386!
1387   lonmax_ins_omp = 130.
1388   call getin('lonmax_ins',lonmax_ins_omp)
1389!
1390!Config Key  = latmin_ins
1391!Config Desc = -20. 
1392!Config Def  = latitude minimale sorties "bilKP_ins"
1393!Config Help =
1394!
1395   latmin_ins_omp = -20.
1396   call getin('latmin_ins',latmin_ins_omp)
1397!
1398!Config Key  = latmax_ins
1399!Config Desc = 20.
1400!Config Def  = latitude maximale sorties "bilKP_ins"
1401!Config Help =
1402!
1403   latmax_ins_omp = 20.
1404   call getin('latmax_ins',latmax_ins_omp)
1405!
1406!Config Key  = ecrit_hf
1407!Config Desc =
1408!Config Def  = 1./8. !toutes les 3h
1409!Config Help =
1410!
1411  ecrit_hf_omp = 1./8.
1412  call getin('ecrit_hf',ecrit_hf_omp)
1413!
1414!Config Key  = ecrit_ins
1415!Config Desc =
1416!Config Def  = 1./48. ! toutes les 1/2 h
1417!Config Help =
1418!
1419  ecrit_ins_omp = 1./48.
1420  call getin('ecrit_ins',ecrit_ins_omp)
1421!
1422!Config Key  = ecrit_day
1423!Config Desc =
1424!Config Def  = 1.0 !tous les jours
1425!Config Help = nombre de jours pour ecriture fichier histday.nc
1426!
1427  ecrit_day_omp = 1.0
1428  call getin('ecrit_day',ecrit_day_omp)
1429!
1430!Config Key  = ecrit_mth
1431!Config Desc =
1432!Config Def  = 30. !tous les 30jours (1 fois par mois)
1433!Config Help =
1434!
1435  ecrit_mth_omp = 30.
1436  call getin('ecrit_mth',ecrit_mth_omp)
1437!
1438!Config Key  = ecrit_tra
1439!Config Desc =
1440!Config Def  = 30. !tous les 30jours (1 fois par mois)
1441!Config Help =
1442!
1443  ecrit_tra_omp = 0.
1444  call getin('ecrit_tra',ecrit_tra_omp)
1445!
1446!Config Key  = ecrit_reg
1447!Config Desc =
1448!Config Def  = 0.25  !4 fois par jour
1449!Config Help =
1450!
1451  ecrit_reg_omp = 0.25   !4 fois par jour
1452  call getin('ecrit_reg',ecrit_reg_omp)
1453!
1454!
1455!
1456! PARAMETRES CDRAG
1457!
1458!Config Key  = f_cdrag_ter
1459!Config Desc =
1460!Config Def  = 0.8
1461!Config Help =
1462!
1463  f_cdrag_ter_omp = 0.8
1464  call getin('f_cdrag_ter',f_cdrag_ter_omp)
1465!
1466!Config Key  = f_cdrag_oce
1467!Config Desc =
1468!Config Def  = 0.8
1469!Config Help =
1470!
1471  f_cdrag_oce_omp = 0.8
1472  call getin('f_cdrag_oce',f_cdrag_oce_omp)
1473!
1474! RUGORO
1475!Config Key  = f_rugoro
1476!Config Desc =
1477!Config Def  = 0.
1478!Config Help =
1479!
1480  f_rugoro_omp = 0.
1481  call getin('f_rugoro',f_rugoro_omp)
1482
1483! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
1484!
1485!Config Key  = supcrit1
1486!Config Desc =
1487!Config Def  = .540
1488!Config Help =
1489!
1490  supcrit1_omp = .540
1491  call getin('supcrit1',supcrit1_omp)
1492
1493!
1494!Config Key  = supcrit2
1495!Config Desc =
1496!Config Def  = .600
1497!Config Help =
1498!
1499  supcrit2_omp = .600
1500  call getin('supcrit2',supcrit2_omp)
1501
1502!
1503! PARAMETERS FOR THE MIXING DISTRIBUTION
1504! iflag_mix: 0=OLD,
1505!            1=NEW (JYG),           
1506!            2=NEW + conv. depth inhib. by tropos. dryness
1507! '2' is NOT operationnal and should not be used.
1508!
1509!Config Key  = iflag_mix
1510!Config Desc =
1511!Config Def  = 1
1512!Config Help =
1513!
1514  iflag_mix_omp = 1
1515  call getin('iflag_mix',iflag_mix_omp)
1516
1517!
1518!Config Key  = scut
1519!Config Desc =
1520!Config Def  = 0.95
1521!Config Help =
1522!
1523  scut_omp = 0.95
1524  call getin('scut',scut_omp)
1525
1526!
1527!Config Key  = qqa1
1528!Config Desc =
1529!Config Def  = 1.0
1530!Config Help =
1531!
1532  qqa1_omp = 1.0
1533  call getin('qqa1',qqa1_omp)
1534
1535!
1536!Config Key  = qqa2
1537!Config Desc =
1538!Config Def  = 0.0
1539!Config Help =
1540!
1541  qqa2_omp = 0.0
1542  call getin('qqa2',qqa2_omp)
1543
1544!
1545!Config Key  = gammas
1546!Config Desc =
1547!Config Def  = 0.05
1548!Config Help =
1549!
1550  gammas_omp = 0.05
1551  call getin('gammas',gammas_omp)
1552
1553!
1554!Config Key  = Fmax
1555!Config Desc =
1556!Config Def  = 0.65
1557!Config Help =
1558!
1559  Fmax_omp = 0.65
1560  call getin('Fmax',Fmax_omp)
1561
1562!
1563!Config Key  = alphas 
1564!Config Desc =
1565!Config Def  = -5.
1566!Config Help =
1567!
1568  alphas_omp = -5.
1569  call getin('alphas',alphas_omp)
1570
1571!Config key = ok_strato
1572!Config  Desc = activation de la version strato
1573!Config  Def  = .FALSE.
1574!Config  Help = active la version stratosphérique de LMDZ de F. Lott
1575
1576  ok_strato_omp=.FALSE.
1577  CALL getin('ok_strato',ok_strato_omp)
1578     
1579!Config  key = ok_hines
1580!Config  Desc = activation de la parametrisation de hines
1581!Config  Def  = .FALSE.
1582!Config  Help = Clefs controlant la parametrization de Hines
1583!               Et la sponge layer (Runs Stratospheriques)
1584
1585  ok_hines_omp=.FALSE.
1586  CALL getin('ok_hines',ok_hines_omp)
1587
1588!Config Key  = OK_LES                                               
1589!Config Desc = Pour des sorties LES                                 
1590!Config Def  = .false.                                             
1591!Config Help = Pour creer le fichier histLES contenant les sorties 
1592!              LES                                                 
1593!                                                                   
1594  ok_LES_omp = .false.                                             
1595  call getin('OK_LES', ok_LES_omp)                                 
1596
1597!Config Key  = callstats                                               
1598!Config Desc = Pour des sorties callstats                                 
1599!Config Def  = .false.                                             
1600!Config Help = Pour creer le fichier stats contenant les sorties 
1601!              stats                                                 
1602!                                                                   
1603  callstats_omp = .false.                                             
1604  call getin('callstats', callstats_omp)                                 
1605!
1606!Config Key  = ecrit_LES
1607!Config Desc = Frequence d'ecriture des resultats du LES en nombre de jours;
1608!              par defaut 1., i.e. 1 jour
1609!Config Def  = 1./8.
1610!Config Help = ...
1611!
1612!
1613  ecrit_LES_omp = 1./8.
1614  call getin('ecrit_LES', ecrit_LES_omp)
1615!
1616  read_climoz = 0 ! default value
1617  call getin('read_climoz', read_climoz)
1618
1619  carbon_cycle_tr_omp=.FALSE.
1620  CALL getin('carbon_cycle_tr',carbon_cycle_tr_omp)
1621
1622  carbon_cycle_cpl_omp=.FALSE.
1623  CALL getin('carbon_cycle_cpl',carbon_cycle_cpl_omp)
1624
1625!$OMP END MASTER
1626!$OMP BARRIER
1627
1628    R_ecc = R_ecc_omp
1629    R_peri = R_peri_omp
1630    R_incl = R_incl_omp
1631    solaire = solaire_omp
1632    co2_ppm = co2_ppm_omp
1633    RCO2 = RCO2_omp
1634    CH4_ppb = CH4_ppb_omp
1635    RCH4 = RCH4_omp
1636    N2O_ppb = N2O_ppb_omp
1637    RN2O = RN2O_omp
1638    CFC11_ppt = CFC11_ppt_omp
1639    RCFC11 = RCFC11_omp
1640    CFC12_ppt = CFC12_ppt_omp
1641    RCFC12 = RCFC12_omp
1642    RCO2_act = RCO2
1643    RCH4_act = RCH4
1644    RN2O_act = RN2O
1645    RCFC11_act = RCFC11
1646    RCFC12_act = RCFC12
1647    ok_4xCO2atm = ok_4xCO2atm_omp
1648    RCO2_per = RCO2_per_omp
1649    RCH4_per = RCH4_per_omp
1650    RN2O_per = RN2O_per_omp
1651    RCFC11_per = RCFC11_per_omp
1652    RCFC12_per = RCFC12_per_omp
1653   
1654    cycle_diurne = cycle_diurne_omp
1655    soil_model = soil_model_omp
1656    new_oliq = new_oliq_omp
1657    ok_orodr = ok_orodr_omp
1658    ok_orolf = ok_orolf_omp
1659    ok_limitvrai = ok_limitvrai_omp
1660    nbapp_rad = nbapp_rad_omp
1661    iflag_con = iflag_con_omp
1662    iflag_ener_conserv = iflag_ener_conserv_omp
1663
1664    epmax = epmax_omp
1665    ok_adj_ema = ok_adj_ema_omp
1666    iflag_clw = iflag_clw_omp
1667    cld_lc_lsc = cld_lc_lsc_omp
1668    cld_lc_con = cld_lc_con_omp
1669    cld_tau_lsc = cld_tau_lsc_omp
1670    cld_tau_con = cld_tau_con_omp
1671    ffallv_lsc = ffallv_lsc_omp
1672    ffallv_con = ffallv_con_omp
1673    coef_eva = coef_eva_omp
1674    reevap_ice = reevap_ice_omp
1675    iflag_pdf = iflag_pdf_omp
1676    solarlong0 = solarlong0_omp
1677    qsol0 = qsol0_omp
1678    inertie_sol = inertie_sol_omp
1679    inertie_ice = inertie_ice_omp
1680    inertie_sno = inertie_sno_omp
1681    rad_froid = rad_froid_omp
1682    rad_chau1 = rad_chau1_omp
1683    rad_chau2 = rad_chau2_omp
1684    t_glace_min = t_glace_min_omp
1685    t_glace_max = t_glace_max_omp
1686    rei_min = rei_min_omp
1687    rei_max = rei_max_omp
1688    top_height = top_height_omp
1689    overlap = overlap_omp
1690    cdmmax = cdmmax_omp
1691    cdhmax = cdhmax_omp
1692    ksta = ksta_omp
1693    ksta_ter = ksta_ter_omp
1694    ok_kzmin = ok_kzmin_omp
1695    fmagic = fmagic_omp
1696    pmagic = pmagic_omp
1697    iflag_pbl = iflag_pbl_omp
1698    lev_histhf = lev_histhf_omp
1699    lev_histday = lev_histday_omp
1700    lev_histmth = lev_histmth_omp
1701    lev_histins = lev_histins_omp
1702    lev_histLES = lev_histLES_omp
1703    lev_histdayNMC = lev_histdayNMC_omp
1704    ok_histNMC(:) = ok_histNMC_omp(:)
1705    freq_outNMC(:) = freq_outNMC_omp(:)
1706    freq_calNMC(:) = freq_calNMC_omp(:)
1707
1708    type_ocean = type_ocean_omp
1709    version_ocean = version_ocean_omp
1710
1711    ok_veget=.true.
1712    type_veget=type_veget_omp
1713    if (type_veget=='n' .or. type_veget=='bucket' .or. type_veget=='betaclim' ) &
1714         & then
1715         ok_veget=.false.
1716    endif
1717
1718    ok_newmicro = ok_newmicro_omp
1719    ok_journe = ok_journe_omp
1720    ok_hf = ok_hf_omp
1721    ok_mensuel = ok_mensuel_omp
1722    ok_instan = ok_instan_omp
1723    freq_ISCCP = freq_ISCCP_omp
1724    ecrit_ISCCP = ecrit_ISCCP_omp
1725    freq_COSP = freq_COSP_omp
1726    ok_ade = ok_ade_omp
1727    ok_aie = ok_aie_omp
1728    ok_cdnc = ok_cdnc_omp
1729    aerosol_couple = aerosol_couple_omp
1730    flag_aerosol=flag_aerosol_omp
1731    flag_aerosol_strat=flag_aerosol_strat_omp
1732    new_aod=new_aod_omp
1733    aer_type = aer_type_omp
1734    bl95_b0 = bl95_b0_omp
1735    bl95_b1 = bl95_b1_omp
1736    fact_cldcon = fact_cldcon_omp
1737    facttemps = facttemps_omp
1738    ratqsbas = ratqsbas_omp
1739    ratqshaut = ratqshaut_omp
1740    tau_ratqs = tau_ratqs_omp
1741
1742    iflag_radia = iflag_radia_omp
1743    iflag_rrtm = iflag_rrtm_omp
1744    iflag_cldcon = iflag_cldcon_omp
1745    iflag_ratqs = iflag_ratqs_omp
1746    ip_ebil_phy = ip_ebil_phy_omp
1747    iflag_thermals = iflag_thermals_omp
1748    iflag_thermals_ed = iflag_thermals_ed_omp
1749    iflag_thermals_optflux = iflag_thermals_optflux_omp
1750    nsplit_thermals = nsplit_thermals_omp
1751    tau_thermals = tau_thermals_omp
1752    alp_bl_k = alp_bl_k_omp
1753!!! nrlmd le 10/04/2012
1754    iflag_trig_bl = iflag_trig_bl_omp
1755    s_trig = s_trig_omp
1756    tau_trig_shallow = tau_trig_shallow_omp
1757    tau_trig_deep = tau_trig_deep_omp
1758    iflag_clos_bl = iflag_clos_bl_omp
1759!!! fin nrlmd le 10/04/2012
1760    iflag_coupl = iflag_coupl_omp
1761    iflag_clos = iflag_clos_omp
1762    iflag_wake = iflag_wake_omp
1763    alp_offset = alp_offset_omp
1764    iflag_cvl_sigd = iflag_cvl_sigd_omp
1765    type_run = type_run_omp
1766    ok_isccp = ok_isccp_omp
1767    ok_cosp = ok_cosp_omp
1768    ok_mensuelCOSP = ok_mensuelCOSP_omp
1769    ok_journeCOSP = ok_journeCOSP_omp
1770    ok_hfCOSP = ok_hfCOSP_omp
1771    seuil_inversion=seuil_inversion_omp
1772    lonmin_ins = lonmin_ins_omp
1773    lonmax_ins = lonmax_ins_omp
1774    latmin_ins = latmin_ins_omp
1775    latmax_ins = latmax_ins_omp
1776    ecrit_hf   = ecrit_hf_omp
1777    ecrit_ins   = ecrit_ins_omp
1778    ecrit_day = ecrit_day_omp
1779    ecrit_mth = ecrit_mth_omp
1780    ecrit_tra = ecrit_tra_omp
1781    ecrit_reg = ecrit_reg_omp
1782    cvl_corr = cvl_corr_omp
1783    ok_lic_melt = ok_lic_melt_omp
1784    f_cdrag_ter=f_cdrag_ter_omp
1785    f_cdrag_oce=f_cdrag_oce_omp
1786    f_rugoro=f_rugoro_omp
1787    supcrit1 = supcrit1_omp
1788    supcrit2 = supcrit2_omp
1789    iflag_mix = iflag_mix_omp
1790    scut = scut_omp
1791    qqa1 = qqa1_omp
1792    qqa2 = qqa2_omp
1793    gammas = gammas_omp
1794    Fmax = Fmax_omp
1795    alphas = alphas_omp
1796    ok_strato = ok_strato_omp
1797    ok_hines = ok_hines_omp
1798    ok_LES = ok_LES_omp
1799    callstats = callstats_omp
1800    ecrit_LES = ecrit_LES_omp
1801    carbon_cycle_tr = carbon_cycle_tr_omp
1802    carbon_cycle_cpl = carbon_cycle_cpl_omp
1803
1804! Test of coherence between type_ocean and version_ocean
1805    IF (type_ocean=='couple' .AND. (version_ocean/='opa8' .AND. version_ocean/='nemo') ) THEN
1806       WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid in coupled configuration'
1807       CALL abort_gcm('conf_phys','version_ocean not valid',1)
1808    END IF
1809
1810    IF (type_ocean=='slab' .AND. version_ocean=='xxxxxx') THEN
1811       version_ocean='sicOBS'
1812    ELSE IF (type_ocean=='slab' .AND. version_ocean/='sicOBS') THEN
1813       WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean'
1814       CALL abort_gcm('conf_phys','version_ocean not valid',1)
1815    END IF
1816
1817! Test sur new_aod. Ce flag permet de retrouver les resultats de l'AR4
1818! il n'est utilisable que lors du couplage avec le SO4 seul
1819    IF (ok_ade .OR. ok_aie) THEN
1820       IF ( flag_aerosol .EQ. 0 ) THEN
1821          CALL abort_gcm('conf_phys','flag_aerosol=0 not compatible avec ok_ade ou ok_aie=.TRUE.',1)
1822       END IF
1823       IF ( .NOT. new_aod .AND.  flag_aerosol .NE. 1) THEN
1824          CALL abort_gcm('conf_phys','new_aod=.FALSE. not compatible avec flag_aerosol=1',1)
1825       END IF
1826    END IF
1827
1828! ok_cdnc must be set to y if ok_aie is activated
1829    IF (ok_aie .AND. .NOT. ok_cdnc) THEN
1830       CALL abort_gcm('conf_phys', 'ok_cdnc must be set to y if ok_aie is activated',1)
1831    ENDIF
1832
1833!$OMP MASTER
1834
1835  write(lunout,*)' ##############################################'
1836  write(lunout,*)' Configuration des parametres de la physique: '
1837  write(lunout,*)' Type ocean = ', type_ocean
1838  write(lunout,*)' Version ocean = ', version_ocean
1839  write(lunout,*)' Config veget = ', ok_veget,type_veget
1840  write(lunout,*)' Sortie journaliere = ', ok_journe
1841  write(lunout,*)' Sortie haute frequence = ', ok_hf
1842  write(lunout,*)' Sortie mensuelle = ', ok_mensuel
1843  write(lunout,*)' Sortie instantanee = ', ok_instan
1844  write(lunout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
1845  write(lunout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
1846  write(lunout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
1847  write(lunout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
1848  write(lunout,*)' Excentricite = ',R_ecc
1849  write(lunout,*)' Equinoxe = ',R_peri
1850  write(lunout,*)' Inclinaison =',R_incl
1851  write(lunout,*)' Constante solaire =',solaire
1852  write(lunout,*)' co2_ppm =',co2_ppm
1853  write(lunout,*)' RCO2_act = ',RCO2_act
1854  write(lunout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
1855  write(lunout,*)' N2O_ppb =',N2O_ppb,' RN2O_act=  ',RN2O_act
1856  write(lunout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act=  ',RCFC11_act
1857  write(lunout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act=  ',RCFC12_act
1858  write(lunout,*)' ok_4xCO2atm=',ok_4xCO2atm
1859  write(lunout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
1860  write(lunout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
1861  write(lunout,*)' RCFC12_per = ',RCFC12_per
1862  write(lunout,*)' cvl_corr=', cvl_corr
1863  write(lunout,*)'ok_lic_melt=', ok_lic_melt
1864  write(lunout,*)'cycle_diurne=',cycle_diurne
1865  write(lunout,*)'soil_model=',soil_model
1866  write(lunout,*)'new_oliq=',new_oliq
1867  write(lunout,*)'ok_orodr=',ok_orodr
1868  write(lunout,*)'ok_orolf=',ok_orolf
1869  write(lunout,*)'ok_limitvrai=',ok_limitvrai
1870  write(lunout,*)'nbapp_rad=',nbapp_rad
1871  write(lunout,*)'iflag_con=',iflag_con
1872  write(lunout,*)'iflag_ener_conserv=',iflag_ener_conserv
1873  write(lunout,*)' epmax = ', epmax
1874  write(lunout,*)' ok_adj_ema = ', ok_adj_ema
1875  write(lunout,*)' iflag_clw = ', iflag_clw
1876  write(lunout,*)' cld_lc_lsc = ', cld_lc_lsc
1877  write(lunout,*)' cld_lc_con = ', cld_lc_con
1878  write(lunout,*)' cld_tau_lsc = ', cld_tau_lsc
1879  write(lunout,*)' cld_tau_con = ', cld_tau_con
1880  write(lunout,*)' ffallv_lsc = ', ffallv_lsc
1881  write(lunout,*)' ffallv_con = ', ffallv_con
1882  write(lunout,*)' coef_eva = ', coef_eva
1883  write(lunout,*)' reevap_ice = ', reevap_ice
1884  write(lunout,*)' iflag_pdf = ', iflag_pdf
1885  write(lunout,*)' iflag_cldcon = ', iflag_cldcon
1886  write(lunout,*)' iflag_radia = ', iflag_radia
1887  write(lunout,*)' iflag_rrtm = ', iflag_rrtm
1888  write(lunout,*)' iflag_ratqs = ', iflag_ratqs
1889  write(lunout,*)' seuil_inversion = ', seuil_inversion
1890  write(lunout,*)' fact_cldcon = ', fact_cldcon
1891  write(lunout,*)' facttemps = ', facttemps
1892  write(lunout,*)' ok_newmicro = ',ok_newmicro
1893  write(lunout,*)' ratqsbas = ',ratqsbas
1894  write(lunout,*)' ratqshaut = ',ratqshaut
1895  write(lunout,*)' tau_ratqs = ',tau_ratqs
1896  write(lunout,*)' top_height = ',top_height
1897  write(lunout,*)' rad_froid = ',rad_froid
1898  write(lunout,*)' rad_chau1 = ',rad_chau1
1899  write(lunout,*)' rad_chau2 = ',rad_chau2
1900  write(lunout,*)' t_glace_min = ',t_glace_min
1901  write(lunout,*)' t_glace_max = ',t_glace_max
1902  write(lunout,*)' rei_min = ',rei_min
1903  write(lunout,*)' rei_max = ',rei_max
1904  write(lunout,*)' overlap = ',overlap
1905  write(lunout,*)' cdmmax = ',cdmmax
1906  write(lunout,*)' cdhmax = ',cdhmax
1907  write(lunout,*)' ksta = ',ksta
1908  write(lunout,*)' ksta_ter = ',ksta_ter
1909  write(lunout,*)' ok_kzmin = ',ok_kzmin
1910  write(lunout,*)' fmagic = ',fmagic
1911  write(lunout,*)' pmagic = ',pmagic
1912  write(lunout,*)' ok_ade = ',ok_ade
1913  write(lunout,*)' ok_aie = ',ok_aie
1914  write(lunout,*)' aerosol_couple = ', aerosol_couple
1915  write(lunout,*)' flag_aerosol = ', flag_aerosol
1916  write(lunout,*)' flag_aerosol_strat = ', flag_aerosol_strat
1917  write(lunout,*)' new_aod = ', new_aod
1918  write(lunout,*)' aer_type = ',aer_type
1919  write(lunout,*)' bl95_b0 = ',bl95_b0
1920  write(lunout,*)' bl95_b1 = ',bl95_b1
1921  write(lunout,*)' lev_histhf = ',lev_histhf
1922  write(lunout,*)' lev_histday = ',lev_histday
1923  write(lunout,*)' lev_histmth = ',lev_histmth
1924  write(lunout,*)' lev_histins = ',lev_histins
1925  write(lunout,*)' lev_histLES = ',lev_histLES
1926  write(lunout,*)' lev_histdayNMC = ',lev_histdayNMC
1927  write(lunout,*)' ok_histNMC = ',ok_histNMC
1928  write(lunout,*)' freq_outNMC = ',freq_outNMC
1929  write(lunout,*)' freq_calNMC = ',freq_calNMC
1930  write(lunout,*)' iflag_pbl = ', iflag_pbl
1931  write(lunout,*)' iflag_thermals = ', iflag_thermals
1932  write(lunout,*)' iflag_thermals_ed = ', iflag_thermals_ed
1933  write(lunout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux
1934  write(lunout,*)' iflag_clos = ', iflag_clos
1935  write(lunout,*)' type_run = ',type_run
1936  write(lunout,*)' ok_isccp = ',ok_isccp
1937  write(lunout,*)' ok_cosp = ',ok_cosp
1938  write(lunout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP
1939  write(lunout,*)' ok_journeCOSP = ',ok_journeCOSP
1940  write(lunout,*)' ok_hfCOSP =',ok_hfCOSP
1941  write(lunout,*)' solarlong0 = ', solarlong0
1942  write(lunout,*)' qsol0 = ', qsol0
1943  write(lunout,*)' inertie_sol = ', inertie_sol
1944  write(lunout,*)' inertie_ice = ', inertie_ice
1945  write(lunout,*)' inertie_sno = ', inertie_sno
1946  write(lunout,*)' f_cdrag_ter = ',f_cdrag_ter
1947  write(lunout,*)' f_cdrag_oce = ',f_cdrag_oce
1948  write(lunout,*)' f_rugoro = ',f_rugoro
1949  write(lunout,*)' supcrit1 = ', supcrit1
1950  write(lunout,*)' supcrit2 = ', supcrit2
1951  write(lunout,*)' iflag_mix = ', iflag_mix
1952  write(lunout,*)' scut = ', scut
1953  write(lunout,*)' qqa1 = ', qqa1
1954  write(lunout,*)' qqa2 = ', qqa2
1955  write(lunout,*)' gammas = ', gammas
1956  write(lunout,*)' Fmax = ', Fmax
1957  write(lunout,*)' alphas = ', alphas
1958  write(lunout,*)' iflag_wake = ', iflag_wake
1959  write(lunout,*)' alp_offset = ', alp_offset
1960!!! nrlmd le 10/04/2012
1961  write(lunout,*)' iflag_trig_bl = ', iflag_trig_bl
1962  write(lunout,*)' s_trig = ', s_trig
1963  write(lunout,*)' tau_trig_shallow = ', tau_trig_shallow
1964  write(lunout,*)' tau_trig_deep = ', tau_trig_deep
1965  write(lunout,*)' iflag_clos_bl = ', iflag_clos_bl
1966!!! fin nrlmd le 10/04/2012
1967
1968  write(lunout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
1969   lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
1970  write(lunout,*)' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
1971   ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
1972
1973  write(lunout,*) 'ok_strato = ', ok_strato
1974  write(lunout,*) 'ok_hines = ',  ok_hines
1975  write(lunout,*) 'read_climoz = ', read_climoz
1976  write(lunout,*) 'carbon_cycle_tr = ', carbon_cycle_tr
1977  write(lunout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl
1978 
1979!$OMP END MASTER
1980
1981  return
1982 
1983  end subroutine conf_phys
1984
1985end module conf_phys_m
1986!
1987!#################################################################
1988!
1989
1990   subroutine conf_interface(tau_calv)
1991
1992   use IOIPSL
1993   implicit none
1994   include "iniprint.h"
1995! Configuration de l'interace atm/surf
1996!
1997! tau_calv:    temps de relaxation pour la fonte des glaciers
1998
1999  REAL          :: tau_calv
2000  REAL,SAVE     :: tau_calv_omp
2001
2002!
2003!Config Key  = tau_calv
2004!Config Desc = temps de relaxation pour fonte des glaciers en jours
2005!Config Def  = 1 an
2006!Config Help =
2007!
2008  tau_calv_omp = 360.*10.
2009!$OMP MASTER
2010  call getin('tau_calv',tau_calv_omp)
2011!$OMP END MASTER
2012!$OMP BARRIER
2013
2014  tau_calv=tau_calv_omp
2015 
2016!$OMP MASTER
2017  write(lunout,*)' ##############################################'
2018  WRITE(lunout,*)' Configuration de l''interface atm/surfaces  : '
2019  WRITE(lunout,*)' tau_calv = ',tau_calv
2020!$OMP END MASTER
2021
2022  return
2023
2024  end subroutine conf_interface
Note: See TracBrowser for help on using the repository browser.