Changeset 4013 for LMDZ6/branches/Ocean_skin/libf/phylmd/conf_phys_m.F90
- Timestamp:
- Nov 19, 2021, 4:58:59 PM (3 years ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
- Property svn:mergeinfo changed
-
LMDZ6/branches/Ocean_skin/libf/phylmd/conf_phys_m.F90
r3798 r4013 17 17 iflag_cld_th, & 18 18 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, & 21 21 flag_bc_internal_mixture, bl95_b0, bl95_b1,& 22 22 read_climoz, & … … 27 27 USE phys_cal_mod 28 28 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 29 31 USE mod_grid_phy_lmdz, ONLY: klon_glo 30 32 USE print_control_mod, ONLY: lunout … … 65 67 ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc 66 68 ! 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) 67 70 ! 68 71 … … 77 80 INTEGER :: flag_aerosol 78 81 INTEGER :: flag_aerosol_strat 82 INTEGER :: flag_volc_surfstrat 79 83 LOGICAL :: flag_aer_feedback 80 84 LOGICAL :: flag_bc_internal_mixture … … 88 92 CHARACTER (len = 8), SAVE :: aer_type_omp 89 93 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 94 101 LOGICAL, SAVE :: ok_newmicro_omp 95 102 LOGICAL, SAVE :: ok_all_xml_omp … … 102 109 INTEGER, SAVE :: flag_aerosol_omp 103 110 INTEGER, SAVE :: flag_aerosol_strat_omp 111 INTEGER, SAVE :: flag_volc_surfstrat_omp 104 112 LOGICAL, SAVE :: flag_aer_feedback_omp 105 113 LOGICAL, SAVE :: flag_bc_internal_mixture_omp … … 174 182 INTEGER,SAVE :: iflag_cloudth_vert_omp 175 183 INTEGER,SAVE :: iflag_rain_incloud_vol_omp 184 INTEGER,SAVE :: iflag_vice_omp 176 185 REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp 177 186 REAL,SAVE :: t_glace_min_omp, t_glace_max_omp 178 187 REAL,SAVE :: exposant_glace_omp 188 INTEGER,SAVE :: iflag_gammasat_omp, iflag_mpc_bl_omp 179 189 REAL,SAVE :: rei_min_omp, rei_max_omp 180 INTEGER,SAVE :: iflag_sic_omp 190 INTEGER,SAVE :: iflag_sic_omp, iflag_inertie_omp 181 191 REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp 182 192 REAL,SAVE :: inertie_lic_omp … … 237 247 LOGICAL, SAVE :: carbon_cycle_rad_omp 238 248 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 239 253 LOGICAL, SAVE :: adjust_tropopause_omp 240 254 LOGICAL, SAVE :: ok_daily_climoz_omp 255 LOGICAL, SAVE :: ok_new_lscp_omp 256 LOGICAL, SAVE :: ok_icefra_lscp_omp 241 257 242 258 INTEGER, INTENT(OUT):: read_climoz ! read ozone climatology, OpenMP shared … … 327 343 ! Martin et Etienne 328 344 !Config Key = landice_opt 329 !Config Desc = which landice snow model (BULK, SISVATor INLANDSIS)345 !Config Desc = which landice snow model (BULK, or INLANDSIS) 330 346 !Config Def = 0 331 347 landice_opt_omp = 0 … … 334 350 335 351 !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 336 359 !Config Key = iflag_tsurf_inlandsis 337 360 !Config Desc = which method to calculate tsurf in INLANDSIS 338 361 !Config Def = 0 339 iflag_tsurf_inlandsis_omp = 0362 iflag_tsurf_inlandsis_omp = 1 340 363 CALL getin('iflag_tsurf_inlandsis', iflag_tsurf_inlandsis_omp) 341 364 365 342 366 !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 355 373 356 374 !Etienne 357 375 !Config Key = SnoMod 358 376 !Config Desc = activation of snow modules in inlandsis 359 !Config Def = 1377 !Config Def = .TRUE. 360 378 SnoMod_omp = .TRUE. 361 379 CALL getin('SnoMod', SnoMod_omp) … … 364 382 !Config Key = BloMod 365 383 !Config Desc = activation of blowing snow in inlandsis 366 !Config Def = 1384 !Config Def = .FALSE. 367 385 BloMod_omp = .FALSE. 368 386 CALL getin('BloMod', BloMod_omp) … … 371 389 !Config Key = ok_outfor 372 390 !Config Desc = activation of output ascii file in inlandsis 373 !Config Def = 1374 ok_outfor_omp = .FALSE.391 !Config Def = .FALSE. 392 ok_outfor_omp = .FALSE. 375 393 CALL getin('ok_outfor', ok_outfor_omp) 376 394 377 395 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) 378 479 379 480 !================================================================== … … 457 558 ok_volcan_omp = .FALSE. 458 559 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) 459 570 460 571 ! … … 1231 1342 CALL getin('iflag_sic',iflag_sic_omp) 1232 1343 ! 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 ! 1233 1352 !Config Key = inertie_sic 1234 1353 !Config Desc = … … 1318 1437 1319 1438 ! 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 ! 1320 1460 !Config Key = iflag_t_glace 1321 1461 !Config Desc = … … 1343 1483 iflag_rain_incloud_vol_omp = 0 1344 1484 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 1345 1496 1346 1497 ! … … 2198 2349 !Config Help = .FALSE. ensure much fewer (no calendar dependency) 2199 2350 ! 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 2201 2373 ecrit_LES_omp = 1./8. 2202 2374 CALL getin('ecrit_LES', ecrit_LES_omp) … … 2214 2386 CALL getin('carbon_cycle_rad',carbon_cycle_rad_omp) 2215 2387 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 2217 2400 ! level_coupling_esm : level of coupling of the biogeochemical fields between LMDZ, ORCHIDEE and NEMO 2218 2401 ! Definitions of level_coupling_esm in physiq.def … … 2227 2410 level_coupling_esm_omp=0 ! default value 2228 2411 CALL getin('level_coupling_esm',level_coupling_esm_omp) 2229 ! << PC2230 2412 2231 2413 !$OMP END MASTER … … 2291 2473 albsno0 = albsno0_omp 2292 2474 iflag_sic = iflag_sic_omp 2475 iflag_inertie = iflag_inertie_omp 2293 2476 inertie_sol = inertie_sol_omp 2294 2477 inertie_sic = inertie_sic_omp … … 2301 2484 t_glace_max = t_glace_max_omp 2302 2485 exposant_glace = exposant_glace_omp 2486 iflag_gammasat=iflag_gammasat_omp 2487 iflag_mpc_bl=iflag_mpc_bl_omp 2303 2488 iflag_t_glace = iflag_t_glace_omp 2304 2489 iflag_cloudth_vert=iflag_cloudth_vert_omp 2305 2490 iflag_rain_incloud_vol=iflag_rain_incloud_vol_omp 2491 iflag_vice=iflag_vice_omp 2306 2492 iflag_ice_thermo = iflag_ice_thermo_omp 2307 2493 rei_min = rei_min_omp … … 2344 2530 ok_veget=.FALSE. 2345 2531 ENDIF 2346 ! SISVAT andINLANDSIS2532 ! INLANDSIS 2347 2533 !================================================= 2348 2534 landice_opt = landice_opt_omp 2349 2535 iflag_tsurf_inlandsis = iflag_tsurf_inlandsis_omp 2350 iflag_ albzenith = iflag_albzenith_omp2351 n_dtis=n_dtis_omp2536 iflag_temp_inlandsis = iflag_temp_inlandsis_omp 2537 iflag_albcalc = iflag_albcalc_omp 2352 2538 SnoMod=SnoMod_omp 2353 2539 BloMod=BloMod_omp 2354 2540 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 2355 2552 !================================================= 2356 2553 ok_all_xml = ok_all_xml_omp … … 2370 2567 ok_cdnc = ok_cdnc_omp 2371 2568 ok_volcan = ok_volcan_omp 2569 flag_volc_surfstrat = flag_volc_surfstrat_omp 2372 2570 aerosol_couple = aerosol_couple_omp 2373 2571 chemistry_couple = chemistry_couple_omp 2374 flag_aerosol =flag_aerosol_omp2375 flag_aerosol_strat =flag_aerosol_strat_omp2376 flag_aer_feedback =flag_aer_feedback_omp2572 flag_aerosol = flag_aerosol_omp 2573 flag_aerosol_strat = flag_aerosol_strat_omp 2574 flag_aer_feedback = flag_aer_feedback_omp 2377 2575 flag_bc_internal_mixture=flag_bc_internal_mixture_omp 2378 2576 aer_type = aer_type_omp … … 2493 2691 carbon_cycle_rad = carbon_cycle_rad_omp 2494 2692 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 2495 2699 2496 2700 ! Test of coherence between type_ocean and version_ocean … … 2515 2719 ENDIF 2516 2720 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 2517 2726 IF (NSW.NE.2.AND.NSW.NE.4.AND.NSW.NE.6) THEN 2518 2727 WRITE(lunout,*) ' ERROR iflag_rrtm=1 and NSW<>2,4,6 not possible' … … 2590 2799 CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with flag_aerosol=6',1) 2591 2800 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 2592 2809 2593 2810 ! Test on carbon cycle … … 2694 2911 WRITE(lunout,*) ' t_glace_max = ',t_glace_max 2695 2912 WRITE(lunout,*) ' exposant_glace = ',exposant_glace 2913 WRITE(lunout,*) ' iflag_gammasat = ',iflag_gammasat 2914 WRITE(lunout,*) ' iflag_mpc_bl = ',iflag_mpc_bl 2696 2915 WRITE(lunout,*) ' iflag_t_glace = ',iflag_t_glace 2697 2916 WRITE(lunout,*) ' iflag_cloudth_vert = ',iflag_cloudth_vert 2698 2917 WRITE(lunout,*) ' iflag_rain_incloud_vol = ',iflag_rain_incloud_vol 2918 WRITE(lunout,*) ' iflag_vice = ',iflag_vice 2699 2919 WRITE(lunout,*) ' iflag_ice_thermo = ',iflag_ice_thermo 2700 2920 WRITE(lunout,*) ' rei_min = ',rei_min … … 2712 2932 WRITE(lunout,*) ' ok_ade = ',ok_ade 2713 2933 WRITE(lunout,*) ' ok_volcan = ',ok_volcan 2934 WRITE(lunout,*) ' flag_volc_surfstrat = ',flag_volc_surfstrat 2714 2935 WRITE(lunout,*) ' ok_aie = ',ok_aie 2715 2936 WRITE(lunout,*) ' ok_alw = ',ok_alw … … 2757 2978 WRITE(lunout,*) ' albsno0 = ', albsno0 2758 2979 WRITE(lunout,*) ' iflag_sic = ', iflag_sic 2980 WRITE(lunout,*) ' iflag_inertie = ', iflag_inertie 2759 2981 WRITE(lunout,*) ' inertie_sol = ', inertie_sol 2760 2982 WRITE(lunout,*) ' inertie_sic = ', inertie_sic … … 2807 3029 WRITE(lunout,*) ' adjust_tropopause = ', adjust_tropopause 2808 3030 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 2809 3033 WRITE(lunout,*) ' read_climoz = ', read_climoz 2810 3034 WRITE(lunout,*) ' carbon_cycle_tr = ', carbon_cycle_tr … … 2812 3036 WRITE(lunout,*) ' carbon_cycle_rad = ', carbon_cycle_rad 2813 3037 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 2814 3042 WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis 2815 WRITE(lunout,*) ' iflag_ albzenith = ', iflag_albzenith2816 WRITE(lunout,*) ' n_dtis = ', n_dtis3043 WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis 3044 WRITE(lunout,*) ' iflag_albcalc = ', iflag_albcalc 2817 3045 WRITE(lunout,*) ' SnoMod = ', SnoMod 2818 3046 WRITE(lunout,*) ' BloMod = ', BloMod 2819 3047 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 2821 3059 2822 3060 !$OMP END MASTER
Note: See TracChangeset
for help on using the changeset viewer.