Ignore:
Timestamp:
Nov 19, 2021, 4:58:59 PM (3 years ago)
Author:
lguez
Message:

Sync latest trunk changes to Ocean_skin

Location:
LMDZ6/branches/Ocean_skin
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Ocean_skin

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

    r3798 r4013  
    1717       iflag_cld_th, &
    1818       iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
    19        ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, aerosol_couple, chemistry_couple, &
    20        flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
     19       ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat, aerosol_couple, &
     20       chemistry_couple, flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
    2121       flag_bc_internal_mixture, bl95_b0, bl95_b1,&
    2222       read_climoz, &
    … …  
    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
    2931    USE mod_grid_phy_lmdz, ONLY: klon_glo
    3032    USE print_control_mod, ONLY: lunout
    … …  
    6567    ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
    6668    ! ok_volcan: activate volcanic diags (SW heat & LW cool rate, SW & LW flux)
     69    ! flag_volc_surfstrat: VolMIP flag, activate forcing surface cooling rate (=1), strato heating rate (=2) or nothing (=0, default)
    6770    !
    6871
    … …  
    7780    INTEGER              :: flag_aerosol
    7881    INTEGER              :: flag_aerosol_strat
     82    INTEGER              :: flag_volc_surfstrat
    7983    LOGICAL              :: flag_aer_feedback
    8084    LOGICAL              :: flag_bc_internal_mixture
    … …  
    8892    CHARACTER (len = 8), SAVE  :: aer_type_omp
    8993    INTEGER, SAVE       :: landice_opt_omp
    90     INTEGER, SAVE       :: n_dtis_omp
    91     INTEGER, SAVE       :: iflag_tsurf_inlandsis_omp
    92     INTEGER, SAVE       :: iflag_albzenith_omp   
    93     LOGICAL, SAVE       :: SnoMod_omp,BloMod_omp,ok_outfor_omp
     94    INTEGER, SAVE       :: iflag_tsurf_inlandsis_omp,iflag_temp_inlandsis_omp
     95    INTEGER, SAVE       :: iflag_albcalc_omp,iflag_z0m_snow_omp   
     96    LOGICAL, SAVE       :: SnoMod_omp,BloMod_omp,ok_outfor_omp,ok_zsn_ii_omp
     97    LOGICAL, SAVE       :: discret_xf_omp,opt_runoff_ac_omp 
     98    LOGICAL, SAVE       :: is_ok_slush_omp,is_ok_z0h_rn_omp,is_ok_density_kotlyakov_omp
     99    REAL, SAVE          :: prescribed_z0m_snow_omp,correc_alb_omp
     100    REAL, SAVE          :: buf_sph_pol_omp,buf_siz_pol_omp
    94101    LOGICAL, SAVE       :: ok_newmicro_omp
    95102    LOGICAL, SAVE       :: ok_all_xml_omp
    … …  
    102109    INTEGER, SAVE       :: flag_aerosol_omp
    103110    INTEGER, SAVE       :: flag_aerosol_strat_omp
     111    INTEGER, SAVE       :: flag_volc_surfstrat_omp
    104112    LOGICAL, SAVE       :: flag_aer_feedback_omp
    105113    LOGICAL, SAVE       :: flag_bc_internal_mixture_omp
    … …  
    174182    INTEGER,SAVE :: iflag_cloudth_vert_omp
    175183    INTEGER,SAVE :: iflag_rain_incloud_vol_omp
     184    INTEGER,SAVE :: iflag_vice_omp
    176185    REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
    177186    REAL,SAVE :: t_glace_min_omp, t_glace_max_omp
    178187    REAL,SAVE :: exposant_glace_omp
     188    INTEGER,SAVE :: iflag_gammasat_omp, iflag_mpc_bl_omp
    179189    REAL,SAVE :: rei_min_omp, rei_max_omp
    180     INTEGER,SAVE :: iflag_sic_omp
     190    INTEGER,SAVE :: iflag_sic_omp, iflag_inertie_omp
    181191    REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
    182192    REAL,SAVE :: inertie_lic_omp
    … …  
    237247    LOGICAL, SAVE :: carbon_cycle_rad_omp
    238248    INTEGER, SAVE :: level_coupling_esm_omp
     249    LOGICAL, SAVE :: read_fco2_ocean_cor_omp
     250    REAL, SAVE    :: var_fco2_ocean_cor_omp
     251    LOGICAL, SAVE :: read_fco2_land_cor_omp
     252    REAL, SAVE    :: var_fco2_land_cor_omp
    239253    LOGICAL, SAVE :: adjust_tropopause_omp
    240254    LOGICAL, SAVE :: ok_daily_climoz_omp
     255    LOGICAL, SAVE :: ok_new_lscp_omp
     256    LOGICAL, SAVE :: ok_icefra_lscp_omp
    241257
    242258    INTEGER, INTENT(OUT):: read_climoz ! read ozone climatology, OpenMP shared
    … …  
    327343    ! Martin et Etienne
    328344    !Config Key  = landice_opt
    329     !Config Desc = which landice snow model (BULK, SISVAT or INLANDSIS)
     345    !Config Desc = which landice snow model (BULK, or INLANDSIS)
    330346    !Config Def  = 0
    331347    landice_opt_omp = 0
    … …  
    334350
    335351    !Etienne
     352    !Config Key  = iflag_temp_inlandsis
     353    !Config Desc = which method to calculate temp within the soil in INLANDSIS
     354    !Config Def  = 0
     355    iflag_temp_inlandsis_omp = 0
     356    CALL getin('iflag_temp_inlandsis', iflag_temp_inlandsis_omp)
     357
     358    !Etienne
    336359    !Config Key  = iflag_tsurf_inlandsis
    337360    !Config Desc = which method to calculate tsurf in INLANDSIS
    338361    !Config Def  = 0
    339     iflag_tsurf_inlandsis_omp = 0
     362    iflag_tsurf_inlandsis_omp = 1
    340363    CALL getin('iflag_tsurf_inlandsis', iflag_tsurf_inlandsis_omp)
    341364
     365
    342366    !Etienne
    343     !Config Key  = iflag_albzenith
    344     !Config Desc = method to account for albedo sensitivity to solar zenith angle
    345     !Config Def  = 0
    346     iflag_albzenith_omp = 0
    347     CALL getin('iflag_albzenith', iflag_albzenith_omp)
    348 
    349     !Etienne
    350     !Config Key  = n_dtis
    351     !Config Desc = number of subtimesteps for INLANDSIS
    352     !Config Def  = 1
    353     n_dtis_omp = 1
    354     CALL getin('n_dtis', n_dtis_omp)
     367    !Config Key  = iflag_albcalc
     368    !Config Desc = method to calculate snow albedo in INLANDSIS
     369    !Config Def  = 0
     370    iflag_albcalc_omp = 0
     371    CALL getin('iflag_albcalc', iflag_albcalc_omp)
     372
    355373
    356374    !Etienne
    357375    !Config Key  = SnoMod
    358376    !Config Desc = activation of snow modules in inlandsis
    359     !Config Def  = 1
     377    !Config Def  = .TRUE.
    360378    SnoMod_omp = .TRUE.
    361379    CALL getin('SnoMod', SnoMod_omp)
    … …  
    364382    !Config Key  = BloMod
    365383    !Config Desc = activation of blowing snow in inlandsis
    366     !Config Def  = 1
     384    !Config Def  = .FALSE.
    367385    BloMod_omp = .FALSE.
    368386    CALL getin('BloMod', BloMod_omp)
    … …  
    371389    !Config Key  = ok_outfor
    372390    !Config Desc = activation of output ascii file in inlandsis
    373     !Config Def  = 1
    374    ok_outfor_omp = .FALSE.
     391    !Config Def  = .FALSE.
     392    ok_outfor_omp = .FALSE.
    375393    CALL getin('ok_outfor', ok_outfor_omp)
    376394
    377395
     396    !Etienne
     397    !Config Key  = ok_sn_ii
     398    !Config Desc = activation of ice/snow detection
     399    !Config Def  = .TRUE.
     400    ok_zsn_ii_omp = .TRUE.
     401    CALL getin('ok_zsn_ii', ok_zsn_ii_omp)
     402
     403
     404    !Etienne
     405    !Config Key  = discret_xf
     406    !Config Desc = snow discretization following XF
     407    !Config Def  = .TRUE.
     408    discret_xf_omp = .TRUE.
     409    CALL getin('discret_xf', discret_xf_omp)
     410
     411
     412    !Etienne
     413    !Config Key  = is_ok_slush
     414    !Config Desc = activation of the slush option
     415    !Config Def  = .TRUE.
     416    is_ok_slush_omp = .TRUE.
     417    CALL getin('is_ok_slush', is_ok_slush_omp)
     418
     419    !Etienne
     420    !Config Key  = opt_runoff_ac
     421    !Config Desc = option runoff AC
     422    !Config Def  = .TRUE.
     423    opt_runoff_ac_omp = .TRUE.
     424    CALL getin('opt_runoff_ac', opt_runoff_ac_omp)
     425
     426    !Etienne
     427    !Config Key  = is_ok_z0h_rn
     428    !Config Desc = z0h calculation following RN method
     429    !Config Def  = .TRUE.
     430    is_ok_z0h_rn_omp = .TRUE.
     431    CALL getin('is_ok_z0h_rn', is_ok_z0h_rn_omp)
     432
     433
     434    !Etienne
     435    !Config Key  = is_ok_density_kotlyakov
     436    !Config Desc = snow density calculation following kotlyakov
     437    !Config Def  = .FALSE.
     438    is_ok_density_kotlyakov_omp = .FALSE.
     439    CALL getin('is_ok_density_kotlyakov', is_ok_density_kotlyakov_omp)
     440
     441
     442    !Etienne
     443    !Config Key  = prescribed_z0m_snow
     444    !Config Desc = prescribed snow z0m
     445    !Config Def  = 0.005
     446    prescribed_z0m_snow_omp = 0.005
     447    CALL getin('prescribed_z0m_snow', prescribed_z0m_snow_omp)
     448
     449
     450    !Etienne
     451    !Config Key  = iflag_z0m_snow
     452    !Config Desc = method to calculate snow z0m
     453    !Config Def  = 0
     454    iflag_z0m_snow_omp = 0
     455    CALL getin('iflag_z0m_snow', iflag_z0m_snow_omp)
     456
     457
     458    !Etienne
     459    !Config Key  = correc_alb
     460    !Config Desc = correction term for albedo
     461    !Config Def  = 1.01
     462    correc_alb_omp=1.01
     463    CALL getin('correc_alb', correc_alb_omp)
     464
     465
     466    !Etienne
     467    !Config Key  = buf_sph_pol
     468    !Config Desc = sphericity of buffer layer in polar regions
     469    !Config Def  = 99.
     470    buf_sph_pol_omp=99.
     471    CALL getin('buf_sph_pol', buf_sph_pol_omp)
     472
     473    !Etienne
     474    !Config Key  = buf_siz_pol
     475    !Config Desc = grain size of buffer layer in polar regions in e-4m
     476    !Config Def  = 4.
     477    buf_siz_pol_omp=4.
     478    CALL getin('buf_siz_pol', buf_siz_pol_omp)
    378479
    379480    !==================================================================
    … …  
    457558    ok_volcan_omp = .FALSE.
    458559    CALL getin('ok_volcan', ok_volcan_omp)
     560
     561    !
     562    !Config Key  = flag_volc_surfstrat
     563    !Config Desc = impose cooling rate at the surface (=1),
     564    !              heating rate in the strato (=2), or nothing (=0)
     565    !Config Def  = 0
     566    !Config Help = Used in radlwsw_m.F
     567    !
     568    flag_volc_surfstrat_omp = 0 ! NL: SURFSTRAT
     569    CALL getin('flag_volc_surfstrat', flag_volc_surfstrat_omp)
    459570
    460571    !
    … …  
    12311342    CALL getin('iflag_sic',iflag_sic_omp)
    12321343    !
     1344    !Config Key  = iflag_inertie
     1345    !Config Desc =
     1346    !Config Def  = 0
     1347    !Config Help =
     1348    !
     1349    iflag_inertie_omp = 0
     1350    CALL getin('iflag_inertie',iflag_inertie_omp)
     1351    !
    12331352    !Config Key  = inertie_sic
    12341353    !Config Desc = 
    … …  
    13181437
    13191438    !
     1439    !Config Key  = iflag_gammasat
     1440    !Config Desc = 
     1441    !Config Def  = 0
     1442    !Config Help =
     1443    !
     1444    iflag_gammasat_omp=0
     1445    CALL getin('iflag_gammasat',iflag_gammasat_omp)
     1446
     1447
     1448    !
     1449    !Config Key  = iflag_mpc_bl
     1450    !Config Desc =
     1451    !Config Def  = 0
     1452    !Config Help =
     1453    !
     1454    iflag_mpc_bl_omp=0
     1455    CALL getin('iflag_mpc_bl',iflag_mpc_bl_omp)
     1456
     1457
     1458
     1459    !
    13201460    !Config Key  = iflag_t_glace
    13211461    !Config Desc = 
    … …  
    13431483    iflag_rain_incloud_vol_omp = 0
    13441484    CALL getin('iflag_rain_incloud_vol',iflag_rain_incloud_vol_omp)
     1485
     1486    !
     1487    !Config Key  = iflag_vice
     1488    !Config Desc = 
     1489    !Config Def  = 0
     1490    !Config Help =
     1491    !
     1492    iflag_vice_omp = 0
     1493    CALL getin('iflag_vice',iflag_vice_omp)
     1494
     1495
    13451496
    13461497    !
    … …  
    21982349    !Config Help = .FALSE. ensure much fewer (no calendar dependency)
    21992350    !  and lighter monthly climoz files, inetrpolated in time at gcm run time.
    2200     !
     2351   
     2352    ok_new_lscp_omp = .FALSE.
     2353    CALL getin('ok_new_lscp', ok_new_lscp_omp)
     2354    !
     2355    !Config Key  = ok_new_lscp_omp
     2356    !Config Desc = new cloud scheme ith ice and mixed phase (Etienne and JB)
     2357    !Config Def  = .FALSE.
     2358    !Config Help = ...
     2359
     2360
     2361
     2362    ok_icefra_lscp_omp = .FALSE.
     2363    CALL getin('ok_icefra_lscp', ok_icefra_lscp_omp)
     2364    !
     2365    !Config Key  = ok_icefra_lscp_omp
     2366    !Config Desc = ice fraction in radiation from lscp
     2367    !Config Def  = .FALSE.
     2368    !Config Help = ...
     2369
     2370
     2371
     2372
    22012373    ecrit_LES_omp = 1./8.
    22022374    CALL getin('ecrit_LES', ecrit_LES_omp)
    … …  
    22142386    CALL getin('carbon_cycle_rad',carbon_cycle_rad_omp)
    22152387
    2216     ! >> PC
     2388    read_fco2_ocean_cor_omp=.FALSE.
     2389    CALL getin('read_fco2_ocean_cor',read_fco2_ocean_cor_omp)
     2390
     2391    var_fco2_ocean_cor_omp=0. ! default value
     2392    CALL getin('var_fco2_ocean_cor',var_fco2_ocean_cor_omp)
     2393
     2394    read_fco2_land_cor_omp=.FALSE.
     2395    CALL getin('read_fco2_land_cor',read_fco2_land_cor_omp)
     2396
     2397    var_fco2_land_cor_omp=0. ! default value
     2398    CALL getin('var_fco2_land_cor',var_fco2_land_cor_omp)
     2399
    22172400    ! level_coupling_esm : level of coupling of the biogeochemical fields between LMDZ, ORCHIDEE and NEMO
    22182401    ! Definitions of level_coupling_esm in physiq.def
    … …  
    22272410    level_coupling_esm_omp=0 ! default value
    22282411    CALL getin('level_coupling_esm',level_coupling_esm_omp)
    2229     ! << PC
    22302412
    22312413    !$OMP END MASTER
    … …  
    22912473    albsno0 = albsno0_omp
    22922474    iflag_sic = iflag_sic_omp
     2475    iflag_inertie = iflag_inertie_omp
    22932476    inertie_sol = inertie_sol_omp
    22942477    inertie_sic = inertie_sic_omp
    … …  
    23012484    t_glace_max = t_glace_max_omp
    23022485    exposant_glace = exposant_glace_omp
     2486    iflag_gammasat=iflag_gammasat_omp
     2487    iflag_mpc_bl=iflag_mpc_bl_omp
    23032488    iflag_t_glace = iflag_t_glace_omp
    23042489    iflag_cloudth_vert=iflag_cloudth_vert_omp
    23052490    iflag_rain_incloud_vol=iflag_rain_incloud_vol_omp
     2491    iflag_vice=iflag_vice_omp
    23062492    iflag_ice_thermo = iflag_ice_thermo_omp
    23072493    rei_min = rei_min_omp
    … …  
    23442530       ok_veget=.FALSE.
    23452531    ENDIF
    2346     ! SISVAT and INLANDSIS
     2532    ! INLANDSIS
    23472533    !=================================================
    23482534    landice_opt = landice_opt_omp
    23492535    iflag_tsurf_inlandsis = iflag_tsurf_inlandsis_omp
    2350     iflag_albzenith = iflag_albzenith_omp
    2351     n_dtis=n_dtis_omp
     2536    iflag_temp_inlandsis = iflag_temp_inlandsis_omp
     2537    iflag_albcalc = iflag_albcalc_omp
    23522538    SnoMod=SnoMod_omp
    23532539    BloMod=BloMod_omp
    23542540    ok_outfor=ok_outfor_omp
     2541    is_ok_slush=is_ok_slush_omp
     2542    opt_runoff_ac=opt_runoff_ac_omp
     2543    is_ok_z0h_rn=is_ok_z0h_rn_omp
     2544    is_ok_density_kotlyakov=is_ok_density_kotlyakov_omp
     2545    prescribed_z0m_snow=prescribed_z0m_snow_omp
     2546    correc_alb=correc_alb_omp
     2547    iflag_z0m_snow=iflag_z0m_snow_omp
     2548    ok_zsn_ii=ok_zsn_ii_omp
     2549    discret_xf=discret_xf_omp
     2550    buf_sph_pol=buf_sph_pol_omp
     2551    buf_siz_pol=buf_siz_pol_omp
    23552552    !=================================================
    23562553    ok_all_xml = ok_all_xml_omp
    … …  
    23702567    ok_cdnc = ok_cdnc_omp
    23712568    ok_volcan = ok_volcan_omp
     2569    flag_volc_surfstrat = flag_volc_surfstrat_omp
    23722570    aerosol_couple = aerosol_couple_omp
    23732571    chemistry_couple = chemistry_couple_omp
    2374     flag_aerosol=flag_aerosol_omp
    2375     flag_aerosol_strat=flag_aerosol_strat_omp
    2376     flag_aer_feedback=flag_aer_feedback_omp
     2572    flag_aerosol = flag_aerosol_omp
     2573    flag_aerosol_strat = flag_aerosol_strat_omp
     2574    flag_aer_feedback = flag_aer_feedback_omp
    23772575    flag_bc_internal_mixture=flag_bc_internal_mixture_omp
    23782576    aer_type = aer_type_omp
    … …  
    24932691    carbon_cycle_rad = carbon_cycle_rad_omp
    24942692    level_coupling_esm = level_coupling_esm_omp
     2693    ok_new_lscp = ok_new_lscp_omp
     2694    ok_icefra_lscp=ok_icefra_lscp_omp
     2695    read_fco2_ocean_cor = read_fco2_ocean_cor_omp
     2696    var_fco2_ocean_cor = var_fco2_ocean_cor_omp
     2697    read_fco2_land_cor = read_fco2_land_cor_omp
     2698    var_fco2_land_cor = var_fco2_land_cor_omp
    24952699
    24962700    ! Test of coherence between type_ocean and version_ocean
    … …  
    25152719      ENDIF
    25162720    ELSE IF (iflag_rrtm .EQ. 1) THEN
     2721      IF (NSW.NE.2.AND.NSW.NE.4.AND.NSW.NE.6) THEN
     2722        WRITE(lunout,*) ' ERROR iflag_rrtm=1 and NSW<>2,4,6 not possible'
     2723        CALL abort_physic('conf_phys','choice NSW not valid',1)
     2724      ENDIF
     2725   ELSE IF (iflag_rrtm .EQ. 2) THEN
    25172726      IF (NSW.NE.2.AND.NSW.NE.4.AND.NSW.NE.6) THEN
    25182727        WRITE(lunout,*) ' ERROR iflag_rrtm=1 and NSW<>2,4,6 not possible'
    … …  
    25902799       CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with flag_aerosol=6',1)
    25912800    ENDIF
     2801
     2802    ! test sur flag_volc_surfstrat
     2803    IF (flag_volc_surfstrat.LT.0.OR.flag_volc_surfstrat.GT.2) THEN
     2804       CALL abort_physic('conf_phys', 'flag_volc_surfstrat can only be 0 1 or 2',1)
     2805    ENDIF
     2806    IF ((.NOT.ok_volcan.OR..NOT.ok_ade.OR..NOT.ok_aie).AND.flag_volc_surfstrat.GT.0) THEN
     2807       CALL abort_physic('conf_phys', 'ok_ade, ok_aie, ok_volcan need to be activated if flag_volc_surfstrat is 1 or 2',1)
     2808    ENDIF
    25922809
    25932810    ! Test on carbon cycle
    … …  
    26942911    WRITE(lunout,*) ' t_glace_max = ',t_glace_max
    26952912    WRITE(lunout,*) ' exposant_glace = ',exposant_glace
     2913    WRITE(lunout,*) ' iflag_gammasat = ',iflag_gammasat
     2914    WRITE(lunout,*) ' iflag_mpc_bl = ',iflag_mpc_bl
    26962915    WRITE(lunout,*) ' iflag_t_glace = ',iflag_t_glace
    26972916    WRITE(lunout,*) ' iflag_cloudth_vert = ',iflag_cloudth_vert
    26982917    WRITE(lunout,*) ' iflag_rain_incloud_vol = ',iflag_rain_incloud_vol
     2918    WRITE(lunout,*) ' iflag_vice = ',iflag_vice
    26992919    WRITE(lunout,*) ' iflag_ice_thermo = ',iflag_ice_thermo
    27002920    WRITE(lunout,*) ' rei_min = ',rei_min
    … …  
    27122932    WRITE(lunout,*) ' ok_ade = ',ok_ade
    27132933    WRITE(lunout,*) ' ok_volcan = ',ok_volcan
     2934    WRITE(lunout,*) ' flag_volc_surfstrat = ',flag_volc_surfstrat
    27142935    WRITE(lunout,*) ' ok_aie = ',ok_aie
    27152936    WRITE(lunout,*) ' ok_alw = ',ok_alw
    … …  
    27572978    WRITE(lunout,*) ' albsno0 = ', albsno0
    27582979    WRITE(lunout,*) ' iflag_sic = ', iflag_sic
     2980    WRITE(lunout,*) ' iflag_inertie = ', iflag_inertie
    27592981    WRITE(lunout,*) ' inertie_sol = ', inertie_sol
    27602982    WRITE(lunout,*) ' inertie_sic = ', inertie_sic
    … …  
    28073029    WRITE(lunout,*) ' adjust_tropopause = ', adjust_tropopause
    28083030    WRITE(lunout,*) ' ok_daily_climoz = ',ok_daily_climoz
     3031    WRITE(lunout,*) ' ok_new_lscp = ', ok_new_lscp
     3032    WRITE(lunout,*) ' ok_icefra_lscp = ', ok_icefra_lscp
    28093033    WRITE(lunout,*) ' read_climoz = ', read_climoz
    28103034    WRITE(lunout,*) ' carbon_cycle_tr = ', carbon_cycle_tr
    … …  
    28123036    WRITE(lunout,*) ' carbon_cycle_rad = ', carbon_cycle_rad
    28133037    WRITE(lunout,*) ' level_coupling_esm = ', level_coupling_esm
     3038    WRITE(lunout,*) ' read_fco2_ocean_cor = ', read_fco2_ocean_cor
     3039    WRITE(lunout,*) ' var_fco2_ocean_cor = ', var_fco2_ocean_cor
     3040    WRITE(lunout,*) ' read_fco2_land_cor = ', read_fco2_land_cor
     3041    WRITE(lunout,*) ' var_fco2_land_cor = ', var_fco2_land_cor
    28143042    WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
    2815     WRITE(lunout,*) ' iflag_albzenith = ', iflag_albzenith
    2816     WRITE(lunout,*) ' n_dtis = ', n_dtis
     3043    WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis
     3044    WRITE(lunout,*) ' iflag_albcalc = ', iflag_albcalc
    28173045    WRITE(lunout,*) ' SnoMod = ', SnoMod
    28183046    WRITE(lunout,*) ' BloMod = ', BloMod
    28193047    WRITE(lunout,*) ' ok_outfor = ', ok_outfor
    2820 
     3048    WRITE(lunout,*) ' is_ok_slush = ', is_ok_slush
     3049    WRITE(lunout,*) ' opt_runoff_ac = ', opt_runoff_ac
     3050    WRITE(lunout,*) ' is_ok_z0h_rn = ', is_ok_z0h_rn
     3051    WRITE(lunout,*) ' is_ok_density_kotlyakov = ', is_ok_density_kotlyakov
     3052    WRITE(lunout,*) ' prescribed_z0m_snow = ', prescribed_z0m_snow
     3053    WRITE(lunout,*) ' iflag_z0m_snow = ', iflag_z0m_snow
     3054    WRITE(lunout,*) ' ok_zsn_ii = ', ok_zsn_ii
     3055    WRITE(lunout,*) ' discret_xf = ', discret_xf
     3056    WRITE(lunout,*) ' correc_alb= ', correc_alb
     3057    WRITE(lunout,*) ' buf_sph_pol = ', buf_sph_pol
     3058    WRITE(lunout,*) ' buf_siz_pol= ', buf_siz_pol
    28213059
    28223060    !$OMP END MASTER
Note: See TracChangeset for help on using the changeset viewer.