Changeset 3130 for trunk/LMDZ.MARS/libf/phymars/albedocaps.F90

Timestamp:

Nov 20, 2023, 1:25:31 PM (13 months ago)

Author:

jbclement

Message:

PEM:
The perennial co2 ice is now taken into account with co2 frost (qsurf) to compute the tendency and to make the update + Rework of how co2 frost is converted to perennial co2 ice at the end of the PEM run + Correction of the value of 'threshold_co2_frost2perennial' to correspond to 10 m + Perennial co2 ice is now handled outside 'paleoclimate' in "phyetat0_mod.F90" of the Mars PCM + Some cleanings.

/!\ Commit for the PEM management of co2 ice before a rework of ice management in the PEM!
JBC

File:

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

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_perenialco2
+USE paleoclimate_mod, ONLY: paleoclimate,albedo_perennialco2
 
 implicit none
@@ -21 +21 @@
 integer,intent(in) :: ngrid
 real,intent(in) :: piceco2(ngrid) ! amount of CO2 ice on the surface (kg/m2)
-real,intent(inout) :: piceco2_peren(ngrid) ! amount of perenial co2 ice (kg/m^2)
+real,intent(inout) :: piceco2_peren(ngrid) ! amount of perennial co2 ice (kg/m^2)
 real,intent(out) :: psolaralb(ngrid,2) ! albedo of the surface
 real,intent(out) :: emisref(ngrid) ! emissivity of the surface
@@ -109 +109 @@
       if(paleoclimate) then
          if((piceco2_peren(ig).gt.0.).and.(piceco2(ig).lt.piceco2_peren(ig))) then
-             psolaralb(ig,1) = albedo_perenialco2
-             psolaralb(ig,2) = albedo_perenialco2
+             psolaralb(ig,1) = albedo_perennialco2
+             psolaralb(ig,2) = albedo_perennialco2
              piceco2_peren(ig) = piceco2(ig)
          endif