Ignore:
Timestamp:
Mar 24, 2004, 9:59:48 AM (20 years ago)
Author:
lmdzadmin
Message:

Ajout/modif sorties 3D niveaux NMC
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmthNMC.h

    r500 r508  
    44       IF(region.EQ.'3d') THEN
    55         zsto = dtime
     6cIM 150304 BEG
    67         zout = dtime * ecrit_mth
     8c        zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
     9cIM 150304 END
     10cIM 120304 BEG
     11c       zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))
     12c       zsto = zout
     13c       zsto = dtime
     14cIM 120304 END
    715       ELSEIF(region.EQ.'1d') THEN
    816         zsto = dtime
     
    4048
    4149c Champs interpolles sur des niveaux de pression du NMC
    42          CALL histdef(nid_nmc, "u1000",
    43      . "Zonal wind 1000mb",
    44      .                "m/s",
    45      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    46      .                "ave(X)", zsto,zout)
    47 
    48          CALL histdef(nid_nmc, "v1000",
    49      .                "Meridional wind 1000mb", "m/s",
    50      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    51      .                "ave(X)", zsto,zout)
    52 
    53          CALL histdef(nid_nmc, "u925", "Zonal wind 925mb",
    54      .                "m/s",
    55      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    56      .                "ave(X)", zsto,zout)
    57 
    58          CALL histdef(nid_nmc, "v925", "Meridional wind 925mb",
    59      .                "m/s",
    60      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    61      .                "ave(X)", zsto,zout)
    62 
    63          CALL histdef(nid_nmc, "u850", "Zonal wind 850mb",
    64      .                "m/s",
    65      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    66      .                "ave(X)", zsto,zout)
    67 
    68          CALL histdef(nid_nmc, "v850", "Meridional wind 850mb",
    69      .                "m/s",
     50c IMIMIM 110304 BEG
     51c
     52         DO k=1, nlevSTD
     53c
     54         bb=clevSTD(k)
     55c
     56         IF(k.GE.2) THEN
     57          aa=clevSTD(k)
     58          bb=aa(1:lnblnk1(aa))
     59         ENDIF
     60c
     61          CALL histdef(nid_nmc, "t"//bb,
     62     .                 "Temperature"//bb//"mb","K",
    7063     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    7164     .                "ave(X)", zsto,zout)
    7265c
    73          CALL histdef(nid_nmc, "u700", "Zonal wind 700mb",
    74      .                "m/s",
    75      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    76      .                "ave(X)", zsto,zout)
    77 
    78          CALL histdef(nid_nmc, "v700", "Meridional wind 700mb",
    79      .                "m/s",
     66         CALL histdef(nid_nmc, "phi"//bb,
     67     .                "Geopotential"//bb//"mb", "m2/s2",
    8068     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    8169     .                "ave(X)", zsto,zout)
    8270c
    83          CALL histdef(nid_nmc, "u600", "Zonal wind 600mb",
    84      .                "m/s",
    85      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    86      .                "ave(X)", zsto,zout)
    87 
    88          CALL histdef(nid_nmc, "v600", "Meridional wind 600mb",
    89      .                "m/s",
     71          CALL histdef(nid_nmc, "q"//bb,
     72     .                 "Specific humidity"//bb//"mb","kg/kg",
    9073     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    9174     .                "ave(X)", zsto,zout)
    9275c
    93          CALL histdef(nid_nmc, "u500", "Zonal wind 500mb",
    94      .                "m/s",
    95      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    96      .                "ave(X)", zsto,zout)
    97 
    98          CALL histdef(nid_nmc, "v500", "Meridional wind 500mb",
    99      .                "m/s",
     76         CALL histdef(nid_nmc, "rh"//bb,
     77     .                 "Relative humidity"//bb//"mb", "%",
    10078     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    10179     .                "ave(X)", zsto,zout)
    10280c
    103          CALL histdef(nid_nmc, "u400", "Zonal wind 400mb",
    104      .                "m/s",
    105      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    106      .                "ave(X)", zsto,zout)
    107 
    108          CALL histdef(nid_nmc, "v400", "Meridional wind 400mb",
    109      .                "m/s",
     81          CALL histdef(nid_nmc, "u"//bb,
     82     .                 "Zonal wind"//bb//"mb","K",
    11083     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    11184     .                "ave(X)", zsto,zout)
    11285c
    113          CALL histdef(nid_nmc, "u300", "Zonal wind 300mb",
    114      .                "m/s",
    115      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    116      .                "ave(X)", zsto,zout)
    117 
    118          CALL histdef(nid_nmc, "v300", "Meridional wind 300mb",
    119      .                "m/s",
     86          CALL histdef(nid_nmc, "v"//bb,
     87     .                 "Meridional wind"//bb//"mb","K",
    12088     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    12189     .                "ave(X)", zsto,zout)
    12290c
    123          CALL histdef(nid_nmc, "u250", "Zonal wind 250mb",
    124      .                "m/s",
    125      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    126      .                "ave(X)", zsto,zout)
     91         ENDDO !nlevSTD
     92c IMIMIM 110304 END
    12793
    128          CALL histdef(nid_nmc, "v250", "Meridional wind 250mb",
    129      .                "m/s",
    130      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    131      .                "ave(X)", zsto,zout)
    132 c
    133          CALL histdef(nid_nmc, "u200", "Zonal wind 200mb",
    134      .                "m/s",
    135      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    136      .                "ave(X)", zsto,zout)
    137 
    138          CALL histdef(nid_nmc, "v200", "Meridional wind 200mb",
    139      .                "m/s",
    140      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    141      .                "ave(X)", zsto,zout)
    142 c
    143          CALL histdef(nid_nmc, "u150", "Zonal wind 150mb",
    144      .                "m/s",
    145      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    146      .                "ave(X)", zsto,zout)
    147 
    148          CALL histdef(nid_nmc, "v150", "Meridional wind 150mb",
    149      .                "m/s",
    150      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    151      .                "ave(X)", zsto,zout)
    152 c
    153          CALL histdef(nid_nmc, "u100", "Zonal wind 100mb",
    154      .                "m/s",
    155      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    156      .                "ave(X)", zsto,zout)
    157 
    158          CALL histdef(nid_nmc, "v100", "Meridional wind 100mb",
    159      .                "m/s",
    160      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    161      .                "ave(X)", zsto,zout)
    162 c
    163          CALL histdef(nid_nmc, "u70", "Zonal wind 70mb",
    164      .                "m/s",
    165      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    166      .                "ave(X)", zsto,zout)
    167 
    168          CALL histdef(nid_nmc, "v70", "Meridional wind 70mb",
    169      .                "m/s",
    170      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    171      .                "ave(X)", zsto,zout)
    172 c
    173          CALL histdef(nid_nmc, "u50", "Zonal wind 50mb",
    174      .                "m/s",
    175      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    176      .                "ave(X)", zsto,zout)
    177 
    178          CALL histdef(nid_nmc, "v50", "Meridional wind 50mb",
    179      .                "m/s",
    180      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    181      .                "ave(X)", zsto,zout)
    182 c
    183          CALL histdef(nid_nmc, "u30", "Zonal wind 30mb",
    184      .                "m/s",
    185      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    186      .                "ave(X)", zsto,zout)
    187 
    188          CALL histdef(nid_nmc, "v30", "Meridional wind 30mb",
    189      .                "m/s",
    190      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    191      .                "ave(X)", zsto,zout)
    192 c
    193          CALL histdef(nid_nmc, "u20", "Zonal wind 20mb",
    194      .                "m/s",
    195      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    196      .                "ave(X)", zsto,zout)
    197 
    198          CALL histdef(nid_nmc, "v20", "Meridional wind 20mb",
    199      .                "m/s",
    200      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    201      .                "ave(X)", zsto,zout)
    202 c
    203          CALL histdef(nid_nmc, "u10", "Zonal wind 10mb",
    204      .                "m/s",
    205      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    206      .                "ave(X)", zsto,zout)
    207 
    208          CALL histdef(nid_nmc, "v10", "Meridional wind 10mb",
    209      .                "m/s",
    210      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    211      .                "ave(X)", zsto,zout)
    212 c
    21394         CALL histend(nid_nmc)
    21495c
Note: See TracChangeset for help on using the changeset viewer.