Changeset 1889 for trunk/LMDZ.TITAN/libf/phytitan
- Timestamp:
- Jan 9, 2018, 12:26:53 PM (7 years ago)
- Location:
- trunk/LMDZ.TITAN/libf/phytitan
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/comchem_h.F90
r1887 r1889 19 19 !$OMP_THREADPRIVATE(zlay_kim) 20 20 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 21 27 22 28 ! Allocatable arrays for start2archive -
trunk/LMDZ.TITAN/libf/phytitan/iostart.F90
r1871 r1889 944 944 NF90_PUT_ATT, NF90_NOERR, nf90_strerror, & 945 945 nf90_inq_dimid, nf90_inquire_dimension, NF90_INQ_VARID 946 USE comchem_h, only: nlaykim_up 946 947 USE comsoil_h, only: nsoilmx 947 948 USE mod_phys_lmdz_para, only: is_master … … 997 998 ! We know it is an "mlayer" kind of 1D array 998 999 idim1d=idim3 1000 ELSEIF (var_size==nlaykim_up) THEN 1001 ! We know it is an "preskim" kind of 1D array 1002 idim1d=idim8 999 1003 ELSE 1000 1004 PRINT *, "put_var_rgen error : wrong dimension" -
trunk/LMDZ.TITAN/libf/phytitan/phyredem.F90
r1789 r1889 9 9 alb,ith,pzmea,pzstd,pzsig,pzgam,pzthe) 10 10 ! create physics restart file and write time-independent variables 11 use comchem_h, only: preskim 11 12 use comsoil_h, only: volcapa, mlayer 12 13 use geometry_mod, only: cell_area … … 104 105 call put_var("soildepth","Soil mid-layer depth",mlayer) 105 106 107 ! Write the mid-layer upper chemistry pressure 108 call put_var("preskim","Upper chemistry mid-layer pressure",preskim) 109 106 110 ! Write longitudes 107 111 call put_field("longitude","Longitudes of physics grid",lonfi) -
trunk/LMDZ.TITAN/libf/phytitan/tabfi_mod.F90
r1871 r1889 55 55 emissiv 56 56 use comsoil_h, only: volcapa 57 use comchem_h, only: nlaykim_up58 57 use iostart, only: get_var 59 58 use mod_phys_lmdz_para, only: is_parallel … … 150 149 dtemisice(:)=0 !time scale for snow metamorphism 151 150 volcapa=1000000 ! volumetric heat capacity of subsurface 152 ! chemistry153 nlaykim_up=70 ! size of vertical grid for upper chemistry154 151 155 152 ELSE … … 207 204 ! soil properties 208 205 volcapa = tab_cntrl(tab0+35) ! volumetric heat capacity 209 ! chemistry210 nlaykim_up = nint(tab_cntrl(tab0+40)) ! size of vertical grid for upper chemistry211 206 !----------------------------------------------------------------------- 212 207 ! Save some constants for later use (as routine arguments) … … 261 256 262 257 write(*,5) '(35) volcapa',tab_cntrl(tab0+35),volcapa 263 264 write(*,5) '(40) nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)265 258 266 259 write(*,*) … … 571 564 572 565 write(*,5) '(35) volcapa',tab_cntrl(tab0+35),volcapa 573 574 write(*,5) '(40) nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)575 566 576 567 write(*,*)
Note: See TracChangeset
for help on using the changeset viewer.