Changeset 3586 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jan 20, 2025, 10:31:40 AM (8 days ago)

Author:

jbclement

Message:

Mars PCM:

• Albedo is now set properly in every situation given the presence of CO2/H2O ice and frost + water albedo is done solely in "physiq".
• Small update to make the display of code version status/diff clearer.

JBC

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• albedocaps.F90 (modified) (3 diffs)
• co2condens_mod.F (modified) (1 diff)
• dyn1d/init_testphys1d_mod.F90 (modified) (1 diff)
• dyn1d/testphys1d.F90 (modified) (1 diff)
• physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/albedocaps.F90

@@ -6 +6 @@
 use ioipsl_getin_p_mod, only: getin_p
 use geometry_mod, only: latitude ! grid point latitudes (rad)
-use surfdat_h, only: TESicealbedo, TESice_Ncoef, TESice_Scoef,       &
-                     emisice, albedice, watercaptag, albedo_h2o_cap, &
-                     emissiv, albedodat, albedo_perennialco2
+use surfdat_h, only: TESicealbedo, TESice_Ncoef, TESice_Scoef, &
+                     emisice, emissiv, albedice, albedodat, albedo_perennialco2
 USE mod_phys_lmdz_transfert_para, ONLY: bcast
 USE mod_phys_lmdz_para, ONLY: is_master
@@ -96 +95 @@
         icap = 1 ! Northern hemisphere
     endif
-
-    if (piceco2(ig) > 0) then ! CO2 frost
+    
+    ! Initialization with bare ground
+    ! set emissivity of surface to be bare ground emissivity
+    emisref(ig) = emissiv
+    ! set the surface albedo to bare ground albedo
+    psolaralb(ig,:) = albedodat(ig)
+
+    ! CO2 perennial ice
+    if (paleoclimate .and. piceco2_peren(ig) > 0.) then
+        psolaralb(ig,:) = albedo_perennialco2(icap)
+        emisref(ig) = emisice(icap)
+    endif
+    ! CO2 frost
+    if (piceco2(ig) > 0.) then
         ! set emissivity of surface to be the ice emissivity
         emisref(ig) = emisice(icap)
@@ -105 +116 @@
             psolaralb(ig,2) = psolaralb(ig,1)
         else
-            psolaralb(ig,1) = albedice(icap)
-            psolaralb(ig,2) = albedice(icap)
+            psolaralb(ig,:) = albedice(icap)
         endif
-    else if (paleoclimate .and. piceco2_peren(ig) > 0.) then ! CO2 perennial ice
-        psolaralb(ig,1) = albedo_perennialco2(icap)
-        psolaralb(ig,2) = albedo_perennialco2(icap)
-        emisref(ig) = emisice(icap)
-    else if (watercaptag(ig) .and. water) then ! H2O ice in polar caps
-        ! there is a water ice cap: set the surface albedo to the water ice one
-        ! to do : emissivity
-        emisref(ig) = 1
-        psolaralb(ig,:) = albedo_h2o_cap
-    else ! Bare ground
-        ! set emissivity of surface to be bare ground emissivity
-        emisref(ig) = emissiv
-        ! set the surface albedo to bare ground albedo
-        psolaralb(ig,:) = albedodat(ig)
-    endif ! of if (piceco2(ig).gt.0)
+    endif
 enddo ! of ig=1,ngrid
 

trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F

@@ -612 +612 @@
         alb_tmp = psolaralb(:,:,islope)
         emisref_tmp = 0.
-        perennial_co2ice_tmp =  perennial_co2ice(:,islope) 
+        perennial_co2ice_tmp =  perennial_co2ice(:,islope)
         CALL albedocaps(zls,ngrid,piceco2_tmp,perennial_co2ice_tmp,
      &                  alb_tmp,emisref_tmp)
         perennial_co2ice(:,islope) = perennial_co2ice_tmp
-        psolaralb(:,1,islope) =  alb_tmp(:,1)
-        psolaralb(:,2,islope) =  alb_tmp(:,2)
+        psolaralb(:,:,islope) =  alb_tmp(:,:)
         emisref(:,islope) = emisref_tmp
       ENDDO

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -140 +140 @@
 obliquit = 25.2           ! Obliquity (deg) ~25.2
 eccentric = 0.0934        ! Eccentricity (0.0934)
+Lsperi = 251.             ! Solar longitude of perihelion
 
 ! Planetary Boundary Layer and Turbulence parameters

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

@@ -6 +6 @@
 use phyredem,            only: physdem0, physdem1
 use watersat_mod,        only: watersat
-use tracer_mod,          only: igcm_h2o_vap, igcm_h2o_ice, igcm_co2, noms
+use tracer_mod,          only: igcm_h2o_vap, igcm_h2o_ice, noms
 use comcstfi_h,          only: pi, g, rcp, cpp
-use time_phylmdz_mod, only: daysec
+use time_phylmdz_mod,    only: daysec
 use dimradmars_mod,      only: tauvis, totcloudfrac, albedo
 use dust_param_mod,      only: tauscaling

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -2366 +2366 @@
 c              ALWAYS PLACE these lines after co2condens !!!
 c       -------------------------------------------------------------
-         do ig=1,ngrid
+         do ig = 1,ngrid
           do islope = 1,nslope
-           if ((qsurf(ig,igcm_co2,islope).eq.0).and.
-     &        (qsurf(ig,igcm_h2o_ice,islope)
-     &       .gt.frost_albedo_threshold)) then
-             if ((watercaptag(ig)).and.(cst_cap_albedo)) then
-               albedo(ig,1,islope) = albedo_h2o_cap
-               albedo(ig,2,islope) = albedo_h2o_cap
-             else
-               albedo(ig,1,islope) = albedo_h2o_frost
-               albedo(ig,2,islope) = albedo_h2o_frost
-             endif !((watercaptag(ig)).and.(cst_cap_albedo)) then
-c              write(*,*) "frost thickness", qsurf(ig,igcm_h2o_ice)
-c              write(*,*) "physiq.F frost :"
-c     &        ,latitude(ig)*180./pi, longitude(ig)*180./pi
+           if (abs(qsurf(ig,igcm_co2,islope)) < 1.e-10) then ! No CO2 frost
+
+             if (qsurf(ig,igcm_h2o_ice,islope) > frost_albedo_threshold)
+     & then ! There is H2O frost
+               if (cst_cap_albedo .and. watercaptag(ig) .and.
+     & abs(perennial_co2ice(ig,islope)) < 1.e-10) then ! Water cap remains unchanged by water frost deposition and no CO2 perennial ice
+                 albedo(ig,:,islope) = albedo_h2o_cap
+                 emis(ig,islope) = 1.
+               else
+                 albedo(ig,:,islope) = albedo_h2o_frost
+                 emis(ig,islope) = 1.
+               endif
+             else ! No H2O frost
+               if (abs(perennial_co2ice(ig,islope)) < 1.e-10 .and.
+     & watercaptag(ig)) then ! No CO2 perennial ice but there is water cap
+                 albedo(ig,:,islope) = albedo_h2o_cap
+                 emis(ig,islope) = 1.
+               endif
+             endif
+
            endif
           enddo ! islope