Ignore:
Timestamp:
Dec 19, 2002, 5:46:39 PM (22 years ago)
Author:
lmdzadmin
Message:

Convergence avec la version de Ionela dec 2002

YOMCST.? : suppression RI0 (IM)
albedo.F : facteur 1.2 sur le nouveau calcul (IM)
clesphys.h : rajout de différentes ctes (concentration des gaz) (IM)
clmain.F : separation des flux LW, SW (JLD)

remplace qsurf par yqsol (IM)

conf_phys.F90 : rajout de différentes ctes (gaz + orbite) (IM)
convect3.F : DPINV+SIGD*0.5*(EVAP(1)+EVAP(2)) (SBL)
cv3_routines.F:
cvparam3.h : compatibilite avec conema3 TEMPORAIRE (FH)
phyetat0.F : lecture de co2_ppm et solaire pour tests de coherence
phyredem.F : co2_ppm et solaire passé en common
physiq.F : separation flux LW, SW

rajout diagnostiques (slp, w500)
suppression iflag_con = 4
clwcon0=qcondc (FH)
position dU "ENDIF ! ok_cvl"

radlwsw.F : passage des concentrations gaz dans un common (IM)

PEMIS(i) = 1.0 (JLD pour cohérence ORCHIDEE)

stdlevvar.F90 :
suphec.F : suppression init. des ctes orbitales (IM)

