Changeset 4264


Ignore:
Timestamp:
Sep 21, 2022, 9:43:30 AM (20 months ago)
Author:
dcugnet
Message:

Exceptions for REPROBUS with old conventions for HNO3 in the start.nc file were wrongly placed.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90

    r4263 r4264  
    159159    var = tracers(iq)%name
    160160    oldVar = new2oldName(var)
    161     !--------------------------------------------------------------------------------------------------------------------------
    162     IF(NF90_INQ_VARID(fID, var, vID) == NF90_NoErr) THEN                                 !=== REGULAR CASE
    163       CALL get_var2(var,q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
    164     !--------------------------------------------------------------------------------------------------------------------------
    165     ELSE IF(NF90_INQ_VARID(fID, oldVar, vID) == NF90_NoErr) THEN                         !=== OLD NAME
    166       CALL msg('Tracer <'//TRIM(var)//'> is missing => initialized to <'//TRIM(oldVar)//'>', modname)
    167       CALL get_var2(oldVar, q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
    168     !--------------------------------------------------------------------------------------------------------------------------
    169 #ifdef INCA
    170     ELSE IF(NF90_INQ_VARID(fID, 'OX',   vID) == NF90_NoErr .AND. var == 'O3') THEN       !=== INCA: OX INSTEAD OF O3
    171       CALL msg('Tracer <O3> is missing => initialized to <OX>', modname)
    172       CALL get_var2( 'OX' , q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
    173 #endif
    174     !--------------------------------------------------------------------------------------------------------------------------
    175161#ifdef REPROBUS
    176     ELSE IF(lOldHNO3 .AND. var == 'HNO3') THEN                                           !=== REPROBUS: HNO3_g INSTEAD OF HNO3
     162    IF(lOldHNO3 .AND. var == 'HNO3') THEN                                                !=== REPROBUS: HNO3_g INSTEAD OF HNO3
    177163      CALL msg('Tracer <HNO3_g> is used for <HNO3> (old convention)', modname)
    178164      CALL get_var2( 'HNO3_g', q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
     
    181167      CALL get_var2( 'HNO3', q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
    182168#endif
     169    !--------------------------------------------------------------------------------------------------------------------------
     170    ELSE IF(NF90_INQ_VARID(fID, var, vID) == NF90_NoErr) THEN                            !=== REGULAR CASE
     171      CALL get_var2(var,q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
     172    !--------------------------------------------------------------------------------------------------------------------------
     173    ELSE IF(NF90_INQ_VARID(fID, oldVar, vID) == NF90_NoErr) THEN                         !=== OLD NAME
     174      CALL msg('Tracer <'//TRIM(var)//'> is missing => initialized to <'//TRIM(oldVar)//'>', modname)
     175      CALL get_var2(oldVar, q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
     176    !--------------------------------------------------------------------------------------------------------------------------
     177#ifdef INCA
     178    ELSE IF(NF90_INQ_VARID(fID, 'OX',   vID) == NF90_NoErr .AND. var == 'O3') THEN       !=== INCA: OX INSTEAD OF O3
     179      CALL msg('Tracer <O3> is missing => initialized to <OX>', modname)
     180      CALL get_var2( 'OX' , q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:)
     181#endif
     182    !--------------------------------------------------------------------------------------------------------------------------
    183183    ELSE IF(tracers(iq)%iso_iGroup == iH2O .AND. niso > 0) THEN                          !=== WATER ISOTOPES
    184184      iName    = tracers(iq)%iso_iName
Note: See TracChangeset for help on using the changeset viewer.