Changeset 4106


Ignore:
Timestamp:
Mar 22, 2022, 11:35:36 AM (2 years ago)
Author:
acozic
Message:

Add condition to write some variables in output (without inca, or only with new physics)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r4089 r4106  
    942942       CALL histwrite_phy(o_SWupTOAclr, zx_tmp_fi2d)
    943943
    944        IF (vars_defined) THEN
    945           zx_tmp_fi2d(:) = swupc0(:,klevp1)*swradcorr(:)
    946        ENDIF
    947        CALL histwrite_phy(o_SWupTOAcleanclr, zx_tmp_fi2d)
     944       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     945          IF (vars_defined) THEN
     946             zx_tmp_fi2d(:) = swupc0(:,klevp1)*swradcorr(:)
     947          ENDIF
     948          CALL histwrite_phy(o_SWupTOAcleanclr, zx_tmp_fi2d)
     949       ENDIF
    948950
    949951       IF (vars_defined) THEN
    … …  
    10121014       CALL histwrite_phy(o_SWupSFCclr, zx_tmp_fi2d)
    10131015
    1014        IF (vars_defined) THEN
    1015           zx_tmp_fi2d(:) = swupc0(:,1)*swradcorr(:)
    1016        ENDIF
    1017        CALL histwrite_phy(o_SWupSFCcleanclr, zx_tmp_fi2d)
     1016       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1017          IF (vars_defined) THEN
     1018             zx_tmp_fi2d(:) = swupc0(:,1)*swradcorr(:)
     1019          ENDIF
     1020          CALL histwrite_phy(o_SWupSFCcleanclr, zx_tmp_fi2d)
     1021       ENDIF
    10181022
    10191023       IF (vars_defined) THEN
    … …  
    10271031       CALL histwrite_phy(o_SWdnSFCclr, zx_tmp_fi2d)
    10281032
    1029        IF (vars_defined) THEN
    1030           zx_tmp_fi2d(:) = swdnc0(:,1)*swradcorr(:)
    1031        ENDIF
    1032        CALL histwrite_phy(o_SWdnSFCcleanclr, zx_tmp_fi2d)
     1033       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1034          IF (vars_defined) THEN
     1035             zx_tmp_fi2d(:) = swdnc0(:,1)*swradcorr(:)
     1036          ENDIF
     1037          CALL histwrite_phy(o_SWdnSFCcleanclr, zx_tmp_fi2d)
     1038       ENDIF
    10331039
    10341040       CALL histwrite_phy(o_fdiffSWdnSFC, solswfdiff)
    … …  
    10471053       CALL histwrite_phy(o_LWdnSFCclr, sollwdownclr)
    10481054
    1049        IF (vars_defined) THEN
    1050           zx_tmp_fi2d(:) = lwupc0(:,klevp1)
    1051        ENDIF
    1052        CALL histwrite_phy(o_LWupTOAcleanclr, zx_tmp_fi2d)
    1053        IF (vars_defined) THEN
    1054           zx_tmp_fi2d(:) = -1.*lwdnc0(:,1)
    1055        ENDIF
    1056        CALL histwrite_phy(o_LWdnSFCcleanclr, zx_tmp_fi2d)
     1055       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1056          IF (vars_defined) THEN
     1057             zx_tmp_fi2d(:) = lwupc0(:,klevp1)
     1058          ENDIF
     1059          CALL histwrite_phy(o_LWupTOAcleanclr, zx_tmp_fi2d)
     1060       ENDIF
     1061
     1062       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1063          IF (vars_defined) THEN
     1064             zx_tmp_fi2d(:) = -1.*lwdnc0(:,1)
     1065          ENDIF
     1066          CALL histwrite_phy(o_LWdnSFCcleanclr, zx_tmp_fi2d)
     1067       ENDIF
    10571068
    10581069       CALL histwrite_phy(o_bils, bils)
    … …  
    15581569!This is warranted by treating INCA aerosols as offline aerosols
    15591570       IF (flag_aerosol.GT.0) THEN
    1560           CALL histwrite_phy(o_od443aer, od443aer)
    1561           CALL histwrite_phy(o_od550aer, od550aer)
    1562           CALL histwrite_phy(o_od865aer, od865aer)
    1563           CALL histwrite_phy(o_abs550aer, abs550aer)
    1564           CALL histwrite_phy(o_od550lt1aer, od550lt1aer)
    1565           CALL histwrite_phy(o_sconcso4, sconcso4)
    1566           CALL histwrite_phy(o_sconcno3, sconcno3)
    1567           CALL histwrite_phy(o_sconcoa, sconcoa)
    1568           CALL histwrite_phy(o_sconcbc, sconcbc)
    1569           CALL histwrite_phy(o_sconcss, sconcss)
    1570           CALL histwrite_phy(o_sconcdust, sconcdust)
    1571           CALL histwrite_phy(o_concso4, concso4)
    1572           CALL histwrite_phy(o_concno3, concno3)
    1573           CALL histwrite_phy(o_concoa, concoa)
    1574           CALL histwrite_phy(o_concbc, concbc)
    1575           CALL histwrite_phy(o_concss, concss)
    1576           CALL histwrite_phy(o_concdust, concdust)
    1577           CALL histwrite_phy(o_loadso4, loadso4)
    1578           CALL histwrite_phy(o_loadoa, loadoa)
    1579           CALL histwrite_phy(o_loadbc, loadbc)
    1580           CALL histwrite_phy(o_loadss, loadss)
    1581           CALL histwrite_phy(o_loaddust, loaddust)
    1582           CALL histwrite_phy(o_loadno3, loadno3)
    1583           CALL histwrite_phy(o_dryod550aer, dryod550aer)
    1584           DO naero = 1, naero_tot-1
    1585              CALL histwrite_phy(o_drytausumaero(naero),drytausum_aero(:,naero))
    1586           END DO
     1571          IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1572
     1573             CALL histwrite_phy(o_od443aer, od443aer)
     1574             CALL histwrite_phy(o_od550aer, od550aer)
     1575             CALL histwrite_phy(o_od865aer, od865aer)
     1576             CALL histwrite_phy(o_abs550aer, abs550aer)
     1577             CALL histwrite_phy(o_od550lt1aer, od550lt1aer)
     1578             CALL histwrite_phy(o_sconcso4, sconcso4)
     1579             CALL histwrite_phy(o_sconcno3, sconcno3)
     1580             CALL histwrite_phy(o_sconcoa, sconcoa)
     1581             CALL histwrite_phy(o_sconcbc, sconcbc)
     1582             CALL histwrite_phy(o_sconcss, sconcss)
     1583             CALL histwrite_phy(o_sconcdust, sconcdust)
     1584             CALL histwrite_phy(o_concso4, concso4)
     1585             CALL histwrite_phy(o_concno3, concno3)
     1586             CALL histwrite_phy(o_concoa, concoa)
     1587             CALL histwrite_phy(o_concbc, concbc)
     1588             CALL histwrite_phy(o_concss, concss)
     1589             CALL histwrite_phy(o_concdust, concdust)
     1590             CALL histwrite_phy(o_loadso4, loadso4)
     1591             CALL histwrite_phy(o_loadoa, loadoa)
     1592             CALL histwrite_phy(o_loadbc, loadbc)
     1593             CALL histwrite_phy(o_loadss, loadss)
     1594             CALL histwrite_phy(o_loaddust, loaddust)
     1595             CALL histwrite_phy(o_loadno3, loadno3)
     1596             CALL histwrite_phy(o_dryod550aer, dryod550aer)
     1597             DO naero = 1, naero_tot-1
     1598                CALL histwrite_phy(o_drytausumaero(naero),drytausum_aero(:,naero))
     1599             END DO
     1600          ENDIF
    15871601       ENDIF
    15881602       !--STRAT AER
    … …  
    16781692          CALL histwrite_phy(o_solswad0, zx_tmp_fi2d)
    16791693         
    1680           CALL histwrite_phy(o_toplwad, toplwad_aero)
    1681           CALL histwrite_phy(o_toplwad0, toplwad0_aero)
    1682           CALL histwrite_phy(o_sollwad, sollwad_aero)
    1683           CALL histwrite_phy(o_sollwad0, sollwad0_aero)
     1694          IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1695
     1696             CALL histwrite_phy(o_toplwad, toplwad_aero)
     1697             CALL histwrite_phy(o_toplwad0, toplwad0_aero)
     1698             CALL histwrite_phy(o_sollwad, sollwad_aero)
     1699             CALL histwrite_phy(o_sollwad0, sollwad0_aero)
     1700          ENDIF
    16841701          !====MS forcing diagnostics
    16851702          !ym warning : topsw_aero, solsw_aero, topsw0_aero, solsw0_aero are not defined by model
    … …  
    17521769       ! Champs 3D:
    17531770       IF (ok_ade .OR. ok_aie) then
    1754           CALL histwrite_phy(o_ec550aer, ec550aer)
     1771          IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     1772             CALL histwrite_phy(o_ec550aer, ec550aer)
     1773          ENDIF
    17551774       ENDIF
    17561775       CALL histwrite_phy(o_lwcon, flwc)
    … …  
    22092228       CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1)
    22102229
    2211        IF (vars_defined) THEN
    2212          DO k=1, klevp1
    2213            zx_tmp_fi3d1(:,k)=swupc0(:,k)*swradcorr(:)
    2214          ENDDO
    2215        ENDIF
    2216        CALL histwrite_phy(o_rsucsaf, zx_tmp_fi3d1)
     2230       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     2231          IF (vars_defined) THEN
     2232             DO k=1, klevp1
     2233                zx_tmp_fi3d1(:,k)=swupc0(:,k)*swradcorr(:)
     2234             ENDDO
     2235          ENDIF
     2236          CALL histwrite_phy(o_rsucsaf, zx_tmp_fi3d1)
     2237       ENDIF
    22172238
    22182239       IF (vars_defined) THEN
    … …  
    22232244       CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
    22242245
    2225 
    2226        IF (vars_defined) THEN
    2227          DO k=1, klevp1
    2228            zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
    2229          ENDDO
    2230        ENDIF
    2231        CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
     2246       IF ((type_trac/='inca') .OR. (config_inca=='aeNP')) THEN
     2247          IF (vars_defined) THEN
     2248             DO k=1, klevp1
     2249                zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
     2250             ENDDO
     2251          ENDIF
     2252          CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
     2253       ENDIF
    22322254
    22332255       CALL histwrite_phy(o_rlu, lwup)
Note: See TracChangeset for help on using the changeset viewer.