Ignore:
Timestamp:
May 17, 2021, 3:35:58 PM (3 years ago)
Author:
evignon
Message:

Commission de la nouvelle interface LMDZ-SISVAT
la prochaine commission consistera a supprimer l'ancien repertoire sisvat
et a faire un peu de nettoyage.

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r3876 r3900  
    2727    USE phys_cal_mod
    2828    USE carbon_cycle_mod,  ONLY: carbon_cycle_tr, carbon_cycle_cpl, carbon_cycle_rad, level_coupling_esm
    29     USE carbon_cycle_mod,  ONLY: read_fco2_ocean_cor,var_fco2_ocean_cor
    30     USE carbon_cycle_mod,  ONLY: read_fco2_land_cor,var_fco2_land_cor
    3129    USE mod_grid_phy_lmdz, ONLY: klon_glo
    3230    USE print_control_mod, ONLY: lunout
    … …  
    9088    CHARACTER (len = 8), SAVE  :: aer_type_omp
    9189    INTEGER, SAVE       :: landice_opt_omp
    92     INTEGER, SAVE       :: n_dtis_omp
    93     INTEGER, SAVE       :: iflag_tsurf_inlandsis_omp
    94     INTEGER, SAVE       :: iflag_albzenith_omp   
    95     LOGICAL, SAVE       :: SnoMod_omp,BloMod_omp,ok_outfor_omp
     90    INTEGER, SAVE       :: iflag_tsurf_inlandsis_omp,iflag_temp_inlandsis_omp
     91    INTEGER, SAVE       :: iflag_albcalc_omp,iflag_z0m_snow_omp   
     92    LOGICAL, SAVE       :: SnoMod_omp,BloMod_omp,ok_outfor_omp,ok_zsn_ii_omp
     93    LOGICAL, SAVE       :: discret_xf_omp,opt_runoff_ac_omp 
     94    LOGICAL, SAVE       :: is_ok_slush_omp,is_ok_z0h_rn_omp,is_ok_density_kotlyakov_omp
     95    REAL, SAVE          :: prescribed_z0m_snow_omp,correc_alb_omp
     96    REAL, SAVE          :: buf_sph_pol_omp,buf_siz_pol_omp
    9697    LOGICAL, SAVE       :: ok_newmicro_omp
    9798    LOGICAL, SAVE       :: ok_all_xml_omp
    … …  
    239240    LOGICAL, SAVE :: carbon_cycle_rad_omp
    240241    INTEGER, SAVE :: level_coupling_esm_omp
    241     LOGICAL, SAVE :: read_fco2_ocean_cor_omp
    242     REAL, SAVE    :: var_fco2_ocean_cor_omp
    243     LOGICAL, SAVE :: read_fco2_land_cor_omp
    244     REAL, SAVE    :: var_fco2_land_cor_omp
    245242    LOGICAL, SAVE :: adjust_tropopause_omp
    246243    LOGICAL, SAVE :: ok_daily_climoz_omp
    … …  
    340337
    341338    !Etienne
     339    !Config Key  = iflag_temp_inlandsis
     340    !Config Desc = which method to calculate temp within the soil in INLANDSIS
     341    !Config Def  = 0
     342    iflag_temp_inlandsis_omp = 0
     343    CALL getin('iflag_temp_inlandsis', iflag_temp_inlandsis_omp)
     344
     345    !Etienne
    342346    !Config Key  = iflag_tsurf_inlandsis
    343347    !Config Desc = which method to calculate tsurf in INLANDSIS
    344348    !Config Def  = 0
    345     iflag_tsurf_inlandsis_omp = 0
     349    iflag_tsurf_inlandsis_omp = 1
    346350    CALL getin('iflag_tsurf_inlandsis', iflag_tsurf_inlandsis_omp)
    347351
     352
    348353    !Etienne
    349     !Config Key  = iflag_albzenith
    350     !Config Desc = method to account for albedo sensitivity to solar zenith angle
    351     !Config Def  = 0
    352     iflag_albzenith_omp = 0
    353     CALL getin('iflag_albzenith', iflag_albzenith_omp)
    354 
    355     !Etienne
    356     !Config Key  = n_dtis
    357     !Config Desc = number of subtimesteps for INLANDSIS
    358     !Config Def  = 1
    359     n_dtis_omp = 1
    360     CALL getin('n_dtis', n_dtis_omp)
     354    !Config Key  = iflag_albcalc
     355    !Config Desc = method to calculate snow albedo in INLANDSIS
     356    !Config Def  = 0
     357    iflag_albcalc_omp = 0
     358    CALL getin('iflag_albcalc', iflag_albcalc_omp)
     359
    361360
    362361    !Etienne
    363362    !Config Key  = SnoMod
    364363    !Config Desc = activation of snow modules in inlandsis
    365     !Config Def  = 1
     364    !Config Def  = .TRUE.
    366365    SnoMod_omp = .TRUE.
    367366    CALL getin('SnoMod', SnoMod_omp)
    … …  
    370369    !Config Key  = BloMod
    371370    !Config Desc = activation of blowing snow in inlandsis
    372     !Config Def  = 1
     371    !Config Def  = .FALSE.
    373372    BloMod_omp = .FALSE.
    374373    CALL getin('BloMod', BloMod_omp)
    … …  
    377376    !Config Key  = ok_outfor
    378377    !Config Desc = activation of output ascii file in inlandsis
    379     !Config Def  = 1
     378    !Config Def  = .FALSE.
    380379    ok_outfor_omp = .FALSE.
    381380    CALL getin('ok_outfor', ok_outfor_omp)
     381
     382
     383    !Etienne
     384    !Config Key  = ok_sn_ii
     385    !Config Desc = activation of ice/snow detection
     386    !Config Def  = .TRUE.
     387    ok_zsn_ii = .TRUE.
     388    CALL getin('ok_zsn_ii', ok_zsn_ii_omp)
     389
     390
     391    !Etienne
     392    !Config Key  = discret_xf
     393    !Config Desc = snow discretization following XF
     394    !Config Def  = .TRUE.
     395    discret_xf = .TRUE.
     396    CALL getin('discret_xf', discret_xf_omp)
     397
     398
     399    !Etienne
     400    !Config Key  = is_ok_slush
     401    !Config Desc = activation of the slush option
     402    !Config Def  = .TRUE.
     403    is_ok_slush_omp = .TRUE.
     404    CALL getin('is_ok_slush', is_ok_slush_omp)
     405
     406    !Etienne
     407    !Config Key  = opt_runoff_ac
     408    !Config Desc = option runoff AC
     409    !Config Def  = .TRUE.
     410    opt_runoff_ac_omp = .TRUE.
     411    CALL getin('opt_runoff_ac', opt_runoff_ac_omp)
     412
     413    !Etienne
     414    !Config Key  = is_ok_z0h_rn
     415    !Config Desc = z0h calculation following RN method
     416    !Config Def  = .TRUE.
     417    is_ok_z0h_rn = .TRUE.
     418    CALL getin('is_ok_z0h_rn', is_ok_z0h_rn_omp)
     419
     420
     421    !Etienne
     422    !Config Key  = is_ok_density_kotlyakov
     423    !Config Desc = snow density calculation following kotlyakov
     424    !Config Def  = .FALSE.
     425    is_ok_density_kotlyakov = .FALSE.
     426    CALL getin('is_ok_density_kotlyakov', is_ok_density_kotlyakov_omp)
     427
     428
     429    !Etienne
     430    !Config Key  = prescribed_z0m_snow
     431    !Config Desc = prescribed snow z0m
     432    !Config Def  = 0.005
     433    prescribed_z0m_snow = 0.005
     434    CALL getin('prescribed_z0m_snow', prescribed_z0m_snow_omp)
     435
     436
     437    !Etienne
     438    !Config Key  = iflag_z0m_snow
     439    !Config Desc = method to calculate snow z0m
     440    !Config Def  = 0
     441    iflag_z0m_snow = 0
     442    CALL getin('iflag_z0m_snow', iflag_z0m_snow_omp)
     443
     444
     445    !Etienne
     446    !Config Key  = correc_alb
     447    !Config Desc = correction term for albedo
     448    !Config Def  = 1.01
     449    correc_alb_omp=1.01
     450    CALL getin('correc_alb', correc_alb_omp)
     451
     452
     453    !Etienne
     454    !Config Key  = buf_sph_pol
     455    !Config Desc = sphericity of buffer layer in polar regions
     456    !Config Def  = 99.
     457    buf_sph_pol_omp=99.
     458    CALL getin('buf_sph_pol', buf_sph_pol_omp)
     459
     460    !Etienne
     461    !Config Key  = buf_siz_pol
     462    !Config Desc = grain size of buffer layer in polar regions in e-4m
     463    !Config Def  = 4.
     464    buf_siz_pol_omp=4.
     465    CALL getin('buf_siz_pol', buf_siz_pol_omp)
    382466
    383467    !==================================================================
    … …  
    11401224    CALL getin('coefw_cld_cv',coefw_cld_cv_omp)
    11411225
     1226
     1227
     1228
    11421229    !
    11431230    !Config Key  = iflag_pdf
    … …  
    11481235    iflag_pdf_omp = 0
    11491236    CALL getin('iflag_pdf',iflag_pdf_omp)
    1150 
    11511237    !
    11521238    !Config Key  = fact_cldcon
    … …  
    11751261    ok_newmicro_omp = .TRUE.
    11761262    CALL getin('ok_newmicro',ok_newmicro_omp)
    1177 
    11781263    !
    11791264    !Config Key  = ratqsbas
    … …  
    11841269    ratqsbas_omp = 0.01
    11851270    CALL getin('ratqsbas',ratqsbas_omp)
    1186 
    11871271    !
    11881272    !Config Key  = ratqshaut
    … …  
    22182302    CALL getin('carbon_cycle_rad',carbon_cycle_rad_omp)
    22192303
    2220     read_fco2_ocean_cor_omp=.FALSE.
    2221     CALL getin('read_fco2_ocean_cor',read_fco2_ocean_cor_omp)
    2222 
    2223     var_fco2_ocean_cor_omp=0. ! default value
    2224     CALL getin('var_fco2_ocean_cor',var_fco2_ocean_cor_omp)
    2225 
    2226     read_fco2_land_cor_omp=.FALSE.
    2227     CALL getin('read_fco2_land_cor',read_fco2_land_cor_omp)
    2228 
    2229     var_fco2_land_cor_omp=0. ! default value
    2230     CALL getin('var_fco2_land_cor',var_fco2_land_cor_omp)
    2231 
     2304    ! >> PC
    22322305    ! level_coupling_esm : level of coupling of the biogeochemical fields between LMDZ, ORCHIDEE and NEMO
    22332306    ! Definitions of level_coupling_esm in physiq.def
    … …  
    22422315    level_coupling_esm_omp=0 ! default value
    22432316    CALL getin('level_coupling_esm',level_coupling_esm_omp)
     2317    ! << PC
    22442318
    22452319    !$OMP END MASTER
    … …  
    23582432       ok_veget=.FALSE.
    23592433    ENDIF
    2360     ! SISVAT and INLANDSIS
     2434    ! INLANDSIS
    23612435    !=================================================
    23622436    landice_opt = landice_opt_omp
    23632437    iflag_tsurf_inlandsis = iflag_tsurf_inlandsis_omp
    2364     iflag_albzenith = iflag_albzenith_omp
    2365     n_dtis=n_dtis_omp
     2438    iflag_temp_inlandsis = iflag_temp_inlandsis_omp
     2439    iflag_albcalc = iflag_albcalc_omp
    23662440    SnoMod=SnoMod_omp
    23672441    BloMod=BloMod_omp
    23682442    ok_outfor=ok_outfor_omp
     2443    is_ok_slush=is_ok_slush_omp
     2444    opt_runoff_ac=opt_runoff_ac_omp
     2445    is_ok_z0h_rn=is_ok_z0h_rn_omp
     2446    is_ok_density_kotlyakov=is_ok_density_kotlyakov_omp
     2447    prescribed_z0m_snow=prescribed_z0m_snow_omp
     2448    correc_alb=correc_alb_omp
     2449    iflag_z0m_snow=iflag_z0m_snow_omp
     2450    ok_zsn_ii=ok_zsn_ii_omp
     2451    discret_xf=discret_xf_omp
     2452    buf_sph_pol=buf_sph_pol_omp
     2453    buf_siz_pol=buf_siz_pol_omp
    23692454    !=================================================
    23702455    ok_all_xml = ok_all_xml_omp
    … …  
    25072592    carbon_cycle_rad = carbon_cycle_rad_omp
    25082593    level_coupling_esm = level_coupling_esm_omp
    2509     read_fco2_ocean_cor = read_fco2_ocean_cor_omp
    2510     var_fco2_ocean_cor = var_fco2_ocean_cor_omp
    2511     read_fco2_land_cor = read_fco2_land_cor_omp
    2512     var_fco2_land_cor = var_fco2_land_cor_omp
    25132594
    25142595    ! Test of coherence between type_ocean and version_ocean
    … …  
    28312912    WRITE(lunout,*) ' level_coupling_esm = ', level_coupling_esm
    28322913    WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
    2833     WRITE(lunout,*) ' iflag_albzenith = ', iflag_albzenith
    2834     WRITE(lunout,*) ' n_dtis = ', n_dtis
     2914    WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis
     2915    WRITE(lunout,*) ' iflag_albcalc = ', iflag_albcalc
    28352916    WRITE(lunout,*) ' SnoMod = ', SnoMod
    28362917    WRITE(lunout,*) ' BloMod = ', BloMod
    28372918    WRITE(lunout,*) ' ok_outfor = ', ok_outfor
    2838 
     2919    WRITE(lunout,*) ' is_ok_slush = ', is_ok_slush
     2920    WRITE(lunout,*) ' opt_runoff_ac = ', opt_runoff_ac
     2921    WRITE(lunout,*) ' is_ok_z0h_rn = ', is_ok_z0h_rn
     2922    WRITE(lunout,*) ' is_ok_density_kotlyakov = ', is_ok_density_kotlyakov
     2923    WRITE(lunout,*) ' prescribed_z0m_snow = ', prescribed_z0m_snow
     2924    WRITE(lunout,*) ' iflag_z0m_snow = ', iflag_z0m_snow
     2925    WRITE(lunout,*) ' ok_zsn_ii = ', ok_zsn_ii
     2926    WRITE(lunout,*) ' discret_xf = ', discret_xf
     2927    WRITE(lunout,*) ' correc_alb= ', correc_alb
     2928    WRITE(lunout,*) ' buf_sph_pol = ', buf_sph_pol
     2929    WRITE(lunout,*) ' buf_siz_pol= ', buf_siz_pol
    28392930
    28402931    !$OMP END MASTER
Note: See TracChangeset for help on using the changeset viewer.