Ignore:
Timestamp:
Feb 8, 2024, 2:50:21 PM (11 months ago)
Author:
jbclement
Message:

Mars PCM:
Modification of r3200 for the 1D model: the 'CO2cond_ps' coefficient (now defined in "co2condens_mod.F") is applied not only to update the surface pressure but also all the tracers and the winds from now on. It is done in "physiq_mod.F" just after the call to 'co2condens'.
JBC

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.COMMON/libf/evolution/pem.F90

    r3206 r3207  
    6969use compute_soiltemp_mod,       only: compute_tsoil_pem
    7070use writediagpem_mod,           only: writediagpem, writediagsoilpem
     71use co2condens_mod,             only: CO2cond_ps
    7172
    7273#ifndef CPP_STD
     
    136137real, dimension(:,:,:), allocatable :: q             ! champs advectes
    137138real, dimension(ip1jmp1)            :: ps            ! pression au sol
    138 real                                :: CO2cond_ps    ! Coefficient to control the surface pressure change
    139139real, dimension(:),     allocatable :: ps_start_PCM  ! (ngrid) surface pressure
    140140real, dimension(:,:),   allocatable :: ps_timeseries ! (ngrid x timelen) instantaneous surface pressure
     
    296296    allocate(q(ip1jmp1,llm,nqtot))
    297297    allocate(longitude(ngrid),latitude(ngrid),cell_area(ngrid))
    298     CO2cond_ps = 1.
    299298#else
    300299    allocate(q(1,llm,nqtot))
     
    316315    call init_testphys1d('start1D_evol.txt','startfi_evol.nc',therestart1D,therestartfi,ngrid,nlayer,610.,nq,q, &
    317316                         time_0,ps(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,     &
    318                          play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,CO2cond_ps)
     317                         play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
    319318    ps(2) = ps(1)
    320319    nsplit_phys = 1
Note: See TracChangeset for help on using the changeset viewer.