Changeset 4250


Ignore:
Timestamp:
Sep 16, 2022, 11:32:38 AM (20 months ago)
Author:
oboucher
Message:

added co2_ppm0 in def files (for ESM version only)
cleaned up the radiation inputs

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r4245 r4250  
    137137    INTEGER, SAVE       :: iflag_gusts_omp,iflag_z0_oce_omp
    138138
    139     ! Local
    140     REAL                 :: zzz
    141 
    142139    REAL :: seuil_inversion
    143140    REAL,SAVE :: seuil_inversion_omp
    … …  
    162159
    163160    REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp
    164     REAL,SAVE      :: solaire_omp_init
     161    REAL,SAVE :: solaire_omp_init
    165162    LOGICAL,SAVE :: ok_suntime_rrtm_omp
    166163    REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp
     164    REAL,SAVE :: co2_ppm0_omp
    167165    REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp
    168166    REAL,SAVE :: N2O_ppb_omp, RN2O_omp, N2O_ppb_per_omp, RN2O_per_omp
    … …  
    746744    !
    747745    !Config Key  = co2_ppm
    748     !Config Desc = concentration du gaz carbonique en ppmv
     746    !Config Desc = concentration du CO2 en ppmv
    749747    !Config Def  = 348.
    750748    !Config Help =
    751749    !               
    752     !
    753750    !valeur AMIP II
    754751    co2_ppm_omp = 348.
    755752    CALL getin('co2_ppm', co2_ppm_omp)
    756753    !
    757     !Config Key  = RCO2
    758     !Config Desc = Concentration du CO2
    759     !Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
    760     !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
    761762    !Config Help =
    762763    !               
    763     ! RCO2 = 5.286789092164308E-04
    764     !ancienne valeur
    765     RCO2_omp = co2_ppm_omp * 1.0e-06 * RMCO2 / RMD ! pour co2_ppm=348.
    766 
    767     !  CALL getin('RCO2', RCO2)
     764    co2_ppm0_omp = 284.32
     765    CALL getin('co2_ppm0', co2_ppm0_omp)
    768766    !
    769767    !Config Key  = RCH4
    … …  
    772770    !Config Help =
    773771    !               
    774     !
    775     !valeur AMIP II
    776     !OK  RCH4 = 1.65E-06* 16.043/28.97
    777     ! RCH4 = 9.137366240938903E-07
    778     !
    779     !ancienne valeur
    780     ! RCH4 = 1.72E-06* 16.043/28.97
    781     !OK CALL getin('RCH4', RCH4)
    782     zzz = 1650.
    783     CALL getin('CH4_ppb', zzz)
    784     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
    785775    RCH4_omp = CH4_ppb_omp * 1.0E-09 * RMCH4 / RMD
    786776    !
    … …  
    790780    !Config Help =
    791781    !               
    792     !
    793     !valeur AMIP II
    794     !OK  RN2O = 306.E-09* 44.013/28.97
    795     ! RN2O = 4.648939592682085E-07
    796     !
    797     !ancienne valeur
    798     ! RN2O = 310.E-09* 44.013/28.97
    799     !OK  CALL getin('RN2O', RN2O)
    800     zzz=306.
    801     CALL getin('N2O_ppb', zzz)
    802     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
    803785    RN2O_omp = N2O_ppb_omp * 1.0E-09 * RMN2O / RMD
    804786    !
    … …  
    808790    !Config Help =
    809791    !               
    810     !
    811     !OK RCFC11 = 280.E-12* 137.3686/28.97
    812     zzz = 280.
    813     CALL getin('CFC11_ppt',zzz)
    814     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
    815795    RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * RMCFC11 / RMD
    816     ! RCFC11 = 1.327690990680013E-09
    817     !OK CALL getin('RCFC11', RCFC11)
    818796    !
    819797    !Config Key  = RCFC12
    … …  
    822800    !Config Help =
    823801    !               
    824     !
    825     !OK RCFC12 = 484.E-12* 120.9140/28.97
    826     zzz = 484.
    827     CALL getin('CFC12_ppt',zzz)
    828     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
    829805    RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * RMCFC12 / RMD
    830     ! RCFC12 = 2.020102726958923E-09
    831     !OK CALL getin('RCFC12', RCFC12)
    832 
    833     !ajout CFMIP begin
     806
    834807    !
    835808    !Config Key  = co2_ppm_per
    836     !Config Desc = concentration du co2_ppm_per
     809    !Config Desc = concentration du CO2 perturbé en ppmv (CFMIP)
    837810    !Config Def  = 348.
    838811    !Config Help =
    … …  
    840813    co2_ppm_per_omp = co2_ppm_omp
    841814    CALL getin('co2_ppm_per', co2_ppm_per_omp)
    842     !
    843     !Config Key  = RCO2_per
    844     !Config Desc = Concentration du CO2_per
    845     !Config Def  = co2_ppm_per * 1.0e-06  * 44.011/28.97
    846     !Config Def  = 348. * 1.0e-06  * 44.011/28.97
    847     !Config Help =
    848     !               
     815    !conversion en rapport de mélange massique
    849816    RCO2_per_omp = co2_ppm_per_omp * 1.0e-06 * RMCO2 / RMD
    850 
    851     !Config Key  = ok_4xCO2atm
    852     !Config Desc = Calcul ou non effet radiatif 4xco2
    853     !Config Def  = .FALSE.
    854     !Config Help =
    855817
    856818    !Config Key  = RCH4_per
    … …  
    859821    !Config Help =
    860822    !               
    861     zzz = CH4_ppb_omp
    862     CALL getin('CH4_ppb_per', zzz)
    863     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
    864826    RCH4_per_omp = CH4_ppb_per_omp * 1.0E-09 * RMCH4 / RMD
    865827    !
    … …  
    869831    !Config Help =
    870832    !               
    871     zzz = N2O_ppb_omp
    872     CALL getin('N2O_ppb_per', zzz)
    873     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
    874836    RN2O_per_omp = N2O_ppb_per_omp * 1.0E-09 * RMN2O / RMD
    875837    !
    … …  
    879841    !Config Help =
    880842    !               
    881     zzz = CFC11_ppt_omp
    882     CALL getin('CFC11_ppt_per',zzz)
    883     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
    884846    RCFC11_per_omp=CFC11_ppt_per_omp* 1.0E-12 * RMCFC11 / RMD
    885847    !
    … …  
    889851    !Config Help =
    890852    !               
    891     zzz = CFC12_ppt_omp
    892     CALL getin('CFC12_ppt_per',zzz)
    893     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
    894856    RCFC12_per_omp = CFC12_ppt_per_omp * 1.0E-12 * RMCFC12 / RMD
    895     !ajout CFMIP end
    896857
    897858    !
    … …  
    24392400    ok_suntime_rrtm = ok_suntime_rrtm_omp
    24402401    co2_ppm = co2_ppm_omp
     2402    co2_ppm0 = co2_ppm0_omp
    24412403    RCO2 = RCO2_omp
    24422404    CH4_ppb = CH4_ppb_omp
    … …  
    28702832    WRITE(lunout,*) ' ok_suntime_rrtm =',ok_suntime_rrtm
    28712833    WRITE(lunout,*) ' co2_ppm =',co2_ppm
     2834    WRITE(lunout,*) ' co2_ppm0 =',co2_ppm0
    28722835    WRITE(lunout,*) ' RCO2_act = ',RCO2_act
    28732836    WRITE(lunout,*) ' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
Note: See TracChangeset for help on using the changeset viewer.