Ignore:
Timestamp:
Jan 9, 2018, 12:26:53 PM (7 years ago)
Author:
jvatant
Message:

Making chemistry handling more flexible - Step 2
+ Added the calculation of the pressure grid in newstart
using Vervack profile in gr_kim_vervack routine
+ Next step : regridding !
--JVO

Location:
trunk/LMDZ.TITAN/libf/phytitan
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.TITAN/libf/phytitan/comchem_h.F90

    r1887 r1889  
    1919!$OMP_THREADPRIVATE(zlay_kim)
    2020
     21
     22  ! Variable and allocatables for regriding chemistry newstart
     23  ! ----------------------------------------------------------
     24   
     25  INTEGER :: nlaykimold ! Number of upper atm. layers for chemistry in the start_archive file
     26  REAL, ALLOCATABLE, DIMENSION(:) :: preskimold ! Pressure grid of upper chemistry in the start_archive file
    2127
    2228  ! Allocatable arrays for start2archive
  • trunk/LMDZ.TITAN/libf/phytitan/iostart.F90

    r1871 r1889  
    944944                    NF90_PUT_ATT, NF90_NOERR, nf90_strerror, &
    945945                    nf90_inq_dimid, nf90_inquire_dimension, NF90_INQ_VARID
     946  USE comchem_h, only: nlaykim_up
    946947  USE comsoil_h, only: nsoilmx
    947948  USE mod_phys_lmdz_para, only: is_master
     
    997998        ! We know it is an  "mlayer" kind of 1D array
    998999        idim1d=idim3
     1000      ELSEIF (var_size==nlaykim_up) THEN
     1001        ! We know it is an  "preskim" kind of 1D array
     1002        idim1d=idim8
    9991003      ELSE
    10001004        PRINT *, "put_var_rgen error : wrong dimension"
  • trunk/LMDZ.TITAN/libf/phytitan/phyredem.F90

    r1789 r1889  
    99                         alb,ith,pzmea,pzstd,pzsig,pzgam,pzthe)
    1010! create physics restart file and write time-independent variables
     11  use comchem_h, only: preskim
    1112  use comsoil_h, only: volcapa, mlayer
    1213  use geometry_mod, only: cell_area
     
    104105  call put_var("soildepth","Soil mid-layer depth",mlayer)
    105106 
     107  ! Write the mid-layer upper chemistry pressure
     108  call put_var("preskim","Upper chemistry mid-layer pressure",preskim)
     109 
    106110  ! Write longitudes
    107111  call put_field("longitude","Longitudes of physics grid",lonfi)
  • trunk/LMDZ.TITAN/libf/phytitan/tabfi_mod.F90

    r1871 r1889  
    5555                           emissiv
    5656      use comsoil_h, only: volcapa
    57       use comchem_h, only: nlaykim_up
    5857      use iostart, only: get_var
    5958      use mod_phys_lmdz_para, only: is_parallel
     
    150149        dtemisice(:)=0 !time scale for snow metamorphism
    151150        volcapa=1000000 ! volumetric heat capacity of subsurface
    152 ! chemistry
    153         nlaykim_up=70 ! size of vertical grid for upper chemistry
    154151       
    155152      ELSE
     
    207204! soil properties
    208205      volcapa = tab_cntrl(tab0+35) ! volumetric heat capacity
    209 ! chemistry
    210       nlaykim_up = nint(tab_cntrl(tab0+40)) ! size of vertical grid for upper chemistry
    211206!-----------------------------------------------------------------------
    212207!       Save some constants for later use (as routine arguments)
     
    261256
    262257      write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
    263 
    264       write(*,5) '(40)     nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)
    265258
    266259      write(*,*)
     
    571564 
    572565      write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
    573      
    574       write(*,5) '(40)     nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)
    575566
    576567      write(*,*) 
Note: See TracChangeset for help on using the changeset viewer.