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

Timestamp:

Jan 16, 2025, 9:50:10 AM (6 months ago)

Author:

jbclement

Message:

Mars PCM:

• Correction of r3581: the key-words to get 'albedo_perennialco2(1:2)" from the "callphys.def" are now 'albedo_perennialco2_north' and 'albedo_perennialco2_south'.
• Moving 'albedo_perennialco2' from the 'paleoclimate' module to the 'surfdat_h' module.

JBC

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• albedocaps.F90 (modified) (1 diff)
• conf_phys.F (modified) (4 diffs)
• paleoclimate_mod.F90 (modified) (1 diff)
• surfdat_h.F90 (modified) (7 diffs)

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, &
+use surfdat_h, only: TESicealbedo, TESice_Ncoef, TESice_Scoef,       &
                      emisice, albedice, watercaptag, albedo_h2o_cap, &
-                     emissiv, albedodat
+                     emissiv, albedodat, albedo_perennialco2
 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
 
 implicit none

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

@@ -20 +20 @@
      &                       nuice_ref,nuiceco2_ref
       use surfdat_h, only: albedo_h2o_cap,albedo_h2o_frost,
+     &                     albedo_perennialco2,
      &                     frost_albedo_threshold, inert_h2o_ice,
      &                     frost_metam_threshold,old_wsublimation_scheme
@@ -37 +38 @@
       use aeropacity_mod, only: iddist, topdustref
       USE mod_phys_lmdz_transfert_para, ONLY: bcast
-      USE paleoclimate_mod,ONLY: paleoclimate,albedo_perennialco2,
+      USE paleoclimate_mod,ONLY: paleoclimate,
      &                           lag_layer, include_waterbuoyancy 
       use microphys_h, only: mteta
@@ -56 +57 @@
       
       real :: ecritphy ! to check that this obsolete flag is no longer used...
-      real :: value_albed
+      real :: albedo_perennialco2_north, albedo_perennialco2_south
  
       CHARACTER ch1*12
@@ -385 +386 @@
 
          write(*,*)"Albedo for perennial CO2 ice (northern hemisphere)?"
-         albedo_perennialco2(1) = 0.6 ! default value
-         call getin_p("albedo_perennialco2_north",value_albed)
-         albedo_perennialco2(1) = value_albed
+         albedo_perennialco2_north = 0.6 ! default value
+         call getin_p("albedo_perennialco2_north",albedo_perennialco2_
+     &   north)
+         albedo_perennialco2(1) = albedo_perennialco2_north
          write(*,*)"albedo_perennialco2(1) = ",albedo_perennialco2(1)
          
          write(*,*)"Albedo for perennial CO2 ice (southern hemisphere)?"
-         albedo_perennialco2(2) = 0.85 ! default value
-         call getin_p("albedo_perennialco2_south",value_albed)
-         albedo_perennialco2(2) = value_albed
+         albedo_perennialco2_south = 0.85 ! default value
+         call getin_p("albedo_perennialco2_south",albedo_perennialco2_
+     &   south)
+         albedo_perennialco2(2) = albedo_perennialco2_south
          write(*,*)"albedo_perennialco2(2) = ",albedo_perennialco2(2)
 

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

@@ -17 +17 @@
     real,    allocatable, dimension(:,:) :: lag_co2_ice           ! Thickness of the lag before CO2 ice [m]
     real,    allocatable, dimension(:,:) :: d_coef                ! Diffusion coefficent
-    real, dimension(2)                   :: albedo_perennialco2   ! Albedo for perennial co2 ice [1]
     logical                              :: lag_layer             ! Does lag layer is present?
     logical                              :: 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)
+!$OMP THREADPRIVATE(h2o_ice_depth,lag_co2_ice,d_coef,lag_layer,include_waterbuoyancy)
 
 !=======================================================================

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

@@ -2 +2 @@
 
   ! arrays are allocated in conf_phys