nouvelles E/S (ini_hist..., write_hist...)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histday.h

    r418 r433  
    3434c
    3535C
    36       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
    37       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
    38       CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    39 c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
    40 c    .              'Surface Temperature','K')
    41 C
    42       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
    43       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    44       CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    45 c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
    46 c    .              'Surface Temperature','K')
    47 C
    48       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
    49       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    50       CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    51 c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
    52 c    .              'Surface Temperature','K')
    53 C
    54       zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
    55       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    56       CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    57 c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
    58 c    .              'Surface Temperature','K')
    59 C
    60 cccIM
     36C
    6137      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    6238      CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    6339c
    64       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    65       CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
    66      .               iim*jjmp1,ndex2d)
    67 c
    68       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
    69       CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
    70      .               iim*jjmp1,ndex2d)
    71 c
    72       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    73       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    74       CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
    75      .               ndex2d)
    76 c
    77       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    78       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    79       CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
    80      .               iim*jjmp1,ndex2d)
    81 c
    82       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
    83       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    84       CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
    85      .               iim*jjmp1,ndex2d)
    86 c
    87       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
    88       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    89       CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    90      .               ndex2d)
    91 c
    92       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
    93       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    94       CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
    95      .               ndex2d)
    96 c
    97       zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
    98       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    99       CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
    100      .               ndex2d)
    101 c
    10240      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
    10341      CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    … …  
    11149     .               ndex2d)
    11250c
    113       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
    114       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    115       CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
    116      .               ndex2d)
    117 c
    118       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
    119       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    120       CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
    121      .               ndex2d)
    122 c
    123       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
    124       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    125       CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    126      .               ndex2d)
    127 c
    128       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
    129       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    130       CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
    131      .               ndex2d)
    132 c
    133       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
    134       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    135       CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
    136      .               ndex2d)
    137 c
    138       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
    139       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    140       CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
    141      .               ndex2d)
    142 c
    143       zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
    144       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    145       CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
    146      .               ndex2d)
    147 C
    148       zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
    149       CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
    150       CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
    151      .               ndex2d)
    152 C
    15351      DO i = 1, klon
    15452         zx_tmp_fi2d(i) = paprs(i,1)
    … …  
    231129c
    232130c
     131cXXX      DO i = 1, klon
     132cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
     133cXXX      ENDDO
     134cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     135cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     136c
     137      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
     138      CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     139c     call writephy(fid_day,prof2d_av,'cldl',cldl,
     140c    .              'Low-level cloudiness','-')
     141c
     142      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
     143      CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     144c     call writephy(fid_day,prof2d_av,'cldm',cldm,
     145c    .              'Mid-level cloudiness','-')
     146c
     147      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
     148      CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     149c     call writephy(fid_day,prof2d_av,'cldh',cldh,
     150c    .              'High-level cloudiness','-')
     151c
     152      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
     153      CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     154c     call writephy(fid_day,prof2d_av,'cldt',cldt,
     155c    .              'Total cloudiness','-')
     156c
     157      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
     158      CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     159c     call writephy(fid_day,prof2d_av,'cldq',cldq,
     160c    .              'Cloud liquid water path','-')
     161c
     162      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
     163      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     164      CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
     165     .                               iim*jjmp1,ndex2d)
     166c
     167      zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
     168      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     169      CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
     170     .                               iim*jjmp1,ndex2d)
     171c
     172      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
     173      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     174      CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
     175     .                               iim*jjmp1,ndex2d)
     176c
     177      zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
     178      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     179      CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
     180     .                               iim*jjmp1,ndex2d)
     181
     182
     183c   Ecriture de champs dynamiques sur des niveaux de pression
     184
     185      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u850,zx_tmp_2d)
     186      CALL histwrite(nid_day,"u850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     187
     188      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v850,zx_tmp_2d)
     189      CALL histwrite(nid_day,"v850",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     190
     191      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u500,zx_tmp_2d)
     192      CALL histwrite(nid_day,"u500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     193
     194      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v500,zx_tmp_2d)
     195      CALL histwrite(nid_day,"v500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     196
     197      CALL gr_fi_ecrit(1, klon,iim,jjmp1, u200,zx_tmp_2d)
     198      CALL histwrite(nid_day,"u200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     199
     200      CALL gr_fi_ecrit(1, klon,iim,jjmp1, v200,zx_tmp_2d)
     201      CALL histwrite(nid_day,"v200",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     202
     203      CALL gr_fi_ecrit(1, klon,iim,jjmp1, phi500,zx_tmp_2d)
     204      CALL histwrite(nid_day,"phi500",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     205
     206      CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
     207      CALL histwrite(nid_day,"slp",itap,zx_tmp_2d,iim*jjmp1,ndex2d)
     208
     209
     210c
     211c=================================================================
     212c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
     213c=================================================================
     214      if ( lev_histday.gt.1 ) then
     215
     216c   En attendant un eventuel debugage.
     217
     218      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     219      CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
     220     .               iim*jjmp1,ndex2d)
     221c
     222      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
     223      CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
     224     .               iim*jjmp1,ndex2d)
     225c
     226c=================================================================
     227c=================================================================
     228c=================================================================
     229c   ECRITURE DES CHAMPS SUR LES SOUS SURFACES
     230c=================================================================
     231
     232      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
     233      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
     234      CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     235c     call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter),
     236c    .              'Surface Temperature','K')
     237C
     238      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
     239      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     240      CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     241c     call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic),
     242c    .              'Surface Temperature','K')
     243C
     244      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
     245      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     246      CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     247c     call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce),
     248c    .              'Surface Temperature','K')
     249C
     250      zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
     251      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
     252      CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     253c     call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic),
     254c    .              'Surface Temperature','K')
     255C
     256cccIM
     257      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     258      CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
     259     .               ndex2d)
     260c
     261      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
     262      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     263      CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
     264     .               iim*jjmp1,ndex2d)
     265c
     266      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
     267      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
     268      CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
     269     .               iim*jjmp1,ndex2d)
     270c
     271      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
     272      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     273      CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
     274     .               ndex2d)
     275c
     276      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
     277      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     278      CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
     279     .               ndex2d)
     280c
     281      zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
     282      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     283      CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
     284     .               ndex2d)
     285c
     286      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
     287      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     288      CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
     289     .               ndex2d)
     290c
     291      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
     292      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     293      CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
     294     .               ndex2d)
     295c
     296      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
     297      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     298      CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
     299     .               ndex2d)
     300c
     301      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
     302      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     303      CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
     304     .               ndex2d)
     305c
     306      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
     307      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     308      CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
     309     .               ndex2d)
     310c
     311      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
     312      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     313      CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
     314     .               ndex2d)
     315c
     316      zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
     317      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     318      CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
     319     .               ndex2d)
     320C
     321      zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
     322      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
     323      CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
     324     .               ndex2d)
     325C
    233326      DO nsrf = 1, nbsrf
    234327CXXX
    … …  
    301394C
    302395      END DO 
    303 C
    304 cXXX      DO i = 1, klon
    305 cXXX         zx_tmp_fi2d(i) = pctsrf(i,is_sic)
    306 cXXX      ENDDO
    307 cXXX      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    308 cXXX      CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    309 c
    310       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
    311       CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    312 c     call writephy(fid_day,prof2d_av,'cldl',cldl,
    313 c    .              'Low-level cloudiness','-')
    314 c
    315       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
    316       CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    317 c     call writephy(fid_day,prof2d_av,'cldm',cldm,
    318 c    .              'Mid-level cloudiness','-')
    319 c
    320       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
    321       CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    322 c     call writephy(fid_day,prof2d_av,'cldh',cldh,
    323 c    .              'High-level cloudiness','-')
    324 c
    325       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
    326       CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    327 c     call writephy(fid_day,prof2d_av,'cldt',cldt,
    328 c    .              'Total cloudiness','-')
    329 c
    330       CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
    331       CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    332 c     call writephy(fid_day,prof2d_av,'cldq',cldq,
    333 c    .              'Cloud liquid water path','-')
     396c=================================================================
     397c   FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
     398c=================================================================
     399
     400c=================================================================
     401c   ECRITURE DES CHAMPS 3D
     402c=================================================================
     403
    334404c
    335405c Champs 3D:
    … …  
    380450c    .              'Air pressure','Pa')
    381451cccIM
    382       zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, klevp1)
    383       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    384       CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
    385      .                               iim*jjmp1,ndex2d)
    386 c
    387       zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1)
    388       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    389       CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
    390      .                               iim*jjmp1,ndex2d)
    391 c
    392       zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, klevp1)
    393       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    394       CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
    395      .                               iim*jjmp1,ndex2d)
    396 c
    397       zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1)
    398       CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
    399       CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
    400      .                               iim*jjmp1,ndex2d)
    401 c
     452
     453c=================================================================
     454c   FIN ECRITURE DES CHAMPS 3D
     455c=================================================================
     456
     457c=================================================================
     458c=================================================================
     459      endif !  lev_histday.gt.1
     460c=================================================================
     461
     462
    402463      if (ok_sync) then
    403464c       call writephy_sync(fid_day)
    404465        call histsync(nid_day)
    405466      endif
     467
    406468      ENDIF
Note: See TracChangeset for help on using the changeset viewer.