Ignore:
Timestamp:
Mar 29, 2023, 3:14:27 PM (18 months ago)
Author:
lguez
Message:

Sync latest trunk changes to branch LMDZ_ECRad

Location:
LMDZ6/branches/LMDZ_ECRad
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/LMDZ_ECRad

  • LMDZ6/branches/LMDZ_ECRad/libf/phylmd/conf_phys_m.F90

    r4114 r4482  
    133133    REAL, SAVE          :: f_rugoro_omp   , z0min_omp
    134134    REAL, SAVE          :: z0m_seaice_omp,z0h_seaice_omp
     135    REAL, SAVE          :: z0m_landice_omp,z0h_landice_omp
    135136    REAL, SAVE          :: min_wind_speed_omp,f_gust_wk_omp,f_gust_bl_omp,f_qsat_oce_omp, f_z0qh_oce_omp
    136137    INTEGER, SAVE       :: iflag_gusts_omp,iflag_z0_oce_omp
    137 
    138     ! Local
    139     REAL                 :: zzz
    140138
    141139    REAL :: seuil_inversion
    … …  
    161159
    162160    REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp
    163     REAL,SAVE      :: solaire_omp_init
     161    REAL,SAVE :: solaire_omp_init
    164162    LOGICAL,SAVE :: ok_suntime_rrtm_omp
    165163    REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp
     164    REAL,SAVE :: co2_ppm0_omp
    166165    REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp
    167166    REAL,SAVE :: N2O_ppb_omp, RN2O_omp, N2O_ppb_per_omp, RN2O_per_omp
    … …  
    186185    REAL,SAVE :: t_glace_min_omp, t_glace_max_omp
    187186    REAL,SAVE :: exposant_glace_omp
    188     INTEGER,SAVE :: iflag_gammasat_omp, iflag_mpc_bl_omp
     187    INTEGER,SAVE :: iflag_gammasat_omp
    189188    REAL,SAVE :: rei_min_omp, rei_max_omp
    190189    INTEGER,SAVE :: iflag_sic_omp, iflag_inertie_omp
    … …  
    199198    REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp,f_ri_cd_min_omp
    200199    LOGICAL,SAVE :: ok_kzmin_omp
    201     REAL, SAVE   :: pbl_lmixmin_alpha_omp
    202200    REAL, SAVE ::  fmagic_omp, pmagic_omp
    203201    INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
    … …  
    227225    LOGICAL, SAVE :: ok_lic_cond_omp
    228226    !
     227    REAL, SAVE    :: zrel_oro_t_omp, zstd_orodr_t_omp
     228    REAL, SAVE    :: zpmm_orodr_t_omp, zpmm_orolf_t_omp
    229229    INTEGER, SAVE :: iflag_cycle_diurne_omp
    230230    LOGICAL, SAVE :: soil_model_omp,new_oliq_omp
    … …  
    255255    LOGICAL, SAVE :: ok_new_lscp_omp
    256256    LOGICAL, SAVE :: ok_icefra_lscp_omp
    257     LOGICAL, SAVE :: ok_radliq_snow_omp
    258257
    259258
    … …  
    745744    !
    746745    !Config Key  = co2_ppm
    747     !Config Desc = concentration du gaz carbonique en ppmv
     746    !Config Desc = concentration du CO2 en ppmv
    748747    !Config Def  = 348.
    749748    !Config Help =
    750749    !               
    751     !
    752750    !valeur AMIP II
    753751    co2_ppm_omp = 348.
    754752    CALL getin('co2_ppm', co2_ppm_omp)
    755753    !
    756     !Config Key  = RCO2
    757     !Config Desc = Concentration du CO2
    758     !Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
    759     !Config Def  = 348. * 1.0e-06  * 44.011/28.97
     754    !conversion en rapport de mélange massique
     755    RCO2_omp = co2_ppm_omp * 1.0e-06 * RMCO2 / RMD
     756
     757    !
     758    !Config Key  = co2_ppm0
     759    !Config Desc = concentration initiale du CO2 en ppmv pour la version ESM avec CO2 interactif dans le cas
     760    !              où cette concentration de figure pas dans l'état de redémarrage de la physique
     761    !Config Def  = 284.32
    760762    !Config Help =
    761763    !               
    762     ! RCO2 = 5.286789092164308E-04
    763     !ancienne valeur
    764     RCO2_omp = co2_ppm_omp * 1.0e-06 * RMCO2 / RMD ! pour co2_ppm=348.
    765 
    766     !  CALL getin('RCO2', RCO2)
     764    co2_ppm0_omp = 284.32
     765    CALL getin('co2_ppm0', co2_ppm0_omp)
    767766    !
    768767    !Config Key  = RCH4
    … …  
    771770    !Config Help =
    772771    !               
    773     !
    774     !valeur AMIP II
    775     !OK  RCH4 = 1.65E-06* 16.043/28.97
    776     ! RCH4 = 9.137366240938903E-07
    777     !
    778     !ancienne valeur
    779     ! RCH4 = 1.72E-06* 16.043/28.97
    780     !OK CALL getin('RCH4', RCH4)
    781     zzz = 1650.
    782     CALL getin('CH4_ppb', zzz)
    783     CH4_ppb_omp = zzz
     772    CH4_ppb_omp = 1650.
     773    CALL getin('CH4_ppb', CH4_ppb_omp)
     774    !conversion en rapport de mélange massique
    784775    RCH4_omp = CH4_ppb_omp * 1.0E-09 * RMCH4 / RMD
    785776    !
    … …  
    789780    !Config Help =
    790781    !               
    791     !
    792     !valeur AMIP II
    793     !OK  RN2O = 306.E-09* 44.013/28.97
    794     ! RN2O = 4.648939592682085E-07
    795     !
    796     !ancienne valeur
    797     ! RN2O = 310.E-09* 44.013/28.97
    798     !OK  CALL getin('RN2O', RN2O)
    799     zzz=306.
    800     CALL getin('N2O_ppb', zzz)
    801     N2O_ppb_omp = zzz
     782    N2O_ppb_omp = 306.
     783    CALL getin('N2O_ppb', N2O_ppb_omp)
     784    !conversion en rapport de mélange massique
    802785    RN2O_omp = N2O_ppb_omp * 1.0E-09 * RMN2O / RMD
    803786    !
    … …  
    807790    !Config Help =
    808791    !               
    809     !
    810     !OK RCFC11 = 280.E-12* 137.3686/28.97
    811     zzz = 280.
    812     CALL getin('CFC11_ppt',zzz)
    813     CFC11_ppt_omp = zzz
     792    CFC11_ppt_omp = 280.
     793    CALL getin('CFC11_ppt',CFC11_ppt_omp)
     794    !conversion en rapport de mélange massique
    814795    RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * RMCFC11 / RMD
    815     ! RCFC11 = 1.327690990680013E-09
    816     !OK CALL getin('RCFC11', RCFC11)
    817796    !
    818797    !Config Key  = RCFC12
    … …  
    821800    !Config Help =
    822801    !               
    823     !
    824     !OK RCFC12 = 484.E-12* 120.9140/28.97
    825     zzz = 484.
    826     CALL getin('CFC12_ppt',zzz)
    827     CFC12_ppt_omp = zzz
     802    CFC12_ppt_omp = 484.
     803    CALL getin('CFC12_ppt',CFC12_ppt_omp)
     804    !conversion en rapport de mélange massique
    828805    RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * RMCFC12 / RMD
    829     ! RCFC12 = 2.020102726958923E-09
    830     !OK CALL getin('RCFC12', RCFC12)
    831 
    832     !ajout CFMIP begin
     806
    833807    !
    834808    !Config Key  = co2_ppm_per
    835     !Config Desc = concentration du co2_ppm_per
     809    !Config Desc = concentration du CO2 perturbé en ppmv (CFMIP)
    836810    !Config Def  = 348.
    837811    !Config Help =
    … …  
    839813    co2_ppm_per_omp = co2_ppm_omp
    840814    CALL getin('co2_ppm_per', co2_ppm_per_omp)
    841     !
    842     !Config Key  = RCO2_per
    843     !Config Desc = Concentration du CO2_per
    844     !Config Def  = co2_ppm_per * 1.0e-06  * 44.011/28.97
    845     !Config Def  = 348. * 1.0e-06  * 44.011/28.97
    846     !Config Help =
    847     !               
     815    !conversion en rapport de mélange massique
    848816    RCO2_per_omp = co2_ppm_per_omp * 1.0e-06 * RMCO2 / RMD
    849 
    850     !Config Key  = ok_4xCO2atm
    851     !Config Desc = Calcul ou non effet radiatif 4xco2
    852     !Config Def  = .FALSE.
    853     !Config Help =
    854817
    855818    !Config Key  = RCH4_per
    … …  
    858821    !Config Help =
    859822    !               
    860     zzz = CH4_ppb_omp
    861     CALL getin('CH4_ppb_per', zzz)
    862     CH4_ppb_per_omp = zzz
     823    CH4_ppb_per_omp = CH4_ppb_omp
     824    CALL getin('CH4_ppb_per', CH4_ppb_per_omp)
     825    !conversion en rapport de mélange massique
    863826    RCH4_per_omp = CH4_ppb_per_omp * 1.0E-09 * RMCH4 / RMD
    864827    !
    … …  
    868831    !Config Help =
    869832    !               
    870     zzz = N2O_ppb_omp
    871     CALL getin('N2O_ppb_per', zzz)
    872     N2O_ppb_per_omp = zzz
     833    N2O_ppb_per_omp = N2O_ppb_omp
     834    CALL getin('N2O_ppb_per', N2O_ppb_per_omp)
     835    !conversion en rapport de mélange massique
    873836    RN2O_per_omp = N2O_ppb_per_omp * 1.0E-09 * RMN2O / RMD
    874837    !
    … …  
    878841    !Config Help =
    879842    !               
    880     zzz = CFC11_ppt_omp
    881     CALL getin('CFC11_ppt_per',zzz)
    882     CFC11_ppt_per_omp = zzz
     843    CFC11_ppt_per_omp = CFC11_ppt_omp
     844    CALL getin('CFC11_ppt_per',CFC11_ppt_per_omp)
     845    !conversion en rapport de mélange massique
    883846    RCFC11_per_omp=CFC11_ppt_per_omp* 1.0E-12 * RMCFC11 / RMD
    884847    !
    … …  
    888851    !Config Help =
    889852    !               
    890     zzz = CFC12_ppt_omp
    891     CALL getin('CFC12_ppt_per',zzz)
    892     CFC12_ppt_per_omp = zzz
     853    CFC12_ppt_per_omp = CFC12_ppt_omp
     854    CALL getin('CFC12_ppt_per',CFC12_ppt_per_omp)
     855    !conversion en rapport de mélange massique
    893856    RCFC12_per_omp = CFC12_ppt_per_omp * 1.0E-12 * RMCFC12 / RMD
    894     !ajout CFMIP end
    895857
    896858    !
    … …  
    937899    ok_orolf_omp = .TRUE.
    938900    CALL getin('ok_orolf', ok_orolf_omp)
     901
     902
     903    !Config  Key  =  zrel_oro_t
     904    !Config  Desc = zrel_oro_t
     905    !Config  Def  = 9999.
     906    !Config  Help = Connais pas !
     907    zrel_oro_t_omp = 9999.
     908    CALL getin('zrel_oro_t', zrel_oro_t_omp)
     909
     910    !Config  Key  =  zstd_orodr_t
     911    !Config  Desc = zstd_orodr_t
     912    !Config  Def  = 0.
     913    !Config  Help = Connais pas !
     914    zstd_orodr_t_omp = 10.
     915    CALL getin('zstd_orodr_t', zstd_orodr_t_omp)
     916
     917
     918    !Config  Key  =  zpmm_orodr_t
     919    !Config  Desc = zpmm_orodr_t
     920    !Config  Def  = 0.
     921    !Config  Help = Connais pas !
     922    zpmm_orodr_t_omp = 100.
     923    CALL getin('zpmm_orodr_t', zpmm_orodr_t_omp)
     924
     925
     926    !Config  Key  =  zpmm_orolf_t
     927    !Config  Desc = zpmm_orolf_t
     928    !Config  Def  = 0.
     929    !Config  Help = Connais pas !
     930    zpmm_orolf_t_omp = 100.
     931    CALL getin('zpmm_orolf_t', zpmm_orolf_t_omp)
    939932
    940933    !Config  Key  = ok_limitvrai
    … …  
    14571450
    14581451    !
    1459     !Config Key  = iflag_mpc_bl
    1460     !Config Desc =
    1461     !Config Def  = 0
    1462     !Config Help =
    1463     !
    1464     iflag_mpc_bl_omp=0
    1465     CALL getin('iflag_mpc_bl',iflag_mpc_bl_omp)
    1466 
    1467 
    1468 
    1469     !
    14701452    !Config Key  = iflag_t_glace
    14711453    !Config Desc = 
    … …  
    16351617    ok_kzmin_omp = .TRUE.
    16361618    CALL getin('ok_kzmin',ok_kzmin_omp)
    1637 
    1638     pbl_lmixmin_alpha_omp=0.0
    1639     CALL getin('pbl_lmixmin_alpha',pbl_lmixmin_alpha_omp)
    16401619
    16411620    !
    … …  
    21232102    z0m_seaice_omp = 0.002 ; CALL getin('z0m_seaice',z0m_seaice_omp)
    21242103    z0h_seaice_omp = 0.002 ; CALL getin('z0h_seaice',z0h_seaice_omp)
     2104
     2105
     2106    z0m_landice_omp = 0.001 ; CALL getin('z0m_landice',z0m_landice_omp)
     2107    z0h_landice_omp = 0.001 ; CALL getin('z0h_landice',z0h_landice_omp)
    21252108
    21262109    f_rugoro_omp = 0.
    … …  
    23762359    !Config Help = ...
    23772360
    2378 
    2379     ok_radliq_snow_omp = .FALSE.
    2380     CALL getin('ok_radliq_snow', ok_radliq_snow_omp)
    2381     !
    2382     !Config Key  = ok_radliq_snow_omp
    2383     !Config Desc = take into account snowfall for radiation
    2384     !Config Def  = .FALSE.
    2385     !Config Help = ...
    23862361
    23872362    ecrit_LES_omp = 1./8.
    … …  
    24342409    ok_suntime_rrtm = ok_suntime_rrtm_omp
    24352410    co2_ppm = co2_ppm_omp
     2411    co2_ppm0 = co2_ppm0_omp
    24362412    RCO2 = RCO2_omp
    24372413    CH4_ppb = CH4_ppb_omp
    … …  
    24592435    ok_orodr = ok_orodr_omp
    24602436    ok_orolf = ok_orolf_omp
     2437    zrel_oro_t=zrel_oro_t_omp
     2438    zstd_orodr_t=zstd_orodr_t_omp
     2439    zpmm_orodr_t=zpmm_orodr_t_omp
     2440    zpmm_orolf_t=zpmm_orolf_t_omp
    24612441    ok_limitvrai = ok_limitvrai_omp
    24622442    nbapp_rad = nbapp_rad_omp
    … …  
    25002480    exposant_glace = exposant_glace_omp
    25012481    iflag_gammasat=iflag_gammasat_omp
    2502     iflag_mpc_bl=iflag_mpc_bl_omp
    25032482    iflag_t_glace = iflag_t_glace_omp
    25042483    iflag_cloudth_vert=iflag_cloudth_vert_omp
    … …  
    25202499    f_ri_cd_min = f_ri_cd_min_omp
    25212500    ok_kzmin = ok_kzmin_omp
    2522     pbl_lmixmin_alpha=pbl_lmixmin_alpha_omp
    25232501    fmagic = fmagic_omp
    25242502    pmagic = pmagic_omp
    … …  
    26662644    z0m_seaice=z0m_seaice_omp
    26672645    z0h_seaice=z0h_seaice_omp
     2646    z0m_landice=z0m_landice_omp
     2647    z0h_landice=z0h_landice_omp
    26682648
    26692649    f_rugoro=f_rugoro_omp
    … …  
    27082688    ok_new_lscp = ok_new_lscp_omp
    27092689    ok_icefra_lscp=ok_icefra_lscp_omp
    2710     ok_radliq_snow=ok_radliq_snow_omp
    27112690    read_fco2_ocean_cor = read_fco2_ocean_cor_omp
    27122691    var_fco2_ocean_cor = var_fco2_ocean_cor_omp
    … …  
    28632842    WRITE(lunout,*) ' ok_suntime_rrtm =',ok_suntime_rrtm
    28642843    WRITE(lunout,*) ' co2_ppm =',co2_ppm
     2844    WRITE(lunout,*) ' co2_ppm0 =',co2_ppm0
    28652845    WRITE(lunout,*) ' RCO2_act = ',RCO2_act
    28662846    WRITE(lunout,*) ' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
    … …  
    28812861    WRITE(lunout,*) ' ok_orodr=',ok_orodr
    28822862    WRITE(lunout,*) ' ok_orolf=',ok_orolf
     2863    WRITE(lunout,*) ' zrel_oro_t=',zrel_oro_t
     2864    WRITE(lunout,*) ' zstd_orodr_t=',zstd_orodr_t
     2865    WRITE(lunout,*) ' zpmm_orodr_t=',zpmm_orodr_t
     2866    WRITE(lunout,*) ' zpmm_orolf_t=',zpmm_orolf_t
    28832867    WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai
    28842868    WRITE(lunout,*) ' nbapp_rad=',nbapp_rad
    … …  
    29292913    WRITE(lunout,*) ' exposant_glace = ',exposant_glace
    29302914    WRITE(lunout,*) ' iflag_gammasat = ',iflag_gammasat
    2931     WRITE(lunout,*) ' iflag_mpc_bl = ',iflag_mpc_bl
    29322915    WRITE(lunout,*) ' iflag_t_glace = ',iflag_t_glace
    29332916    WRITE(lunout,*) ' iflag_cloudth_vert = ',iflag_cloudth_vert
    … …  
    29482931    WRITE(lunout,*) ' f_ri_cd_min = ',f_ri_cd_min
    29492932    WRITE(lunout,*) ' ok_kzmin = ',ok_kzmin
    2950     WRITE(lunout,*) ' pbl_lmixmin_alpha = ',pbl_lmixmin_alpha
    29512933    WRITE(lunout,*) ' fmagic = ',fmagic
    29522934    WRITE(lunout,*) ' pmagic = ',pmagic
    … …  
    30483030    WRITE(lunout,*) ' ok_new_lscp = ', ok_new_lscp
    30493031    WRITE(lunout,*) ' ok_icefra_lscp = ', ok_icefra_lscp
    3050     WRITE(lunout,*) ' ok_radliq_snow = ', ok_radliq_snow
    30513032    WRITE(lunout,*) ' read_climoz = ', read_climoz
    30523033    WRITE(lunout,*) ' carbon_cycle_tr = ', carbon_cycle_tr
Note: See TracChangeset for help on using the changeset viewer.