Changeset 3343 for trunk/LMDZ.MARS

Timestamp:

May 23, 2024, 3:05:55 PM (6 months ago)

Author:

llange

Message:

Mars PCM
Small correction of how qsurf of co2 is initialized in the south pole.
Some correction for the albedo of co2 ice when using paleoclimate: a uniform value (albedo_co2_cap) is used for both hemisphere for CO2 frost, and a different (higher) value is used for perennial ice.
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/albedocaps.F90 (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/paleoclimate_mod.F90 (modified) (1 diff)
• libf/phymars/phyetat0_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4646 +4646 @@
 == 21/05/2024 == JBC
 Addition of the paleoclimate variable 'd_coef' in the writing of the "restartfi.nc" file.
+
+== 23/05/2024 == LL
+Small correction of how qsurf of co2 is initialized in the south pole.
+Some correction for the albedo of co2 ice when using paleoclimate: a uniform value (albedo_co2_cap) is used for both hemisphere for CO2 frost, and a different (higher) value is used for perennial ice. 

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

@@ -11 +11 @@
 USE mod_phys_lmdz_transfert_para, ONLY: bcast
 USE mod_phys_lmdz_para, ONLY: is_master
-USE paleoclimate_mod, ONLY: paleoclimate,albedo_perennialco2
+USE paleoclimate_mod, ONLY: paleoclimate,albedo_perennialco2,albedo_co2_cap
 
 implicit none
@@ -105 +105 @@
       psolaralb(ig,2)=psolaralb(ig,1)
     else
-      psolaralb(ig,:)=albedice(icap)
-      if (paleoclimate .and. piceco2_peren(ig) > 0. .and. abs(piceco2(ig)) < 1.e-10) psolaralb(ig,:) = albedo_perennialco2
-    endif
+      if(paleoclimate) then
+            psolaralb(ig,1) = albedo_co2_cap
+            psolaralb(ig,2) = albedo_co2_cap
+      else
+            psolaralb(ig,1) = albedice(icap)
+            psolaralb(ig,2) = albedice(icap)
+      endif
+    endif
+  else if(paleoclimate .and. piceco2_peren(ig) > 0.) then
+     psolaralb(ig,1) = albedo_perennialco2
+     psolaralb(ig,2) = albedo_perennialco2
   else if (watercaptag(ig) .and. water) then
     ! there is a water ice cap: set the surface albedo to the water ice one

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

@@ -38 +38 @@
       USE mod_phys_lmdz_transfert_para, ONLY: bcast
       USE paleoclimate_mod,ONLY: paleoclimate,albedo_perennialco2,
+     &                           albedo_co2_cap,
      &                           lag_layer, include_waterbuoyancy 
       use microphys_h, only: mteta
@@ -375 +376 @@
          call getin_p("albedo_perennialco2",albedo_perennialco2)
          write(*,*)"albedo_perennialco2 = ",albedo_perennialco2
+
+         write(*,*)"Albedo for seasonal CO2 ice?"
+         albedo_co2_cap = 0.6 ! default value
+         call getin_p("albedo_co2_cap",albedo_co2_cap)
+         write(*,*)"albedo_co2_cap = ",albedo_co2_cap
 
          write(*,*)"Include water buoyancy effect??"

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

@@ -21 +21 @@
     logical, save                              :: lag_layer             ! Does lag layer is present?
     logical, save                              :: include_waterbuoyancy ! Include the effect of water buoyancy when computing the sublimation of water ice ?
-!$OMP THREADPRIVATE(h2o_ice_depth,lag_co2_ice,d_coef,albedo_perennialco2,lag_layer,include_waterbuoyancy)
+    real,    save                              :: albedo_co2_cap
+!$OMP THREADPRIVATE(h2o_ice_depth,lag_co2_ice,d_coef,albedo_perennialco2,lag_layer,include_waterbuoyancy,albedo_co2_cap)
 
 !=======================================================================

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

@@ -842 +842 @@
         if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
             perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
-            qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) + perennial_co2ice(ngrid,:) ! perennial ice + frost
+            qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) 
         endif
     endif ! not found
@@ -852 +852 @@
     if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
         perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
-        qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) + perennial_co2ice(ngrid,:) ! perennial ice + frost
+        qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:)
     endif
   endif !startphy_file