/[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

trunk/libf/phylmd/conf_phys.f90 revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC trunk/phylmd/conf_phys.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC
# Line 1  Line 1 
1  module conf_phys_m  module conf_phys_m
2    
   ! This module is clean: no C preprocessor directive, no include line.  
   
3    implicit none    implicit none
4    
5    integer iflag_pbl    integer:: iflag_pbl = 1 ! for the planetary boundary layer
6      REAL:: rad_chau1 = 13., rad_chau2 = 9.
7    
8  contains  contains
9    
10    subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &    subroutine conf_phys
11         fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &  
12         ratqsbas,ratqshaut,if_ebil, &      ! From phylmd/conf_phys.F90, version 1.7 2005/07/05 07:21:23
13         ok_ade, ok_aie, &  
14         bl95_b0, bl95_b1,&      ! Configuration de la "physique" de LMDZ.
15         iflag_thermals,nsplit_thermals)  
16        USE clesphys, ONLY: bug_ozone, cdhmax, cdmmax, cfc11_ppt, cfc12_ppt, &
17      ! From phylmd/conf_phys.F90,v 1.7 2005/07/05 07:21:23           ch4_ppb, co2_ppm, ecrit_day, ecrit_hf, ecrit_hf2mth, ecrit_ins, &
18             ecrit_mth, ecrit_reg, ecrit_tra, ksta, ksta_ter, latmax_ins, &
19      use IOIPSL, only: getin           latmin_ins, lev_histday, lev_histhf, lev_histmth, lonmax_ins, &
20      use clesphys, only: solaire, co2_ppm, rco2, ch4_ppb, rch4, n2o_ppb, rn2o, &           lonmin_ins, n2o_ppb, ok_isccp, ok_kzmin, ok_regdyn, overlap, rcfc11, &
21           cfc11_ppt, rcfc11, cfc12_ppt, rcfc12, top_height, overlap, cdmmax, &           rcfc12, rch4, rco2, rn2o, solaire, top_height, type_run
22           cdhmax, ksta, ksta_ter, ok_kzmin, lev_histhf, lev_histday, &      use clesphys2, only: read_clesphys2
23           lev_histmth, type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins, &      USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &
24             cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice
25        USE conema3_m, ONLY: epmax, iflag_clw, ok_adj_ema
26        use unit_nml_m, only: unit_nml
27        USE yomcst, ONLY: r_ecc, r_incl, r_peri
28    
29        namelist /conf_phys_nml/ R_ecc, R_peri, R_incl, solaire, co2_ppm, &
30             CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt, epmax, ok_adj_ema, &
31             iflag_clw, cld_lc_lsc, cld_lc_con, cld_tau_lsc, cld_tau_con, &
32             ffallv_lsc, ffallv_con, coef_eva, reevap_ice, iflag_pdf, &
33             top_height, overlap, cdmmax, cdhmax, ksta, ksta_ter, &
34             ok_kzmin, iflag_pbl, lev_histhf, lev_histday, lev_histmth, &
35             type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins, &
36           latmin_ins, latmax_ins, ecrit_ins, ecrit_hf, ecrit_hf2mth, &           latmin_ins, latmax_ins, ecrit_ins, ecrit_hf, ecrit_hf2mth, &
37           ecrit_day, ecrit_mth, ecrit_tra, ecrit_reg, bug_ozone           ecrit_day, ecrit_mth, ecrit_tra, ecrit_reg, bug_ozone
     use YOMCST  
     use conema3_m  
     use comfisrtilp  
     use nuagecom  
   
     !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12  
   
     ! Configuration de la "physique" de LMDZ a l'aide de la fonction  
     ! GETIN de IOIPSL  
   
     ! ocean:      type d'ocean (force, slab, couple)  
     ! ok_veget:   type de modele de vegetation  
     ! ok_journe:  sorties journalieres  
     ! ok_mensuel: sorties mensuelles  
     ! ok_instan:  sorties instantanees  
     ! ok_ade, ok_aie: apply or not aerosol direct and indirect effects  
     ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc  
   
     ! Sortie:  
     character(len=*), intent(out):: ocean  
     logical              :: ok_veget, ok_newmicro  
     logical              :: ok_journe, ok_mensuel, ok_instan          
     LOGICAL              :: ok_ade, ok_aie  
     REAL                 :: bl95_b0, bl95_b1  
     real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut  
     integer              :: iflag_cldcon, if_ebil  
   
     ! Local  
     integer              :: numout = 6  
     real                 :: zzz  
