Changeset 3328 for trunk/LMDZ.COMMON/libf/evolution
- Timestamp:
- May 15, 2024, 6:51:20 PM (7 months ago)
- Location:
- trunk/LMDZ.COMMON/libf/evolution
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/changelog.txt
r3327 r3328 307 307 - Correction of the computation for the change of sublimating ice necessary to the stopping criteria (the mistake was introduced in r3149) + simplification of the algorithm + renaming of variables with more explicit names; 308 308 - Cleaning of "soil_thermalproperties_mod.F90". 309 310 == 15/05/2024 == LL 311 - Remove a wrong line which modified the initial surface sublimating during 312 the main time loop 313 - Cleaning the previous cleaning commit. 314 -
trunk/LMDZ.COMMON/libf/evolution/ice_table_mod.F90
r3327 r3328 178 178 do ig = 1,ngrid 179 179 do islope = 1,nslope 180 print*, ig, islope, ice_table_depth(ig,islope)181 180 call compute_Tice_pem(nsoil_PEM,tsoil(ig,:,islope),tsurf(ig,islope),ice_table_depth(ig,islope),Tice(ig,islope)) 182 181 rho(ig,islope) = -3.5353e-4*Tice(ig,islope)**2 + 0.0351* Tice(ig,islope) + 933.5030 ! Rottgers, 2012 -
trunk/LMDZ.COMMON/libf/evolution/pem.F90
r3327 r3328 223 223 real :: totmassh2o_adsorbded ! Total mass of H2O that is exchanged because of adsorption / desoprtion over the planets [kg] 224 224 logical :: bool_sublim ! logical to check if there is sublimation or not 225 logical, dimension(:,:), allocatable :: co2ice_has_disappeared ! logical to check if a co2ice reservoir has already disappeared at a previous timestep 225 226 real, dimension(:,:), allocatable :: porefillingice_thickness_prev_iter ! ngrid x nslope: Thickness of the ice table at the previous iteration [m] 226 227 real, dimension(:), allocatable :: delta_h2o_icetablesublim ! ngrid x Total mass of the H2O that has sublimated / condenses from the ice table [kg] … … 496 497 allocate(watersoil_density_PEM_timeseries(ngrid,nsoilmx_PEM,nslope,timelen)) 497 498 allocate(delta_co2_adsorbded(ngrid)) 499 allocate(co2ice_has_disappeared(ngrid,nslope)) 498 500 allocate(porefillingice_thickness_prev_iter(ngrid,nslope)) 499 501 allocate(delta_h2o_icetablesublim(ngrid)) … … 564 566 ini_co2ice_sublim = .false. 565 567 ini_h2oice_sublim = .false. 568 co2ice_has_disappeared(:,:) = .false. 566 569 do i = 1,ngrid 567 570 Total_surface = Total_surface + cell_area(i) … … 789 792 do ig = 1,ngrid 790 793 do islope = 1,nslope 791 if (co2_ice_ini(ig,islope) > 0. .and. co2_ice(ig,islope) < 1.e-10) then ! co2 ice disappeared, look for closest point without co2ice 794 if (co2_ice_ini(ig,islope) > 0. .and. co2_ice(ig,islope) < 1.e-10 .and. (.not.(co2ice_has_disappeared(ig,islope)))) then ! co2 ice disappeared, look for closest point without co2ice 795 co2ice_has_disappeared(ig,islope) = .true. 792 796 if (latitude_deg(ig) > 0) then 793 797 do ig_loop = ig,ngrid … … 813 817 enddo 814 818 endif 815 co2_ice_ini(ig,islope) = 0.816 819 if ((co2_ice(ig,islope) < 1.e-10) .and. (h2o_ice(ig,islope) > frost_albedo_threshold)) then 817 820 albedo(ig,1,islope) = albedo_h2o_frost … … 832 835 enddo 833 836 tsurf_ave(ig,islope) = ave/timelen 834 ! set the surface albedo to be the ice albedo835 ! if (latitude_deg(ig) > 0) then836 ! icap = 1837 ! else838 ! icap = 2839 ! endif840 ! albedo(ig,1,islope) = albedice(icap)841 ! albedo(ig,2,islope) = albedice(icap)842 ! emis(ig,islope) = emisice(icap)843 837 endif 844 838 enddo
Note: See TracChangeset
for help on using the changeset viewer.