source: LMDZ5/trunk/libf/phylmd/conf_phys_m.F90 @ 1724

Last change on this file since 1724 was 1724, checked in by musat, 11 years ago

Ajout coefficient beta (=Evap/Evap? potentielle) fonction de la latitude
et de la saison (voir these AbdulKhadre? Traore).

Frederic Hourdin
IM

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