Changeset 2918


Ignore:
Timestamp:
Mar 21, 2023, 7:29:38 PM (22 months ago)
Author:
llange
Message:

PEM
Implementing Eran's idea: introduce a flag in the run_PEM.def to compute CO2 glacier flows or not
LL

Location:
trunk/LMDZ.COMMON/libf/evolution
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/evolution/co2glaciers_mod.F90

    r2856 r2918  
    11      module co2glaciers_mod
    22        implicit none
     3        LOGICAL  co2glaciersflow ! True by default, to compute co2 ice flow. Read in  pem.def
    34
    45!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  • trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90

    r2895 r2918  
    1818  USE comsoil_h_pem, only: soil_pem,fluxgeo,water_reservoir_nom,depth_breccia,depth_bedrock,reg_thprop_dependp
    1919  USE adsorption_mod,only: adsorption_pem
     20  USE co2glaciers_mod, only: co2glaciersflow
     21
    2022  CHARACTER(len=20),parameter :: modname ='conf_pem'
    2123
    22 !PEM parameter
     24!PEM parameters
     25
     26
    2327
    2428  year_bp_ini=0.
    … …  
    4246  Max_iter_pem=99999999
    4347  CALL getin('Max_iter_pem', Max_iter_pem)
     48
     49  co2glaciersflow = .true.
     50  CALL getin('co2glaciersflow', co2glaciersflow)
    4451
    4552  soil_pem=.true.

  • trunk/LMDZ.COMMON/libf/evolution/pem.F90

    r2897 r2918  
    8585      use conf_pem_mod, only: conf_pem
    8686      use pemredem, only:  pemdem0,pemdem1
    87       use co2glaciers_mod,only: co2glaciers_evol
     87      use co2glaciers_mod,only: co2glaciers_evol,co2glaciersflow
    8888      use criterion_pem_stop_mod,only: criterion_waterice_stop,criterion_co2_stop
    8989!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! SOIL
    … …  
    855855!------------------------
    856856
    857       print *, "Co2 glacier flows"
    858 
    859     call co2glaciers_evol(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,vmr_co2_pem_phys,ps_timeseries,&
     857      print *, "CO2 glacier flows"
     858
     859      if (co2glaciersflow) then
     860       call co2glaciers_evol(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,vmr_co2_pem_phys,ps_phys_timeseries,&
    860861                         global_ave_press_GCM,global_ave_press_new,co2ice_slope,flag_co2flow,flag_co2flow_mesh)
     862      endif
     863
     864
    861865
    862866     DO islope=1, nslope
Note: See TracChangeset for help on using the changeset viewer.