Changeset 3550 for LMDZ6


Ignore:
Timestamp:
Aug 7, 2019, 2:31:31 PM (5 years ago)
Author:
oboucher
Message:

reordering loops for tracer output

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r3523 r3550  
    22342234!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    22352235       IF (iflag_phytrac == 1 ) then
    2236        IF (nqtot.GE.nqo+1) THEN
    2237           DO iq=nqo+1, nqtot
    2238             IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN
     2236         IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN
     2237           DO iq=nqo+1, nqtot
    22392238             !--3D fields
    22402239             CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
     
    22632262             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
    22642263#endif
    2265             ENDIF   !--type_trac
     2264           ENDDO !--iq
     2265         ENDIF   !--type_trac
    22662266!
    2267             IF (type_trac == 'co2i') THEN
     2267         IF (type_trac == 'co2i') THEN
     2268           DO iq=nqo+1, nqtot
    22682269             !--3D fields
    22692270             CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
     
    22802281             ENDIF
    22812282             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
    2282              !--CO2 net fluxes
    2283              CALL histwrite_phy(o_flx_co2_land,  fco2_land)
    2284              CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean)
    2285              CALL histwrite_phy(o_flx_co2_ff,    fco2_ff)
    2286              CALL histwrite_phy(o_flx_co2_bb,    fco2_bb)
    2287             ENDIF !--type_trac co2i
    2288 
    2289           ENDDO !--iq
    2290        ENDIF !--nqtot > nqo+1
    2291 
    2292        IF (type_trac == 'repr') THEN
     2283           ENDDO !--iq
     2284           !--CO2 net fluxes
     2285           CALL histwrite_phy(o_flx_co2_land,  fco2_land)
     2286           CALL histwrite_phy(o_flx_co2_ocean, fco2_ocean)
     2287           CALL histwrite_phy(o_flx_co2_ff,    fco2_ff)
     2288           CALL histwrite_phy(o_flx_co2_bb,    fco2_bb)
     2289         ENDIF !--type_trac co2i
     2290
     2291         IF (type_trac == 'repr') THEN
    22932292#ifdef REPROBUS
    22942293           DO iq=1,nbnas
     
    22962295           ENDDO
    22972296#endif
    2298        ENDIF
     2297         ENDIF
    22992298
    23002299       ENDIF   !(iflag_phytrac==1)
Note: See TracChangeset for help on using the changeset viewer.