Ignore:
Timestamp:
Sep 4, 2018, 1:39:46 PM (6 years ago)
Author:
oboucher
Message:

Adding the reading of flag level_coupling_esm needed for the ESM version
Cleaning up the conf_phys_m.F90 routine with no further change

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

    r3378 r3384  
    2626    USE surface_data
    2727    USE phys_cal_mod
    28     USE carbon_cycle_mod,  ONLY: carbon_cycle_tr, carbon_cycle_cpl
     28    USE carbon_cycle_mod,  ONLY: carbon_cycle_tr, carbon_cycle_cpl, level_coupling_esm
    2929    USE mod_grid_phy_lmdz, ONLY: klon_glo
    3030    USE print_control_mod, ONLY: lunout
    3131
    32     include "conema3.h"
    33     include "fisrtilp.h"
    34     include "nuage.h"
    35     include "YOMCST.h"
    36     include "YOMCST2.h"
    37 
    38     include "thermcell.h"
    39 
     32    INCLUDE "conema3.h"
     33    INCLUDE "fisrtilp.h"
     34    INCLUDE "nuage.h"
     35    INCLUDE "YOMCST.h"
     36    INCLUDE "YOMCST2.h"
     37    INCLUDE "thermcell.h"
    4038
    4139    !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
    42     include "clesphys.h"
    43     include "compbl.h"
    44     include "comsoil.h"
    45     include "YOEGWD.h"
     40    INCLUDE "clesphys.h"
     41    INCLUDE "compbl.h"
     42    INCLUDE "comsoil.h"
     43    INCLUDE "YOEGWD.h"
    4644    !
    4745    ! Configuration de la "physique" de LMDZ a l'aide de la fonction
    … …  
    4947    !
    5048    ! LF 05/2001
    51     !
    52 
    5349    !
    5450    ! type_ocean:      type d'ocean (force, slab, couple)
    … …  
    6763    ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
    6864    !
    69 
    7065
    7166    ! Sortie:
    … …  
    108103    REAL,SAVE           :: tau_cld_cv_omp, coefw_cld_cv_omp
    109104    INTEGER, SAVE       :: iflag_cld_cv_omp
    110 
    111105
    112106    REAL, SAVE          :: ratqshaut_omp
    … …  
    201195    INTEGER, SAVE :: levout_histNMC_omp(3)
    202196    LOGICAL, SAVE :: ok_histNMC_omp(3)
    203     REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3)
     197    REAL, SAVE    :: freq_outNMC_omp(3), freq_calNMC_omp(3)
    204198    CHARACTER*4, SAVE :: type_run_omp
    205     LOGICAL,SAVE :: ok_cosp_omp, ok_airs_omp
    206     LOGICAL,SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
    207     REAL,SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
    208     REAL,SAVE :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp
    209     REAL,SAVE :: ecrit_ins_omp
    210     REAL,SAVE :: ecrit_LES_omp
    211     REAL,SAVE :: ecrit_tra_omp
    212     REAL,SAVE :: cvl_comp_threshold_omp
    213     REAL,SAVE :: cvl_sig2feed_omp
    214     REAL,SAVE :: cvl_corr_omp
    215     LOGICAL,SAVE :: ok_lic_melt_omp
    216     LOGICAL,SAVE :: ok_lic_cond_omp
    217     !
    218     INTEGER,SAVE :: iflag_cycle_diurne_omp
    219     LOGICAL,SAVE :: soil_model_omp,new_oliq_omp
    220     LOGICAL,SAVE :: ok_orodr_omp, ok_orolf_omp, ok_limitvrai_omp
     199    LOGICAL, SAVE :: ok_cosp_omp, ok_airs_omp
     200    LOGICAL, SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
     201    REAL, SAVE    :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
     202    REAL, SAVE    :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp
     203    REAL, SAVE    :: ecrit_ins_omp
     204    REAL, SAVE    :: ecrit_LES_omp
     205    REAL, SAVE    :: ecrit_tra_omp
     206    REAL, SAVE    :: cvl_comp_threshold_omp
     207    REAL, SAVE    :: cvl_sig2feed_omp
     208    REAL, SAVE    :: cvl_corr_omp
     209    LOGICAL, SAVE :: ok_lic_melt_omp
     210    LOGICAL, SAVE :: ok_lic_cond_omp
     211    !
     212    INTEGER, SAVE :: iflag_cycle_diurne_omp
     213    LOGICAL, SAVE :: soil_model_omp,new_oliq_omp
     214    LOGICAL, SAVE :: ok_orodr_omp, ok_orolf_omp, ok_limitvrai_omp
    221215    INTEGER, SAVE :: nbapp_rad_omp, iflag_con_omp
    222216    INTEGER, SAVE :: nbapp_cv_omp, nbapp_wk_omp
    … …  
    231225    REAL, SAVE    :: sso_gkdrag_omp,sso_grahil_omp,sso_grcrit_omp
    232226    REAL, SAVE    :: sso_gfrcri_omp,sso_gkwake_omp,sso_gklift_omp
    233     LOGICAL,SAVE  :: ok_qch4_omp
    234     LOGICAL,SAVE  :: carbon_cycle_tr_omp
    235     LOGICAL,SAVE  :: carbon_cycle_cpl_omp
    236     LOGICAL,SAVE  :: adjust_tropopause_omp
    237     LOGICAL,SAVE  :: ok_daily_climoz_omp
     227    LOGICAL, SAVE :: ok_qch4_omp
     228    LOGICAL, SAVE :: carbon_cycle_tr_omp
     229    LOGICAL, SAVE :: carbon_cycle_cpl_omp
     230    INTEGER, SAVE :: level_coupling_esm_omp
     231    LOGICAL, SAVE :: adjust_tropopause_omp
     232    LOGICAL, SAVE :: ok_daily_climoz_omp
    238233
    239234    INTEGER, INTENT(OUT):: read_climoz ! read ozone climatology, OpenMP shared
    … …  
    21482143    carbon_cycle_cpl_omp=.FALSE.
    21492144    CALL getin('carbon_cycle_cpl',carbon_cycle_cpl_omp)
     2145
     2146    ! >> PC
     2147    ! level_coupling_esm : level of coupling of the biogeochemical fields between LMDZ, ORCHIDEE and NEMO
     2148    ! Definitions of level_coupling_esm in physiq.def
     2149    ! level_coupling_esm = 0  ! No field exchange between LMDZ and ORCHIDEE models
     2150    !                         ! No field exchange between LMDZ and NEMO
     2151    ! level_coupling_esm = 1  ! Field exchange between LMDZ and ORCHIDEE models
     2152    !                         ! No field exchange between LMDZ and NEMO models
     2153    ! level_coupling_esm = 2  ! No field exchange between LMDZ and ORCHIDEE models
     2154    !                         ! Field exchange between LMDZ and NEMO models
     2155    ! level_coupling_esm = 3  ! Field exchange between LMDZ and ORCHIDEE models
     2156    !                         ! Field exchange between LMDZ and NEMO models
     2157    level_coupling_esm_omp=0 ! default value
     2158    CALL getin('level_coupling_esm',level_coupling_esm_omp)
     2159    ! << PC
    21502160
    21512161    !$OMP END MASTER
    … …  
    24042414    carbon_cycle_tr = carbon_cycle_tr_omp
    24052415    carbon_cycle_cpl = carbon_cycle_cpl_omp
     2416    level_coupling_esm = level_coupling_esm_omp
    24062417
    24072418    ! Test of coherence between type_ocean and version_ocean
    … …  
    25112522       WRITE(lunout,*)'ifl_pbltree is now changed to zero'
    25122523       ifl_pbltree=0
    2513     END IF
     2524    ENDIF
    25142525
    25152526    !$OMP MASTER
    25162527
    2517     write(lunout,*)' ##############################################'
    2518     write(lunout,*)' Configuration des parametres de la physique: '
    2519     write(lunout,*)' Type ocean = ', type_ocean
    2520     write(lunout,*)' Version ocean = ', version_ocean
    2521     write(lunout,*)' Config veget = ', ok_veget,type_veget
    2522     write(lunout,*)' Snow model SISVAT : ok_snow = ', ok_snow
    2523     write(lunout,*)' Config xml pour XIOS : ok_all_xml = ', ok_all_xml
    2524     write(lunout,*)' Sortie journaliere = ', ok_journe
    2525     write(lunout,*)' Sortie haute frequence = ', ok_hf
    2526     write(lunout,*)' Sortie mensuelle = ', ok_mensuel
    2527     write(lunout,*)' Sortie instantanee = ', ok_instan
    2528     write(lunout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
    2529     write(lunout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
    2530     write(lunout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
    2531     write(lunout,*)' Frequence appel simulateur AIRS, freq_AIRS =', freq_AIRS
    2532     write(lunout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
    2533     write(lunout,*)' Excentricite = ',R_ecc
    2534     write(lunout,*)' Equinoxe = ',R_peri
    2535     write(lunout,*)' Inclinaison =',R_incl
    2536     write(lunout,*)' Constante solaire =',solaire
    2537     write(lunout,*)' ok_suntime_rrtm =',ok_suntime_rrtm
    2538     write(lunout,*)' co2_ppm =',co2_ppm
    2539     write(lunout,*)' RCO2_act = ',RCO2_act
    2540     write(lunout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
    2541     write(lunout,*)' N2O_ppb =',N2O_ppb,' RN2O_act=  ',RN2O_act
    2542     write(lunout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act=  ',RCFC11_act
    2543     write(lunout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act=  ',RCFC12_act
    2544     write(lunout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
    2545     write(lunout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
    2546     write(lunout,*)' RCFC12_per = ',RCFC12_per
    2547     write(lunout,*)' cvl_comp_threshold=', cvl_comp_threshold
    2548     write(lunout,*)' cvl_sig2feed=', cvl_sig2feed
    2549     write(lunout,*)' cvl_corr=', cvl_corr
    2550     write(lunout,*)'ok_lic_melt=', ok_lic_melt
    2551     write(lunout,*)'ok_lic_cond=', ok_lic_cond
    2552     write(lunout,*)'iflag_cycle_diurne=',iflag_cycle_diurne
    2553     write(lunout,*)'soil_model=',soil_model
    2554     write(lunout,*)'new_oliq=',new_oliq
    2555     write(lunout,*)'ok_orodr=',ok_orodr
    2556     write(lunout,*)'ok_orolf=',ok_orolf
    2557     write(lunout,*)'ok_limitvrai=',ok_limitvrai
    2558     write(lunout,*)'nbapp_rad=',nbapp_rad
    2559     write(lunout,*)'iflag_con=',iflag_con
    2560     write(lunout,*)'nbapp_cv=',nbapp_cv
    2561     write(lunout,*)'nbapp_wk=',nbapp_wk
    2562     write(lunout,*)'iflag_ener_conserv=',iflag_ener_conserv
    2563     write(lunout,*)'ok_conserv_q=',ok_conserv_q
    2564     write(lunout,*)'iflag_fisrtilp_qsat=',iflag_fisrtilp_qsat
    2565     write(lunout,*)'iflag_bergeron=',iflag_bergeron
    2566     write(lunout,*)' epmax = ', epmax
    2567     write(lunout,*)' coef_epmax_cape = ', coef_epmax_cape
    2568     write(lunout,*)' ok_adj_ema = ', ok_adj_ema
    2569     write(lunout,*)' iflag_clw = ', iflag_clw
    2570     write(lunout,*)' cld_lc_lsc = ', cld_lc_lsc
    2571     write(lunout,*)' cld_lc_con = ', cld_lc_con
    2572     write(lunout,*)' cld_tau_lsc = ', cld_tau_lsc
    2573     write(lunout,*)' cld_tau_con = ', cld_tau_con
    2574     write(lunout,*)' ffallv_lsc = ', ffallv_lsc
    2575     write(lunout,*)' ffallv_con = ', ffallv_con
    2576     write(lunout,*)' coef_eva = ', coef_eva
    2577     write(lunout,*)' reevap_ice = ', reevap_ice
    2578     write(lunout,*)' iflag_pdf = ', iflag_pdf
    2579     write(lunout,*)' iflag_cld_th = ', iflag_cld_th
    2580     write(lunout,*)' iflag_cld_cv = ', iflag_cld_cv
    2581     write(lunout,*)' tau_cld_cv = ', tau_cld_cv
    2582     write(lunout,*)' coefw_cld_cv = ', coefw_cld_cv
    2583     write(lunout,*)' iflag_radia = ', iflag_radia
    2584     write(lunout,*)' iflag_rrtm = ', iflag_rrtm
    2585     write(lunout,*)' NSW = ', NSW
    2586     write(lunout,*)' iflag_albedo = ', iflag_albedo !albedo SB
    2587     write(lunout,*)' ok_chlorophyll =',ok_chlorophyll ! albedo SB
    2588     write(lunout,*)' iflag_ratqs = ', iflag_ratqs
    2589     write(lunout,*)' seuil_inversion = ', seuil_inversion
    2590     write(lunout,*)' fact_cldcon = ', fact_cldcon
    2591     write(lunout,*)' facttemps = ', facttemps
    2592     write(lunout,*)' ok_newmicro = ',ok_newmicro
    2593     write(lunout,*)' ratqsbas = ',ratqsbas
    2594     write(lunout,*)' ratqshaut = ',ratqshaut
    2595     write(lunout,*)' tau_ratqs = ',tau_ratqs
    2596     write(lunout,*)' top_height = ',top_height
    2597     write(lunout,*)' rad_froid = ',rad_froid
    2598     write(lunout,*)' rad_chau1 = ',rad_chau1
    2599     write(lunout,*)' rad_chau2 = ',rad_chau2
    2600     write(lunout,*)' t_glace_min = ',t_glace_min
    2601     write(lunout,*)' t_glace_max = ',t_glace_max
    2602     write(lunout,*)' exposant_glace = ',exposant_glace
    2603     write(lunout,*)' iflag_t_glace = ',iflag_t_glace
    2604     write(lunout,*)' iflag_cloudth_vert = ',iflag_cloudth_vert
    2605     write(lunout,*)' iflag_rain_incloud_vol = ',iflag_rain_incloud_vol
    2606     write(lunout,*)' iflag_ice_thermo = ',iflag_ice_thermo
    2607     write(lunout,*)' rei_min = ',rei_min
    2608     write(lunout,*)' rei_max = ',rei_max
    2609     write(lunout,*)' overlap = ',overlap
    2610     write(lunout,*)' cdmmax = ',cdmmax
    2611     write(lunout,*)' cdhmax = ',cdhmax
    2612     write(lunout,*)' ksta = ',ksta
    2613     write(lunout,*)' ksta_ter = ',ksta_ter
    2614     write(lunout,*)' f_ri_cd_min = ',f_ri_cd_min
    2615     write(lunout,*)' ok_kzmin = ',ok_kzmin
    2616     write(lunout,*)' pbl_lmixmin_alpha = ',pbl_lmixmin_alpha
    2617     write(lunout,*)' fmagic = ',fmagic
    2618     write(lunout,*)' pmagic = ',pmagic
    2619     write(lunout,*)' ok_ade = ',ok_ade
    2620     write(lunout,*)' ok_aie = ',ok_aie
    2621     write(lunout,*)' ok_alw = ',ok_alw
    2622     write(lunout,*)' aerosol_couple = ', aerosol_couple
    2623     write(lunout,*)' chemistry_couple = ', chemistry_couple
    2624     write(lunout,*)' flag_aerosol = ', flag_aerosol
    2625     write(lunout,*)' flag_aerosol_strat= ', flag_aerosol_strat
    2626     write(lunout,*)' new_aod = ', new_aod
    2627     write(lunout,*)' aer_type = ',aer_type
    2628     write(lunout,*)' bl95_b0 = ',bl95_b0
    2629     write(lunout,*)' bl95_b1 = ',bl95_b1
    2630     write(lunout,*)' lev_histhf = ',lev_histhf
    2631     write(lunout,*)' lev_histday = ',lev_histday
    2632     write(lunout,*)' lev_histmth = ',lev_histmth
    2633     write(lunout,*)' lev_histins = ',lev_histins
    2634     write(lunout,*)' lev_histLES = ',lev_histLES
    2635     write(lunout,*)' lev_histdayNMC = ',lev_histdayNMC
    2636     write(lunout,*)' levout_histNMC = ',levout_histNMC
    2637     write(lunout,*)' ok_histNMC = ',ok_histNMC
    2638     write(lunout,*)' freq_outNMC = ',freq_outNMC
    2639     write(lunout,*)' freq_calNMC = ',freq_calNMC
    2640     write(lunout,*)' iflag_pbl = ', iflag_pbl
     2528    WRITE(lunout,*) ' ##############################################'
     2529    WRITE(lunout,*) ' Configuration des parametres de la physique: '
     2530    WRITE(lunout,*) ' Type ocean = ', type_ocean
     2531    WRITE(lunout,*) ' Version ocean = ', version_ocean
     2532    WRITE(lunout,*) ' Config veget = ', ok_veget,type_veget
     2533    WRITE(lunout,*) ' Snow model SISVAT : ok_snow = ', ok_snow
     2534    WRITE(lunout,*) ' Config xml pour XIOS : ok_all_xml = ', ok_all_xml
     2535    WRITE(lunout,*) ' Sortie journaliere = ', ok_journe
     2536    WRITE(lunout,*) ' Sortie haute frequence = ', ok_hf
     2537    WRITE(lunout,*) ' Sortie mensuelle = ', ok_mensuel
     2538    WRITE(lunout,*) ' Sortie instantanee = ', ok_instan
     2539    WRITE(lunout,*) ' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
     2540    WRITE(lunout,*) ' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
     2541    WRITE(lunout,*) ' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
     2542    WRITE(lunout,*) ' Frequence appel simulateur AIRS, freq_AIRS =', freq_AIRS
     2543    WRITE(lunout,*) ' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
     2544    WRITE(lunout,*) ' Excentricite = ',R_ecc
     2545    WRITE(lunout,*) ' Equinoxe = ',R_peri
     2546    WRITE(lunout,*) ' Inclinaison =',R_incl
     2547    WRITE(lunout,*) ' Constante solaire =',solaire
     2548    WRITE(lunout,*) ' ok_suntime_rrtm =',ok_suntime_rrtm
     2549    WRITE(lunout,*) ' co2_ppm =',co2_ppm
     2550    WRITE(lunout,*) ' RCO2_act = ',RCO2_act
     2551    WRITE(lunout,*) ' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
     2552    WRITE(lunout,*) ' N2O_ppb =',N2O_ppb,' RN2O_act=  ',RN2O_act
     2553    WRITE(lunout,*) ' CFC11_ppt=',CFC11_ppt,' RCFC11_act=  ',RCFC11_act
     2554    WRITE(lunout,*) ' CFC12_ppt=',CFC12_ppt,' RCFC12_act=  ',RCFC12_act
     2555    WRITE(lunout,*) ' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
     2556    WRITE(lunout,*) ' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
     2557    WRITE(lunout,*) ' RCFC12_per = ',RCFC12_per
     2558    WRITE(lunout,*) ' cvl_comp_threshold=', cvl_comp_threshold
     2559    WRITE(lunout,*) ' cvl_sig2feed=', cvl_sig2feed
     2560    WRITE(lunout,*) ' cvl_corr=', cvl_corr
     2561    WRITE(lunout,*) ' ok_lic_melt=', ok_lic_melt
     2562    WRITE(lunout,*) ' ok_lic_cond=', ok_lic_cond
     2563    WRITE(lunout,*) ' iflag_cycle_diurne=',iflag_cycle_diurne
     2564    WRITE(lunout,*) ' soil_model=',soil_model
     2565    WRITE(lunout,*) ' new_oliq=',new_oliq
     2566    WRITE(lunout,*) ' ok_orodr=',ok_orodr
     2567    WRITE(lunout,*) ' ok_orolf=',ok_orolf
     2568    WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai
     2569    WRITE(lunout,*) ' nbapp_rad=',nbapp_rad
     2570    WRITE(lunout,*) ' iflag_con=',iflag_con
     2571    WRITE(lunout,*) ' nbapp_cv=',nbapp_cv
     2572    WRITE(lunout,*) ' nbapp_wk=',nbapp_wk
     2573    WRITE(lunout,*) ' iflag_ener_conserv=',iflag_ener_conserv
     2574    WRITE(lunout,*) ' ok_conserv_q=',ok_conserv_q
     2575    WRITE(lunout,*) ' iflag_fisrtilp_qsat=',iflag_fisrtilp_qsat
     2576    WRITE(lunout,*) ' iflag_bergeron=',iflag_bergeron
     2577    WRITE(lunout,*) ' epmax = ', epmax
     2578    WRITE(lunout,*) ' coef_epmax_cape = ', coef_epmax_cape
     2579    WRITE(lunout,*) ' ok_adj_ema = ', ok_adj_ema
     2580    WRITE(lunout,*) ' iflag_clw = ', iflag_clw
     2581    WRITE(lunout,*) ' cld_lc_lsc = ', cld_lc_lsc
     2582    WRITE(lunout,*) ' cld_lc_con = ', cld_lc_con
     2583    WRITE(lunout,*) ' cld_tau_lsc = ', cld_tau_lsc
     2584    WRITE(lunout,*) ' cld_tau_con = ', cld_tau_con
     2585    WRITE(lunout,*) ' ffallv_lsc = ', ffallv_lsc
     2586    WRITE(lunout,*) ' ffallv_con = ', ffallv_con
     2587    WRITE(lunout,*) ' coef_eva = ', coef_eva
     2588    WRITE(lunout,*) ' reevap_ice = ', reevap_ice
     2589    WRITE(lunout,*) ' iflag_pdf = ', iflag_pdf
     2590    WRITE(lunout,*) ' iflag_cld_th = ', iflag_cld_th
     2591    WRITE(lunout,*) ' iflag_cld_cv = ', iflag_cld_cv
     2592    WRITE(lunout,*) ' tau_cld_cv = ', tau_cld_cv
     2593    WRITE(lunout,*) ' coefw_cld_cv = ', coefw_cld_cv
     2594    WRITE(lunout,*) ' iflag_radia = ', iflag_radia
     2595    WRITE(lunout,*) ' iflag_rrtm = ', iflag_rrtm
     2596    WRITE(lunout,*) ' NSW = ', NSW
     2597    WRITE(lunout,*) ' iflag_albedo = ', iflag_albedo !albedo SB
     2598    WRITE(lunout,*) ' ok_chlorophyll =',ok_chlorophyll ! albedo SB
     2599    WRITE(lunout,*) ' iflag_ratqs = ', iflag_ratqs
     2600    WRITE(lunout,*) ' seuil_inversion = ', seuil_inversion
     2601    WRITE(lunout,*) ' fact_cldcon = ', fact_cldcon
     2602    WRITE(lunout,*) ' facttemps = ', facttemps
     2603    WRITE(lunout,*) ' ok_newmicro = ',ok_newmicro
     2604    WRITE(lunout,*) ' ratqsbas = ',ratqsbas
     2605    WRITE(lunout,*) ' ratqshaut = ',ratqshaut
     2606    WRITE(lunout,*) ' tau_ratqs = ',tau_ratqs
     2607    WRITE(lunout,*) ' top_height = ',top_height
     2608    WRITE(lunout,*) ' rad_froid = ',rad_froid
     2609    WRITE(lunout,*) ' rad_chau1 = ',rad_chau1
     2610    WRITE(lunout,*) ' rad_chau2 = ',rad_chau2
     2611    WRITE(lunout,*) ' t_glace_min = ',t_glace_min
     2612    WRITE(lunout,*) ' t_glace_max = ',t_glace_max
     2613    WRITE(lunout,*) ' exposant_glace = ',exposant_glace
     2614    WRITE(lunout,*) ' iflag_t_glace = ',iflag_t_glace
     2615    WRITE(lunout,*) ' iflag_cloudth_vert = ',iflag_cloudth_vert
     2616    WRITE(lunout,*) ' iflag_rain_incloud_vol = ',iflag_rain_incloud_vol
     2617    WRITE(lunout,*) ' iflag_ice_thermo = ',iflag_ice_thermo
     2618    WRITE(lunout,*) ' rei_min = ',rei_min
     2619    WRITE(lunout,*) ' rei_max = ',rei_max
     2620    WRITE(lunout,*) ' overlap = ',overlap
     2621    WRITE(lunout,*) ' cdmmax = ',cdmmax
     2622    WRITE(lunout,*) ' cdhmax = ',cdhmax
     2623    WRITE(lunout,*) ' ksta = ',ksta
     2624    WRITE(lunout,*) ' ksta_ter = ',ksta_ter
     2625    WRITE(lunout,*) ' f_ri_cd_min = ',f_ri_cd_min
     2626    WRITE(lunout,*) ' ok_kzmin = ',ok_kzmin
     2627    WRITE(lunout,*) ' pbl_lmixmin_alpha = ',pbl_lmixmin_alpha
     2628    WRITE(lunout,*) ' fmagic = ',fmagic
     2629    WRITE(lunout,*) ' pmagic = ',pmagic
     2630    WRITE(lunout,*) ' ok_ade = ',ok_ade
     2631    WRITE(lunout,*) ' ok_aie = ',ok_aie
     2632    WRITE(lunout,*) ' ok_alw = ',ok_alw
     2633    WRITE(lunout,*) ' aerosol_couple = ', aerosol_couple
     2634    WRITE(lunout,*) ' chemistry_couple = ', chemistry_couple
     2635    WRITE(lunout,*) ' flag_aerosol = ', flag_aerosol
     2636    WRITE(lunout,*) ' flag_aerosol_strat= ', flag_aerosol_strat
     2637    WRITE(lunout,*) ' new_aod = ', new_aod
     2638    WRITE(lunout,*) ' aer_type = ',aer_type
     2639    WRITE(lunout,*) ' bl95_b0 = ',bl95_b0
     2640    WRITE(lunout,*) ' bl95_b1 = ',bl95_b1
     2641    WRITE(lunout,*) ' lev_histhf = ',lev_histhf
     2642    WRITE(lunout,*) ' lev_histday = ',lev_histday
     2643    WRITE(lunout,*) ' lev_histmth = ',lev_histmth
     2644    WRITE(lunout,*) ' lev_histins = ',lev_histins
     2645    WRITE(lunout,*) ' lev_histLES = ',lev_histLES
     2646    WRITE(lunout,*) ' lev_histdayNMC = ',lev_histdayNMC
     2647    WRITE(lunout,*) ' levout_histNMC = ',levout_histNMC
     2648    WRITE(lunout,*) ' ok_histNMC = ',ok_histNMC
     2649    WRITE(lunout,*) ' freq_outNMC = ',freq_outNMC
     2650    WRITE(lunout,*) ' freq_calNMC = ',freq_calNMC
     2651    WRITE(lunout,*) ' iflag_pbl = ', iflag_pbl
    26412652!FC
    2642     write(lunout,*)' ifl_pbltree = ', ifl_pbltree
    2643     write(lunout,*)' Cd_frein = ', Cd_frein
    2644     write(lunout,*)' iflag_pbl_split = ', iflag_pbl_split
    2645     write(lunout,*)' iflag_order2_sollw = ', iflag_order2_sollw
    2646     write(lunout,*)' iflag_thermals = ', iflag_thermals
    2647     write(lunout,*)' iflag_thermals_ed = ', iflag_thermals_ed
    2648     write(lunout,*)' fact_thermals_ed_dz = ', fact_thermals_ed_dz
    2649     write(lunout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux
    2650     write(lunout,*)' iflag_thermals_closure = ', iflag_thermals_closure
    2651     write(lunout,*)' iflag_clos = ', iflag_clos
    2652     write(lunout,*)' coef_clos_ls = ', coef_clos_ls
    2653     write(lunout,*)' type_run = ',type_run
    2654     write(lunout,*)' ok_cosp = ',ok_cosp
    2655     write(lunout,*)' ok_airs = ',ok_airs
    2656 
    2657     write(lunout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP
    2658     write(lunout,*)' ok_journeCOSP = ',ok_journeCOSP
    2659     write(lunout,*)' ok_hfCOSP =',ok_hfCOSP
    2660     write(lunout,*)' solarlong0 = ', solarlong0
    2661     write(lunout,*)' qsol0 = ', qsol0
    2662     write(lunout,*)' evap0 = ', evap0
    2663     write(lunout,*)' albsno0 = ', albsno0
    2664     write(lunout,*)' iflag_sic = ', iflag_sic
    2665     write(lunout,*)' inertie_sol = ', inertie_sol
    2666     write(lunout,*)' inertie_sic = ', inertie_sic
    2667     write(lunout,*)' inertie_lic = ', inertie_lic
    2668     write(lunout,*)' inertie_sno = ', inertie_sno
    2669     write(lunout,*)' f_cdrag_ter = ',f_cdrag_ter
    2670     write(lunout,*)' f_cdrag_oce = ',f_cdrag_oce
    2671     write(lunout,*)' f_rugoro = ',f_rugoro
    2672     write(lunout,*)' z0min = ',z0min
    2673     write(lunout,*)' supcrit1 = ', supcrit1
    2674     write(lunout,*)' supcrit2 = ', supcrit2
    2675     write(lunout,*)' iflag_mix = ', iflag_mix
    2676     write(lunout,*)' iflag_mix_adiab = ', iflag_mix_adiab
    2677     write(lunout,*)' scut = ', scut
    2678     write(lunout,*)' qqa1 = ', qqa1
    2679     write(lunout,*)' qqa2 = ', qqa2
    2680     write(lunout,*)' gammas = ', gammas
    2681     write(lunout,*)' Fmax = ', Fmax
    2682     write(lunout,*)' tmax_fonte_cv = ', tmax_fonte_cv
    2683     write(lunout,*)' alphas = ', alphas
    2684     write(lunout,*)' iflag_wake = ', iflag_wake
    2685     write(lunout,*)' alp_offset = ', alp_offset
     2653    WRITE(lunout,*) ' ifl_pbltree = ', ifl_pbltree
     2654    WRITE(lunout,*) ' Cd_frein = ', Cd_frein
     2655    WRITE(lunout,*) ' iflag_pbl_split = ', iflag_pbl_split
     2656    WRITE(lunout,*) ' iflag_order2_sollw = ', iflag_order2_sollw
     2657    WRITE(lunout,*) ' iflag_thermals = ', iflag_thermals
     2658    WRITE(lunout,*) ' iflag_thermals_ed = ', iflag_thermals_ed
     2659    WRITE(lunout,*) ' fact_thermals_ed_dz = ', fact_thermals_ed_dz
     2660    WRITE(lunout,*) ' iflag_thermals_optflux = ', iflag_thermals_optflux
     2661    WRITE(lunout,*) ' iflag_thermals_closure = ', iflag_thermals_closure
     2662    WRITE(lunout,*) ' iflag_clos = ', iflag_clos
     2663    WRITE(lunout,*) ' coef_clos_ls = ', coef_clos_ls
     2664    WRITE(lunout,*) ' type_run = ',type_run
     2665    WRITE(lunout,*) ' ok_cosp = ',ok_cosp
     2666    WRITE(lunout,*) ' ok_airs = ',ok_airs
     2667
     2668    WRITE(lunout,*) ' ok_mensuelCOSP = ',ok_mensuelCOSP
     2669    WRITE(lunout,*) ' ok_journeCOSP = ',ok_journeCOSP
     2670    WRITE(lunout,*) ' ok_hfCOSP =',ok_hfCOSP
     2671    WRITE(lunout,*) ' solarlong0 = ', solarlong0
     2672    WRITE(lunout,*) ' qsol0 = ', qsol0
     2673    WRITE(lunout,*) ' evap0 = ', evap0
     2674    WRITE(lunout,*) ' albsno0 = ', albsno0
     2675    WRITE(lunout,*) ' iflag_sic = ', iflag_sic
     2676    WRITE(lunout,*) ' inertie_sol = ', inertie_sol
     2677    WRITE(lunout,*) ' inertie_sic = ', inertie_sic
     2678    WRITE(lunout,*) ' inertie_lic = ', inertie_lic
     2679    WRITE(lunout,*) ' inertie_sno = ', inertie_sno
     2680    WRITE(lunout,*) ' f_cdrag_ter = ',f_cdrag_ter
     2681    WRITE(lunout,*) ' f_cdrag_oce = ',f_cdrag_oce
     2682    WRITE(lunout,*) ' f_rugoro = ',f_rugoro
     2683    WRITE(lunout,*) ' z0min = ',z0min
     2684    WRITE(lunout,*) ' supcrit1 = ', supcrit1
     2685    WRITE(lunout,*) ' supcrit2 = ', supcrit2
     2686    WRITE(lunout,*) ' iflag_mix = ', iflag_mix
     2687    WRITE(lunout,*) ' iflag_mix_adiab = ', iflag_mix_adiab
     2688    WRITE(lunout,*) ' scut = ', scut
     2689    WRITE(lunout,*) ' qqa1 = ', qqa1
     2690    WRITE(lunout,*) ' qqa2 = ', qqa2
     2691    WRITE(lunout,*) ' gammas = ', gammas
     2692    WRITE(lunout,*) ' Fmax = ', Fmax
     2693    WRITE(lunout,*) ' tmax_fonte_cv = ', tmax_fonte_cv
     2694    WRITE(lunout,*) ' alphas = ', alphas
     2695    WRITE(lunout,*) ' iflag_wake = ', iflag_wake
     2696    WRITE(lunout,*) ' alp_offset = ', alp_offset
    26862697    ! nrlmd le 10/04/2012
    2687     write(lunout,*)' iflag_trig_bl = ', iflag_trig_bl
    2688     write(lunout,*)' s_trig = ', s_trig
    2689     write(lunout,*)' tau_trig_shallow = ', tau_trig_shallow
    2690     write(lunout,*)' tau_trig_deep = ', tau_trig_deep
    2691     write(lunout,*)' iflag_clos_bl = ', iflag_clos_bl
     2698    WRITE(lunout,*) ' iflag_trig_bl = ', iflag_trig_bl
     2699    WRITE(lunout,*) ' s_trig = ', s_trig
     2700    WRITE(lunout,*) ' tau_trig_shallow = ', tau_trig_shallow
     2701    WRITE(lunout,*) ' tau_trig_deep = ', tau_trig_deep
     2702    WRITE(lunout,*) ' iflag_clos_bl = ', iflag_clos_bl
    26922703    ! fin nrlmd le 10/04/2012
    26932704
    2694     write(lunout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
     2705    WRITE(lunout,*) ' lonmin lonmax latmin latmax bilKP_ins =',&
    26952706         lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
    2696     write(lunout,*)' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
     2707    WRITE(lunout,*) ' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
    26972708         ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
    26982709
    2699     write(lunout,*) 'ok_strato = ', ok_strato
    2700     write(lunout,*) 'ok_hines = ',  ok_hines
    2701     write(lunout,*) 'ok_gwd_rando = ',  ok_gwd_rando
    2702     write(lunout,*) 'ok_qch4 = ',  ok_qch4
    2703     write(lunout,*) 'gwd_rando_ruwmax = ', gwd_rando_ruwmax
    2704     write(lunout,*) 'gwd_rando_sat = ', gwd_rando_sat
    2705     write(lunout,*) 'gwd_front_ruwmax = ', gwd_front_ruwmax
    2706     write(lunout,*) 'gwd_front_sat = ', gwd_front_sat
    2707     write(lunout,*) 'SSO gkdrag =',gkdrag
    2708     write(lunout,*) 'SSO grahilo=',grahilo
    2709     write(lunout,*) 'SSO grcrit=',grcrit
    2710     write(lunout,*) 'SSO gfrcrit=',gfrcrit
    2711     write(lunout,*) 'SSO gkwake=',gkwake
    2712     write(lunout,*) 'SSO gklift=',gklift
    2713     write(lunout,*) 'adjust_tropopause = ', adjust_tropopause
    2714     write(lunout,*) 'ok_daily_climoz = ',ok_daily_climoz
    2715     write(lunout,*) 'read_climoz = ', read_climoz
    2716     write(lunout,*) 'carbon_cycle_tr = ', carbon_cycle_tr
    2717     write(lunout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl
     2710    WRITE(lunout,*) ' ok_strato = ', ok_strato
     2711    WRITE(lunout,*) ' ok_hines = ',  ok_hines
     2712    WRITE(lunout,*) ' ok_gwd_rando = ',  ok_gwd_rando
     2713    WRITE(lunout,*) ' ok_qch4 = ',  ok_qch4
     2714    WRITE(lunout,*) ' gwd_rando_ruwmax = ', gwd_rando_ruwmax
     2715    WRITE(lunout,*) ' gwd_rando_sat = ', gwd_rando_sat
     2716    WRITE(lunout,*) ' gwd_front_ruwmax = ', gwd_front_ruwmax
     2717    WRITE(lunout,*) ' gwd_front_sat = ', gwd_front_sat
     2718    WRITE(lunout,*) ' SSO gkdrag =',gkdrag
     2719    WRITE(lunout,*) ' SSO grahilo=',grahilo
     2720    WRITE(lunout,*) ' SSO grcrit=',grcrit
     2721    WRITE(lunout,*) ' SSO gfrcrit=',gfrcrit
     2722    WRITE(lunout,*) ' SSO gkwake=',gkwake
     2723    WRITE(lunout,*) ' SSO gklift=',gklift
     2724    WRITE(lunout,*) ' adjust_tropopause = ', adjust_tropopause
     2725    WRITE(lunout,*) ' ok_daily_climoz = ',ok_daily_climoz
     2726    WRITE(lunout,*) ' read_climoz = ', read_climoz
     2727    WRITE(lunout,*) ' carbon_cycle_tr = ', carbon_cycle_tr
     2728    WRITE(lunout,*) ' carbon_cycle_cpl = ', carbon_cycle_cpl
     2729    WRITE(lunout,*) ' level_coupling_esm = ', level_coupling_esm
    27182730
    27192731    !$OMP END MASTER
    … …  
    27362748  !
    27372749  ! tau_calv:    temps de relaxation pour la fonte des glaciers
    2738 
     2750  !
    27392751  REAL          :: tau_calv
    2740   REAL,SAVE     :: tau_calv_omp
    2741 
     2752  REAL, SAVE    :: tau_calv_omp
    27422753  !
    27432754  !Config Key  = tau_calv
    … …  
    27512762  !$OMP END MASTER
    27522763  !$OMP BARRIER
    2753 
     2764  !
    27542765  tau_calv=tau_calv_omp
    2755 
     2766  !
    27562767  !$OMP MASTER
    2757   write(lunout,*)' ##############################################'
     2768  WRITE(lunout,*)' ##############################################'
    27582769  WRITE(lunout,*)' Configuration de l''interface atm/surfaces  : '
    27592770  WRITE(lunout,*)' tau_calv = ',tau_calv
    27602771  !$OMP END MASTER
    2761 
     2772  !
    27622773  RETURN
    27632774
Note: See TracChangeset for help on using the changeset viewer.