Ignore:
Timestamp:
Aug 17, 2006, 5:41:51 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/V3_test/libf/phylmd/ini_histday.h

    r684 r704  
    1111         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
    1212c
    13          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    14          DO i = 1, iim
    15             zx_lon(i,1) = rlon(i+1)
    16             zx_lon(i,jjmp1) = rlon(i+1)
    17          ENDDO
     13cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
     14cym         DO i = 1, iim
     15cym            zx_lon(i,1) = rlon(i+1)
     16cym            zx_lon(i,jjmp1) = rlon(i+1)
     17cym         ENDDO
    1818         DO ll=1,klev
    1919            znivsig(ll)=float(ll)
    2020         ENDDO
    21          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
    22          CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
    23      .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
    24      .                 nhori, nid_day)
     21cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
     22cym         CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
     23cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     24cym     .                 nhori, nid_day)
     25
     26         CALL histbeg_phy("histday",itau_phy, zjulian, dtime,
     27     .                    nhori, nid_day)
     28
     29         write(*,*)'Journee ', itau_phy, zjulian
    2530         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb",
    2631     .                 klev, presnivs/100., nvert)
     
    2934c
    3035         CALL histdef(nid_day, "phis", "Surface geop. height", "m2/s2",
    31      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     36     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    3237     .                "once", zstophy,zout)
    3338c
    3439         CALL histdef(nid_day, "aire", "Grid area", "-",
    35      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     40     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    3641     .                "once", zstophy,zout)
    3742c
    3843         CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","-",
    39      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     44     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4045     .                "once", zstophy,zout)
    4146c
    4247         CALL histdef(nid_day, "contfracOR","% sfce terre OR", "-",
    43      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     48     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4449     .                "once", zstophy,zout)
    4550c
     
    4752c
    4853         CALL histdef(nid_day, "tsol", "Surface Temperature", "K",
    49      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     54     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    5055     .                "ave(X)", zstophy,zout)
    5156c
    5257         CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
    53      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    54      .                "ave(X)", zstophy,zout)
    55 c
    56 c Champs retires momentannéement en attendant un hypothetique
     58     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
     59     .                "ave(X)", zstophy,zout)
     60c
     61c Champs retires momentannment en attendant un hypothetique
    5762c debugage
    5863
    5964         CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
    6065     .                "K",
    61      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     66     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    6267     .                t2mincels, zstophy,zout)
    6368c
    6469         CALL histdef(nid_day, "t2m_max", "Temp. 2m max.",
    6570     .                "K",
    66      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     71     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    6772     .                t2maxcels, zstophy,zout)
    6873c
    6974         CALL histdef(nid_day, "plul", "Large-scale Precip.",
    7075     .   "kg/(s*m2)",
    71      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     76     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7277     .                "ave(X)", zstophy,zout)
    7378c
    7479         CALL histdef(nid_day, "pluc", "Convective Precip.",
    7580     .   "kg/(s*m2)",
    76      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     81     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7782     .                "ave(X)", zstophy,zout)
    7883c
    7984         CALL histdef(nid_day, "snowl", "Solid Large-scale Precip.",
    8085     .   "kg/(m2*s)",
    81      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     86     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    8287     .                "ave(X)", zstophy,zout)
    8388c
     
    8590         CALL histdef(nid_day, "snowc", "Solid Convective Precip.",
    8691     .   "kg/(m2*s)",
    87      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     92     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    8893     .                "ave(X)", zstophy,zout)
    8994         ENDIF !snowc=0.
    9095c
    9196         CALL histdef(nid_day, "flat", "Latent heat flux", "W/m2",
    92      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     97     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    9398     .                "ave(X)", zstophy,zout)
    9499c
    95100         CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
    96      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     101     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    97102     .                "inst(X)", zstoday,zout)
    98103c
    99104         CALL histdef(nid_day, "q2m", "Specific humidity", "kg/kg",
    100      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     105     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    101106     .                "ave(X)", zstophy,zout)
    102107c
    103108         CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
    104      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     109     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    105110     .                "ave(X)", zstophy,zout)
    106111c
    107112         CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
    108      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     113     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    109114     .                "ave(X)", zstophy,zout)
    110115c
    111116         CALL histdef(nid_day, "wind10m","10-m wind speed","m/s",
    112      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     117     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    113118     .                "ave(X)", zstophy,zout)
    114119c
    115120         CALL histdef(nid_day, "wind10max", "10-m wind speed max.",
    116121     .                "m/s",
    117      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     122     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    118123     .                t2maxcels, zstophy,zout)
    119124c
    120125         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
    121      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     126     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    122127     .                "ave(X)", zstophy,zout)
    123128c
    124129         CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
    125130     .                , "kg/(s*m2)",
    126      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     131     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    127132     .                "ave(X)", zstophy,zout)
    128133c
    129134         CALL histdef(nid_day, "snowf", "Snow fall", "kg/(m2*s)",
    130      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     135     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    131136     .                "ave(X)", zstophy,zout)
    132137c
    133138c        CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
    134 c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     139c    .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    135140c    .                "ave(X)", zstophy,zout)
    136141c
    137142         CALL histdef(nid_day, "evap", "Evaporation", "kg/(m2*s)",
    138      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     143     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    139144     .                "ave(X)", zstophy,zout)
    140145c
    141146         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    142      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     147     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    143148     .                "ave(X)", zstophy,zout)
    144149c
    145150         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    146      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     151     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    147152     .                "ave(X)", zstophy,zout)
    148153c
    149154         CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
    150155     .                "W/m2",
    151      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     156     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    152157     .                "ave(X)", zstophy,zout)
    153158c
    154159         CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
    155      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     160     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    156161     .                "ave(X)", zstophy,zout)
    157162c
    158163         CALL histdef(nid_day, "radsol", "Rayonnement au sol", "W/m2",
    159      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     164     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    160165     .                "ave(X)", zstophy,zout)
    161166c
    162167         CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
    163      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     168     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    164169     .                "ave(X)", zstophy,zout)
    165170c
    166171         CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
    167      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     172     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    168173     .                "ave(X)", zstophy,zout)
    169174c
    170175         CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
    171      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     176     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    172177     .                "ave(X)", zstophy,zout)
    173178c
    174179         CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
    175      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     180     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    176181     .                "ave(X)", zstophy,zout)
    177182c
    178183         CALL histdef(nid_day, "cldq", "Cloud liquid water path",
    179184     .                "kg/m2",
    180      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     185     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    181186     .                "ave(X)", zstophy,zout)
    182187c
    183188         CALL histdef(nid_day, "prw", "Precipitable water", "kg/m2",
    184      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     189     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    185190     .                "ave(X)", zstophy,zout)
    186191c
     
    203208          CALL histdef(nid_day, "u"//bb2,
    204209     .                 "Zonal wind "//bb2//"mb","m/s",
    205      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     210     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    206211     .                "inst(X)", zout,zout)
    207212c
    208213          CALL histdef(nid_day, "v"//bb2,
    209214     .                 "Meridional wind "//bb2//"mb","m/s",
    210      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     215     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    211216     .                "inst(X)", zout,zout)
    212217c
     
    217222         CALL histdef(nid_day, "w"//bb2, "Vertical wind "//bb2//"mb",
    218223     .                "m/s",
    219      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     224     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    220225     .                "inst(X)", zout,zout)
    221226c
    222227         CALL histdef(nid_day,"phi"//bb2, "Geopotential "//bb2//"mb",
    223228     .                "m",
    224      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     229     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    225230     .                "inst(X)", zout,zout)
    226231c
     
    230235c
    231236         CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
    232      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     237     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    233238     .                "ave(X)", zstophy,zout)
    234239c
    235240         CALL histdef(nid_day, "cape_max", "CAPE max.",
    236241     .                "J/kg",
    237      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     242     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    238243     .                capemaxcels, zstophy,zout)
    239244c
    240245         CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
    241      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     246     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    242247     .                "ave(X)", zstophy,zout)
    243248c
    244249         CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
    245      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     250     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    246251     .                "ave(X)", zstophy,zout)
    247252c
    248253         CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
    249      .                iim,jjmp1,nhori, 1,1,1,-99,
     254     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    250255     .                32, "ave(X)", zstophy,zout)
    251256c
     
    257262          CALL histdef(nid_day, "lmt_bils",
    258263     .         "Bilan au sol atmosphere forcee", "W/m2",
    259      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     264     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    260265     .                "ave(X)", zstophy,zout)
    261266c
     
    264269          CALL histdef(nid_day, "slab_bils",
    265270     .                "Bilan au sol slab", "W/m2",
    266      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     271     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    267272     .                "ave(X)", zstophy,zout)
    268273c
    269274          CALL histdef(nid_day, "tslab", "Slab SST ", "K",
    270      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     275     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    271276     .                "ave(X)", zstophy,zout)
    272277c
    273278          CALL histdef(nid_day, "seaice", "Slab seaice", "kg/m2",
    274      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     279     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    275280     .                "ave(X)", zstophy,zout)
    276281c
    277282          CALL histdef(nid_day, "siceh", "Slab seaice height", "m",
    278      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     283     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    279284     .                "ave(X)", zstophy,zout)
    280285c
     
    283288         CALL histdef(nid_day, "fluxo",
    284289     .                "Flux turbulents ocean-atmosphere", "W/m2",
    285      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     290     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    286291     .                "ave(X)", zstophy,zout)
    287292c
    288293         CALL histdef(nid_day, "fluxg",
    289294     .                "Flux ocean-glace de mer", "W/m2",
    290      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     295     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    291296     .                "ave(X)", zstophy,zout)
    292297c
     
    296301c
    297302         CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
    298      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     303     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    299304     .                "ave(X)", zstophy,zout)
    300305c
    301306         CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2/K",
    302      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     307     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    303308     .                "ave(X)", zstophy,zout)
    304309c
     
    314319c
    315320         CALL histdef(nid_day, "temp", "Air temperature", "K",
    316      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     321     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    317322     .                "ave(X)", zstophy,zout)
    318323c
    319324         CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
    320      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     325     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    321326     .                "ave(X)", zstophy,zout)
    322327c
    323328         CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
    324      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     329     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    325330     .                "ave(X)", zstophy,zout)
    326331c
    327332         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
    328      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     333     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    329334     .                "ave(X)", zstophy,zout)
    330335c
    331336         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
    332      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     337     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    333338     .                "ave(X)", zstophy,zout)
    334339c
    335340         CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
    336      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     341     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    337342     .                "ave(X)", zstophy,zout)
    338343c
    339344         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
    340      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     345     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    341346     .                "ave(X)", zstophy,zout)
    342347c
     
    354359c
    355360         CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
    356      .                iim,jjmp1,nhori, 1,1,1,-99,
     361     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    357362     .                32, "ave(X)", zstophy,zout)
    358363c
    359364         CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
    360      .                iim,jjmp1,nhori, 1,1,1,-99,
     365     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    361366     .                32, "ave(X)", zstophy,zout)
    362367c
    363368         CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
    364      .                iim,jjmp1,nhori, 1,1,1,-99,
     369     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    365370     .                32, "ave(X)", zstophy,zout)
    366371c
    367372         CALL histdef(nid_day, "SWupTOAclr",
    368373     .                "SWup clear sky at TOA","W/m2",
    369      .                iim,jjmp1,nhori, 1,1,1,-99,
     374     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    370375     .                32, "ave(X)", zstophy,zout)
    371376
    372377         CALL histdef(nid_day, "SWupSFCclr",
    373378     .                "SWup clear sky at surface","W/m2",
    374      .                iim,jjmp1,nhori, 1,1,1,-99,
     379     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    375380     .                32, "ave(X)", zstophy,zout)
    376381
    377382         CALL histdef(nid_day, "SWdnTOAclr",
    378383     .                "SWdn clear sky at TOA","W/m2",
    379      .                iim,jjmp1,nhori, 1,1,1,-99,
     384     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    380385     .                32, "ave(X)", zstophy,zout)
    381386
    382387         CALL histdef(nid_day, "SWdnSFCclr",
    383388     .                "SWdn clear sky at surface","W/m2",
    384      .                iim,jjmp1,nhori, 1,1,1,-99,
     389     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    385390     .                32, "ave(X)", zstophy,zout)
    386391c
     
    388393         CALL histdef(nid_day, "LWdnSFC",
    389394     .                "LW down at surface","W/m2",
    390      .                iim,jjmp1,nhori, 1,1,1,-99,
     395     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    391396     .                32, "ave(X)", zstophy,zout)
    392397c
    393398         CALL histdef(nid_day, "LWupSFC",
    394399     .                "LW down at surface","W/m2",
    395      .                iim,jjmp1,nhori, 1,1,1,-99,
     400     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    396401     .                32, "ave(X)", zstophy,zout)
    397402c
     
    400405c    .                "Cld top pressure","mb",
    401406     .                "Cld top pressure (cf ISCCP simulator)","Pa",
    402      .                iim,jjmp1,nhori, 1,1,1,-99,
     407     .                iim,jjphy_nb,nhori, 1,1,1,-99,
    403408     .                32, "ave(X)", zstophy,zout)
    404409cIM     ENDIF !(iflag_con.GE.3) THEN
     
    411416c
    412417         CALL histdef(nid_day, "tter", "Surface Temperature ter", "K",
    413      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     418     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    414419     .                "ave(X)", zstophy,zout)
    415420c
    416421         CALL histdef(nid_day, "tlic", "Surface Temperature lic", "K",
    417      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     422     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    418423     .                "ave(X)", zstophy,zout)
    419424c
    420425         CALL histdef(nid_day, "toce", "Surface Temperature oce", "K",
    421      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     426     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    422427     .                "ave(X)", zstophy,zout)
    423428c
    424429         CALL histdef(nid_day, "tsic", "Surface Temperature sic", "K",
    425      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     430     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    426431     .                "ave(X)", zstophy,zout)
    427432c
    428433         CALL histdef(nid_day, "t2mter", "Temp.tere 2m", "K",
    429      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     434     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    430435     .                "ave(X)", zstophy,zout)
    431436c
    432437         CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
    433      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     438     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    434439     .                "ave(X)", zstophy,zout)
    435440c
    436441         CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
    437      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     442     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    438443     .                "ave(X)", zstophy,zout)
    439444c
    440445         CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
    441      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     446     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    442447     .                "ave(X)", zstophy,zout)
    443448c
    444449         CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
    445450     .                "K",
    446      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     451     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    447452     .                t2mincels, zstophy,zout)
    448453c
    449454         CALL histdef(nid_day, "t2mter_max", "Temp.terre 2m max.",
    450455     .                "K",
    451      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     456     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    452457     .                t2maxcels, zstophy,zout)
    453458c
    454459         CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
    455      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     460     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    456461     .                "ave(X)", zstophy,zout)
    457462c
    458463         CALL histdef(nid_day, "u10mlic", "Vent zonal lic 10m", "m/s",
    459      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     464     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    460465     .                "ave(X)", zstophy,zout)
    461466c
    462467         CALL histdef(nid_day, "u10moce", "Vent zonal oce 10m", "m/s",
    463      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     468     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    464469     .                "ave(X)", zstophy,zout)
    465470c
    466471         CALL histdef(nid_day, "u10msic", "Vent zonal sic 10m",
    467      .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
     472     .                 "m/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    468473     .                "ave(X)", zstophy,zout)
    469474c
    470475         CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
    471      .                "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
     476     .                "m/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    472477     .                "ave(X)", zstophy,zout)
    473478c
    474479         CALL histdef(nid_day, "v10mlic", "Vent meridien lic 10m",
    475      .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
     480     .                 "m/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    476481     .                "ave(X)", zstophy,zout)
    477482c
    478483         CALL histdef(nid_day, "v10moce", "Vent meridien oce 10m",
    479      .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
     484     .                 "m/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    480485     .                "ave(X)", zstophy,zout)
    481486c
    482487         CALL histdef(nid_day, "v10msic", "Vent meridien sic 10m",
    483      .                 "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32,
     488     .                 "m/s",iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    484489     .                "ave(X)", zstophy,zout)
    485490c
     
    488493           call histdef(nid_day, "pourc_"//clnsurf(nsrf),
    489494     $         "% "//clnsurf(nsrf), "%", 
    490      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     495     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    491496     $         "ave(X)", zstophy,zout)
    492497C
    493498           call histdef(nid_day, "fract_"//clnsurf(nsrf),
    494499     $         "Fraction "//clnsurf(nsrf), "1", 
    495      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     500     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    496501     $         "ave(X)", zstophy,zout)
    497502C
    498503           call histdef(nid_day, "tsol_"//clnsurf(nsrf),
    499504     $         "Temperature "//clnsurf(nsrf), "K", 
    500      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     505     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    501506     $         "ave(X)", zstophy,zout)
    502507C
    503508           call histdef(nid_day, "sens_"//clnsurf(nsrf),
    504509     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
    505      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     510     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    506511     $         "ave(X)", zstophy,zout)
    507512c
    508513           call histdef(nid_day, "lat_"//clnsurf(nsrf),
    509514     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
    510      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     515     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    511516     $         "ave(X)", zstophy,zout)
    512517C
    513518           call histdef(nid_day, "taux_"//clnsurf(nsrf),
    514519     $         "Zonal wind stress "//clnsurf(nsrf),"Pa",
    515      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     520     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    516521     $         "ave(X)", zstophy,zout)
    517522
    518523           call histdef(nid_day, "tauy_"//clnsurf(nsrf),
    519524     $         "Meridional wind stress "//clnsurf(nsrf), "Pa", 
    520      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     525     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    521526     $         "ave(X)", zstophy,zout)
    522527C
    523528           call histdef(nid_day, "albe_"//clnsurf(nsrf),
    524529     $         "Albedo surf. "//clnsurf(nsrf), "-", 
    525      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     530     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    526531     $         "ave(X)", zstophy,zout)
    527532C
    528533           call histdef(nid_day, "rugs_"//clnsurf(nsrf),
    529534     $         "Rugosite "//clnsurf(nsrf), "m", 
    530      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     535     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    531536     $         "ave(X)", zstophy,zout)
    532537C
     
    545550         CALL histdef(nid_day, "tops0", "CS Solar rad. at TOA",
    546551     .                "W/m2",
    547      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     552     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    548553     .                "ave(X)", zstophy,zout)
    549554c
    550555         CALL histdef(nid_day, "topl0", "CS IR rad. at TOA",
    551556     .                "W/m2",
    552      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     557     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    553558     .                "ave(X)", zstophy,zout)
    554559c
    555560         CALL histdef(nid_day, "sols0", "CS Net Solar rad. at surf.",
    556561     .                "W/m2",
    557      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     562     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    558563     .                "ave(X)", zstophy,zout)
    559564c
    560565         CALL histdef(nid_day, "soll0", "CS Net IR rad. at surface",
    561566     .                "W/m2",
    562      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     567     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    563568     .                "ave(X)", zstophy,zout)
    564569c
    565570         CALL histdef(nid_day, "rneb", "Cloud fraction", "-",
    566      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     571     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    567572     .                "ave(X)", zstophy,zout)
    568573c
    569574         CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction"
    570575     .                , "-",
    571      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     576     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    572577     .                "ave(X)", zstophy,zout)
    573578c
    574579         CALL histdef(nid_day,"lwcon","Cloud liquid water content",
    575580     .                "kg/kg",
    576      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     581     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    577582     .                "ave(X)", zstophy,zout)
    578583c
    579584        CALL histdef(nid_day,"iwcon","Cloud ice water content","kg/kg",
    580      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     585     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    581586     .                "ave(X)", zstophy,zout)
    582587c
    583588         CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
    584      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     589     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    585590     .                "ave(X)", zstophy,zout)
    586591c
    587592         CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
    588      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     593     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    589594     .                "ave(X)", zstophy,zout)
    590595c
    591596        CALL histdef(nid_day,"meantaucld",
    592597     .                "ISCCP mean cloud optical thickness","1",
    593      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     598     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    594599     .                "ave(X)", zstophy,zout)
    595600c
    596601        CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
    597      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     602     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    598603     .                "ave(X)", zstophy,zout)
    599604c
    600605        CALL histdef(nid_day,"cldemi","Cloud optical emissivity","1",
    601      .                iim,jjmp1,nhori, klev,1,klev, nvert, 32,
     606     .                iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    602607     .                "ave(X)", zstophy,zout)
    603608c
     
    605610c
    606611         CALL histdef(nid_day, "dtcon", "Convection dT", "K/s",
    607      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     612     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    608613     .                "ave(X)", zstophy,zout)
    609614c
    610615         CALL histdef(nid_day, "dqcon", "Convection dQ", "(kg/kg)/s",
    611      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     616     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    612617     .                "ave(X)", zstophy,zout)
    613618c
    614619         CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
    615      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     620     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    616621     .                "ave(X)", zstophy,zout)
    617622c
    618623         CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
    619      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     624     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    620625     .                "ave(X)", zstophy,zout)
    621626c
    622627         CALL histdef(nid_day, "tpot",
    623628     .                "Surface air potential temperature", "K",
    624      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     629     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    625630     .                "ave(X)", zstophy,zout)
    626631c
    627632         CALL histdef(nid_day, "tpote",
    628633     .        "Surface air equivalent potential temperature", "K",
    629      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     634     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    630635     .                "ave(X)", zstophy,zout)
    631636c
Note: See TracChangeset for help on using the changeset viewer.