Changeset 4041 for trunk/LMDZ.MARS

Timestamp:

Feb 5, 2026, 9:30:17 AM (2 weeks ago)

Author:

jbclement

Message:

Mars PCM:
Deletion of depreciated keys to set water ice/frost albedo.
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/conf_phys.F (modified) (3 diffs)
• libf/phymars/tabfi.F (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -5028 +5028 @@
 == 04/02/2026 == EM
 Minor updates/fixes on "xvik" to comply with picky gfortran compiler.
- 
+
+== 05/02/2026 == JBC
+Deletion of depreciated keys to set water ice/frost albedo.

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

@@ -934 +934 @@
          write(*,*) " caps = ",caps
 
-! JN : now separated between albedo_h2o_cap and 
-! albedo_h2o_frost. Retrocompatible with old
-! callphys.def with albedo_h2o_ice
-         write(*,*) "water ice albedo ? Old settings use ",
-     &              "albedo_h2o_ice, new settings use ", 
-     &              "albedo_h2o_cap and albedo_h2o_frost " 
+! Water ice albedo
+         write(*,*) "water cap albedo ?"
          albedo_h2o_cap=0.35
-         albedo_h2o_frost=0.35
-         call getin_p("albedo_h2o_ice",albedo_h2o_cap)
-         albedo_h2o_frost=albedo_h2o_cap
          call getin_p("albedo_h2o_cap",albedo_h2o_cap)
          write(*,*) " albedo_h2o_cap = ",albedo_h2o_cap
+         write(*,*) "water frist albedo ?"
+         albedo_h2o_frost=0.35
          call getin_p("albedo_h2o_frost",albedo_h2o_frost)
          write(*,*) " albedo_h2o_frost = ",albedo_h2o_frost
-
-! Northern polar cap albedo (JN 2021)
-         write(*,*)"Watercaptag albedo is unchanged by water frost",
+         write(*,*) "frost thickness threshold for albedo ?"
+         frost_albedo_threshold=0.005 ! 5.4 mic (i.e 0.005 kg.m-2)
+         call getin_p("frost_albedo_threshold",
+     &    frost_albedo_threshold)
+         write(*,*) " frost_albedo_threshold = ",
+     &            frost_albedo_threshold
+         write(*,*)"Water cap albedo is unchanged by water frost",
      &              " deposition (default is false)"
          cst_cap_albedo=.false. ! default value
@@ -972 +971 @@
      &            frost_metam_threshold
 
-
 ! inert_h2o_ice
          write(*,*) "water ice thermal inertia ?"
@@ -978 +976 @@
          call getin_p("inert_h2o_ice",inert_h2o_ice)
          write(*,*) " inert_h2o_ice = ",inert_h2o_ice
-! frost_albedo_threshold
-         write(*,*) "frost thickness threshold for albedo ?"
-         frost_albedo_threshold=0.005 ! 5.4 mic (i.e 0.005 kg.m-2)
-         call getin_p("frost_albedo_threshold",
-     &    frost_albedo_threshold)
-         write(*,*) " frost_albedo_threshold = ",
-     &            frost_albedo_threshold
 
 ! TMP: old_wsublimation_scheme

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

@@ -108 +108 @@
 c Reference pressure
 c-------------------------------------
-c     pressrf = 670.            ! Pression de reference (Pa) ~650
+c     pressrf = 670. ! Pression de reference (Pa) ~650
 
 c Infos about Mars for the dynamics and physics
 c----------------------------------------------------------
-      rad=3397200.          ! radius of Mars (m)  ~3397200 m
-      daysec=88775.         ! length of a sol (s)  ~88775 s
-      omeg=4.*asin(1.)/(daysec)       ! rotation rate  (rad.s-1)
-      g=3.72                ! gravity (m.s-2) ~3.72
-      mugaz=43.49           ! Molar mass of the atmosphere (g.mol-1) ~43.49
-      rcp=.256793         ! = r/cp  ~0.256793
+      rad=3397200.              ! radius of Mars (m)  ~3397200 m
+      daysec=88775.             ! length of a sol (s)  ~88775 s
+      omeg=4.*asin(1.)/(daysec) ! rotation rate  (rad.s-1)
+      g=3.72                    ! gravity (m.s-2) ~3.72
+      mugaz=43.49               ! Molar mass of the atmosphere (g.mol-1) ~43.49
+      rcp=.256793               ! = r/cp  ~0.256793
 
 c Informations about Mars, only for physics
 c-----------------------------------------------------
-      year_day = 669.       !Modif FH: length of year (sols) ~668.6
-      periheli = 206.66         ! min. Sun-Mars distance (Mkm) ~206.66
-      aphelie = 249.22          ! max. Sun-Mars distance (Mkm) ~249.22
-      peri_day =  485.    ! date of perihelion (sols since N. spring)
-      obliquit = 25.19         ! Obliquity of the planet (deg) ~25.19
+      year_day = 669.   ! Modif FH: length of year (sols) ~668.6
+      periheli = 206.66 ! min. Sun-Mars distance (Mkm) ~206.66
+      aphelie = 249.22  ! max. Sun-Mars distance (Mkm) ~249.22
+      peri_day =  485.  ! date of perihelion (sols since N. spring)
+      obliquit = 25.19  ! Obliquity of the planet (deg) ~25.19
 
 c Boundary layer and turbulence
 c----------------------------
-      z0_default =  1.e-2       ! surface roughness (m) ~0.01
-      emin_turb = 1.e-6         ! minimal energy ~1.e-8
-      lmixmin = 30              ! mixing length ~100
+      z0_default =  1.e-2 ! surface roughness (m) ~0.01
+      emin_turb = 1.e-6   ! minimal energy ~1.e-8
+      lmixmin = 30        ! mixing length ~100
 
 c Optical properties of polar caps and ground emissivity
 c-----------------------------------------------------
-      emissiv=.95               ! Emissivity of martian soil ~.95
-      emisice(1)=0.95           ! Emissivity of northern cap
-      emisice(2)=0.95           ! Emissivity of southern cap
-      albedice(1)=0.65          ! Albedo of northern cap
-      albedice(2)=0.65          ! Albedo of southern cap
-      iceradius(1) = 100.e-6    ! mean scat radius of CO2 snow (north)
-      iceradius(2) = 100.e-6    ! mean scat radius of CO2 snow (south)
-      dtemisice(1) = 0.4   ! time scale for snow metamorphism (north)
-      dtemisice(2) = 0.4   ! time scale for snow metamorphism (south)
+      emissiv=.95            ! Emissivity of martian soil ~.95
+      emisice(1)=0.95        ! Emissivity of northern cap
+      emisice(2)=0.95        ! Emissivity of southern cap
+      albedice(1)=0.65       ! Albedo of northern cap
+      albedice(2)=0.65       ! Albedo of southern cap
+      iceradius(1) = 100.e-6 ! mean scat radius of CO2 snow (north)
+      iceradius(2) = 100.e-6 ! mean scat radius of CO2 snow (south)
+      dtemisice(1) = 0.4     ! time scale for snow metamorphism (north)
+      dtemisice(2) = 0.4     ! time scale for snow metamorphism (south)
 
 c dust aerosol properties
 c---------------------------------
-      tauvis= 0.2          ! mean visible optical depth
+      tauvis= 0.2 ! mean visible optical depth
 
 c  Ancien code radiatif (non utilise avec le code d'apres 03/96)