Ignore:
Timestamp:
Oct 31, 2017, 12:22:52 PM (7 years ago)
Author:
jvatant
Message:

Adapt modifs related to tracers to newstart and start2archive
+ nqtot in physics -> nqtot_p. Else conflict with nqtot from dyn
TBD : use of nqtot_p everywhere in physics instead of nq ...
--JVO

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

Legend:

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

    r1647 r1815  
    466466  USE mod_grid_phy_lmdz, only: klon_glo
    467467  USE dimphy, only: klev, klevp1
    468   USE tracer_h, only: nqtot
     468  USE tracer_h, only: nqtot_p
    469469  USE comsoil_h, only: nsoilmx
    470470
     
    532532      ENDIF
    533533     
    534       if (nqtot>0) then
     534      if (nqtot_p>0) then
    535535        ! 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)
    537537        IF (ierr/=NF90_NOERR) THEN
    538538          write(*,*)'open_restartphy: problem defining number_of_advected_fields dimension '
  • trunk/LMDZ.TITAN/libf/phytitan/tracer_h.F90

    r1814 r1815  
    2121  IMPLICIT NONE
    2222
    23   INTEGER, SAVE :: nqtot  = 0 !! Total number of tracers
     23  INTEGER, SAVE :: nqtot_p  = 0 !! Total number of physical tracers
    2424  INTEGER, SAVE :: nmicro = 0 !! Number of microphysics tracers.
    2525  INTEGER, SAVE :: nice   = 0 !! Number of microphysics ice tracers (subset of nmicro).
    2626  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)
    2828
    2929  INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: chimi_indx  !! Indexes of all chemical species tracers
     
    9797    verb = .true. ; IF (PRESENT(talk_to_me)) verb = talk_to_me
    9898
    99     ! nqtot could be used everywhere in the physic :)
    100     nqtot=nq
     99    ! nqtot_p could be used everywhere in the physic :)
     100    nqtot_p=nq
    101101
    102102    IF (.NOT.ALLOCATED(noms)) ALLOCATE(noms(nq))
     
    153153    ! Compute number of chemical species:
    154154    !   simply assume that all other tracers ARE chemical species
    155     nchimi = nqtot - nmicro
     155    nchimi = nqtot_p - nmicro
    156156
    157157    ! Titan chemistry requires exactly 44 tracers:
Note: See TracChangeset for help on using the changeset viewer.