-  real,save,allocatable :: albedodat(:) ! albedo of bare ground
-  real,save,allocatable :: phisfi(:) ! geopotential at ground level
-  real, dimension(2) :: albedice ! default albedo for ice (1:Northern hemisphere 2:Southern hemisphere)
-  real, dimension(2) :: emisice ! ice emissivity (1:Northern hemisphere 2:Southern hemisphere)
-  real,save :: emissiv ! emissivity of bare ground
-  logical,save :: TESicealbedo ! use TES ice cap albedoes (if set to .true.)
-  logical,save,allocatable :: watercaptag(:) ! flag for water ice surface
-  real, save, allocatable :: dryness(:)
+  real, allocatable    :: albedodat(:) ! albedo of bare ground
+  real, allocatable    :: phisfi(:) ! geopotential at ground level
+  real, dimension(2)   :: albedice ! default albedo for ice (1:Northern hemisphere 2:Southern hemisphere)
+  real, dimension(2)   :: emisice ! ice emissivity (1:Northern hemisphere 2:Southern hemisphere)
+  real                 :: emissiv ! emissivity of bare ground
+  logical              :: TESicealbedo ! use TES ice cap albedoes (if set to .true.)
+  logical, allocatable :: watercaptag(:) ! flag for water ice surface
+  real,    allocatable :: dryness(:)
 
 !$OMP THREADPRIVATE(albedodat, phisfi,albedice,emisice,emissiv,TESicealbedo,     &
 !$OMP                watercaptag,dryness)
-      
-  logical,save :: temptag !temp tag for water caps
+
+  logical :: temptag !temp tag for water caps
 
 !$OMP THREADPRIVATE(temptag)
-      
-  real,save :: albedo_h2o_cap ! water cap albedo
-  real,save :: albedo_h2o_frost ! water frost albedo
-  real,save :: inert_h2o_ice ! water ice thermal inertia
-  real,save :: frost_albedo_threshold ! water frost thickness on the ground (kg.m^-2, ie mm)
-  real,save :: frost_metam_threshold ! water frost threshold before conversion to ice (kg.m^-2, ie mm)
-  real,save :: TESice_Ncoef ! coefficient for TES ice albedo in Northern hemisphere
-  real,save :: TESice_Scoef ! coefficient for TES ice albedo in Southern hemisphere
-  real,save :: iceradius(2) , dtemisice(2)
-  real,save,allocatable :: zmea(:),zstd(:),zsig(:),zgam(:),zthe(:)
-  real,save,allocatable :: hmons(:),summit(:),base(:)
-  real,save,allocatable :: z0(:) ! surface roughness length (m)
-  real,save :: z0_default ! default (constant over planet) surface roughness (m)
 
