Ignore:
Timestamp:
Jun 4, 2010, 6:56:18 PM (14 years ago)
Author:
musat
Message:

Last corrections for CMIP5:

  • Add O3 at standard level files histmthNMC.nc
  • Add positive attribute "down" for vertical axes for all output files
  • Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
  • Correct units for vertical axes : mb instead of hPa
  • Add mass flux at the bottom of clouds
  • Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
  • Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
  • Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
  • Fields albs, albslw become alb1, alb2 in output files
  • Correct title for rugs_* fields
  • Correct units for pbase and ptop are Pa (not mb)
  • Correct ndayrain field

FH/JLD/JYG/MAF/IM

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r1374 r1398  
    2727     .                 nhori, nid_mthnmc)
    2828c
    29          CALL histvert(nid_mthnmc, "plev", "pressure", "Pa",
    30      .                 nlevSTD, rlevSTD, nvert)
    31 ccc
     29         CALL histvert(nid_mthnmc, "plev", "pressure", "mb",
     30     .                 nlevSTD, rlevSTD/100., nvert,"down")
     31c
     32cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
     33cIM             le time_counter et les bounds
     34cIM
    3235ccc Champs 3D interpolles sur des niveaux de pression du NMC
    3336ccc
    … …  
    4447     .                 "Valeurs non-definies","-",
    4548     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    46      .                "inst(X)", zout,zout)
     49     .                "ave(X)", zout,zout)
    4750c
    4851          CALL histdef(nid_mthnmc, "ta",
    4952     .                 "Air temperature","K",
    5053     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    51      .                "inst(X)", zout,zout)
     54     .                "ave(X)", zout,zout)
    5255c
    5356         CALL histdef(nid_mthnmc, "zg",
    5457     .                "Geopotential height", "m",
    5558     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    56      .                "inst(X)", zout,zout)
     59     .                "ave(X)", zout,zout)
    5760c
    5861          CALL histdef(nid_mthnmc, "hus",
    5962     .                 "Specific humidity","1",
    6063     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    61      .                "inst(X)", zout,zout)
     64     .                "ave(X)", zout,zout)
    6265c
    6366         CALL histdef(nid_mthnmc, "hur",
    6467     .                 "Relative humidity", "%",
    6568     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    66      .                "inst(X)", zout,zout)
     69     .                "ave(X)", zout,zout)
    6770c
    6871          CALL histdef(nid_mthnmc, "ua",
    6972     .                 "Eastward wind","m s-1",
    7073     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    71      .                "inst(X)", zout,zout)
     74     .                "ave(X)", zout,zout)
    7275c
    7376          CALL histdef(nid_mthnmc, "va",
    7477     .                 "Northward wind","m s-1",
    7578     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    76      .                "inst(X)", zout,zout)
     79     .                "ave(X)", zout,zout)
    7780c
    7881          CALL histdef(nid_mthnmc, "wap",
    7982     .                 "Lagrangian tendency of air pressure","Pa s-1",
    8083     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    81      .                "inst(X)", zout,zout)
     84     .                "ave(X)", zout,zout)
    8285c
    8386          CALL histdef(nid_mthnmc, "psbg",
    8487     .         "Pressure sfce below ground","%",
    8588     .         iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    86      .         "inst(X)", zout,zout)
     89     .         "ave(X)", zout,zout)
    8790c
    8891          CALL histdef(nid_mthnmc, "uv",
    8992     .         "uv ",
    9093     .         "m2/s2",iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    91      .         "inst(X)", zout,zout)
     94     .         "ave(X)", zout,zout)
    9295c
    9396          CALL histdef(nid_mthnmc, "vq",
    … …  
    9598     .         "m/s * (kg/kg)",iim,jj_nb,nhori,
    9699     .          nlevSTD,1,nlevSTD, nvert, 32,
    97      .         "inst(X)", zout,zout)
     100     .         "ave(X)", zout,zout)
    98101c
    99102          CALL histdef(nid_mthnmc, "vT",
    … …  
    101104     .         "mK/s",iim,jj_nb,nhori,
    102105     .          nlevSTD,1,nlevSTD, nvert, 32,
    103      .         "inst(X)", zout,zout)
     106     .         "ave(X)", zout,zout)
    104107c
    105108          CALL histdef(nid_mthnmc, "wq",
    … …  
    107110     .         "(Pa/s)*(kg/kg)",iim,jj_nb,nhori,
    108111     .          nlevSTD,1,nlevSTD, nvert, 32,
    109      .         "inst(X)", zout,zout)
     112     .         "ave(X)", zout,zout)
    110113c
    111114          CALL histdef(nid_mthnmc, "vphi",
    … …  
    113116     .         "m2/s",iim,jj_nb,nhori,
    114117     .          nlevSTD,1,nlevSTD, nvert, 32,
    115      .         "inst(X)", zout,zout)
     118     .         "ave(X)", zout,zout)
    116119c
    117120          CALL histdef(nid_mthnmc, "wT",
    … …  
    119122     .         "K*Pa/s",iim,jj_nb,nhori,
    120123     .          nlevSTD,1,nlevSTD, nvert, 32,
    121      .         "inst(X)", zout,zout)
     124     .         "ave(X)", zout,zout)
    122125c
    123126          CALL histdef(nid_mthnmc, "uxu",
    … …  
    125128     .         "m2/s2",iim,jj_nb,nhori,
    126129     .          nlevSTD,1,nlevSTD, nvert, 32,
    127      .         "inst(X)", zout,zout)
     130     .         "ave(X)", zout,zout)
    128131c
    129132          CALL histdef(nid_mthnmc, "vxv",
    … …  
    131134     .         "m2/s2",iim,jj_nb,nhori,
    132135     .          nlevSTD,1,nlevSTD, nvert, 32,
    133      .         "inst(X)", zout,zout)
     136     .         "ave(X)", zout,zout)
    134137c
    135138          CALL histdef(nid_mthnmc, "TxT",
    … …  
    137140     .         "K2",iim,jj_nb,nhori,
    138141     .          nlevSTD,1,nlevSTD, nvert, 32,
    139      .         "inst(X)", zout,zout)
     142     .         "ave(X)", zout,zout)
    140143c
    141 
     144          CALL histdef(nid_mthnmc, "tro3",
     145     .         "Ozone mole fraction",
     146     .         "1e-9",iim,jj_nb,nhori,
     147     .          nlevSTD,1,nlevSTD, nvert, 32,
     148     .         "ave(X)", zout,zout)
     149c
     150          if (read_climoz == 2) THEN
     151           CALL histdef(nid_mthnmc, "tro3_daylight",
     152     .         "Daylight ozone mole fraction",
     153     .         "1e-9",iim,jj_nb,nhori,
     154     .          nlevSTD,1,nlevSTD, nvert, 32,
     155     .         "ave(X)", zout,zout)
     156          endif
     157c
    142158         CALL histend(nid_mthnmc)
    143159c
Note: See TracChangeset for help on using the changeset viewer.