Changeset 1419 for LMDZ4


Ignore:
Timestamp:
Jul 22, 2010, 5:19:06 PM (14 years ago)
Author:
musat
Message:

Modify levels for some outputs to diminsh the CMIP5 storage (JLD)
Add Hines GWD tendencies duhin, dvhin (F.Lott)
IM

Location:
LMDZ4/branches/LMDZ4_AR5/libf/phylmd
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/ini_histhfNMC.h

    r1400 r1419  
    8888     .                "ave(X)", zout,zout)
    8989c
     90          IF (1.EQ.0) THEN
     91c
    9092          CALL histdef(nid_hfnmc, "psbg",
    9193     .         "Pressure sfce below ground","%",
     
    146148     .         "ave(X)", zout,zout)
    147149c
     150          ENDIF !(1.EQ.0) THEN
    148151         CALL histend(nid_hfnmc)
    149152c
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_mod.F90

    r1405 r1419  
    6666  type(ctrl_out),save :: o_t2m_min      = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_min')
    6767  type(ctrl_out),save :: o_t2m_max      = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_max')
    68   type(ctrl_out),save,dimension(4) :: o_t2m_srf      = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_ter'), &
    69                                                  ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_lic'), &
    70                                                  ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_oce'), &
    71                                                  ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_sic') /)
     68  type(ctrl_out),save,dimension(4) :: o_t2m_srf      = (/ ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_ter'), &
     69                                                 ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_lic'), &
     70                                                 ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_oce'), &
     71                                                 ctrl_out((/ 10, 6, 10, 10, 10 /),'t2m_sic') /)
    7272
    7373  type(ctrl_out),save :: o_wind10m      = ctrl_out((/ 1, 1, 1, 10, 10 /),'wind10m')
     
    8080  type(ctrl_out),save :: o_qsurf        = ctrl_out((/ 1, 10, 10, 10, 10 /),'qsurf')
    8181
    82   type(ctrl_out),save,dimension(4) :: o_u10m_srf     = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_ter'), &
    83                                               ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_lic'), &
    84                                               ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_oce'), &
    85                                               ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_sic') /)
    86 
    87   type(ctrl_out),save,dimension(4) :: o_v10m_srf     = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_ter'), &
    88                                               ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_lic'), &
    89                                               ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_oce'), &
    90                                               ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_sic') /)
     82  type(ctrl_out),save,dimension(4) :: o_u10m_srf     = (/ ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_ter'), &
     83                                              ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_lic'), &
     84                                              ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_oce'), &
     85                                              ctrl_out((/ 10, 6, 10, 10, 10 /),'u10m_sic') /)
     86
     87  type(ctrl_out),save,dimension(4) :: o_v10m_srf     = (/ ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_ter'), &
     88                                              ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_lic'), &
     89                                              ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_oce'), &
     90                                              ctrl_out((/ 10, 6, 10, 10, 10 /),'v10m_sic') /)
    9191
    9292  type(ctrl_out),save :: o_qsol         = ctrl_out((/ 1, 10, 10, 10, 10 /),'qsol')
     
    9999  type(ctrl_out),save :: o_snow         = ctrl_out((/ 1, 1, 10, 10, 5 /),'snow')
    100100  type(ctrl_out),save :: o_evap         = ctrl_out((/ 1, 1, 10, 10, 10 /),'evap')
    101   type(ctrl_out),save,dimension(4) :: o_evap_srf     = (/ ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_ter'), &
    102                                            ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_lic'), &
    103                                            ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_oce'), &
    104                                            ctrl_out((/ 1, 1, 10, 10, 10 /),'evap_sic') /)
     101  type(ctrl_out),save,dimension(4) :: o_evap_srf     = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_ter'), &
     102                                           ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_lic'), &
     103                                           ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_oce'), &
     104                                           ctrl_out((/ 1, 6, 10, 10, 10 /),'evap_sic') /)
    105105  type(ctrl_out),save :: o_msnow       = ctrl_out((/ 1, 10, 10, 10, 10 /),'msnow')
    106106  type(ctrl_out),save :: o_fsnow       = ctrl_out((/ 1, 10, 10, 10, 10 /),'fsnow')
     
    135135  type(ctrl_out),save :: o_soll         = ctrl_out((/ 1, 1, 10, 10, 10 /),'soll')
    136136  type(ctrl_out),save :: o_soll0        = ctrl_out((/ 1, 5, 10, 10, 10 /),'soll0')
    137   type(ctrl_out),save :: o_radsol       = ctrl_out((/ 1, 1, 10, 10, 10 /),'radsol')
     137  type(ctrl_out),save :: o_radsol       = ctrl_out((/ 1, 7, 10, 10, 10 /),'radsol')
    138138  type(ctrl_out),save :: o_SWupSFC      = ctrl_out((/ 1, 4, 10, 10, 5 /),'SWupSFC')
    139139  type(ctrl_out),save :: o_SWupSFCclr   = ctrl_out((/ 1, 4, 10, 10, 5 /),'SWupSFCclr')
     
    153153  type(ctrl_out),save :: o_taux         = ctrl_out((/ 1, 10, 10, 10, 10 /),'taux')
    154154  type(ctrl_out),save :: o_tauy         = ctrl_out((/ 1, 10, 10, 10, 10 /),'tauy')
    155   type(ctrl_out),save,dimension(4) :: o_taux_srf     = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_ter'), &
    156                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_lic'), &
    157                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_oce'), &
    158                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'taux_sic') /)
    159 
    160   type(ctrl_out),save,dimension(4) :: o_tauy_srf     = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_ter'), &
    161                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_lic'), &
    162                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_oce'), &
    163                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'tauy_sic') /)
    164 
    165 
    166   type(ctrl_out),save,dimension(4) :: o_pourc_srf    = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_ter'), &
    167                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_lic'), &
    168                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_oce'), &
    169                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'pourc_sic') /)     
    170 
    171   type(ctrl_out),save,dimension(4) :: o_fract_srf    = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_ter'), &
    172                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_lic'), &
    173                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_oce'), &
    174                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'fract_sic') /)
    175 
    176   type(ctrl_out),save,dimension(4) :: o_tsol_srf     = (/ ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_ter'), &
    177                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_lic'), &
    178                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_oce'), &
    179                                                  ctrl_out((/ 1, 4, 10, 10, 10 /),'tsol_sic') /)
    180 
    181   type(ctrl_out),save,dimension(4) :: o_sens_srf     = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_ter'), &
    182                                                  ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_lic'), &
    183                                                  ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_oce'), &
    184                                                  ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_sic') /)
    185 
    186   type(ctrl_out),save,dimension(4) :: o_lat_srf      = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_ter'), &
    187                                                  ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_lic'), &
    188                                                  ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_oce'), &
    189                                                  ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_sic') /)
     155  type(ctrl_out),save,dimension(4) :: o_taux_srf     = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_ter'), &
     156                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_lic'), &
     157                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_oce'), &
     158                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'taux_sic') /)
     159
     160  type(ctrl_out),save,dimension(4) :: o_tauy_srf     = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_ter'), &
     161                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_lic'), &
     162                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_oce'), &
     163                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'tauy_sic') /)
     164
     165
     166  type(ctrl_out),save,dimension(4) :: o_pourc_srf    = (/ ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_ter'), &
     167                                                 ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_lic'), &
     168                                                 ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_oce'), &
     169                                                 ctrl_out((/ 1, 7, 10, 10, 10 /),'pourc_sic') /)     
     170
     171  type(ctrl_out),save,dimension(4) :: o_fract_srf    = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_ter'), &
     172                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_lic'), &
     173                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_oce'), &
     174                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'fract_sic') /)
     175
     176  type(ctrl_out),save,dimension(4) :: o_tsol_srf     = (/ ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_ter'), &
     177                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_lic'), &
     178                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_oce'), &
     179                                                 ctrl_out((/ 1, 6, 10, 10, 10 /),'tsol_sic') /)
     180
     181  type(ctrl_out),save,dimension(4) :: o_sens_srf     = (/ ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_ter'), &
     182                                                 ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_lic'), &
     183                                                 ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_oce'), &
     184                                                 ctrl_out((/ 1, 6, 10, 7, 10 /),'sens_sic') /)
     185
     186  type(ctrl_out),save,dimension(4) :: o_lat_srf      = (/ ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_ter'), &
     187                                                 ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_lic'), &
     188                                                 ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_oce'), &
     189                                                 ctrl_out((/ 1, 6, 10, 7, 10 /),'lat_sic') /)
    190190
    191191  type(ctrl_out),save,dimension(4) :: o_flw_srf      = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_ter'), &
     
    211211
    212212  type(ctrl_out),save :: o_cdrm         = ctrl_out((/ 1, 10, 10, 10, 10 /),'cdrm')
    213   type(ctrl_out),save :: o_cdrh         = ctrl_out((/ 1, 10, 10, 1, 10 /),'cdrh')
     213  type(ctrl_out),save :: o_cdrh         = ctrl_out((/ 1, 10, 10, 7, 10 /),'cdrh')
    214214  type(ctrl_out),save :: o_cldl         = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldl')
    215215  type(ctrl_out),save :: o_cldm         = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldm')
     
    257257! Champs interpolles sur des niveaux de pression ??? a faire correctement
    258258                                             
    259   type(ctrl_out),save,dimension(6) :: o_uSTDlevs     = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'u850'), &
    260                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'u700'), &
    261                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'u500'), &
    262                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'u200'), &
    263                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'u50'), &
    264                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'u10') /)
     259  type(ctrl_out),save,dimension(6) :: o_uSTDlevs     = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'u850'), &
     260                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'u700'), &
     261                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'u500'), &
     262                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'u200'), &
     263                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'u50'), &
     264                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'u10') /)
    265265                                                     
    266266
    267   type(ctrl_out),save,dimension(6) :: o_vSTDlevs     = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'v850'), &
    268                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'v700'), &
    269                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'v500'), &
    270                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'v200'), &
    271                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'v50'), &
    272                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'v10') /)
    273 
    274   type(ctrl_out),save,dimension(6) :: o_wSTDlevs     = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'w850'), &
    275                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'w700'), &
    276                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'w500'), &
    277                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'w200'), &
    278                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'w50'), &
    279                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'w10') /)
    280 
    281   type(ctrl_out),save,dimension(6) :: o_tSTDlevs     = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'t850'), &
    282                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'t700'), &
    283                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'t500'), &
    284                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'t200'), &
    285                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'t50'), &
    286                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'t10') /)
    287 
    288   type(ctrl_out),save,dimension(6) :: o_qSTDlevs     = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'q850'), &
    289                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'q700'), &
    290                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'q500'), &
    291                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'q200'), &
    292                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'q50'), &
    293                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'q10') /)
    294 
    295   type(ctrl_out),save,dimension(6) :: o_zSTDlevs   = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'z850'), &
    296                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'z700'), &
    297                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'z500'), &
    298                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'z200'), &
    299                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'z50'), &
    300                                                      ctrl_out((/ 1, 1, 3, 10, 10 /),'z10') /)
     267  type(ctrl_out),save,dimension(6) :: o_vSTDlevs     = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'v850'), &
     268                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'v700'), &
     269                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'v500'), &
     270                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'v200'), &
     271                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'v50'), &
     272                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'v10') /)
     273
     274  type(ctrl_out),save,dimension(6) :: o_wSTDlevs     = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'w850'), &
     275                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'w700'), &
     276                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'w500'), &
     277                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'w200'), &
     278                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'w50'), &
     279                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'w10') /)
     280
     281  type(ctrl_out),save,dimension(6) :: o_tSTDlevs     = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'t850'), &
     282                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'t700'), &
     283                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'t500'), &
     284                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'t200'), &
     285                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'t50'), &
     286                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'t10') /)
     287
     288  type(ctrl_out),save,dimension(6) :: o_qSTDlevs     = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'q850'), &
     289                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'q700'), &
     290                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'q500'), &
     291                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'q200'), &
     292                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'q50'), &
     293                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'q10') /)
     294
     295  type(ctrl_out),save,dimension(6) :: o_zSTDlevs   = (/ ctrl_out((/ 1, 7, 7, 10, 10 /),'z850'), &
     296                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'z700'), &
     297                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'z500'), &
     298                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'z200'), &
     299                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'z50'), &
     300                                                     ctrl_out((/ 1, 7, 7, 10, 10 /),'z10') /)
    301301
    302302
     
    419419  type(ctrl_out),save :: o_theta        = ctrl_out((/ 2, 3, 4, 10, 10 /),'theta')
    420420  type(ctrl_out),save :: o_ovap         = ctrl_out((/ 2, 3, 4, 10, 10 /),'ovap')
    421   type(ctrl_out),save :: o_ovapinit     = ctrl_out((/ 2, 3, 10, 10, 10 /),'ovapinit')
     421  type(ctrl_out),save :: o_ovapinit     = ctrl_out((/ 2, 10, 10, 10, 10 /),'ovapinit')
    422422  type(ctrl_out),save :: o_wvapp        = ctrl_out((/ 2, 10, 10, 10, 10 /),'wvapp')
    423423  type(ctrl_out),save :: o_geop         = ctrl_out((/ 2, 3, 10, 10, 10 /),'geop')
     
    449449!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    450450
    451   type(ctrl_out),save,dimension(4) :: o_albe_srf     = (/ ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_ter'), &
    452                                                      ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_lic'), &
    453                                                      ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_oce'), &
    454                                                      ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_sic') /)
     451  type(ctrl_out),save,dimension(4) :: o_albe_srf     = (/ ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_ter'), &
     452                                                     ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_lic'), &
     453                                                     ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_oce'), &
     454                                                     ctrl_out((/ 3, 7, 10, 7, 10 /),'albe_sic') /)
    455455
    456456  type(ctrl_out),save,dimension(4) :: o_ages_srf     = (/ ctrl_out((/ 10, 10, 10, 10, 10 /),'ages_ter'), &
     
    459459                                                     ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_sic') /)
    460460
    461   type(ctrl_out),save,dimension(4) :: o_rugs_srf     = (/ ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_ter'), &
    462                                                      ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_lic'), &
    463                                                      ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_oce'), &
    464                                                      ctrl_out((/ 3, 4, 10, 10, 10 /),'rugs_sic') /)
     461  type(ctrl_out),save,dimension(4) :: o_rugs_srf     = (/ ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_ter'), &
     462                                                     ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_lic'), &
     463                                                     ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_oce'), &
     464                                                     ctrl_out((/ 3, 6, 10, 10, 10 /),'rugs_sic') /)
    465465
    466466  type(ctrl_out),save :: o_alb1         = ctrl_out((/ 3, 10, 10, 10, 10 /),'alb1')
     
    477477  type(ctrl_out),save :: o_dudyn        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dudyn')  !AXC
    478478  type(ctrl_out),save :: o_dvdyn        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvdyn')  !AXC
    479   type(ctrl_out),save :: o_dtcon        = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtcon')
     479  type(ctrl_out),save :: o_dtcon        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtcon')
    480480  type(ctrl_out),save :: o_ducon        = ctrl_out((/ 4, 10, 10, 10, 10 /),'ducon')
    481   type(ctrl_out),save :: o_dqcon        = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqcon')
     481  type(ctrl_out),save :: o_dqcon        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqcon')
    482482  type(ctrl_out),save :: o_dtwak        = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtwak')
    483483  type(ctrl_out),save :: o_dqwak        = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqwak')
     
    524524  type(ctrl_out),save :: o_dulif        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dulif')
    525525  type(ctrl_out),save :: o_dvlif        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvlif')
     526  type(ctrl_out),save :: o_duhin        = ctrl_out((/ 4, 10, 10, 10, 10 /),'duhin')
     527  type(ctrl_out),save :: o_dvhin        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvhin')
    526528
    527529! Attention a refaire correctement
     
    12241226     ENDIF
    12251227
     1228     IF (ok_hines) then
     1229 CALL histdef3d(iff,o_duhin%flag,o_duhin%name, "Hines GWD dU", "m/s2")
     1230 CALL histdef3d(iff,o_dvhin%flag,o_dvhin%name, "Hines GWD dV", "m/s2")
     1231     ENDIF
     1232
    12261233
    12271234!Attention : sorties uniquement pour traceurs 3 et 4
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phys_output_write.h

    r1405 r1419  
    15091509        ENDIF
    15101510
     1511       IF (ok_hines) THEN
     1512       IF (o_duhin%flag(iff)<=lev_files(iff)) THEN
     1513      zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys
     1514      CALL histwrite_phy(nid_files(iff),o_duhin%name,itau_w,zx_tmp_fi3d)
     1515       ENDIF
     1516
     1517        IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN
     1518      zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys
     1519      CALL histwrite_phy(nid_files(iff),o_dvhin%name,itau_w,zx_tmp_fi3d)
     1520       ENDIF
     1521        ENDIF
    15111522
    15121523        IF (nqtot.GE.3 .AND. o_trac(1)%flag(iff)<=lev_files(iff)) THEN
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histhfNMC.h

    r1398 r1419  
    3030       CALL histwrite_phy(nid_hfnmc,"wap",itau_w,wwriteSTD3)
    3131c
     32       IF (1.EQ.0) THEN
     33c
    3234       DO k=1, nlevSTD
    3335        DO i=1, klon
     
    6062       CALL histwrite_phy(nid_hfnmc,"TxT",itau_w,T2sumSTD(:,:,3))
    6163c
    62 c     ENDIF !type_run
     64       ENDIF !(1.EQ.0) THEN
    6365c
    6466      if (ok_sync) then
Note: See TracChangeset for help on using the changeset viewer.