-  LOGICAL, SAVE :: old_wsublimation_scheme    ! TEMPORARY : TO USE THE OLD WATER SUBLIMATION SCHEME (i.e., using Cd instead of Ch), true by default
+  real              :: albedo_h2o_cap ! water cap albedo
+  real              :: albedo_h2o_frost ! water frost albedo
+  real              :: inert_h2o_ice ! water ice thermal inertia
+  real              :: frost_albedo_threshold ! water frost thickness on the ground (kg.m^-2, ie mm)
+  real              :: frost_metam_threshold ! water frost threshold before conversion to ice (kg.m^-2, ie mm)
+  real              :: TESice_Ncoef ! coefficient for TES ice albedo in Northern hemisphere
+  real              :: TESice_Scoef ! coefficient for TES ice albedo in Southern hemisphere
+  real              :: iceradius(2) , dtemisice(2)
+  real, allocatable :: zmea(:),zstd(:),zsig(:),zgam(:),zthe(:)
+  real, allocatable :: hmons(:),summit(:),base(:)
+  real, allocatable :: z0(:) ! surface roughness length (m)
+  real              :: z0_default ! default (constant over planet) surface roughness (m)
+
+  logical :: old_wsublimation_scheme    ! TEMPORARY : TO USE THE OLD WATER SUBLIMATION SCHEME (i.e., using Cd instead of Ch), true by default
 
 !$OMP THREADPRIVATE(albedo_h2o_cap,albedo_h2o_frost,inert_h2o_ice,               &
@@ -40 +40 @@
 
   !! mountain top dust flows
-  REAL,SAVE,ALLOCATABLE :: alpha_hmons(:) ! sub-grid scale mountain mesh fraction
-  REAL,SAVE,ALLOCATABLE :: hsummit(:) ! mountain height above the GCM surface
-  LOGICAL,SAVE,ALLOCATABLE :: contains_mons(:) ! is there a mountain in the grid mesh ?
-          
+  real,    allocatable :: alpha_hmons(:) ! sub-grid scale mountain mesh fraction
+  real,    allocatable :: hsummit(:) ! mountain height above the GCM surface
+  logical, allocatable :: contains_mons(:) ! is there a mountain in the grid mesh ?
+
 !$OMP THREADPRIVATE(alpha_hmons,hsummit,contains_mons)
 
   !! variables
-  REAL,SAVE,ALLOCATABLE :: tsurf(:,:)   ! Surface temperature (K)
-  REAL,SAVE,ALLOCATABLE :: emis(:,:)    ! Thermal IR surface emissivity
-  REAL,SAVE,ALLOCATABLE :: capcal(:,:) ! surface heat capacity (J m-2 K-1)
-  REAL,SAVE,ALLOCATABLE :: fluxgrd(:,:) ! surface conduction flux (W.m-2)
-  REAL,ALLOCATABLE,SAVE :: qsurf(:,:,:) ! tracer on surface (e.g. kg.m-2)
-  REAL,SAVE,ALLOCATABLE :: watercap(:,:) ! Surface water ice (kg.m-2)
-  REAL,SAVE,ALLOCATABLE :: perennial_co2ice(:,:) ! Perennial CO2 ice (kg.m-2)
+  real, allocatable  :: tsurf(:,:)   ! Surface temperature (K)
+  real, allocatable  :: emis(:,:)    ! Thermal IR surface emissivity
+  real, allocatable  :: capcal(:,:) ! surface heat capacity (J m-2 K-1)
+  real, allocatable  :: fluxgrd(:,:) ! surface conduction flux (W.m-2)
+  real, allocatable  :: qsurf(:,:,:) ! tracer on surface (e.g. kg.m-2)
+  real, allocatable  :: watercap(:,:) ! Surface water ice (kg.m-2)
+  real, allocatable  :: perennial_co2ice(:,:) ! Perennial CO2 ice (kg.m-2)
+  real, dimension(2) :: albedo_perennialco2   ! Albedo for perennial co2 ice (-)
 !$OMP THREADPRIVATE(tsurf,emis,capcal,fluxgrd,qsurf,watercap,perennial_co2ice)
 
@@ -59 +60 @@
 
   subroutine ini_surfdat_h(ngrid,nq,nslope)
-  
+
   implicit none
   integer,intent(in) :: ngrid ! number of atmospheric columns
-  integer,intent(in) :: nq ! number of tracers  
-  integer,intent(in) :: nslope ! number of sub-grid scale slope  
+  integer,intent(in) :: nq ! number of tracers
+  integer,intent(in) :: nslope ! number of sub-grid scale slope
+
     allocate(albedodat(ngrid))
     allocate(phisfi(ngrid))
@@ -87 +89 @@
     allocate(hsummit(ngrid))
     allocate(contains_mons(ngrid))
-       
+
   end subroutine ini_surfdat_h
 
@@ -118 +120 @@
     if (allocated(hsummit))          deallocate(hsummit)
     if (allocated(contains_mons))    deallocate(contains_mons)
-    
+
   end subroutine end_surfdat_h
 
   subroutine ini_surfdat_h_slope_var(ngrid,nq,nslope)
-  
+
   implicit none
   integer,intent(in) :: ngrid ! number of atmospheric columns
-  integer,intent(in) :: nq ! number of tracers  
-  integer,intent(in) :: nslope ! number of sub-grid scale slope  
+  integer,intent(in) :: nq ! number of tracers
+  integer,intent(in) :: nslope ! number of sub-grid scale slope
     allocate(qsurf(ngrid,nq,nslope))
     allocate(tsurf(ngrid,nslope))
@@ -134 +136 @@
     allocate(capcal(ngrid,nslope))
     allocate(fluxgrd(ngrid,nslope))
-       
+
   end subroutine ini_surfdat_h_slope_var
 
@@ -148 +150 @@
     if (allocated(capcal))           deallocate(capcal)
     if (allocated(fluxgrd))          deallocate(fluxgrd)
-    
+
   end subroutine end_surfdat_h_slope_var