/[lmdze]/trunk/phylmd/conf_phys.f
ViewVC logotype

Diff of /trunk/phylmd/conf_phys.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

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