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

Last change on this file since 701 was 699, checked in by Laurent Fairhead, 19 years ago

En fait le bug sur l'ozone n'en etait pas un! on revient a la situation precedente MPL/JLD
LF

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