Ignore:
Timestamp:
Jan 15, 2025, 6:35:44 PM (2 weeks ago)
Author:
jbclement
Message:

Mars PCM:

  • Deletion of the variable 'albedo_co2_cap' which was redundant with 'albedo_perennialco2'.
  • 'albedo_perennialco2' is now a vector of dimension 2 to account for different albedos of CO2 ice between the northern and southern hemispheres (by default 0.6 and 0.85 respectively).

JBC

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/conf_phys.F

    r3468 r3581  
    3838      USE mod_phys_lmdz_transfert_para, ONLY: bcast
    3939      USE paleoclimate_mod,ONLY: paleoclimate,albedo_perennialco2,
    40      &                           albedo_co2_cap,
    4140     &                           lag_layer, include_waterbuoyancy
    4241      use microphys_h, only: mteta
     
    5756     
    5857      real :: ecritphy ! to check that this obsolete flag is no longer used...
     58      real :: value_albed
    5959 
    6060      CHARACTER ch1*12
     
    384384         write(*,*)"paleoclimate = ",paleoclimate
    385385
    386          write(*,*)"Albedo for perennial CO2 ice?"
    387          albedo_perennialco2 = 0.85 ! default value
    388          call getin_p("albedo_perennialco2",albedo_perennialco2)
    389          write(*,*)"albedo_perennialco2 = ",albedo_perennialco2
    390 
    391          write(*,*)"Albedo for seasonal CO2 ice?"
    392          albedo_co2_cap = 0.6 ! default value
    393          call getin_p("albedo_co2_cap",albedo_co2_cap)
    394          write(*,*)"albedo_co2_cap = ",albedo_co2_cap
     386         write(*,*)"Albedo for perennial CO2 ice (northern hemisphere)?"
     387         albedo_perennialco2(1) = 0.6 ! default value
     388         call getin_p("albedo_perennialco2_north",value_albed)
     389         albedo_perennialco2(1) = value_albed
     390         write(*,*)"albedo_perennialco2(1) = ",albedo_perennialco2(1)
     391         
     392         write(*,*)"Albedo for perennial CO2 ice (southern hemisphere)?"
     393         albedo_perennialco2(2) = 0.85 ! default value
     394         call getin_p("albedo_perennialco2_south",value_albed)
     395         albedo_perennialco2(2) = value_albed
     396         write(*,*)"albedo_perennialco2(2) = ",albedo_perennialco2(2)
    395397
    396398         write(*,*)"Include water buoyancy effect??"
Note: See TracChangeset for help on using the changeset viewer.