Changeset 897 for LMDZ4


Ignore:
Timestamp:
Feb 6, 2008, 6:48:54 PM (16 years ago)
Author:
Laurent Fairhead
Message:

Compatibilite des sorties avec ENSEMBLES SD
LF

Location:
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histday.h

    r892 r897  
    192192     .                "ave(X)", zstophy,zout)
    193193c
     194         CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
     195     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     196     .                "ave(X)", zstophy,zout)
     197c
     198         CALL histdef(nid_day, "cape_max", "CAPE max.",
     199     .                "J/kg",
     200     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     201     .                capemaxcels, zstophy,zout)
     202c
     203         CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
     204     .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
     205     .                "ave(X)", zstophy,zout)
     206c
     207         CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
     208     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     209     .                "ave(X)", zstophy,zout)
     210c
     211         CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
     212     .                iim,jj_nb,nhori, 1,1,1,-99,
     213     .                32, "ave(X)", zstophy,zout)
     214c
    194215c  Champs dynamiques sur niveaux de pression
    195216c
    … …  
    236257         ENDDO !nlevSTD
    237258c
    238          CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
    239      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    240      .                "ave(X)", zstophy,zout)
    241 c
    242          CALL histdef(nid_day, "cape_max", "CAPE max.",
    243      .                "J/kg",
    244      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    245      .                capemaxcels, zstophy,zout)
    246 c
    247          CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
    248      .                "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32,
    249      .                "ave(X)", zstophy,zout)
    250 c
    251          CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
    252      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    253      .                "ave(X)", zstophy,zout)
    254 c
    255          CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
    256      .                iim,jj_nb,nhori, 1,1,1,-99,
    257      .                32, "ave(X)", zstophy,zout)
     259
    258260c
    259261cIM diagnostiques flux ocean-atm ou ocean-glace de mer
    … …  
    310312     .                "ave(X)", zstophy,zout)
    311313c
    312       ENDIF !lev_histday.GE.2
    313 c=================================================================
    314       IF(lev_histday.GE.3) THEN
    315 c=================================================================
    316 c   INITIALISATION DES CHAMPS 3D
    317 c=================================================================
    318 c=================================================================
    319 c=================================================================
    320 c Champs 3D:
    321 c
    322          CALL histdef(nid_day, "temp", "Air temperature", "K",
    323      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    324      .                "ave(X)", zstophy,zout)
    325 c
    326          CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
    327      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    328      .                "ave(X)", zstophy,zout)
    329 c
    330          CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
    331      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    332      .                "ave(X)", zstophy,zout)
    333 c
    334          CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
    335      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    336      .                "ave(X)", zstophy,zout)
    337 c
    338          CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
    339      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    340      .                "ave(X)", zstophy,zout)
    341 c
    342          CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s",
    343      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    344      .                "ave(X)", zstophy,zout)
    345 c
    346          CALL histdef(nid_day, "pres", "Air pressure", "Pa",
    347      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    348      .                "ave(X)", zstophy,zout)
    349 c
    350 c=================================================================
    351 c   FIN INITIALISATION DES CHAMPS 3D
    352 c=================================================================
    353       ENDIF !  lev_histday.GE.3
    354 c
    355       IF (lev_histday.GE.4) THEN
    356 c=================================================================
    357 c
    358 c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
    359 c
    360 c=================================================================
     314         CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
     315     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     316     .                "ave(X)", zstophy,zout)
    361317c
    362318         CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
    … …  
    376332     .                iim,jj_nb,nhori, 1,1,1,-99,
    377333     .                32, "ave(X)", zstophy,zout)
    378 
     334c
    379335         CALL histdef(nid_day, "SWupSFCclr",
    380336     .                "SWup clear sky at surface","W/m2",
    381337     .                iim,jj_nb,nhori, 1,1,1,-99,
    382338     .                32, "ave(X)", zstophy,zout)
    383 
     339c
    384340         CALL histdef(nid_day, "SWdnTOAclr",
    385341     .                "SWdn clear sky at TOA","W/m2",
    386342     .                iim,jj_nb,nhori, 1,1,1,-99,
    387343     .                32, "ave(X)", zstophy,zout)
    388 
     344c
    389345         CALL histdef(nid_day, "SWdnSFCclr",
    390346     .                "SWdn clear sky at surface","W/m2",
    … …  
    399355c
    400356         CALL histdef(nid_day, "LWupSFC",
    401      .                "LW down at surface","W/m2",
    402      .                iim,jj_nb,nhori, 1,1,1,-99,
    403      .                32, "ave(X)", zstophy,zout)
     357     .                "LW up at surface","W/m2",
     358     .                iim,jj_nb,nhori, 1,1,1,-99,
     359     .                32, "ave(X)", zstophy,zout)
     360c
     361         CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
     362     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     363     .                "ave(X)", zstophy,zout)
     364c
     365         CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
     366     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     367     .                "ave(X)", zstophy,zout)
     368c
     369         CALL histdef(nid_day, "taux","Zonal wind stress","Pa",
     370     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     371     .                "ave(X)", zstophy,zout)
     372c
     373         CALL histdef(nid_day, "tauy","Meridional wind stress ", "Pa",
     374     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     375     .                "ave(X)", zstophy,zout)
    404376c
    405377cIM     IF (iflag_con.GE.3) THEN
    … …  
    413385cIM     ENDIF !(iflag_con.GE.3) THEN
    414386c
    415 c=================================================================
    416 c=================================================================
     387      ENDIF !lev_histday.GE.2
     388c=================================================================
     389      IF(lev_histday.GE.3) THEN
     390c=================================================================
     391c   INITIALISATION DES CHAMPS 3D
     392c=================================================================
     393c=================================================================
     394c=================================================================
     395c Champs 3D:
     396c
     397         CALL histdef(nid_day, "temp", "Air temperature", "K",
     398     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     399     .                "ave(X)", zstophy,zout)
     400c
     401         CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg",
     402     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     403     .                "ave(X)", zstophy,zout)
     404c
     405         CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2",
     406     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     407     .                "ave(X)", zstophy,zout)
     408c
     409         CALL histdef(nid_day, "vitu", "Zonal wind", "m/s",
     410     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     411     .                "ave(X)", zstophy,zout)
     412c
     413         CALL histdef(nid_day, "vitv", "Meridional wind", "m/s",
     414     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     415     .                "ave(X)", zstophy,zout)
     416c
     417         CALL histdef(nid_day, "vitw", "Vertical wind", "m/s",
     418     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     419     .                "ave(X)", zstophy,zout)
     420c
     421         CALL histdef(nid_day, "pres", "Air pressure", "Pa",
     422     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     423     .                "ave(X)", zstophy,zout)
     424c
     425c=================================================================
     426c   FIN INITIALISATION DES CHAMPS 3D
     427c=================================================================
     428      ENDIF !  lev_histday.GE.3
     429c
     430      IF (lev_histday.GE.4) THEN
    417431c=================================================================
    418432c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
    … …  
    590604     .                "ave(X)", zstophy,zout)
    591605c
    592          CALL histdef(nid_day,"lwp","Cloud water path","kg/m2",
    593      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    594      .                "ave(X)", zstophy,zout)
    595 c
    596          CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2",
    597      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    598      .                "ave(X)", zstophy,zout)
    599 c
     606        IF (ok_isccp) THEN
    600607        CALL histdef(nid_day,"meantaucld",
    601608     .                "ISCCP mean cloud optical thickness","1",
    602609     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    603610     .                "ave(X)", zstophy,zout)
     611        ENDIF
    604612c
    605613        CALL histdef(nid_day,"cldtau","Cloud optical thickness","1",
    … …  
    621629     .                "ave(X)", zstophy,zout)
    622630c
    623          CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",
    624      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    625      .                "ave(X)", zstophy,zout)
    626 c
    627631         CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%",
    628632     .                iim,jj_nb,nhori, 1,1,1, -99, 32,

  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histhf.h

    r776 r897  
    9292         DO k=1, nlevSTD
    9393c
    94 c        bb=clevSTD(k)
    95 c
    96          IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
    97          IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
     94         bb2=clevSTD(k)
    9895c
    9996c        IF(k.GE.2) THEN
    … …  
    10299c        ENDIF
    103100c
    104           IF(bb2.EQ."500") THEN
    105           CALL histdef(nid_hf, "phi"//bb2,
     101         bb2=clevSTD(k)
     102         IF(bb2.EQ."500") THEN
     103         CALL histdef(nid_hf, "phi"//bb2,
    106104     .                "Geopotential "//bb2//"mb", "m",
    107105     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    108106     .                "inst(X)", zstohf,zout)
    109           ENDIF !bb2.EQ."500"
     107         ENDIF !bb2.EQ."500"
    110108         ENDDO !k=1, nlevSTD
    111109c
    … …  
    139137c
    140138         DO k=1, nlevSTD
     139c
     140c        bb=clevSTD(k)
     141c
     142c        IF(k.GE.2) THEN
     143c         IF(k.GE.1.AND.k.LE.10) bb2=clevSTD(k)
     144c         IF(k.GE.11.AND.k.LE.17) bb3=clevSTD(k)
     145c         aa=clevSTD(k)
     146c         bb=aa(1:lnblnk1(aa))
     147c        ENDIF
     148c
     149         bb2=clevSTD(k)
     150         IF(bb2.EQ."1000".OR.bb2.EQ."925".OR.
     151     $      bb2.EQ."850".OR.bb2.EQ."700".OR.
     152     $      bb2.EQ."500".OR.bb2.EQ."300".OR.
     153     $      bb2.EQ."250".OR.bb2.EQ."200".OR.
     154     $      bb2.EQ."150".OR.bb2.EQ."100".OR.
     155     $      bb2.EQ."50".OR.bb2.EQ."30".OR.
     156     $      bb2.EQ."10")THEN
     157c
     158          CALL histdef(nid_hf, "t"//bb2,
     159     .                 "Temperature "//bb2//"mb","K",
     160     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     161     .                "inst(X)", zstohf,zout)
     162c
     163          IF(bb2.NE."500") THEN
     164          CALL histdef(nid_hf, "phi"//bb2,
     165     .                "Geopotential "//bb2//"mb", "m",
     166     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     167     .                "inst(X)", zstohf,zout)
     168          ENDIF
     169c
     170          IF(bb2.GE."200") THEN
     171          CALL histdef(nid_hf, "q"//bb2,
     172     .                 "Specific humidity "//bb2//"mb","kg/kg",
     173     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     174     .                "inst(X)", zstohf,zout)
     175          ENDIF
     176c
     177          CALL histdef(nid_hf, "u"//bb2,
     178     .                 "Zonal wind "//bb2//"mb","m/s",
     179     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     180     .                "inst(X)", zstohf,zout)
     181c
     182          CALL histdef(nid_hf, "v"//bb2,
     183     .                 "Meridional wind "//bb2//"mb","m/s",
     184     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     185     .                "inst(X)", zstohf,zout)
     186c
     187          CALL histdef(nid_hf, "w"//bb2,
     188     .                 "Vertical motion "//bb2//"mb","Pa/s",
     189     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
     190     .                "inst(X)", zstohf,zstohf)
     191         ENDIF
     192c
     193         ENDDO !nlevSTD
     194c
     195      ENDIF !lev_histhf.GE.3
     196c
     197      IF(lev_histhf.GE.4) THEN
     198c
     199cIM sorties hf 3D
     200c
     201#undef histhf3d
     202#define histhf3d
     203#ifdef histhf3d
     204#include "ini_histhf3d.h"
     205#endif
     206      ENDIF !lev_histhf.GE.4
     207c
     208      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
     209c
     210      IF(lev_histhf.GE.1) THEN
     211c
     212        DO k=1, nlevSTD
     213c
    141214c
    142215c        bb=clevSTD(k)
    … …  
    149222c        ENDIF
    150223c
    151          IF(bb2.EQ."850".or.bb2.EQ."700".OR.
    152      $      bb2.EQ."500".or.bb2.EQ."200") THEN
    153 c
    154           CALL histdef(nid_hf, "t"//bb2,
    155      .                 "Temperature "//bb2//"mb","K",
    156      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    157      .                "inst(X)", zstohf,zout)
    158 c
    159           IF(bb2.NE."500") THEN
    160           CALL histdef(nid_hf, "phi"//bb2,
    161      .                "Geopotential "//bb2//"mb", "m",
    162      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    163      .                "inst(X)", zstohf,zout)
    164           ENDIF
    165 c
    166           CALL histdef(nid_hf, "q"//bb2,
    167      .                 "Specific humidity "//bb2//"mb","kg/kg",
    168      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    169      .                "inst(X)", zstohf,zout)
     224         IF(bb2.EQ."850".or.bb2.EQ."200") THEN
    170225c
    171226          CALL histdef(nid_hf, "u"//bb2,
    172227     .                 "Zonal wind "//bb2//"mb","m/s",
    173      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    174      .                "inst(X)", zstohf,zout)
     228     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
     229     .         "inst(X)", zstohf,zstohf)
    175230c
    176231          CALL histdef(nid_hf, "v"//bb2,
    177232     .                 "Meridional wind "//bb2//"mb","m/s",
    178      .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    179      .                "inst(X)", zstohf,zout)
    180 c
    181          ENDIF !IF(bb2.EQ."850".or.bb2.EQ."700".OR.
    182 c
    183          ENDDO !nlevSTD
    184 c
    185       ENDIF !lev_histhf.GE.3
    186 c
    187       IF(lev_histhf.GE.4) THEN
    188 c
    189 cIM sorties hf 3D
    190 c
    191 #undef histhf3d
    192 #define histhf3d
    193 #ifdef histhf3d
    194 #include "ini_histhf3d.h"
    195 #endif
    196       ENDIF !lev_histhf.GE.4
    197 c
    198       ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
    199 c
    200       IF(lev_histhf.GE.1) THEN
    201 c
    202         DO k=1, nlevSTD
    203 c
     233     .         iim,jj_nb,nhori, 1,1,1, -99, 32,
     234     .         "inst(X)", zstohf,zstohf)
     235c
     236         ENDIF !bb2.EQ."850".or.bb2.EQ."200"
     237        ENDDO !k=1, nlevSTD
     238c
     239        CALL histdef(nid_hf, "topl",
     240     .       "Outgoing longwave radiation", "W/m2",
     241     .       iim,jj_nb,nhori, 1,1,1, -99, 32,
     242     .       "ave(X)", zstophy,zout)
     243c
     244        CALL histdef(nid_hf, "precip","Total precipitation rate",
     245     .       "kg/(s*m2)",iim,jj_nb,nhori, 1,1,1, -99, 32,
     246     .       "ave(X)", zstophy,zout)
     247c
     248        CALL histdef(nid_hf, "slp",
     249     .       "Mean sea-level pressure", "Pa",
     250     .       iim,jj_nb,nhori, 1,1,1, -99, 32,
     251     .       "inst(X)", zstohf,zout)
     252c
     253cIM 280504 BEG
     254      ENDIF !(lev_histhf.GE.1) THEN
     255c
     256      IF(lev_histhf.GE.2) THEN
     257c
     258        DO k=1, nlevSTD
    204259c
    205260c        bb=clevSTD(k)
    … …  
    212267c        ENDIF
    213268c
    214          IF(bb2.EQ."850".or.bb2.EQ."200") THEN
    215 c
    216           CALL histdef(nid_hf, "u"//bb2,
    217      .                 "Zonal wind "//bb2//"mb","m/s",
    218      .         iim,jj_nb,nhori, 1,1,1, -99, 32,
    219      .         "inst(X)", zstohf,zstohf)
    220 c
    221           CALL histdef(nid_hf, "v"//bb2,
    222      .                 "Meridional wind "//bb2//"mb","m/s",
    223      .         iim,jj_nb,nhori, 1,1,1, -99, 32,
    224      .         "inst(X)", zstohf,zstohf)
    225 c
    226          ENDIF !bb2.EQ."850".or.bb2.EQ."200"
    227         ENDDO !k=1, nlevSTD
    228 c
    229         CALL histdef(nid_hf, "topl",
    230      .       "Outgoing longwave radiation", "W/m2",
    231      .       iim,jj_nb,nhori, 1,1,1, -99, 32,
    232      .       "ave(X)", zstophy,zout)
    233 c
    234         CALL histdef(nid_hf, "precip","Total precipitation rate",
    235      .       "kg/(s*m2)",iim,jj_nb,nhori, 1,1,1, -99, 32,
    236      .       "ave(X)", zstophy,zout)
    237 c
    238         CALL histdef(nid_hf, "slp",
    239      .       "Mean sea-level pressure", "Pa",
    240      .       iim,jj_nb,nhori, 1,1,1, -99, 32,
    241      .       "inst(X)", zstohf,zout)
    242 c
    243 cIM 280504 BEG
    244       ENDIF !(lev_histhf.GE.1) THEN
    245 c
    246       IF(lev_histhf.GE.2) THEN
    247 c
    248         DO k=1, nlevSTD
    249 c
    250 c        bb=clevSTD(k)
    251 c
    252 c        IF(k.GE.2) THEN
    253          IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
    254          IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
    255 c         aa=clevSTD(k)
    256 c         bb=aa(1:lnblnk1(aa))
    257 c        ENDIF
    258 c
    259269         IF(bb2.EQ."850".OR.bb2.EQ."500") THEN
    260270c

  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h

    r896 r897  
    292292     .                "ave(X)", zstophy,zout)
    293293c
    294 c
    295294c Effets des aerosols
    296295c
    … …  
    12861285     .                "ave(X)", zstophy,zout)
    12871286c
    1288 c          CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
     1287c Effets des aerosols
     1288c
     1289         IF (ok_ade.OR.ok_aie) THEN
     1290             CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",
     1291     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
     1292     .                    "ave(X)", zsto,zout)
     1293c
     1294             CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",
     1295     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
     1296     .                    "ave(X)", zsto,zout)
     1297c
     1298             CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",
     1299     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
     1300     .                    "ave(X)", zsto,zout)
     1301c
     1302             CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",
     1303     .                    iim,jj_nb,nhori, 1,1,1, -99, 32,
     1304     .                    "ave(X)", zsto,zout)
     1305         ENDIF
     1306c
     1307c        CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa",
    12891308c    .                iim,jj_nb,nhori, 1,1,1, -99, 32,
    12901309c    .                "ave(X)", zstophy,zout)

  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histday.h

    r892 r897  
    284284      CALL histwrite_phy(nid_day,"fder",itau_w,fder)
    285285c
    286       ENDIF !lev_histday.GE.2
    287 c
    288       IF(lev_histday.GE.3) THEN
    289 c=================================================================
    290 c   ECRITURE DES CHAMPS 3D
    291 c=================================================================
    292 
    293 c
    294 c Champs 3D:
    295 c
    296 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
    297       CALL histwrite_phy(nid_day,"temp",itau_w,t_seri)
    298 c
    299 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
    300       CALL histwrite_phy(nid_day,"ovap",itau_w,qx(:,:,ivap))
    301 c
    302 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
    303       CALL histwrite_phy(nid_day,"geop",itau_w,zphi)
    304 c
    305 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
    306       CALL histwrite_phy(nid_day,"vitu",itau_w,u_seri)
    307 c
    308 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
    309       CALL histwrite_phy(nid_day,"vitv",itau_w,v_seri)
    310 c
    311 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
    312       CALL histwrite_phy(nid_day,"vitw",itau_w,omega)
    313 c
    314 cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
    315       CALL histwrite_phy(nid_day,"pres",itau_w,pplay)
    316 c
    317       ENDIF !lev_histday.GE.3
    318 c=================================================================
    319 c   FIN ECRITURE DES CHAMPS 3D
    320 c=================================================================
    321       IF(lev_histday.ge.4) THEN
    322 c=================================================================
    323 c
    324 c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
    325 c
    326 c=================================================================
     286      DO i=1, klon
     287       zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
     288      ENDDO
     289c
     290cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1, zx_tmp_2d,zx_tmp_2d)
     291      CALL histwrite_phy(nid_day,"rh2m",itau_w,zx_tmp_fi2d)
    327292      zx_tmp_fi2d(1 : klon) = swup( 1 : klon, klevp1 )
    328293cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    … …  
    360325cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    361326      CALL histwrite_phy(nid_day,"LWupSFC",itau_w,zx_tmp_fi2d)
     327c
     328      zx_tmp_fi2d(1:klon) = flwp(1:klon)
     329cym     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     330      CALL histwrite_phy(nid_day,"lwp",itau_w,zx_tmp_fi2d)
     331c
     332      zx_tmp_fi2d(1:klon) = fiwp(1:klon)
     333cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     334      CALL histwrite_phy(nid_day,"iwp",itau_w,zx_tmp_fi2d)
     335c
     336      zx_tmp_fi2d(1 : klon) = zxfluxu(1 : klon, 1)
     337cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     338      CALL histwrite_phy(nid_day,"taux",itau_w,zx_tmp_fi2d)
     339c
     340      zx_tmp_fi2d(1 : klon) = zxfluxv(1 : klon, 1)
     341cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
     342      CALL histwrite_phy(nid_day,"tauy",itau_w,zx_tmp_fi2d)
    362343c
    363344cIM   IF (iflag_con.GE.3) THEN
    … …  
    369350cIM   ENDIF
    370351c
     352      ENDIF !lev_histday.GE.2
     353c
     354      IF(lev_histday.GE.3) THEN
     355c=================================================================
     356c   ECRITURE DES CHAMPS 3D
     357c=================================================================
     358c
     359c Champs 3D:
     360c
     361cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
     362      CALL histwrite_phy(nid_day,"temp",itau_w,t_seri)
     363c
     364cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
     365      CALL histwrite_phy(nid_day,"ovap",itau_w,qx(:,:,ivap))
     366c
     367cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
     368      CALL histwrite_phy(nid_day,"geop",itau_w,zphi)
     369c
     370cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
     371      CALL histwrite_phy(nid_day,"vitu",itau_w,u_seri)
     372c
     373cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
     374      CALL histwrite_phy(nid_day,"vitv",itau_w,v_seri)
     375c
     376cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
     377      CALL histwrite_phy(nid_day,"vitw",itau_w,omega)
     378c
     379cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
     380      CALL histwrite_phy(nid_day,"pres",itau_w,pplay)
     381c
     382      ENDIF !lev_histday.GE.3
     383c=================================================================
     384c   FIN ECRITURE DES CHAMPS 3D
     385c=================================================================
     386      IF(lev_histday.ge.4) THEN
     387c=================================================================
     388c
     389c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
     390c
     391c=================================================================
    371392      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
    372393cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
    … …  
    522543      CALL histwrite_phy(nid_day,"iwcon",itau_w,fiwc)
    523544c
    524       zx_tmp_fi2d(1:klon) = flwp(1:klon)
    525 cym     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    526       CALL histwrite_phy(nid_day,"lwp",itau_w,zx_tmp_fi2d)
    527 c
    528       zx_tmp_fi2d(1:klon) = fiwp(1:klon)
    529 cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    530       CALL histwrite_phy(nid_day,"iwp",itau_w,zx_tmp_fi2d)
    531 c
    532545cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1, meantaucld,zx_tmp_2d)
     546      IF (ok_isccp) THEN
    533547      CALL histwrite_phy(nid_day,"meantaucld",itau_w,meantaucld)
     548      ENDIF
    534549c
    535550cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau,zx_tmp_3d)
    … …  
    548563cym      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
    549564      CALL histwrite_phy(nid_day,"dqcon",itau_w,zx_tmp_fi3d)
    550 c
    551       DO i=1, klon
    552        zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
    553       ENDDO
    554 c
    555 cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1, zx_tmp_2d,zx_tmp_2d)
    556       CALL histwrite_phy(nid_day,"rh2m",itau_w,zx_tmp_fi2d)
    557565c
    558566cym      CALL gr_fi_ecrit(1,klon,iim,jjmp1, qsat2m,zx_tmp_2d)

  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histhf.h

    r782 r897  
    117117c
    118118c        IF(k.GE.2) THEN
    119          IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
    120          IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
     119c         IF(k.GE.1.AND.k.LE.10) bb2=clevSTD(k)
     120c         IF(k.GE.11.AND.k.LE.17) bb3=clevSTD(k)
    121121c         aa=clevSTD(k)
    122122c         bb=aa(1:lnblnk1(aa))
    123123c        ENDIF
    124124c
    125       IF(bb2.EQ."850".OR.bb2.EQ."700".OR.
    126      $   bb2.EQ."500".OR.bb2.EQ."200") THEN
     125      bb2=clevSTD(k)
     126      IF(bb2.EQ."1000".OR.bb2.EQ."925".OR.
     127     $      bb2.EQ."850".OR.bb2.EQ."700".OR.
     128     $      bb2.EQ."500".OR.bb2.EQ."300".OR.
     129     $      bb2.EQ."250".OR.bb2.EQ."200".OR.
     130     $      bb2.EQ."150".OR.bb2.EQ."100".OR.
     131     $      bb2.EQ."50".OR.bb2.EQ."30".OR.
     132     $      bb2.EQ."10")THEN
    127133c
    128134cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
    … …  
    135141      ENDIF !bb2.NE."500"
    136142c
     143      IF(bb2.GE."200") THEN
     144c
    137145cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
    138146      CALL histwrite_phy(nid_hf,"q"//bb2,itau_w, qlevSTD(:,k))
     147      ENDIF
    139148c
    140149cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
    … …  
    144153      CALL histwrite_phy(nid_hf,"v"//bb2,itau_w, vlevSTD(:,k))
    145154c
    146       ENDIF !bb2.EQ."850".OR.bb2.EQ."700"
     155cym      CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlevSTD(:,k),zx_tmp_2d)
     156      CALL histwrite_phy(nid_hf,"w"//bb2,itau_w, wlevSTD(:,k))
     157c
     158      ENDIF
    147159      ENDDO !nlevENS
    148160c

  • LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h

    r896 r897  
    321321        CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w,
    322322     $      zx_tmp_fi2d)
    323        ENDIF !nsrf.GT.2
     323       ENDIF !nsrf.GE.2
    324324C
    325325        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
    … …  
    12861286      CALL histwrite_phy(nid_mth,"fder",itau_w,fder)
    12871287c
     1288c Effets des aerosols
     1289c
     1290      IF (ok_ade.OR.ok_aie) THEN
     1291        zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon)
     1292        CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d)
     1293c
     1294        zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon)
     1295        CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d)
     1296c
     1297        zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon)
     1298        CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d)
     1299c
     1300        zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon)
     1301        CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d)
     1302      ENDIF
     1303c
    12881304c
    12891305c      DO i = 1, klon
Note: See TracChangeset for help on using the changeset viewer.