38    
39      integer :: iflag_thermals,nsplit_thermals      namelist /nuagecom/ rad_chau1, rad_chau2
40    
41      !-----------------------------------------------------------      !-----------------------------------------------------------
42    
43      print *, "Call sequence information: conf_phys"      print *, "Call sequence information: conf_phys"
44        call read_clesphys2
45    
46      !Config Key  = OCEAN      R_ecc = 0.016715 ! AMIP II
47      !Config Desc = Type d'ocean      R_peri = 102.7 ! AMIP II
48      !Config Def  = force      R_incl = 23.441 ! AMIP II
49      !Config Help = Type d'ocean utilise: force ou slab mais pas couple      solaire = 1365. ! AMIP II
50      !      co2_ppm = 348. ! AMIP II
51      ocean = 'force '      CH4_ppb = 1650.
52      call getin('OCEAN', ocean)      N2O_ppb = 306.
53      !      CFC11_ppt = 280.
54      !Config Key  = VEGET      CFC12_ppt = 484.
     !Config Desc = Type de modele de vegetation  
     !Config Def  = .false.  
     !Config Help = Type de modele de vegetation utilise  
     !  
     ok_veget = .false.  
     call getin('VEGET', ok_veget)  
     !  
     !Config Key  = OK_journe  
     !Config Desc = Pour des sorties journalieres  
     !Config Def  = .false.  
     !Config Help = Pour creer le fichier histday contenant les sorties  
     !              journalieres  
     !  
     ok_journe = .false.  
     call getin('OK_journe', ok_journe)  
     !  
     !Config Key  = OK_mensuel  
     !Config Desc = Pour des sorties mensuelles  
     !Config Def  = .true.  
     !Config Help = Pour creer le fichier histmth contenant les sorties  
     !              mensuelles  
     !  
     ok_mensuel = .true.  
     call getin('OK_mensuel', ok_mensuel)  
     !  
     !Config Key  = OK_instan  
     !Config Desc = Pour des sorties instantanees  
     !Config Def  = .false.  
     !Config Help = Pour creer le fichier histins contenant les sorties  
     !              instantanees  
     !  
     ok_instan = .false.  
     call getin('OK_instan', ok_instan)  
     !  
     !Config Key  = ok_ade  
     !Config Desc = Aerosol direct effect or not?  
     !Config Def  = .false.  
     !Config Help = Used in radlwsw.F  
     !  
     ok_ade = .false.  
     call getin('ok_ade', ok_ade)  
   
     !  
     !Config Key  = ok_aie  
     !Config Desc = Aerosol indirect effect or not?  
     !Config Def  = .false.  
     !Config Help = Used in nuage.F and radlwsw.F  
     !  
     ok_aie = .false.  
     call getin('ok_aie', ok_aie)  
   
     !  
     !Config Key  = bl95_b0  
     !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)  
     !Config Def  = .false.  
     !Config Help = Used in nuage.F  
     !  
     bl95_b0 = 2.  
     call getin('bl95_b0', bl95_b0)  
   
     !Config Key  = bl95_b1  
     !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)  
     !Config Def  = .false.  
     !Config Help = Used in nuage.F  
     !  
     bl95_b1 = 0.2  
     call getin('bl95_b1', bl95_b1)  
   
     !  
     !  
     !Config Key  = if_ebil  
     !Config Desc = Niveau de sortie pour les diags bilan d'energie  
     !Config Def  = 0  
     !Config Help =  
     !                
     !  
     if_ebil = 0  
     call getin('if_ebil', if_ebil)  
     !!  
     !! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG  
     !!  
     !Config Key  = R_ecc  
     !Config Desc = Excentricite  
     !Config Def  = 0.016715  
     !Config Help =  
     !                
     !valeur AMIP II  
     R_ecc = 0.016715  
     call getin('R_ecc', R_ecc)  
     !!  
     !Config Key  = R_peri  
     !Config Desc = Equinoxe  
     !Config Def  =  
     !Config Help =  
     !                
     !  
     !valeur AMIP II  
     R_peri = 102.7  
     call getin('R_peri', R_peri)  
     !!  
     !Config Key  = R_incl  
     !Config Desc = Inclinaison  
     !Config Def  =  
     !Config Help =  
     !                
     !  
     !valeur AMIP II  
     R_incl = 23.441  
     call getin('R_incl', R_incl)  
     !!  
     !Config Key  = solaire  
     !Config Desc = Constante solaire en W/m2  
     !Config Def  = 1365.  
     !Config Help =  
     !                
     !  
     !valeur AMIP II  
     solaire = 1365.  
     call getin('solaire', solaire)  
     !!  
     !Config Key  = co2_ppm  
     !Config Desc = concentration du gaz carbonique en ppmv  
     !Config Def  = 348.  
     !Config Help =  
     !                
     !  
     !valeur AMIP II  
     co2_ppm = 348.  
     call getin('co2_ppm', co2_ppm)  
     !!  
     !Config Key  = RCO2  
     !Config Desc = Concentration du CO2  
     !Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97  
     !Config Def  = 348. * 1.0e-06  * 44.011/28.97  
     !Config Help =  
     !                
     ! RCO2 = 5.286789092164308E-04  
     !ancienne valeur  
     RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.  
   
     !!  call getin('RCO2', RCO2)  
     !!  
     !Config Key  = RCH4  
     !Config Desc = Concentration du CH4  
     !Config Def  = 1.65E-06* 16.043/28.97  
     !Config Help =  
     !                
     !  
     !valeur AMIP II  
     !OK  RCH4 = 1.65E-06* 16.043/28.97  
     ! RCH4 = 9.137366240938903E-07  
     !  
     !ancienne valeur  
     ! RCH4 = 1.72E-06* 16.043/28.97  
     !OK call getin('RCH4', RCH4)  
     zzz = 1650.  
     call getin('CH4_ppb', zzz)  
     CH4_ppb = zzz  
     RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97  
     !!  
     !Config Key  = RN2O  
     !Config Desc = Concentration du N2O  
     !Config Def  = 306.E-09* 44.013/28.97  
     !Config Help =  
     !                
     !  
     !valeur AMIP II  
     !OK  RN2O = 306.E-09* 44.013/28.97  
     ! RN2O = 4.648939592682085E-07  
     !  
     !ancienne valeur  
     ! RN2O = 310.E-09* 44.013/28.97  
     !OK  call getin('RN2O', RN2O)  
     zzz=306.  
     call getin('N2O_ppb', zzz)  
     N2O_ppb = zzz  
     RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97  
     !!  
     !Config Key  = RCFC11  
     !Config Desc = Concentration du CFC11  
     !Config Def  = 280.E-12* 137.3686/28.97  
     !Config Help =  
     !                
     !  
     !OK RCFC11 = 280.E-12* 137.3686/28.97  
     zzz = 280.  
     call getin('CFC11_ppt',zzz)  
     CFC11_ppt = zzz  
     RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97  
     ! RCFC11 = 1.327690990680013E-09  
     !OK call getin('RCFC11', RCFC11)  
     !!  
     !Config Key  = RCFC12  
     !Config Desc = Concentration du CFC12  
     !Config Def  = 484.E-12* 120.9140/28.97  
     !Config Help =  
     !                
     !  
     !OK RCFC12 = 484.E-12* 120.9140/28.97  
     zzz = 484.  
     call getin('CFC12_ppt',zzz)  
     CFC12_ppt = zzz  
     RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97  
     ! RCFC12 = 2.020102726958923E-09  
     !OK call getin('RCFC12', RCFC12)  
     !!  
     !! Constante solaire & Parametres orbitaux & taux gaz effet de serre END  
     !!  
     !! KE  
     !  
     !Config Key  = epmax  
     !Config Desc = Efficacite precip  
     !Config Def  = 0.993  
     !Config Help =  
     !  
