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_histdayNMC.h

    r1374 r1398  
    3030c
    3131        IF(lev_histdayNMC.EQ.nlevSTD) THEN
    32          CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
    33      .                 nlevSTD, rlevSTD, nvert)
     32         CALL histvert(nid_daynmc, "plev", "pressure", "mb",
     33     .                 nlevSTD, rlevSTD/100., nvert,"down")
    3434        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
    35          CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
    36      .                 nlevSTD8, rlevSTD8, nvert)
     35         CALL histvert(nid_daynmc, "plev", "pressure", "mb",
     36     .                 nlevSTD8, rlevSTD8/100., nvert,"down")
    3737        ENDIF
    38 ccc
     38c
     39cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
     40cIM             le time_counter et les bounds
     41c
    3942ccc Champs 3D interpolles sur des niveaux de pression du NMC
    4043ccc
    … …  
    5154     .                 "Valeurs non-definies","-",
    5255     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    53      .                "inst(X)", zout,zout)
     56     .                "ave(X)", zout,zout)
    5457c
    5558          CALL histdef(nid_daynmc, "ta",
    5659     .                 "Air temperature","K",
    5760     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    58      .                "inst(X)", zout,zout)
     61     .                "ave(X)", zout,zout)
    5962c
    6063         CALL histdef(nid_daynmc, "zg",
    6164     .                "Geopotential height", "m",
    6265     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    63      .                "inst(X)", zout,zout)
     66     .                "ave(X)", zout,zout)
    6467c
    6568          CALL histdef(nid_daynmc, "hus",
    6669     .                 "Specific humidity","1",
    6770     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    68      .                "inst(X)", zout,zout)
     71     .                "ave(X)", zout,zout)
    6972c
    7073         CALL histdef(nid_daynmc, "hur",
    7174     .                 "Relative humidity", "%",
    7275     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    73      .                "inst(X)", zout,zout)
     76     .                "ave(X)", zout,zout)
    7477c
    7578          CALL histdef(nid_daynmc, "ua",
    7679     .                 "Eastward wind","m s-1",
    7780     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    78      .                "inst(X)", zout,zout)
     81     .                "ave(X)", zout,zout)
    7982c
    8083          CALL histdef(nid_daynmc, "va",
    8184     .                 "Northward wind","m s-1",
    8285     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    83      .                "inst(X)", zout,zout)
     86     .                "ave(X)", zout,zout)
    8487c
    8588          CALL histdef(nid_daynmc, "wap",
    8689     .                 "Lagrangian tendency of air pressure","Pa s-1",
    8790     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
    88      .                "inst(X)", zout,zout)
     91     .                "ave(X)", zout,zout)
    8992c
    9093        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
    … …  
    9396     .                 "Valeurs non-definies","-",
    9497     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    95      .                "inst(X)", zout,zout)
     98     .                "ave(X)", zout,zout)
    9699c
    97100          CALL histdef(nid_daynmc, "ta",
    98101     .                 "Air temperature","K",
    99102     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    100      .                "inst(X)", zout,zout)
     103     .                "ave(X)", zout,zout)
    101104c
    102105         CALL histdef(nid_daynmc, "zg",
    103106     .                "Geopotential height", "m",
    104107     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    105      .                "inst(X)", zout,zout)
     108     .                "ave(X)", zout,zout)
    106109c
    107110          CALL histdef(nid_daynmc, "hus",
    108111     .                 "Specific humidity","1",
    109112     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    110      .                "inst(X)", zout,zout)
     113     .                "ave(X)", zout,zout)
    111114c
    112115         CALL histdef(nid_daynmc, "hur",
    113116     .                 "Relative humidity", "%",
    114117     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    115      .                "inst(X)", zout,zout)
     118     .                "ave(X)", zout,zout)
    116119c
    117120          CALL histdef(nid_daynmc, "ua",
    118121     .                 "Eastward wind","m s-1",
    119122     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    120      .                "inst(X)", zout,zout)
     123     .                "ave(X)", zout,zout)
    121124c
    122125          CALL histdef(nid_daynmc, "va",
    123126     .                 "Northward wind","m s-1",
    124127     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    125      .                "inst(X)", zout,zout)
     128     .                "ave(X)", zout,zout)
    126129c
    127130          CALL histdef(nid_daynmc, "wap",
    128131     .                 "Lagrangian tendency of air pressure","Pa s-1",
    129132     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
    130      .                "inst(X)", zout,zout)
     133     .                "ave(X)", zout,zout)
    131134c
    132135        ENDIF
Note: See TracChangeset for help on using the changeset viewer.