Changeset 4034 for trunk/LMDZ.MARS/libf

Timestamp:

Jan 29, 2026, 6:12:08 PM (4 weeks ago)

Author:

emillour

Message:

Mars PCM:
Update albedocaps.F90 to enable specifying the maximum CO2 ice cap albedo
in callphys.def (rather than have it hard coded to 0.9) as max_icecap_albedo
(default value 0.9 for backward compatibility).
EM

File:

• trunk/LMDZ.MARS/libf/phymars/albedocaps.F90 (modified) (3 diffs)

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

@@ -2 +2 @@
 
 implicit none
+
+real,save :: max_icecap_albedo ! maximum allowed value for albedo
+!$OMP THREADPRIVATE(max_icecap_albedo)
 
 contains
@@ -89 +92 @@
     write(*,*)" albedocaps: TESice_Scoef = ",TESice_Scoef
   endif
+
+  ! load maximum allowed value for cap albedo
+  max_icecap_albedo = 0.9 ! default
+  call getin_p("max_icecap_albedo",max_icecap_albedo)
+  write(*,*)" albedocaps: max_icecap_albedo = ",max_icecap_albedo
 
   call read_TES_icecap_albedo(zls_old,pi,radeg,TESlon,TESlatn,TESlats,TESls,TESalbn,TESalbs)
@@ -525 +533 @@
 endif ! of if (icap.eq.1)
 
-! Make sure that returned albedo is never greater than 0.90
-if (alb.gt.0.90) alb=0.90
+! Make sure that returned albedo is never greater than max_icecap_albedo
+if (alb.gt.max_icecap_albedo) alb=max_icecap_albedo
 
 end subroutine TES_icecap_albedo