55      epmax = .993      epmax = .993
     call getin('epmax', epmax)  
     !  
     !Config Key  = ok_adj_ema  
     !Config Desc =    
     !Config Def  = false  
     !Config Help =  
     !  
56      ok_adj_ema = .false.      ok_adj_ema = .false.
     call getin('ok_adj_ema',ok_adj_ema)  
     !  
     !Config Key  = iflag_clw  
     !Config Desc =    
     !Config Def  = 0  
     !Config Help =  
     !  
57      iflag_clw = 0      iflag_clw = 0
     call getin('iflag_clw',iflag_clw)  
     !  
     !Config Key  = cld_lc_lsc  
     !Config Desc =    
     !Config Def  = 2.6e-4  
     !Config Help =  
     !  
58      cld_lc_lsc = 2.6e-4      cld_lc_lsc = 2.6e-4
     call getin('cld_lc_lsc',cld_lc_lsc)  
     !  
     !Config Key  = cld_lc_con  
     !Config Desc =    
     !Config Def  = 2.6e-4  
     !Config Help =  
     !  
59      cld_lc_con = 2.6e-4      cld_lc_con = 2.6e-4
     call getin('cld_lc_con',cld_lc_con)  
     !  
     !Config Key  = cld_tau_lsc  
     !Config Desc =    
     !Config Def  = 3600.  
     !Config Help =  
     !  
