Changeset 1815 for trunk/LMDZ.TITAN/libf/phytitan
- Timestamp:
- Oct 31, 2017, 12:22:52 PM (7 years ago)
- Location:
- trunk/LMDZ.TITAN/libf/phytitan
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/iostart.F90
r1647 r1815 466 466 USE mod_grid_phy_lmdz, only: klon_glo 467 467 USE dimphy, only: klev, klevp1 468 USE tracer_h, only: nqtot 468 USE tracer_h, only: nqtot_p 469 469 USE comsoil_h, only: nsoilmx 470 470 … … 532 532 ENDIF 533 533 534 if (nqtot >0) then534 if (nqtot_p>0) then 535 535 ! only define a tracer dimension if there are tracers 536 ierr=NF90_DEF_DIM(nid_restart,"number_of_advected_fields",nqtot ,idim5)536 ierr=NF90_DEF_DIM(nid_restart,"number_of_advected_fields",nqtot_p,idim5) 537 537 IF (ierr/=NF90_NOERR) THEN 538 538 write(*,*)'open_restartphy: problem defining number_of_advected_fields dimension ' -
trunk/LMDZ.TITAN/libf/phytitan/tracer_h.F90
r1814 r1815 21 21 IMPLICIT NONE 22 22 23 INTEGER, SAVE :: nqtot = 0 !! Total number oftracers23 INTEGER, SAVE :: nqtot_p = 0 !! Total number of physical tracers 24 24 INTEGER, SAVE :: nmicro = 0 !! Number of microphysics tracers. 25 25 INTEGER, SAVE :: nice = 0 !! Number of microphysics ice tracers (subset of nmicro). 26 26 INTEGER, SAVE :: nchimi = 0 !! Number of chemical (gaz species) tracers. 27 !$OMP THREADPRIVATE(nqtot ,nmicro,nice,nchimi)27 !$OMP THREADPRIVATE(nqtot_p,nmicro,nice,nchimi) 28 28 29 29 INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: chimi_indx !! Indexes of all chemical species tracers … … 97 97 verb = .true. ; IF (PRESENT(talk_to_me)) verb = talk_to_me 98 98 99 ! nqtot could be used everywhere in the physic :)100 nqtot =nq99 ! nqtot_p could be used everywhere in the physic :) 100 nqtot_p=nq 101 101 102 102 IF (.NOT.ALLOCATED(noms)) ALLOCATE(noms(nq)) … … 153 153 ! Compute number of chemical species: 154 154 ! simply assume that all other tracers ARE chemical species 155 nchimi = nqtot - nmicro155 nchimi = nqtot_p - nmicro 156 156 157 157 ! Titan chemistry requires exactly 44 tracers:
Note: See TracChangeset
for help on using the changeset viewer.