Ignore:
Timestamp:
May 22, 2007, 3:55:22 PM (17 years ago)
Author:
lsce
Message:

correction de bug pour le couplage avec le modele inca - ACo

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histrac.h

    r740 r759  
    1717
    1818#ifdef INCA_AER
    19          CALL histbeg("histrac_aer", iim,zx_lon, jjm+1,zx_lat,
     19         CALL histbeg("histrac_aer", iim,zx_lon(:,1), jjm+1,zx_lat(1,:),
    2020     .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
    2121     .                 nhori, nid_tra2)
    … …  
    2525     .                 nhori, nid_tra3)
    2626
    27          call histvert(nid_tra2, "presnivs", "presnivs", "mb",
     27         call histvert(nid_tra2, "presnivs", "Vertical levels", "mb",
    2828     .                 klev, presnivs, nvert)
    2929         call histvert(nid_tra3, "presnivs", "presnivs", "mb",
    … …  
    9292#ifdef INCA
    9393#ifdef INCAINFO
     94
     95#ifdef INCA_CH4
    9496         DO it=1, phtcnt
    9597         WRITE(str2,'(i2.2)') it
    … …  
    111113     .                "ave(X)", zsto,zout)
    112114         ENDDO
    113 
     115#endif
     116#ifdef INCA_NMHC
     117         DO it=1, phtcnt
     118         WRITE(str2,'(i2.2)') it
     119         CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",
     120     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     121     .                "ave(X)", zsto,zout)
     122         ENDDO
     123         DO it=1, hetcnt
     124         WRITE(str2,'(i2.2)') it
     125         CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",
     126     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     127     .                "ave(X)", zsto,zout)
     128         ENDDO
     129         DO it=1, extcnt
     130         WRITE(str2,'(i2.2)') it
     131         CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",
     132     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     133     .                "ave(X)", zsto,zout)
     134         ENDDO
     135#endif
    114136         DO it=1, nfs
    115137         WRITE(str2,'(i2.2)') it
    … …  
    120142
    121143#else
     144
     145#ifdef INCA_NMHC
    122146         CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
    123147     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    … …  
    139163     .                "ave(X)", zsto,zout)
    140164#endif
    141 
    142          DO it=1, grpcnt
    143          CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
    144      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    145      .                "ave(X)", zsto,zout)
    146          ENDDO
     165#ifdef INCA_CH4
     166         CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
     167     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     168     .                "ave(X)", zsto,zout)
     169         CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",
     170     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     171     .                "ave(X)", zsto,zout)
     172         CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",
     173     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     174     .                "ave(X)", zsto,zout)
     175         CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",
     176     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     177     .                "ave(X)", zsto,zout)
     178         CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",
     179     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     180     .                "ave(X)", zsto,zout)
     181         CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",
     182     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     183     .                "ave(X)", zsto,zout)
     184#endif
     185#endif
     186
     187!         DO it=1, grpcnt
     188!         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
     189!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     190!     .                "ave(X)", zsto,zout)
     191!         ENDDO
    147192#endif
    148193
    … …  
    349394#endif
    350395
     396#ifdef INCA_CH4
     397         DO it=1, grpcnt
     398         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
     399     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     400     .                "ave(X)", zsto,zout)
     401         ENDDO
     402
     403#endif
     404
     405
    351406#ifdef INCA_NMHC
    352407         CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
    … …  
    371426     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    372427     .                "ave(X)", zsto,zout)
     428
     429         DO it=1, grpcnt
     430         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
     431     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     432     .                "ave(X)", zsto,zout)
     433         ENDDO
     434
     435
    373436#endif
    374437! Fin ajout Anne
    … …  
    376439         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
    377440     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    378      .                "inst(X)",  zsto,zout)
     441     .                "inst(X)", zout,zout)
    379442
    380443         CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
    381444     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    382      .                "inst(X)",  zsto,zout)
     445     .                "inst(X)",  zout,zout)
    383446         CALL histdef(nid_tra, "psrf1", "nature sol", "-",
    384447     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    385      .                "inst(X)",  zsto,zout)
     448     .                "inst(X)",  zout,zout)
    386449         CALL histdef(nid_tra, "psrf2", "nature sol", "-",
    387450     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    388      .                "inst(X)",  zsto,zout)
     451     .                "inst(X)",  zout,zout)
    389452         CALL histdef(nid_tra, "psrf3", "nature sol", "-",
    390453     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    391      .                "inst(X)",  zsto,zout)
     454     .                "inst(X)",  zout,zout)
    392455         CALL histdef(nid_tra, "psrf4", "nature sol", "-",
    393456     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    394      .                "inst(X)",  zsto,zout)
     457     .                "inst(X)",  zout,zout)
    395458         CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
    396459     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    397      .                "inst(X)",  zsto,zout)
     460     .                "inst(X)",  zout,zout)
    398461         CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
    399462     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    400      .                "inst(X)",  zsto,zout)
     463     .                "inst(X)",  zout,zout)
    401464         CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
    402465     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    403      .                "inst",  zsto,zout)
     466     .                "inst(X)",  zout,zout)
    404467         CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
    405468     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
    406      .                "inst(X)",  zsto,zout)
     469     .                "inst(X)",  zout,zout)
    407470         CALL histdef(nid_tra, "pplay", "flux u mont","-",
    408471     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    409      .                "inst(X)", zsto,zout)
     472     .                "inst(X)", zout,zout)
    410473         CALL histdef(nid_tra, "t", "flux u mont","-",
    411474     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    412      .                "inst(X)", zsto,zout)
     475     .                "inst(X)", zout,zout)
    413476         CALL histdef(nid_tra, "mfu", "flux u mont","-",
    414477     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
Note: See TracChangeset for help on using the changeset viewer.