Changeset 508 for LMDZ.3.3/branches


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

Ajout/modif sorties 3D niveaux NMC
IM

Location:
LMDZ.3.3/branches/rel-LF/libf/phylmd
Files:
2 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
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histmthNMC.h

    r500 r508  
    1515c
    1616c  Champs interpolles sur des niveaux de pression du NMC
     17c110304 BEG
     18      DO k=1, nlevSTD
    1719c
    18       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u1000,zx_tmp_2d)
    19       CALL histwrite(nid_nmc,
    20      ."u1000",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     20         bb=clevSTD(k)
    2121c
    22       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v1000,zx_tmp_2d)
    23       CALL histwrite(nid_nmc,
    24      ."v1000",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     22         IF(k.GE.2) THEN
     23          aa=clevSTD(k)
     24          bb=aa(1:lnblnk1(aa))
     25         ENDIF
    2526c
    26       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u925,zx_tmp_2d)
    27       CALL histwrite(nid_nmc,
    28      ."u925",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     27      CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
     28      CALL histwrite(nid_nmc,"t"//bb,itau_w,zx_tmp_2d,
     29     $       iim*jjmp1,ndex2d)
    2930c
    30       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v925,zx_tmp_2d)
    31       CALL histwrite(nid_nmc,
    32      ."v925",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     31      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
     32      CALL histwrite(nid_nmc,"phi"//bb,itau_w,zx_tmp_2d,
     33     $       iim*jjmp1,ndex2d)
    3334c
    34       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
    35       CALL histwrite(nid_nmc,
    36      ."u850",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     35      CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
     36      CALL histwrite(nid_nmc,"q"//bb,itau_w,zx_tmp_2d,
     37     $       iim*jjmp1,ndex2d)
    3738c
    38       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
    39       CALL histwrite(nid_nmc,"v850",itau_w,
    40      .zx_tmp_2d,iim*jjmp1,ndex2d)
     39      CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlevSTD(:,k),zx_tmp_2d)
     40      CALL histwrite(nid_nmc,"rh"//bb,itau_w,zx_tmp_2d,
     41     $       iim*jjmp1,ndex2d)
    4142c
    42       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u700,zx_tmp_2d)
    43       CALL histwrite(nid_nmc,"u700",itau_w,
    44      .zx_tmp_2d,iim*jjmp1,ndex2d)
     43      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
     44      CALL histwrite(nid_nmc,"u"//bb,itau_w,zx_tmp_2d,
     45     $       iim*jjmp1,ndex2d)
    4546c
    46       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v700,zx_tmp_2d)
    47       CALL histwrite(nid_nmc,"v700",itau_w,
    48      .zx_tmp_2d,iim*jjmp1,ndex2d)
     47      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
     48      CALL histwrite(nid_nmc,"v"//bb,itau_w,zx_tmp_2d,
     49     $       iim*jjmp1,ndex2d)
    4950c
    50       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u600,zx_tmp_2d)
    51       CALL histwrite(nid_nmc,"u600",itau_w,
    52      .zx_tmp_2d,iim*jjmp1,ndex2d)
     51      ENDDO !k=1, nlevSTD
     52ccc
     53c110304 END
    5354c
    54       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v600,zx_tmp_2d)
    55       CALL histwrite(nid_nmc,"v600",itau_w,
    56      .zx_tmp_2d,iim*jjmp1,ndex2d)
    57 c
    58       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
    59       CALL histwrite(nid_nmc,"u500",itau_w,
    60      .zx_tmp_2d,iim*jjmp1,ndex2d)
    61 c
    62       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
    63       CALL histwrite(nid_nmc,"v500",itau_w,
    64      .zx_tmp_2d,iim*jjmp1,ndex2d)
    65 c
    66       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u400,zx_tmp_2d)
    67       CALL histwrite(nid_nmc,"u400",itau_w,
    68      .zx_tmp_2d,iim*jjmp1,ndex2d)
    69 c
    70       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v400,zx_tmp_2d)
    71       CALL histwrite(nid_nmc,"v400",itau_w,
    72      .zx_tmp_2d,iim*jjmp1,ndex2d)
    73 c
    74       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u300,zx_tmp_2d)
    75       CALL histwrite(nid_nmc,"u300",itau_w,
    76      .zx_tmp_2d,iim*jjmp1,ndex2d)
    77 c
    78       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v300,zx_tmp_2d)
    79       CALL histwrite(nid_nmc,"v300",itau_w,
    80      .zx_tmp_2d,iim*jjmp1,ndex2d)
    81 c
    82       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u250,zx_tmp_2d)
    83       CALL histwrite(nid_nmc,"u250",itau_w,
    84      .zx_tmp_2d,iim*jjmp1,ndex2d)
    85 c
    86       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v250,zx_tmp_2d)
    87       CALL histwrite(nid_nmc,"v250",itau_w,
    88      .zx_tmp_2d,iim*jjmp1,ndex2d)
    89 c
    90       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
    91       CALL histwrite(nid_nmc,"u200",itau_w,
    92      .zx_tmp_2d,iim*jjmp1,ndex2d)
    93 c
    94       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
    95       CALL histwrite(nid_nmc,"v200",itau_w,
    96      .zx_tmp_2d,iim*jjmp1,ndex2d)
    97 c
    98       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u150,zx_tmp_2d)
    99       CALL histwrite(nid_nmc,"u150",itau_w,
    100      .zx_tmp_2d,iim*jjmp1,ndex2d)
    101 c
    102       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v150,zx_tmp_2d)
    103       CALL histwrite(nid_nmc,"v150",itau_w,
    104      .zx_tmp_2d,iim*jjmp1,ndex2d)
    105 c
    106       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u100,zx_tmp_2d)
    107       CALL histwrite(nid_nmc,"u100",itau_w,
    108      .zx_tmp_2d,iim*jjmp1,ndex2d)
    109 c
    110       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v100,zx_tmp_2d)
    111       CALL histwrite(nid_nmc,"v100",itau_w,
    112      .zx_tmp_2d,iim*jjmp1,ndex2d)
    113 c
    114       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u70,zx_tmp_2d)
    115       CALL histwrite(nid_nmc,"u70",itau_w,
    116      .zx_tmp_2d,iim*jjmp1,ndex2d)
    117 c
    118       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v70,zx_tmp_2d)
    119       CALL histwrite(nid_nmc,"v70",itau_w,
    120      .zx_tmp_2d,iim*jjmp1,ndex2d)
    121 c
    122       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u50,zx_tmp_2d)
    123       CALL histwrite(nid_nmc,"u50",itau_w,
    124      .zx_tmp_2d,iim*jjmp1,ndex2d)
    125 c
    126       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v50,zx_tmp_2d)
    127       CALL histwrite(nid_nmc,"v50",itau_w,
    128      .zx_tmp_2d,iim*jjmp1,ndex2d)
    129 c
    130       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u30,zx_tmp_2d)
    131       CALL histwrite(nid_nmc,"u30",itau_w,
    132      .zx_tmp_2d,iim*jjmp1,ndex2d)
    133 c
    134       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v30,zx_tmp_2d)
    135       CALL histwrite(nid_nmc,"v30",itau_w,
    136      .zx_tmp_2d,iim*jjmp1,ndex2d)
    137 c
    138       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u20,zx_tmp_2d)
    139       CALL histwrite(nid_nmc,"u20",itau_w,
    140      .zx_tmp_2d,iim*jjmp1,ndex2d)
    141 c
    142       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v20,zx_tmp_2d)
    143       CALL histwrite(nid_nmc,"v20",itau_w,
    144      .zx_tmp_2d,iim*jjmp1,ndex2d)
    145 c
    146       CALL gr_fi_ecrit(1, klon,iim,jjmp1, u10,zx_tmp_2d)
    147       CALL histwrite(nid_nmc,"u10",itau_w,
    148      .zx_tmp_2d,iim*jjmp1,ndex2d)
    149 c
    150       CALL gr_fi_ecrit(1, klon,iim,jjmp1, v10,zx_tmp_2d)
    151       CALL histwrite(nid_nmc,"v10",itau_w,
    152      .zx_tmp_2d,iim*jjmp1,ndex2d)
    153 
    15455      if (ok_sync) then
    15556        call histsync(nid_nmc)
Note: See TracChangeset for help on using the changeset viewer.