Ignore:
Timestamp:
Dec 10, 2009, 10:02:56 AM (15 years ago)
Author:
Laurent Fairhead
Message:

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

Location:
LMDZ4/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk

  • LMDZ4/trunk/libf/phylmd/ini_histrac.h

    r1146 r1279  
    11!
    2 ! $Header$
     2! $Id $
    33!
    4       IF (ecrit_tra>0. .AND. config_inca == 'none') THEN
    5 c$OMP MASTER
    6          CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
    7 c
    8          CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,
    9      .                 nhori, nid_tra)
    10          CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
    11      .                 klev, presnivs, nvert)
     4  IF (ecrit_tra>0. .AND. config_inca == 'none') THEN
     5!$OMP MASTER
     6     CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
     7     CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,nhori, nid_tra)
     8     CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",klev, presnivs, nvert)
    129
     10     zsto = pdtphys
     11     zout = ecrit_tra
     12     CALL histdef(nid_tra, "phis", "Surface geop. height", "-",   &
     13          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
     14     CALL histdef(nid_tra, "aire", "Grid area", "-",              &
     15          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
    1316
     17!TRACEURS
     18!----------------
     19     DO it = 1,nbtr
     20        iiq = niadv(it+2)
    1421
    15          zsto = pdtphys
    16          zout = ecrit_tra
    17 c
    18          CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
    19      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    20      .                "once",  zsto,zout)
    21 c
    22          CALL histdef(nid_tra, "aire", "Grid area", "-",
    23      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    24      .                "once",  zsto,zout)
    25          DO it=1,nbtr
    26 C champ 2D
    27          iq=it+2
    28          iiq=niadv(iq)
    29          CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga",
    30      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    31      .                "ave(X)", zsto,zout)
    32          if (lessivage) THEN
    33          CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq),
    34      .              "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    35      .              "ave(X)", zsto,zout)
    36          endif
     22! CONCENTRATIONS
     23        CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga",    &
     24             iim,jj_nb,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
    3725
    38 c---Ajout Olivia
    39          CALL histdef(nid_tra, "d_tr_th_"//tname(iiq),
    40      .                "tendance thermique"// ttext(iiq), "?",
    41      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    42      .                "ave(X)", zsto,zout)
    43 c
    44          if(iflag_con.GE.2) then
    45          CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq),
    46      .                "tendance convection"// ttext(iiq), "?",
    47      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    48      .                "ave(X)", zsto,zout)
    49              endif !(iflag_con.GE.2) then
    50          CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq),
    51      .                "tendance couche limite"// ttext(iiq), "?",
    52      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    53      .                "ave(X)", zsto,zout)
    54 c---fin Olivia   
     26! TD LESSIVAGE
     27        IF (lessivage .AND. aerosol(it)) THEN
     28           CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), &
     29                "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
     30                "ave(X)", zsto,zout)
     31        END IF
    5532
    56          ENDDO
     33! TD THERMIQUES
     34        IF (iflag_thermals.gt.0) THEN
     35           CALL histdef(nid_tra, "d_tr_th_"//tname(iiq),      &
     36                "tendance thermique"// ttext(iiq), "?",       &
     37                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
     38                "ave(X)", zsto,zout)
     39        ENDIF
    5740
    58          CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
    59      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    60      .                "inst(X)",  zout,zout)
     41! TD CONVECTION
     42        IF (iflag_con.GE.2) THEN
     43           CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq),   &
     44                "tendance convection"// ttext(iiq), "?",   &
     45                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
     46                "ave(X)", zsto,zout)
     47        ENDIF
    6148
    62          CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
    63      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    64      .                "inst(X)",  zout,zout)
    65          CALL histdef(nid_tra, "psrf1", "nature sol", "-",
    66      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    67      .                "inst(X)",  zout,zout)
    68          CALL histdef(nid_tra, "psrf2", "nature sol", "-",
    69      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    70      .                "inst(X)",  zout,zout)
    71          CALL histdef(nid_tra, "psrf3", "nature sol", "-",
    72      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    73      .                "inst(X)",  zout,zout)
    74          CALL histdef(nid_tra, "psrf4", "nature sol", "-",
    75      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    76      .                "inst(X)",  zout,zout)
    77          CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
    78      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    79      .                "inst(X)",  zout,zout)
    80          CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
    81      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    82      .                "inst(X)",  zout,zout)
    83          CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
    84      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    85      .                "inst",  zout,zout)
    86          CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
    87      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    88      .                "inst(X)",  zout,zout)
    89          CALL histdef(nid_tra, "pplay", "flux u mont","-",
    90      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    91      .                "inst(X)", zout,zout)
    92          CALL histdef(nid_tra, "t", "flux u mont","-",
    93      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    94      .                "inst(X)", zout,zout)
    95          CALL histdef(nid_tra, "mfu", "flux u mont","-",
    96      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    97      .                "ave(X)", zsto,zout)
    98          CALL histdef(nid_tra, "mfd", "flux u decen","-",
    99      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    100      .                "ave(X)", zsto,zout)
    101          CALL histdef(nid_tra, "en_u", "flux u mont","-",
    102      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    103      .                "ave(X)", zsto,zout)
    104          CALL histdef(nid_tra, "en_d", "flux u mont","-",
    105      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    106      .                "ave(X)", zsto,zout)
    107          CALL histdef(nid_tra, "de_d", "flux u mont","-",
    108      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    109      .                "ave(X)", zsto,zout)
    110          CALL histdef(nid_tra, "de_u", "flux u decen","-",
    111      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    112      .                "ave(X)", zsto,zout)
    113          CALL histdef(nid_tra, "coefh", "turbulent coef","-",
    114      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    115      .                "ave(X)", zsto,zout)
     49! TD COUCHE-LIMITE
     50        CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq),      &
     51             "tendance couche limite"// ttext(iiq), "?",   &
     52             iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
     53             "ave(X)", zsto,zout)
     54     ENDDO
     55!---------------   
     56!
     57! VENT (niveau 1)
     58     CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",      &
     59          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     60          "inst(X)",  zout,zout)     
     61     CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",      &
     62          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     63          "inst(X)",  zout,zout)
    11664
    117 c
    118          CALL histend(nid_tra)
    119          ndex2d = 0
    120          ndex3d = 0
    121          ndex = 0
    122 c$OMP END MASTER
    123       END IF ! ecrit_tra>0. .AND. config_inca == 'none'
     65! TEMPERATURE DU SOL
     66     CALL histdef(nid_tra, "ftsol1", "temper sol", "-",    &
     67          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     68          "inst(X)",  zout,zout)
     69     CALL histdef(nid_tra, "ftsol2", "temper sol", "-",    &
     70          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     71          "inst(X)",  zout,zout)
     72     CALL histdef(nid_tra, "ftsol3", "temper sol", "-",    &
     73          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     74          "inst",  zout,zout)
     75     CALL histdef(nid_tra, "ftsol4", "temper sol", "-",    &
     76          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     77          "inst(X)",  zout,zout)
     78
     79! NATURE DU SOL
     80     CALL histdef(nid_tra, "psrf1", "nature sol", "-",     &
     81          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     82          "inst(X)",  zout,zout)
     83     CALL histdef(nid_tra, "psrf2", "nature sol", "-",     &
     84          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     85          "inst(X)",  zout,zout)
     86     CALL histdef(nid_tra, "psrf3", "nature sol", "-",     &
     87          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     88          "inst(X)",  zout,zout)
     89     CALL histdef(nid_tra, "psrf4", "nature sol", "-",     &
     90          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
     91          "inst(X)",  zout,zout)
     92! DIVERS
     93     CALL histdef(nid_tra, "pplay", "flux u mont","-",     &
     94          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     95          "inst(X)", zout,zout)
     96     CALL histdef(nid_tra, "t", "flux u mont","-",         &
     97          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     98          "inst(X)", zout,zout)
     99     CALL histdef(nid_tra, "mfu", "flux u mont","-",       &
     100          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     101          "ave(X)", zsto,zout)
     102     CALL histdef(nid_tra, "mfd", "flux u decen","-",      &
     103          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     104          "ave(X)", zsto,zout)
     105     CALL histdef(nid_tra, "en_u", "flux u mont","-",      &
     106          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     107          "ave(X)", zsto,zout)
     108     CALL histdef(nid_tra, "en_d", "flux u mont","-",      &
     109          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     110          "ave(X)", zsto,zout)
     111     CALL histdef(nid_tra, "de_d", "flux u mont","-",      &
     112          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     113          "ave(X)", zsto,zout)
     114     CALL histdef(nid_tra, "de_u", "flux u decen","-",     &
     115          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     116          "ave(X)", zsto,zout)
     117     CALL histdef(nid_tra, "coefh", "turbulent coef","-",  &
     118          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
     119          "ave(X)", zsto,zout)   
     120     
     121     CALL histend(nid_tra)
     122!$OMP END MASTER
     123  END IF ! ecrit_tra>0. .AND. config_inca == 'none'
     124 
Note: See TracChangeset for help on using the changeset viewer.