diff --git a/srcF/sfx_api_inmcm.f90 b/srcF/sfx_api_inmcm.f90 index c1354a6ea39f2d3630ebeca755081349f895b99d..37f45df84dd8157e61717a3022324562e483c356 100644 --- a/srcF/sfx_api_inmcm.f90 +++ b/srcF/sfx_api_inmcm.f90 @@ -42,7 +42,7 @@ contains meteo%depth = depth_inm meteo%lai = lai_inm meteo%surface_type = IVEG_sfx - write(*,*) 'surface_type, IVEG_sfx', meteo%surface_type, IVEG_sfx + !write(*,*) 'surface_type, IVEG_sfx', meteo%surface_type, IVEG_sfx end subroutine inmcm_to_sfx_in_cell ! -------------------------------------------------------------------------------- diff --git a/srcF/sfx_config.f90 b/srcF/sfx_config.f90 index 68df2c56de33194ccfce95cc22bedf691b7eff7b..d90c7b246e6059154780e184727052a2d76a371f 100644 --- a/srcF/sfx_config.f90 +++ b/srcF/sfx_config.f90 @@ -167,8 +167,8 @@ contains dataset%filename = get_dataset_filename(id) dataset%nmax = 0 - dataset%surface = surface_lake - dataset%surface_type = surface_lake + dataset%surface = surface_snow + dataset%surface_type = surface_snow dataset%z0_h = -1.0 dataset%lai = 1.0 dataset%depth = 10.0 @@ -176,6 +176,7 @@ contains if (id == dataset_mosaic) then dataset%h = 6.0 dataset%z0_m = 0.0078 + dataset%surface_type = surface_snow else if (id == dataset_irgason) then dataset%h = 9.2 dataset%z0_m = 0.01 diff --git a/srcF/sfx_sheba_noniterative.f90 b/srcF/sfx_sheba_noniterative.f90 index c7912efee0a4025ac5fbc703185c693e93c08abb..9ca876649ea1acab79c5c68a8ed80fba39a163c2 100644 --- a/srcF/sfx_sheba_noniterative.f90 +++ b/srcF/sfx_sheba_noniterative.f90 @@ -154,7 +154,7 @@ contains h = meteo%h(i), & U = meteo%U(i), dT = meteo%dT(i), Tsemi = meteo%Tsemi(i), dQ = meteo%dQ(i), & z0_m = meteo%z0_m(i), depth=meteo%depth(i), lai=meteo%lai(i), surface_type=meteo%surface_type(i)) - write(*,*) 'get_flux', meteo%surface_type(i) + !write(*,*) 'get_flux', meteo%surface_type(i) call get_surface_fluxes(sfx_cell, meteo_cell, numerics) call push_sfx_data(sfx, sfx_cell, i) @@ -302,7 +302,7 @@ contains ! --- get the fluxes ! ---------------------------------------------------------------------------- - write(*,*) 'surface_type', surface_type + !write(*,*) 'surface_type', surface_type if (Rib > 0.0) then ! --- stable stratification block @@ -311,7 +311,7 @@ contains ! Rib = min(Rib, Rib_max) if (surface_type == surface_snow) then - write(*,*) 'snow', Ri_sn, Rib + ! write(*,*) 'snow', Ri_sn, Rib Rib=Rib+Ri_sn endif diff --git a/srcF/sfx_z0m_all_surface.f90 b/srcF/sfx_z0m_all_surface.f90 index 593f663250ee0905416efcf5129116e14cbad955..183403cfaff2ff31ff8133af9c7c90f0a6a6df86 100644 --- a/srcF/sfx_z0m_all_surface.f90 +++ b/srcF/sfx_z0m_all_surface.f90 @@ -84,7 +84,7 @@ module sfx_z0m_all_surface ! --- define dynamic velocity in neutral conditions u_dyn0 = Uc / c - write(*,*) 'ch', z0_m, u_dyn0 + !write(*,*) 'ch', z0_m, u_dyn0 end subroutine ! -------------------------------------------------------------------------------- @@ -125,7 +125,7 @@ module sfx_z0m_all_surface u_dyn0 = Uc / c - write(*,*) 'ow', z0_m, u_dyn0 + !write(*,*) 'ow', z0_m, u_dyn0 end subroutine ! -------------------------------------------------------------------------------- @@ -204,7 +204,7 @@ subroutine get_dynamic_roughness_map(z0_m, u_dyn0, U, h, z0m_map) z0_m=z0m_map h0_m = h / z0_m u_dyn0 = U * kappa / log(h0_m) - write(*,*) 'map', z0_m, u_dyn0 + !write(*,*) 'map', z0_m, u_dyn0 end subroutine ! -------------------------------------------------------------------------------- diff --git a/srcF/sfx_z0t_all_surface.f90 b/srcF/sfx_z0t_all_surface.f90 index c6581f3b82f83bed460de75bfed918af5a7be352..39a570e3567fd20b8af52870bdc229b697257f07 100644 --- a/srcF/sfx_z0t_all_surface.f90 +++ b/srcF/sfx_z0t_all_surface.f90 @@ -49,7 +49,7 @@ module sfx_z0t_all_surface B = min(B, B_max_land) z0_t = z0_m / exp(B) - write(*,*) 'kl_land', z0_t, B + !write(*,*) 'kl_land', z0_t, B end subroutine ! -------------------------------------------------------------------------------- @@ -75,7 +75,7 @@ module sfx_z0t_all_surface B = min(B, B_max_ocean) z0_t = z0_m / exp(B) - write(*,*) 'kl_water', z0_t, B + !write(*,*) 'kl_water', z0_t, B end subroutine