60      cld_tau_lsc = 3600.      cld_tau_lsc = 3600.
     call getin('cld_tau_lsc',cld_tau_lsc)  
     !  
     !Config Key  = cld_tau_con  
     !Config Desc =    
     !Config Def  = 3600.  
     !Config Help =  
     !  
61      cld_tau_con = 3600.      cld_tau_con = 3600.
     call getin('cld_tau_con',cld_tau_con)  
     !  
     !Config Key  = ffallv_lsc  
     !Config Desc =    
     !Config Def  = 1.  
     !Config Help =  
     !  
62      ffallv_lsc = 1.      ffallv_lsc = 1.
     call getin('ffallv_lsc',ffallv_lsc)  
     !  
     !Config Key  = ffallv_con  
     !Config Desc =    
     !Config Def  = 1.  
     !Config Help =  
     !  
63      ffallv_con = 1.      ffallv_con = 1.
     call getin('ffallv_con',ffallv_con)  
     !  
     !Config Key  = coef_eva  
     !Config Desc =    
     !Config Def  = 2.e-5  
     !Config Help =  
     !  
64      coef_eva = 2.e-5      coef_eva = 2.e-5
     call getin('coef_eva',coef_eva)  
     !  
     !Config Key  = reevap_ice  
     !Config Desc =    
     !Config Def  = .false.  
     !Config Help =  
     !  
65      reevap_ice = .false.      reevap_ice = .false.
     call getin('reevap_ice',reevap_ice)  
     !  
     !Config Key  = iflag_cldcon  
     !Config Desc =    
     !Config Def  = 1  
     !Config Help =  
     !  
     iflag_cldcon = 1  
     call getin('iflag_cldcon',iflag_cldcon)  
   
     !  
     !Config Key  = iflag_pdf  
     !Config Desc =    
     !Config Def  = 0  
     !Config Help =  
     !  
66      iflag_pdf = 0      iflag_pdf = 0
     call getin('iflag_pdf',iflag_pdf)  
     !  
     !Config Key  = fact_cldcon  
     !Config Desc =    
     !Config Def  = 0.375  
     !Config Help =  
     !  
     fact_cldcon = 0.375  
     call getin('fact_cldcon',fact_cldcon)  
   
     !  
     !Config Key  = facttemps  
     !Config Desc =    
     !Config Def  = 1.e-4  
     !Config Help =  
     !  
     facttemps = 1.e-4  
     call getin('facttemps',facttemps)  
   
     !  
     !Config Key  = ok_newmicro  
     !Config Desc =    
     !Config Def  = .true.  
     !Config Help =  
     !  
     ok_newmicro = .true.  
     call getin('ok_newmicro',ok_newmicro)  
     !  
     !Config Key  = ratqsbas  
     !Config Desc =    
     !Config Def  = 0.01  
     !Config Help =  
     !  
     ratqsbas = 0.01  
     call getin('ratqsbas',ratqsbas)  
     !  
     !Config Key  = ratqshaut  
     !Config Desc =    
     !Config Def  = 0.3  
     !Config Help =  
     !  
     ratqshaut = 0.3  
     call getin('ratqshaut',ratqshaut)  
   
     !  
     !Config Key  = rad_froid  
     !Config Desc =    
     !Config Def  = 35.0  
     !Config Help =  
     !  
     rad_froid = 35.0  
     call getin('rad_froid',rad_froid)  
   
     !  
     !Config Key  = rad_chau1  
     !Config Desc =    
     !Config Def  = 13.0  
     !Config Help =  
     !  
     rad_chau1 = 13.0  
     call getin('rad_chau1',rad_chau1)  
   
     !  
     !Config Key  = rad_chau2  
     !Config Desc =    
     !Config Def  = 9.0  
     !Config Help =  
     !  
     rad_chau2 = 9.0  
     call getin('rad_chau2',rad_chau2)  
   
     !  
     !Config Key  = top_height  
     !Config Desc =  
     !Config Def  = 3  
     !Config Help =  
     !  
