source: lmdz_wrf/trunk/WRFV3/lmdz/conf_phys_m.F90 @ 1939

Last change on this file since 1939 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

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.