source: LMDZ4/trunk/libf/phylmd/conf_phys.F90 @ 682

Last change on this file since 682 was 659, checked in by lmdzadmin, 19 years ago

Par defaut ecrit_tra=ecritphy
MAF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 17.7 KB
Line 
1!
2! $Header$
3!
4!
5!
6
7  subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
8 &                     fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &
9 &                     ratqsbas,ratqshaut,if_ebil, &
10 &                     ok_ade, ok_aie, &
11 &                     bl95_b0, bl95_b1,&
12 &                     iflag_thermals,nsplit_thermals)
13
14   use IOIPSL
15   implicit none
16
17#include "conema3.h"
18#include "fisrtilp.inc"
19#include "nuage.h"
20#include "YOMCST.inc"
21!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
22include "clesphys.inc"
23include "compbl.h"
24#include "control.inc"
25!
26! Configuration de la "physique" de LMDZ a l'aide de la fonction
27! GETIN de IOIPSL
28!
29! LF 05/2001
30!
31
32!
33! ocean:      type d'ocean (force, slab, couple)
34! ok_veget:   type de modele de vegetation
35! ok_journe:  sorties journalieres
36! ok_mensuel: sorties mensuelles
37! ok_instan:  sorties instantanees
38! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
39! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
40!
41
42
43! Sortie:
44  character (len = 6)  :: ocean
45  logical              :: ok_veget, ok_newmicro
46  logical              :: ok_journe, ok_mensuel, ok_instan       
47  LOGICAL              :: ok_ade, ok_aie
48  REAL                 :: bl95_b0, bl95_b1
49  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
50  integer              :: iflag_cldcon, if_ebil
51
52! Local
53  integer              :: numout = 6
54  real                 :: zzz
55
56  integer :: iflag_thermals,nsplit_thermals
57!
58!
59!
60
61
62!Config Key  = OCEAN
63!Config Desc = Type d'ocean
64!Config Def  = force
65!Config Help = Type d'ocean utilise: force, slab,couple
66!
67  ocean = 'force '
68  call getin('OCEAN', ocean)
69!
70!Config Key  = VEGET
71!Config Desc = Type de modele de vegetation
72!Config Def  = .false.
73!Config Help = Type de modele de vegetation utilise
74!
75  ok_veget = .false.
76  call getin('VEGET', ok_veget)
77!
78!Config Key  = OK_journe
79!Config Desc = Pour des sorties journalieres
80!Config Def  = .false.
81!Config Help = Pour creer le fichier histday contenant les sorties
82!              journalieres
83!
84  ok_journe = .false.
85  call getin('OK_journe', ok_journe)
86!
87!Config Key  = OK_mensuel
88!Config Desc = Pour des sorties mensuelles
89!Config Def  = .true.
90!Config Help = Pour creer le fichier histmth contenant les sorties
91!              mensuelles
92!
93  ok_mensuel = .true.
94  call getin('OK_mensuel', ok_mensuel)
95!
96!Config Key  = OK_instan
97!Config Desc = Pour des sorties instantanees
98!Config Def  = .false.
99!Config Help = Pour creer le fichier histins contenant les sorties
100!              instantanees
101!
102  ok_instan = .false.
103  call getin('OK_instan', ok_instan)
104!
105!Config Key  = ok_ade
106!Config Desc = Aerosol direct effect or not?
107!Config Def  = .false.
108!Config Help = Used in radlwsw.F
109!
110  ok_ade = .false.
111  call getin('ok_ade', ok_ade)
112
113!
114!Config Key  = ok_aie
115!Config Desc = Aerosol indirect effect or not?
116!Config Def  = .false.
117!Config Help = Used in nuage.F and radlwsw.F
118!
119  ok_aie = .false.
120  call getin('ok_aie', ok_aie)
121
122!
123!Config Key  = bl95_b0
124!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
125!Config Def  = .false.
126!Config Help = Used in nuage.F
127!
128  bl95_b0 = 2.
129  call getin('bl95_b0', bl95_b0)
130
131!Config Key  = bl95_b1
132!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
133!Config Def  = .false.
134!Config Help = Used in nuage.F
135!
136  bl95_b1 = 0.2
137  call getin('bl95_b1', bl95_b1)
138
139!
140!
141!Config Key  = if_ebil
142!Config Desc = Niveau de sortie pour les diags bilan d'energie
143!Config Def  = 0
144!Config Help =
145!               
146!
147  if_ebil = 0
148  call getin('if_ebil', if_ebil)
149!!
150!! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
151!!
152!Config Key  = R_ecc
153!Config Desc = Excentricite
154!Config Def  = 0.016715
155!Config Help =
156!               
157!valeur AMIP II
158  R_ecc = 0.016715
159  call getin('R_ecc', R_ecc)
160!!
161!Config Key  = R_peri
162!Config Desc = Equinoxe
163!Config Def  =
164!Config Help =
165!               
166!
167!valeur AMIP II
168  R_peri = 102.7
169  call getin('R_peri', R_peri)
170!!
171!Config Key  = R_incl
172!Config Desc = Inclinaison
173!Config Def  =
174!Config Help =
175!               
176!
177!valeur AMIP II
178  R_incl = 23.441
179  call getin('R_incl', R_incl)
180!!
181!Config Key  = solaire
182!Config Desc = Constante solaire en W/m2
183!Config Def  = 1365.
184!Config Help =
185!               
186!
187!valeur AMIP II
188  solaire = 1365.
189  call getin('solaire', solaire)
190!!
191!Config Key  = co2_ppm
192!Config Desc = concentration du gaz carbonique en ppmv
193!Config Def  = 348.
194!Config Help =
195!               
196!
197!valeur AMIP II
198  co2_ppm = 348.
199  call getin('co2_ppm', co2_ppm)
200!!
201!Config Key  = RCO2
202!Config Desc = Concentration du CO2
203!Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
204!Config Def  = 348. * 1.0e-06  * 44.011/28.97
205!Config Help =
206!               
207! RCO2 = 5.286789092164308E-04
208!ancienne valeur
209  RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.
210
211!!  call getin('RCO2', RCO2)
212!!
213!Config Key  = RCH4
214!Config Desc = Concentration du CH4
215!Config Def  = 1.65E-06* 16.043/28.97
216!Config Help =
217!               
218!
219!valeur AMIP II
220!OK  RCH4 = 1.65E-06* 16.043/28.97
221! RCH4 = 9.137366240938903E-07
222!
223!ancienne valeur
224! RCH4 = 1.72E-06* 16.043/28.97
225!OK call getin('RCH4', RCH4)
226  zzz = 1650.
227  call getin('CH4_ppb', zzz)
228  CH4_ppb = zzz
229  RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
230!!
231!Config Key  = RN2O
232!Config Desc = Concentration du N2O
233!Config Def  = 306.E-09* 44.013/28.97
234!Config Help =
235!               
236!
237!valeur AMIP II
238!OK  RN2O = 306.E-09* 44.013/28.97
239! RN2O = 4.648939592682085E-07
240!
241!ancienne valeur
242! RN2O = 310.E-09* 44.013/28.97
243!OK  call getin('RN2O', RN2O)
244  zzz=306.
245  call getin('N2O_ppb', zzz)
246  N2O_ppb = zzz
247  RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
248!!
249!Config Key  = RCFC11
250!Config Desc = Concentration du CFC11
251!Config Def  = 280.E-12* 137.3686/28.97
252!Config Help =
253!               
254!
255!OK RCFC11 = 280.E-12* 137.3686/28.97
256  zzz = 280.
257  call getin('CFC11_ppt',zzz)
258  CFC11_ppt = zzz
259  RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
260! RCFC11 = 1.327690990680013E-09
261!OK call getin('RCFC11', RCFC11)
262!!
263!Config Key  = RCFC12
264!Config Desc = Concentration du CFC12
265!Config Def  = 484.E-12* 120.9140/28.97
266!Config Help =
267!               
268!
269!OK RCFC12 = 484.E-12* 120.9140/28.97
270  zzz = 484.
271  call getin('CFC12_ppt',zzz)
272  CFC12_ppt = zzz
273  RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
274! RCFC12 = 2.020102726958923E-09
275!OK call getin('RCFC12', RCFC12)
276!!
277!! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
278!!
279!! KE
280!
281!Config Key  = epmax
282!Config Desc = Efficacite precip
283!Config Def  = 0.993
284!Config Help =
285!
286  epmax = .993
287  call getin('epmax', epmax)
288!
289!Config Key  = ok_adj_ema
290!Config Desc = 
291!Config Def  = false
292!Config Help =
293!
294  ok_adj_ema = .false.
295  call getin('ok_adj_ema',ok_adj_ema)
296!
297!Config Key  = iflag_clw
298!Config Desc = 
299!Config Def  = 0
300!Config Help =
301!
302  iflag_clw = 0
303  call getin('iflag_clw',iflag_clw)
304!
305!Config Key  = cld_lc_lsc
306!Config Desc = 
307!Config Def  = 2.6e-4
308!Config Help =
309!
310  cld_lc_lsc = 2.6e-4
311  call getin('cld_lc_lsc',cld_lc_lsc)
312!
313!Config Key  = cld_lc_con
314!Config Desc = 
315!Config Def  = 2.6e-4
316!Config Help =
317!
318  cld_lc_con = 2.6e-4
319  call getin('cld_lc_con',cld_lc_con)
320!
321!Config Key  = cld_tau_lsc
322!Config Desc = 
323!Config Def  = 3600.
324!Config Help =
325!
326  cld_tau_lsc = 3600.
327  call getin('cld_tau_lsc',cld_tau_lsc)
328!
329!Config Key  = cld_tau_con
330!Config Desc = 
331!Config Def  = 3600.
332!Config Help =
333!
334  cld_tau_con = 3600.
335  call getin('cld_tau_con',cld_tau_con)
336!
337!Config Key  = ffallv_lsc
338!Config Desc = 
339!Config Def  = 1.
340!Config Help =
341!
342  ffallv_lsc = 1.
343  call getin('ffallv_lsc',ffallv_lsc)
344!
345!Config Key  = ffallv_con
346!Config Desc = 
347!Config Def  = 1.
348!Config Help =
349!
350  ffallv_con = 1.
351  call getin('ffallv_con',ffallv_con)
352!
353!Config Key  = coef_eva
354!Config Desc = 
355!Config Def  = 2.e-5
356!Config Help =
357!
358  coef_eva = 2.e-5
359  call getin('coef_eva',coef_eva)
360!
361!Config Key  = reevap_ice
362!Config Desc = 
363!Config Def  = .false.
364!Config Help =
365!
366  reevap_ice = .false.
367  call getin('reevap_ice',reevap_ice)
368!
369!Config Key  = iflag_cldcon
370!Config Desc = 
371!Config Def  = 1
372!Config Help =
373!
374  iflag_cldcon = 1
375  call getin('iflag_cldcon',iflag_cldcon)
376
377!
378!Config Key  = iflag_pdf
379!Config Desc = 
380!Config Def  = 0
381!Config Help =
382!
383  iflag_pdf = 0
384  call getin('iflag_pdf',iflag_pdf)
385!
386!Config Key  = fact_cldcon
387!Config Desc = 
388!Config Def  = 0.375
389!Config Help =
390!
391  fact_cldcon = 0.375
392  call getin('fact_cldcon',fact_cldcon)
393
394!
395!Config Key  = facttemps
396!Config Desc = 
397!Config Def  = 1.e-4
398!Config Help =
399!
400  facttemps = 1.e-4
401  call getin('facttemps',facttemps)
402
403!
404!Config Key  = ok_newmicro
405!Config Desc = 
406!Config Def  = .true.
407!Config Help =
408!
409  ok_newmicro = .true.
410  call getin('ok_newmicro',ok_newmicro)
411!
412!Config Key  = ratqsbas
413!Config Desc = 
414!Config Def  = 0.01
415!Config Help =
416!
417  ratqsbas = 0.01
418  call getin('ratqsbas',ratqsbas)
419!
420!Config Key  = ratqshaut
421!Config Desc = 
422!Config Def  = 0.3
423!Config Help =
424!
425  ratqshaut = 0.3
426  call getin('ratqshaut',ratqshaut)
427
428!
429!Config Key  = rad_froid
430!Config Desc = 
431!Config Def  = 35.0
432!Config Help =
433!
434  rad_froid = 35.0
435  call getin('rad_froid',rad_froid)
436
437!
438!Config Key  = rad_chau1
439!Config Desc = 
440!Config Def  = 13.0
441!Config Help =
442!
443  rad_chau1 = 13.0
444  call getin('rad_chau1',rad_chau1)
445
446!
447!Config Key  = rad_chau2
448!Config Desc = 
449!Config Def  = 9.0
450!Config Help =
451!
452  rad_chau2 = 9.0
453  call getin('rad_chau2',rad_chau2)
454
455!
456!Config Key  = top_height
457!Config Desc =
458!Config Def  = 3
459!Config Help =
460!
461  top_height = 3
462  call getin('top_height',top_height)
463
464!
465!Config Key  = overlap
466!Config Desc =
467!Config Def  = 3
468!Config Help =
469!
470  overlap = 3
471  call getin('overlap',overlap)
472
473
474!
475!
476!Config Key  = cdmmax
477!Config Desc =
478!Config Def  = 1.3E-3
479!Config Help =
480!
481  cdmmax = 1.3E-3
482  call getin('cdmmax',cdmmax)
483
484!
485!Config Key  = cdhmax
486!Config Desc =
487!Config Def  = 1.1E-3
488!Config Help =
489!
490  cdhmax = 1.1E-3
491  call getin('cdhmax',cdhmax)
492
493!261103
494!
495!Config Key  = ksta
496!Config Desc =
497!Config Def  = 1.0e-10
498!Config Help =
499!
500  ksta = 1.0e-10
501  call getin('ksta',ksta)
502
503!
504!Config Key  = ksta_ter
505!Config Desc =
506!Config Def  = 1.0e-10
507!Config Help =
508!
509  ksta_ter = 1.0e-10
510  call getin('ksta_ter',ksta_ter)
511
512!
513!Config Key  = ok_kzmin
514!Config Desc =
515!Config Def  = .true.
516!Config Help =
517!
518  ok_kzmin = .true.
519  call getin('ok_kzmin',ok_kzmin)
520
521!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
522! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
523!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
524
525!Config Key  = iflag_pbl
526!Config Desc =
527!Config Def  = 1
528!Config Help =
529!
530  iflag_pbl = 1
531  call getin('iflag_pbl',iflag_pbl)
532!
533!Config Key  = iflag_thermals
534!Config Desc =
535!Config Def  = 0
536!Config Help =
537!
538  iflag_thermals = 0
539  call getin('iflag_thermals',iflag_thermals)
540!
541!
542!Config Key  = nsplit_thermals
543!Config Desc =
544!Config Def  = 1
545!Config Help =
546!
547  nsplit_thermals = 1
548  call getin('nsplit_thermals',nsplit_thermals)
549
550
551
552!
553!Config Key  = lev_histhf
554!Config Desc =
555!Config Def  = 0
556!Config Help =
557!
558  lev_histhf = 0
559  call getin('lev_histhf',lev_histhf)
560
561!
562!Config Key  = lev_histday
563!Config Desc =
564!Config Def  = 1
565!Config Help =
566!
567  lev_histday = 1
568  call getin('lev_histday',lev_histday)
569
570!
571!Config Key  = lev_histmth
572!Config Desc =
573!Config Def  = 2
574!Config Help =
575!
576  lev_histmth = 2
577  call getin('lev_histmth',lev_histmth)
578
579!
580!Config Key  = type_run
581!Config Desc =
582!Config Def  = 'AMIP' ou 'ENSP'
583!Config Help =
584!
585  type_run = 'AMIP'
586  call getin('type_run',type_run)
587
588!
589!Config Key  = ok_isccp
590!Config Desc =
591!Config Def  = .true.
592!Config Help =
593!
594! ok_isccp = .true.
595  ok_isccp = .false.
596  call getin('ok_isccp',ok_isccp)
597
598!
599!
600!Config Key  = ok_regdyn
601!Config Desc =
602!Config Def  = 'AMIP'
603!Config Help =
604!
605! ok_regdyn = .true.
606  ok_regdyn = .false.
607  call getin('ok_regdyn',ok_regdyn)
608!
609! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone
610! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
611!
612!Config Key  = lonmin_ins
613!Config Desc = 100. 
614!Config Def  = longitude minimale sorties "bilKP_ins"
615!Config Help =
616!
617   lonmin_ins = 100.
618   call getin('lonmin_ins',lonmin_ins)
619!
620!Config Key  = lonmax_ins
621!Config Desc = 130.
622!Config Def  = longitude maximale sorties "bilKP_ins"
623!Config Help =
624!
625   lonmax_ins = 130.
626   call getin('lonmax_ins',lonmax_ins)
627!
628!Config Key  = latmin_ins
629!Config Desc = -20. 
630!Config Def  = latitude minimale sorties "bilKP_ins"
631!Config Help =
632!
633   latmin_ins = -20.
634   call getin('latmin_ins',latmin_ins)
635!
636!Config Key  = latmax_ins
637!Config Desc = 20.
638!Config Def  = latitude maximale sorties "bilKP_ins"
639!Config Help =
640!
641   latmax_ins = 20.
642   call getin('latmax_ins',latmax_ins)
643!
644!Config Key  = ecrit_ins
645!Config Desc =
646!Config Def  = NINT(86400./dtime/48.) !a chaque pas de temps physique
647!Config Help =
648!
649! ecrit_ins = NINT(86400./dtime/48.)
650  ecrit_ins = NINT(86400./48.)
651  call getin('ecrit_ins',ecrit_ins)
652!
653!Config Key  = ecrit_hf
654!Config Desc =
655!Config Def  = NINT(86400./dtime *0.25) !toutes les 6h
656!Config Help =
657!
658! ecrit_hf = NINT(86400./dtime *0.25)
659  ecrit_hf = NINT(86400. *0.25)
660  call getin('ecrit_hf',ecrit_hf)
661!
662!Config Key  = ecrit_hf2mth
663!Config Desc =
664!Config Def  = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h
665!Config Help =
666!
667  ecrit_hf2mth = 4*30
668  call getin('ecrit_hf2mth',ecrit_hf2mth)
669!
670!Config Key  = ecrit_day
671!Config Desc =
672!Config Def  = NINT(86400./dtime *1.0) !tous les jours
673!Config Help =
674!
675! ecrit_day = NINT(86400./dtime *1.0)
676  ecrit_day = NINT(86400. *1.0)
677  call getin('ecrit_day',ecrit_day)
678!
679!Config Key  = ecrit_mth
680!Config Desc =
681!Config Def  = NINT(86400./dtime *ecritphy) !1fois par mois
682!Config Help =
683!
684! ecrit_mth = NINT(86400./dtime *ecritphy)
685  ecrit_mth = NINT(86400. *ecritphy)
686  call getin('ecrit_mth',ecrit_mth)
687!
688!Config Key  = ecrit_tra
689!Config Desc =
690!Config Def  = NINT(86400./dtime * ecritphy) !tous les ecritphy jours
691!Config Help =
692!
693! ecrit_tra = NINT(86400./dtime * ecritphy)
694  ecrit_tra = ecritphy
695  call getin('ecrit_tra',ecrit_tra)
696!
697!Config Key  = ecrit_reg
698!Config Desc =
699!Config Def  = NINT(86400./dtime *0.25)  !4 fois par jour
700!Config Help =
701!
702! ecrit_reg = NINT(86400./dtime *0.25)  !4 fois par jour
703  ecrit_reg = NINT(86400. *0.25)  !4 fois par jour
704  call getin('ecrit_reg',ecrit_reg)
705!
706!
707!
708!Config Key  = bug_ozone
709!Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true
710!Config Def  = false
711!Config Help =
712!
713  bug_ozone = .false.
714  call getin('bug_ozone',bug_ozone)
715
716
717  write(numout,*)' ##############################################'
718  write(numout,*)' Configuration des parametres de la physique: '
719  write(numout,*)' Config ocean = ', ocean
720  write(numout,*)' Config veget = ', ok_veget
721  write(numout,*)' Sortie journaliere = ', ok_journe
722  write(numout,*)' Sortie mensuelle = ', ok_mensuel
723  write(numout,*)' Sortie instantanee = ', ok_instan
724  write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil
725  write(numout,*)' Excentricite = ',R_ecc
726  write(numout,*)' Equinoxe = ',R_peri
727  write(numout,*)' Inclinaison =',R_incl
728  write(numout,*)' Constante solaire =',solaire
729  write(numout,*)' co2_ppm =',co2_ppm
730  write(numout,*)' RCO2 = ',RCO2
731  write(numout,*)' CH4_ppb =',CH4_ppb,' RCH4 = ',RCH4
732  write(numout,*)' N2O_ppb =',N2O_ppb,' RN2O =  ',RN2O
733  write(numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11 =  ',RCFC11
734  write(numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12 =  ',RCFC12
735  write(numout,*)' epmax = ', epmax
736  write(numout,*)' ok_adj_ema = ', ok_adj_ema
737  write(numout,*)' iflag_clw = ', iflag_clw
738  write(numout,*)' cld_lc_lsc = ', cld_lc_lsc
739  write(numout,*)' cld_lc_con = ', cld_lc_con
740  write(numout,*)' cld_tau_lsc = ', cld_tau_lsc
741  write(numout,*)' cld_tau_con = ', cld_tau_con
742  write(numout,*)' ffallv_lsc = ', ffallv_lsc
743  write(numout,*)' ffallv_con = ', ffallv_con
744  write(numout,*)' coef_eva = ', coef_eva
745  write(numout,*)' reevap_ice = ', reevap_ice
746  write(numout,*)' iflag_pdf = ', iflag_pdf
747  write(numout,*)' iflag_cldcon = ', iflag_cldcon
748  write(numout,*)' fact_cldcon = ', fact_cldcon
749  write(numout,*)' facttemps = ', facttemps
750  write(numout,*)' ok_newmicro = ',ok_newmicro
751  write(numout,*)' ratqsbas = ',ratqsbas
752  write(numout,*)' ratqshaut = ',ratqshaut
753  write(numout,*)' top_height = ',top_height
754  write(numout,*)' overlap = ',overlap
755  write(numout,*)' cdmmax = ',cdmmax
756  write(numout,*)' cdhmax = ',cdhmax
757  write(numout,*)' ksta = ',ksta
758  write(numout,*)' ksta_ter = ',ksta_ter
759  write(numout,*)' ok_kzmin = ',ok_kzmin
760  write(numout,*)' ok_ade = ',ok_ade
761  write(numout,*)' ok_aie = ',ok_aie
762  write(numout,*)' bl95_b0 = ',bl95_b0
763  write(numout,*)' bl95_b1 = ',bl95_b1
764  write(numout,*)' lev_histhf = ',lev_histhf
765  write(numout,*)' lev_histday = ',lev_histday
766  write(numout,*)' lev_histmth = ',lev_histmth
767  write(numout,*)' iflag_pbl = ', iflag_pbl
768  write(numout,*)' iflag_thermals = ', iflag_thermals
769  write(numout,*)' type_run = ',type_run
770  write(numout,*)' ok_isccp = ',ok_isccp
771  write(numout,*)' ok_regdyn = ',ok_regdyn
772  write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
773 & lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
774  write(numout,*)' ecrit_ ins, hf, hf2mth, day, mth, reg, tra',&
775 & ecrit_ins, ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra
776  write(numout,*)' bug_ozone = ', bug_ozone
777  return
778
779  end subroutine conf_phys
780
781!
782!#################################################################
783!
784
785   subroutine conf_interface(tau_calv)
786
787   use IOIPSL
788   implicit none
789
790! Configuration de l'interace atm/surf
791!
792! tau_calv:    temps de relaxation pour la fonte des glaciers
793
794  REAL          :: tau_calv
795
796! Local
797  integer              :: numout = 6
798!
799!Config Key  = tau_calv
800!Config Desc = temps de relaxation pour fonte des glaciers en jours
801!Config Def  = 1 an
802!Config Help =
803!
804  tau_calv = 360.*10.
805  call getin('tau_calv',tau_calv)
806
807  write(numout,*)' ##############################################'
808  WRITE(numout,*)' Configuration de l''interface atm/surfaces  : '
809  WRITE(numout,*)' tau_calv = ',tau_calv
810  return
811
812  end subroutine conf_interface
Note: See TracBrowser for help on using the repository browser.