Changeset 3032 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Aug 24, 2023, 4:46:03 PM (16 months ago)

Author:

llange

Message:

MARS PEM
Fix bug with albedo (wrong albedo for surface where co2 ice is) and add the porosity in the computation commited previously (missing)
LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• ice_table_mod.F90 (modified) (2 diffs)
• pem.F90 (modified) (5 diffs)

trunk/LMDZ.COMMON/libf/evolution/ice_table_mod.F90

@@ -149 +149 @@
     use comslope_mod, only: subslope_dist,def_slope_mean           
     use comconst_mod,only: pi       
+    use constants_marspem_mod,only:porosity
     implicit none 
 !   inputs
@@ -201 +202 @@
    do ig = 1,ngrid
       do islope = 1,nslope
-        delta_m_h2o(ig)  = delta_m_h2o(ig) +  rho(ig,islope)*(new_ice_table_thickness(ig,islope) - former_ice_table_thickness(ig,islope)) & ! convention > 0. <=> it condenses  
+        delta_m_h2o(ig)  = delta_m_h2o(ig) +  porosity*rho(ig,islope)*(new_ice_table_thickness(ig,islope) - former_ice_table_thickness(ig,islope)) & ! convention > 0. <=> it condenses  
                            *subslope_dist(ig,islope)/cos(def_slope_mean(islope)*pi/180.)
       enddo

trunk/LMDZ.COMMON/libf/evolution/pem.F90

@@ -65 +65 @@
                                   albedodat, zmea, zstd, zsig, zgam, zthe,     &
                                   hmons, summit, base,albedo_h2o_frost,        &
-                                  frost_albedo_threshold, emissiv, watercaptag, perenial_co2ice
+                                  frost_albedo_threshold, emissiv, watercaptag, perenial_co2ice, &
+                                  emisice, albedice
     use dimradmars_mod,     only: totcloudfrac, albedo
     use dust_param_mod,     only: tauscaling
@@ -72 +73 @@
     use planete_h,          only: aphelie, periheli, year_day, peri_day, obliquit
     use comcstfi_h,         only: r, mugaz
+    use paleoclimate_mod,   only: albedo_perenialco2
 #else
     use tracer_h,           only: noms, igcm_h2o_ice, igcm_co2 ! Tracer names
@@ -245 +247 @@
 
 ! Loop variables
-integer :: i, j, ig0, l, ig, nnq, t, islope, ig_loop, islope_loop, iloop, isoil
+integer :: i, j, ig0, l, ig, nnq, t, islope, ig_loop, islope_loop, iloop, isoil,icap
 
 ! Parallel variables
@@ -843 +845 @@
                 enddo
                 tsurf_ave(ig,islope) = ave/timelen
+                ! set the surface albedo to be the ice albedo
+                if (latitude_deg(ig) > 0) then
+                   icap = 1
+                else
+                   icap = 2
+                endif
+                albedo(ig,1,islope) = albedice(icap)
+                albedo(ig,2,islope) = albedice(icap)
+                emis(ig,islope) = emisice(icap)
             endif
         enddo
@@ -1042 +1053 @@
             perenial_co2ice(ig,islope) = 0.5*qsurf(ig,igcm_co2,islope)
             qsurf(ig,igcm_co2,islope)  = 0.5*qsurf(ig,igcm_co2,islope)
+            albedo(ig,1,islope) = albedo_perenialco2
+            albedo(ig,2,islope) = albedo_perenialco2
         endif
     enddo