67      top_height = 3      top_height = 3
     call getin('top_height',top_height)  
   
     !  
     !Config Key  = overlap  
     !Config Desc =  
     !Config Def  = 3  
     !Config Help =  
     !  
68      overlap = 3      overlap = 3
     call getin('overlap',overlap)  
   
   
     !  
     !  
     !Config Key  = cdmmax  
     !Config Desc =  
     !Config Def  = 1.3E-3  
     !Config Help =  
     !  
69      cdmmax = 1.3E-3      cdmmax = 1.3E-3
     call getin('cdmmax',cdmmax)  
   
     !  
     !Config Key  = cdhmax  
     !Config Desc =  
     !Config Def  = 1.1E-3  
     !Config Help =  
     !  
70      cdhmax = 1.1E-3      cdhmax = 1.1E-3
     call getin('cdhmax',cdhmax)  
   
     !261103  
     !  
     !Config Key  = ksta  
     !Config Desc =  
     !Config Def  = 1.0e-10  
     !Config Help =  
     !  
71      ksta = 1.0e-10      ksta = 1.0e-10
     call getin('ksta',ksta)  
   
     !  
     !Config Key  = ksta_ter  
     !Config Desc =  
     !Config Def  = 1.0e-10  
     !Config Help =  
     !  
72      ksta_ter = 1.0e-10      ksta_ter = 1.0e-10
     call getin('ksta_ter',ksta_ter)  
   
     !  
     !Config Key  = ok_kzmin  
     !Config Desc =  
     !Config Def  = .true.  
     !Config Help =  
     !  
73      ok_kzmin = .true.      ok_kzmin = .true.
     call getin('ok_kzmin',ok_kzmin)  
   
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  
     ! PARAMETER FOR THE PLANETARY BOUNDARY LAYER  
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  
   
     !Config Key  = iflag_pbl  
     !Config Desc =  
     !Config Def  = 1  
     !Config Help =  
     !  
     iflag_pbl = 1  
     call getin('iflag_pbl',iflag_pbl)  
     !  
     !Config Key  = iflag_thermals  
     !Config Desc =  
     !Config Def  = 0  
     !Config Help =  
     !  
     iflag_thermals = 0  
     call getin('iflag_thermals',iflag_thermals)  
     !  
     !  
     !Config Key  = nsplit_thermals  
     !Config Desc =  
     !Config Def  = 1  
     !Config Help =  
     !  
     nsplit_thermals = 1  
     call getin('nsplit_thermals',nsplit_thermals)  
   
   
   
     !  
     !Config Key  = lev_histhf  
     !Config Desc =  
     !Config Def  = 0  
     !Config Help =  
     !  
74      lev_histhf = 0      lev_histhf = 0
     call getin('lev_histhf',lev_histhf)  
   
     !  
     !Config Key  = lev_histday  
     !Config Desc =  
     !Config Def  = 1  
     !Config Help =  
     !  
75      lev_histday = 1      lev_histday = 1
     call getin('lev_histday',lev_histday)  
   
     !  
     !Config Key  = lev_histmth  
     !Config Desc =  
     !Config Def  = 2  
     !Config Help =  
     !  
76      lev_histmth = 2      lev_histmth = 2
     call getin('lev_histmth',lev_histmth)  
   
     !  
     !Config Key  = type_run  
     !Config Desc =  
     !Config Def  = 'AMIP' ou 'ENSP'  
     !Config Help =  
     !  
77      type_run = 'AMIP'      type_run = 'AMIP'
     call getin('type_run',type_run)  
   
     !  
     !Config Key  = ok_isccp  
     !Config Desc =  
     !Config Def  = .true.  
     !Config Help =  
     !  
     ! ok_isccp = .true.  
