Ignore:
Timestamp:
May 4, 2005, 5:11:29 PM (19 years ago)
Author:
Laurent Fairhead
Message:

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r633 r634  
    44         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
    55c
    6          CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon)
    7          DO i = 1, iim
    8             zx_lon(i,1) = xlon(i+1)
    9             zx_lon(i,jjm+1) = xlon(i+1)
    10          ENDDO
    11          CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat)
    12          CALL histbeg("histrac", iim,zx_lon(:,1), jjm+1,zx_lat(1,:),
    13      .                 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,
     6cym         CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon)
     7cym         DO i = 1, iim
     8cym            zx_lon(i,1) = xlon(i+1)
     9cym            zx_lon(i,jjm+1) = xlon(i+1)
     10cym         ENDDO
     11cym         CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat)
     12cym         CALL histbeg("histrac", iim,zx_lon(:,1), jjm+1,zx_lat(1,:),
     13cym     .                 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,
     14cym     .                 nhori, nid_tra)
     15         CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,
    1416     .                 nhori, nid_tra)
    1517         CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
    … …  
    1719
    1820#ifdef INCA_AER
    19          CALL histbeg("histrac_aer", iim,zx_lon, jjm+1,zx_lat,
    20      .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
     21cym         CALL histbeg("histrac_aer", iim,zx_lon, jjm+1,zx_lat,
     22cym     .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
     23cym     .                 nhori, nid_tra2)
     24         CALL histbeg_phy("histrac_aer", itau_phy, zjulian, pdtphys,
    2125     .                 nhori, nid_tra2)
    2226
    23          CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,
    24      .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
     27cym         CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,
     28cym     .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
     29cym     .                 nhori, nid_tra3)
     30         CALL histbeg_phy("histrac_inst",itau_phy, zjulian, pdtphys,
    2531     .                 nhori, nid_tra3)
    2632
    … …  
    4248c
    4349         CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
    44      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     50     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4551     .                "once",  zsto,zout)
    4652c
    4753         CALL histdef(nid_tra, "aire", "Grid area", "-",
    48      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     54     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4955     .                "once",  zsto,zout)
    5056#ifdef INCA
    5157         CALL histdef(nid_tra, "ps", "Surface pressure", "Pa",
    52      .                iim,jjm+1,nhori, 1,1,1,-99, 32,
     58     .                iim,jjphy_nb,nhori, 1,1,1,-99, 32,
    5359     .                "ave(X)", zsto,zout)
    5460
    5561         CALL histdef(nid_tra, "ptrop", "Tropopause pressure", "Pa",
    56      .                iim,jjm+1,nhori, 1,1,1,-99, 32,
     62     .                iim,jjphy_nb,nhori, 1,1,1,-99, 32,
    5763     .                "ave(X)", zsto,zout)
    5864
    5965C   3d FIELDS
    6066         CALL histdef(nid_tra, "temp", "Air temperature", "K",
    61      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     67     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    6268     .                "ave(X)", zsto,zout)
    6369
    6470         CALL histdef(nid_tra, "u", "zonal wind component", "m/s",
    65      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     71     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    6672     .                "ave(X)", zsto,zout)
    6773
    6874         CALL histdef(nid_tra, "v", "zonal wind component", "m/s",
    69      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     75     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    7076     .                "ave(X)", zsto,zout)
    7177
    7278         CALL histdef(nid_tra, "h2o", "Specific Humidity", "MMR",
    73      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     79     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    7480     .                "ave(X)", zsto,zout)
    7581
    7682         CALL histdef(nid_tra, "pmid", "Pressure", "Pa",
    77      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     83     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    7884     .                "ave(X)", zsto,zout)
    7985
    8086         CALL histdef(nid_tra, "pdel", "Delta Pressure", "Pa",
    81      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     87     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    8288     .                "ave(X)", zsto,zout)
    8389
    8490! MS info
    8591         CALL histdef(nid_tra, "airm", "Air mass", "kg",
    86      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     92     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    8793     .                "ave(X)", zsto,zout)
    8894
    … …  
    9298         WRITE(str2,'(i2.2)') it
    9399         CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",
    94      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     100     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    95101     .                "ave(X)", zsto,zout)
    96102         ENDDO
    … …  
    98104         WRITE(str2,'(i2.2)') it
    99105         CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",
    100      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     106     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    101107     .                "ave(X)", zsto,zout)
    102108         ENDDO
    … …  
    105111         WRITE(str2,'(i2.2)') it
    106112         CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",
    107      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     113     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    108114     .                "ave(X)", zsto,zout)
    109115         ENDDO
    … …  
    112118         WRITE(str2,'(i2.2)') it
    113119         CALL histdef(nid_tra, "INV"//str2, "INV"//str2, "CM-3",
    114      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     120     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    115121     .                "ave(X)", zsto,zout)
    116122         ENDDO
    … …  
    118124#else
    119125         CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
    120      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     126     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    121127     .                "ave(X)", zsto,zout)
    122128         CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",
    123      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     129     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    124130     .                "ave(X)", zsto,zout)
    125131         CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",
    126      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     132     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    127133     .                "ave(X)", zsto,zout)
    128134         CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",
    129      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     135     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    130136     .                "ave(X)", zsto,zout)
    131137         CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",
    132      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     138     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    133139     .                "ave(X)", zsto,zout)
    134140         CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",
    135      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     141     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    136142     .                "ave(X)", zsto,zout)
    137143#endif
    … …  
    139145         DO it=1, grpcnt
    140146         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
    141      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     147     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    142148     .                "ave(X)", zsto,zout)
    143149         ENDDO
    … …  
    147153
    148154        CALL histdef(nid_tra2, "scavcoef_st","scavcoef_st", "S-1",
    149      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     155     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    150156     .                "ave(X)", zsto,zout)
    151157        CALL histdef(nid_tra2, "scavcoef_cv","scavcoef_cv", "S-1",
    152      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     158     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    153159     .                "ave(X)", zsto,zout)
    154160        CALL histdef(nid_tra2, "AngstroemComp","AngstroemComp",
    155      .      "angs comp", iim,jjm+1,nhori, 1,1,1, -99, 32,
     161     .      "angs comp", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    156162     .      "ave(X)", zsto,zout)
    157163
    … …  
    164170
    165171         CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),
    166      .           "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     172     .           "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    167173     .           "ave(X)", zsto,zout)
    168174         CALL histdef(nid_tra, "Dep_"//solsym(it), "Dep_"//solsym(it),
    169      .           "cm/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     175     .           "cm/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    170176     .           "ave(X)", zsto,zout)
    171177#ifdef INCA_AER
    172178        IF  ((it .ge. trmx) .and. (it .le. trnx)) then
    173179          CALL histdef(nid_tra2, "Sed_"//solsym(it), "Sed_"//solsym(it),
    174      .      "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     180     .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    175181     .      "ave(X)", zsto,zout)
    176182          CALL histdef(nid_tra2, "Dry_"//solsym(it), "Dry_"//solsym(it),
    177      .      "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     183     .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    178184     .      "ave(X)", zsto,zout)
    179185          CALL histdef(nid_tra2, "Wet_"//solsym(it), "Wet_"//solsym(it),
    180      .      "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     186     .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    181187     .      "ave(X)", zsto,zout)
    182188          CALL histdef(nid_tra2, "WetST_"//solsym(it), "WetST_"//solsym(it),
    183      .      "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     189     .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    184190     .      "ave(X)", zsto,zout)
    185191          CALL histdef(nid_tra2, "WetCV_"//solsym(it), "WetCV_"//solsym(it),
    186      .      "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,
     192     .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    187193     .      "ave(X)", zsto,zout)
    188194          CALL histdef(nid_tra2, "Emi_alt_"//solsym(it), "Emi_alt_"//solsym(it),
    189      .      "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     195     .      "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    190196     .      "ave(X)", zsto,zout)
    191197          CALL histdef(nid_tra2, "Load_"//solsym(it), "Load_"//solsym(it),
    192      .      "kg/m2", iim,jjm+1,nhori, 1,1,1, -99, 32,
     198     .      "kg/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    193199     .      "ave(X)", zsto,zout)
    194200          CALL histdef(nid_tra2, "SConc_"//solsym(it), "SConc_"//solsym(it),
    195      .      "kg/m3", iim,jjm+1,nhori, 1,1,1, -99, 32,
     201     .      "kg/m3", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    196202     .      "ave(X)", zsto,zout)
    197203          do la=1,las
    198204           CALL histdef(nid_tra2, "OD_"//cla(la)//solsym(it),
    199205     .                            "OD_"//cla(la)//solsym(it),
    200      .        "opt. depth", iim,jjm+1,nhori, 1,1,1, -99, 32,
     206     .        "opt. depth", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    201207     .        "ave(X)", zsto,zout)
    202208          enddo
    203209          CALL histdef(nid_tra2, "MD_"//solsym(it), "MD_"//solsym(it),
    204      .      "median diameter", iim,jjm+1,nhori, klev,1,klev, nvert, 32,
     210     .      "median diameter", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    205211     .      "ave(X)", zsto,zout)
    206212          CALL histdef(nid_tra3, "Inst_Load_"//solsym(it),
    207213     .                          "Inst_Load_"//solsym(it),
    208      .      "kg/m2", iim,jjm+1,nhori, 1,1,1, -99, 32,
     214     .      "kg/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    209215     .      "inst(X)", zout,zout)
    210216        endif
    211217#endif
    212218         CALL histdef(nid_tra, solsym(it), solsym(it), "VMR",
    213      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     219     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    214220     .                "ave(X)", zsto,zout)
    215221#else
    … …  
    217223         iiq=niadv(iq)
    218224         CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga",
    219      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     225     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    220226     .                "ave(X)", zsto,zout)
    221227         if (lessivage) THEN
    222228         CALL histdef(nid_tra, "fl"//tnom(iq),"Flux "//ttext(iiq),
    223      .              "U/m2/s",iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     229     .              "U/m2/s",iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    224230     .              "ave(X)", zsto,zout)
    225231         endif
    … …  
    228234         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq),
    229235     .                "tendance thermique"// ttext(iiq), "?",
    230      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     236     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    231237     .                "ave(X)", zsto,zout)
    232238         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq),
    233239     .                "tendance convection"// ttext(iiq), "?",
    234      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     240     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    235241     .                "ave(X)", zsto,zout)
    236242         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq),
    237243     .                "tendance couche limite"// ttext(iiq), "?",
    238      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     244     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    239245     .                "ave(X)", zsto,zout)
    240246c---fin Olivia   
    … …  
    246252!#ifdef INCA_CH4
    247253         CALL histdef(nid_tra, "O3_column", "O3_column",
    248      .           "DU", iim,jjm+1,nhori, 1,1,1, -99, 32,
     254     .           "DU", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    249255     .           "ave(X)", zsto,zout)
    250256         CALL histdef(nid_tra, "CO_column", "CO_column",
    251      .           "10^18 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32,
     257     .           "10^18 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    252258     .           "ave(X)", zsto,zout)
    253259         CALL histdef(nid_tra, "CH4_column", "CH4_column",
    254      .           "10^18 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32,
     260     .           "10^18 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    255261     .           "ave(X)", zsto,zout)
    256262         CALL histdef(nid_tra, "NO2_column", "NO2_column",
    257      .           "10^15 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32,
     263     .           "10^15 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    258264     .           "ave(X)", zsto,zout)
    259265         CALL histdef(nid_tra, "O3_ste", "O3_ste",
    260      .           "CM-2 S-1", iim,jjm+1,nhori, 1,1,1, -99, 32,
     266     .           "CM-2 S-1", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    261267     .           "ave(X)", zsto,zout)
    262268         CALL histdef(nid_tra, "O3_prod", "O3_prod", "CM-3 S-1",
    263      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     269     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    264270     .                "ave(X)", zsto,zout)
    265271         CALL histdef(nid_tra, "O3_loss", "O3_loss", "CM-3 S-1",
    266      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     272     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    267273     .                "ave(X)", zsto,zout)
    268274
    … …  
    281287! for sulfur cycle
    282288         CALL histdef(nid_tra, "SO2_p_dmsoh", "SO2_p_dmsoh", "CM-3 S-1",
    283      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     289     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    284290     .                "ave(X)", zsto,zout)
    285291         CALL histdef(nid_tra, "SO2_p_dmsno3", "SO2_p_dmsno3", "CM-3 S-1",
    286      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     292     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    287293     .                "ave(X)", zsto,zout)
    288294         CALL histdef(nid_tra, "SO2_p_h2soh", "SO2_p_h2soh", "CM-3 S-1",
    289      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     295     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    290296     .                "ave(X)", zsto,zout)
    291297         CALL histdef(nid_tra, "SO2_p_dmsooh", "SO2_p_dmsooh", "CM-3 S-1",
    292      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     298     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    293299     .                "ave(X)", zsto,zout)
    294300         CALL histdef(nid_tra, "DMSO_p_dmsoh", "DMSO_p_dmsoh", "CM-3 S-1",
    295      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     301     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    296302     .                "ave(X)", zsto,zout)
    297303         CALL histdef(nid_tra, "ASMSAM_p_dmsooh", "ASMSAM_p_dmsooh", "CM-3 S-1",
    298      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     304     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    299305     .                "ave(X)", zsto,zout)
    300306         CALL histdef(nid_tra, "ASSO4M_p_so2oh", "ASSO4M_p_so2oh", "CM-3 S-1",
    301      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     307     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    302308     .                "ave(X)", zsto,zout)
    303309         CALL histdef(nid_tra, "ASSO4M_p_so2h2o2", "ASSO4M_p_h2o2","CM-3 S-1",
    304      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     310     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    305311     .                "ave(X)", zsto,zout)
    306312         CALL histdef(nid_tra, "ASSO4M_p_so2o3", "ASSO4M_p_so2o3", "CM-3 S-1",
    307      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     313     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    308314     .                "ave(X)", zsto,zout)
    309315
    310316         CALL histdef(nid_tra, "Wet3D_SO2", "Wet3D_SO2",
    311      .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     317     .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    312318     .           "ave(X)", zsto,zout)
    313319         CALL histdef(nid_tra, "Wet3D_DMS", "Wet3D_DMS",
    314      .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     320     .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    315321     .           "ave(X)", zsto,zout)
    316322         CALL histdef(nid_tra, "Wet3D_HNO3", "Wet3D_HNO3",
    317      .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     323     .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    318324     .           "ave(X)", zsto,zout)
    319325         CALL histdef(nid_tra, "Wet3D_H2O2", "Wet3D_H2O2",
    320      .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     326     .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    321327     .           "ave(X)", zsto,zout)
    322328
    323329         CALL histdef(nid_tra, "PH_HIST", "pH over physics time step", "1",
    324      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     330     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    325331     .                "ave(X)", zsto,zout)
    326332#endif
    … …  
    328334#ifdef INCA_NMHC
    329335         CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
    330      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     336     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    331337     .                "ave(X)", zsto,zout)
    332338         CALL histdef(nid_tra, "CO2_nmhcprod", "CO2_nmhcprod", "CM-3 S-1",
    333      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     339     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    334340     .                "ave(X)", zsto,zout)
    335341         CALL histdef(nid_tra, "CO2_radicalprod", "CO2_radicalprod", "CM-3 S-1",
    336      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     342     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    337343     .                "ave(X)", zsto,zout)
    338344         CALL histdef(nid_tra, "HNO3_prod", "HNO3_prod", "CM-3 S-1",
    339      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     345     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    340346     .                "ave(X)", zsto,zout)
    341347         CALL histdef(nid_tra, "HNO3_loss", "HNO3_loss", "CM-3 S-1",
    342      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     348     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    343349     .                "ave(X)", zsto,zout)
    344350         CALL histdef(nid_tra, "CO_prod", "CO_prod", "CM-3 S-1",
    345      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     351     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    346352     .                "ave(X)", zsto,zout)
    347353         CALL histdef(nid_tra, "CO_loss", "CO_loss", "CM-3 S-1",
    348      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     354     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    349355     .                "ave(X)", zsto,zout)
    350356#endif
    … …  
    352358#else
    353359         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
    354      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     360     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    355361     .                "inst(X)",  zsto,zout)
    356362
    357363         CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
    358      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     364     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    359365     .                "inst(X)",  zsto,zout)
    360366         CALL histdef(nid_tra, "psrf1", "nature sol", "-",
    361      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     367     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    362368     .                "inst(X)",  zsto,zout)
    363369         CALL histdef(nid_tra, "psrf2", "nature sol", "-",
    364      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     370     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    365371     .                "inst(X)",  zsto,zout)
    366372         CALL histdef(nid_tra, "psrf3", "nature sol", "-",
    367      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     373     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    368374     .                "inst(X)",  zsto,zout)
    369375         CALL histdef(nid_tra, "psrf4", "nature sol", "-",
    370      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     376     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    371377     .                "inst(X)",  zsto,zout)
    372378         CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
    373      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     379     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    374380     .                "inst(X)",  zsto,zout)
    375381         CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
    376      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     382     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    377383     .                "inst(X)",  zsto,zout)
    378384         CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
    379      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     385     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    380386     .                "inst",  zsto,zout)
    381387         CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
    382      .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     388     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    383389     .                "inst(X)",  zsto,zout)
    384390         CALL histdef(nid_tra, "pplay", "flux u mont","-",
    385      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     391     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    386392     .                "inst(X)", zsto,zout)
    387393         CALL histdef(nid_tra, "t", "flux u mont","-",
    388      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     394     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    389395     .                "inst(X)", zsto,zout)
    390396         CALL histdef(nid_tra, "mfu", "flux u mont","-",
    391      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     397     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    392398     .                "ave(X)", zsto,zout)
    393399         CALL histdef(nid_tra, "mfd", "flux u decen","-",
    394      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     400     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    395401     .                "ave(X)", zsto,zout)
    396402         CALL histdef(nid_tra, "en_u", "flux u mont","-",
    397      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     403     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    398404     .                "ave(X)", zsto,zout)
    399405         CALL histdef(nid_tra, "en_d", "flux u mont","-",
    400      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     406     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    401407     .                "ave(X)", zsto,zout)
    402408         CALL histdef(nid_tra, "de_d", "flux u mont","-",
    403      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     409     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    404410     .                "ave(X)", zsto,zout)
    405411         CALL histdef(nid_tra, "de_u", "flux u decen","-",
    406      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     412     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    407413     .                "ave(X)", zsto,zout)
    408414         CALL histdef(nid_tra, "coefh", "turbulent coef","-",
    409      .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     415     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    410416     .                "ave(X)", zsto,zout)
    411417
Note: See TracChangeset for help on using the changeset viewer.