Ignore:
Timestamp:
Jan 16, 2025, 8:14:15 PM (13 hours ago)
Author:
dcugnet
Message:

Remove tracers attributes "isAdvected" and "isInPhysics" from infotrac (iadv is enough).
Remove tracers attribute "isAdvected" from infotrac_phy (isInPhysics is now equivalent
to former isInPhysics .AND. iadv > 0

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3d_common/infotrac.f90

    r5475 r5481  
    7878!  | nqChildren  | Number of childs            (1st generation only)    | nqfils      | 1:nqtot                |
    7979!  | iadv        | Advection scheme number                              | iadv        | 1,2,10-20(exc.15,19),30|
    80 !  | isAdvected  | Advected tracers flag (.TRUE. if iadv >= 0)          | /           | nqtrue  .TRUE. values  |
    81 !  | isInPhysics | Tracers not extracted from the main table in physics | /           | nqtottr .TRUE. values  |
    8280!  | iso_iGroup  | Isotopes group index in isotopes(:)                  | /           | 1:nbIso                |
    8381!  | iso_iName   | Isotope  name  index in isotopes(iso_iGroup)%trac(:) | iso_indnum  | 1:niso                 |
     
    260258
    261259!==============================================================================================================================
    262 ! 3) Determine the advection scheme choice for water and tracers "iadv" and the fields long name, isAdvected.
     260! 3) Determine the advection scheme choice for water and tracers "iadv" and the field "long name".
    263261!     iadv = 1    "LMDZ-specific humidity transport" (for H2O vapour)          LMV
    264262!     iadv = 2    backward                           (for H2O liquid)          BAK
     
    278276!==============================================================================================================================
    279277   ALLOCATE(ttr(nqtot))
    280    jq = nqtrue+1; tracers(:)%iadv = -1
     278   jq = nqtrue+1
    281279   DO iq = 1, nqtrue
    282280      t1 => tracers(iq)
     
    289287      IF(iad == -1) CALL abort_gcm(modname, msg1, 1)
    290288
    291       !--- SET FIELDS longName, iadv, isAdvected, isInPhysics
     289      !--- SET FIELDS longName and iadv
    292290      t1%longName   = t1%name; IF(iad > 0) t1%longName=TRIM(t1%name)//descrq(iad)
    293291      t1%iadv       = iad
    294       t1%isAdvected = iad >= 0
    295       t1%isInPhysics= delPhase(t1%gen0Name) /= 'H2O' .OR. t1%component /= 'lmdz' !=== MORE EXCEPTIONS ? CO2i, SURSAT CLOUD H2O
    296292      ttr(iq)       = t1
    297293
     
    307303      ttr(jq+1:jq+nm)%longName    = [ (TRIM(t1%longName)//'-'//TRIM(suff(im)), im=1, nm) ]
    308304      ttr(jq+1:jq+nm)%iadv        = [ (-iad,    im=1, nm) ]
    309       ttr(jq+1:jq+nm)%isAdvected  = [ (.FALSE., im=1, nm) ]
    310305      jq = jq + nm
    311306   END DO
Note: See TracChangeset for help on using the changeset viewer.