78      ok_isccp = .false.      ok_isccp = .false.
     call getin('ok_isccp',ok_isccp)  
   
     !  
     !  
     !Config Key  = ok_regdyn  
     !Config Desc =  
     !Config Def  = 'AMIP'  
     !Config Help =  
     !  
     ! ok_regdyn = .true.  
79      ok_regdyn = .false.      ok_regdyn = .false.
     call getin('ok_regdyn',ok_regdyn)  
     !  
     ! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone  
     ! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"  
     !  
     !Config Key  = lonmin_ins  
     !Config Desc = 100.    
     !Config Def  = longitude minimale sorties "bilKP_ins"  
     !Config Help =  
     !  
80      lonmin_ins = 100.      lonmin_ins = 100.
     call getin('lonmin_ins',lonmin_ins)  
     !  
     !Config Key  = lonmax_ins  
     !Config Desc = 130.  
     !Config Def  = longitude maximale sorties "bilKP_ins"  
     !Config Help =  
     !  
81      lonmax_ins = 130.      lonmax_ins = 130.
     call getin('lonmax_ins',lonmax_ins)  
     !  
     !Config Key  = latmin_ins  
     !Config Desc = -20.    
     !Config Def  = latitude minimale sorties "bilKP_ins"  
     !Config Help =  
     !  
82      latmin_ins = -20.      latmin_ins = -20.
     call getin('latmin_ins',latmin_ins)  
     !  
     !Config Key  = latmax_ins  
     !Config Desc = 20.  
     !Config Def  = latitude maximale sorties "bilKP_ins"  
     !Config Help =  
     !  
83      latmax_ins = 20.      latmax_ins = 20.
     call getin('latmax_ins',latmax_ins)  
     !  
     !Config Key  = ecrit_ins  
     !Config Desc =  
     !Config Def  = NINT(86400./dtime/48.) !a chaque pas de temps physique  
     !Config Help =  
     !  
     ! ecrit_ins = NINT(86400./dtime/48.)  
84      ecrit_ins = NINT(86400./48.)      ecrit_ins = NINT(86400./48.)
     call getin('ecrit_ins',ecrit_ins)  
     !  
     !Config Key  = ecrit_hf  
     !Config Desc =  
     !Config Def  = NINT(86400./dtime *0.25) !toutes les 6h  
     !Config Help =  
     !  
     ! ecrit_hf = NINT(86400./dtime *0.25)  
85      ecrit_hf = NINT(86400. *0.25)      ecrit_hf = NINT(86400. *0.25)
86      call getin('ecrit_hf',ecrit_hf)      ecrit_hf2mth = 4*30 ! ecriture mens. a partir de val. inst. toutes les 6h
     !  
     !Config Key  = ecrit_hf2mth  
     !Config Desc =  
     !Config Def  = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h  
     !Config Help =  
     !  
     ecrit_hf2mth = 4*30  
     call getin('ecrit_hf2mth',ecrit_hf2mth)  
     !  
     !Config Key  = ecrit_day  
87      ecrit_day = 86400      ecrit_day = 86400
     call getin('ecrit_day',ecrit_day)  
     !  
88      ecrit_mth = 86400      ecrit_mth = 86400
     call getin('ecrit_mth',ecrit_mth)  
     !  
89      ecrit_tra = 1      ecrit_tra = 1
90      call getin('ecrit_tra',ecrit_tra)      ecrit_reg = NINT(86400. *0.25) ! 4 fois par jour
     !  
     !Config Key  = ecrit_reg  
     !Config Desc =  
     !Config Def  = NINT(86400./dtime *0.25)  !4 fois par jour  
     !Config Help =  
     !  
     ! ecrit_reg = NINT(86400./dtime *0.25)  !4 fois par jour  
     ecrit_reg = NINT(86400. *0.25)  !4 fois par jour  
     call getin('ecrit_reg',ecrit_reg)  
     !  
     !  
     !  
     !Config Key  = bug_ozone  
     !Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true  
     !Config Def  = false  
     !Config Help =  
     !  
91      bug_ozone = .false.      bug_ozone = .false.
     call getin('bug_ozone',bug_ozone)  
92    
93      write(numout,*)' *********'      print *, "Enter namelist 'conf_phys_nml'."
94      write(numout,*)' Configuration des parametres de la physique: '      read(unit=*, nml=conf_phys_nml)
95      write(numout,*)' Config ocean = ', ocean      write(unit_nml, nml=conf_phys_nml)
96      write(numout,*)' Config veget = ', ok_veget  
97      write(numout,*)' Sortie journaliere = ', ok_journe      RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97
98      write(numout,*)' Sortie mensuelle = ', ok_mensuel      RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
99      write(numout,*)' Sortie instantanee = ', ok_instan      RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
100      write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil      RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
101      write(numout,*)' Excentricite = ',R_ecc      RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
102      write(numout,*)' Equinoxe = ',R_peri  
103      write(numout,*)' Inclinaison =',R_incl      print *, ' RCO2 = ', RCO2
104      write(numout,*)' Constante solaire =',solaire      print *, ' RCH4 = ', RCH4
105      write(numout,*)' co2_ppm =',co2_ppm      print *, ' RN2O = ', RN2O
106      write(numout,*)' RCO2 = ',RCO2      print *, ' RCFC11 = ', RCFC11
107      write(numout,*)' CH4_ppb =',CH4_ppb,' RCH4 = ',RCH4      print *, ' RCFC12 = ', RCFC12
108      write(numout,*)' N2O_ppb =',N2O_ppb,' RN2O =  ',RN2O  
109      write(numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11 =  ',RCFC11      print *, "Enter namelist 'nuagecom'."
110      write(numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12 =  ',RCFC12      read(unit=*, nml=nuagecom)
111      write(numout,*)' epmax = ', epmax      write(unit_nml, nml=nuagecom)
     write(numout,*)' ok_adj_ema = ', ok_adj_ema  
     write(numout,*)' iflag_clw = ', iflag_clw  
     write(numout,*)' cld_lc_lsc = ', cld_lc_lsc  
     write(numout,*)' cld_lc_con = ', cld_lc_con  
     write(numout,*)' cld_tau_lsc = ', cld_tau_lsc  
     write(numout,*)' cld_tau_con = ', cld_tau_con  
     write(numout,*)' ffallv_lsc = ', ffallv_lsc  
     write(numout,*)' ffallv_con = ', ffallv_con  
     write(numout,*)' coef_eva = ', coef_eva  
     write(numout,*)' reevap_ice = ', reevap_ice  
     write(numout,*)' iflag_pdf = ', iflag_pdf  
     write(numout,*)' iflag_cldcon = ', iflag_cldcon  
     write(numout,*)' fact_cldcon = ', fact_cldcon  
     write(numout,*)' facttemps = ', facttemps  
     write(numout,*)' ok_newmicro = ',ok_newmicro  
     write(numout,*)' ratqsbas = ',ratqsbas  
     write(numout,*)' ratqshaut = ',ratqshaut  
     write(numout,*)' top_height = ',top_height  
     write(numout,*)' overlap = ',overlap  
     write(numout,*)' cdmmax = ',cdmmax  
     write(numout,*)' cdhmax = ',cdhmax  
     write(numout,*)' ksta = ',ksta  
     write(numout,*)' ksta_ter = ',ksta_ter  
     write(numout,*)' ok_kzmin = ',ok_kzmin  
     write(numout,*)' ok_ade = ',ok_ade  
     write(numout,*)' ok_aie = ',ok_aie  
     write(numout,*)' bl95_b0 = ',bl95_b0  
     write(numout,*)' bl95_b1 = ',bl95_b1  
     write(numout,*)' lev_histhf = ',lev_histhf  
     write(numout,*)' lev_histday = ',lev_histday  
     write(numout,*)' lev_histmth = ',lev_histmth  
     write(numout,*)' iflag_pbl = ', iflag_pbl  
     write(numout,*)' iflag_thermals = ', iflag_thermals  
     write(numout,*)' type_run = ',type_run  
     write(numout,*)' ok_isccp = ',ok_isccp  
     write(numout,*)' ok_regdyn = ',ok_regdyn  
     write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&  
          lonmin_ins, lonmax_ins, latmin_ins, latmax_ins  
     write(numout,*) 'ecrit_ ins, hf, hf2mth, day, mth, reg, tra', ecrit_ins, &  
          ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra  
     write(numout,*)' bug_ozone = ', bug_ozone  
112    
113    end subroutine conf_phys    end subroutine conf_phys
114    
Legend:
Removed from v.12  
changed lines
  Added in v.82
  ViewVC Help
Powered by ViewVC 1.1.21