Ignore:
Timestamp:
Jul 19, 2002, 11:59:46 AM (22 years ago)
Author:
lmdzadmin
Message:

Changement dans les sorties histoires (rain devient precip pour eviter toute
confusion) + unites, commentaires des fichiers histoire
LF

Location:
LMDZ.3.3/branches/rel-LF/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90

    r386 r395  
    55  subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
    66 &                     fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &
    7  &                     ratqsbas,ratqshaut)
     7 &                     ratqsbas,ratqshaut,if_ebil)
    88
    99   use IOIPSL
     
    3434  logical              :: ok_journe, ok_mensuel, ok_instan
    3535  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
    36   integer              :: iflag_cldcon
     36  integer              :: iflag_cldcon, if_ebil
    3737! Local
    3838  integer              :: numout = 6
     
    8585  ok_instan = .false.
    8686  call getin('OK_instan', ok_instan)
     87!
     88!Config Key  = if_ebil
     89!Config Desc = Niveau de sortie pour les diags bilan d'energie
     90!Config Def  = 0
     91!Config Help =
     92!               
     93!
     94  if_ebil = 0
     95  call getin('if_ebil', if_ebil)
    8796!!
    8897!! KE
     
    282291  write(numout,*)' Sortie mensuelle = ', ok_mensuel
    283292  write(numout,*)' Sortie instantanee = ', ok_instan
     293  write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil
    284294  write(numout,*)' epmax = ', epmax
    285295  write(numout,*)' ok_adj_ema = ', ok_adj_ema
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

    r391 r395  
    645645c======================================================================
    646646      modname = 'physiq'
    647       if_ebil = 2
    648647      IF (if_ebil.ge.1) THEN
    649648        DO i=1,klon
     
    675674         call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel,
    676675     .                  ok_instan, fact_cldcon, facttemps,ok_newmicro,
    677      .                  iflag_cldcon,ratqsbas,ratqshaut)
     676     .                  iflag_cldcon,ratqsbas,ratqshaut, if_ebil)
    678677
    679678         DO k = 2, nvm          ! pas de vegetation
     
    881880     .                "ave(X)", zsto,zout)
    882881c
    883          CALL histdef(nid_day, "rain","Precipitation Totale liq+sol"
     882         CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
    884883     .                , "kg/s",
    885884     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     
    890889     .                "ave(X)", zsto,zout)
    891890c
    892          CALL histdef(nid_day, "snow_cov", "Snow cover", "mm",
     891         CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
    893892     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    894893     .                "ave(X)", zsto,zout)
     
    10951094     .                "ave(X)", zsto,zout)
    10961095c
    1097          CALL histdef(nid_mth, "rain", "Precipitation Totale liq+sol",
     1096         CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol",
    10981097     .                "kg/s",
    10991098     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     
    11121111     .                "ave(X)", zsto,zout)
    11131112c
    1114          CALL histdef(nid_mth, "snow_cov", "Snow cover", "mm",
     1113         CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2",
    11151114     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    11161115     .                "ave(X)", zsto,zout)
     
    13361335c
    13371336         CALL histdef(nid_mth, "rhum", "Relative humidity", "-",
     1337     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     1338     .                "ave(X)", zsto,zout)
     1339c
     1340         CALL histdef(nid_mth, "clwcon", "Cloud Liquid water content"
     1341     .                , "kg/kg",
    13381342     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    13391343     .                "ave(X)", zsto,zout)
     
    15621566     .                "inst(X)", zsto,zout)
    15631567c
    1564          CALL histdef(nid_ins, "rain", "Precipitation Totale liq+sol",
     1568         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
    15651569     .                "kg/s",
    15661570     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     
    15711575     .                "inst(X)", zsto,zout)
    15721576c
    1573          CALL histdef(nid_ins, "snow_cov", "Snow cover", "mm",
     1577         CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
    15741578     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    15751579     .                "inst(X)", zsto,zout)
     
    19011905      DO nsrf = 1, nbsrf
    19021906      DO i = 1, klon
    1903 c$$$         frugs(i,nsrf) = MAX(frugs(i,nsrf),0.001
    1904         frugs(i,nsrf) = MAX(frugs(i,nsrf),0.000015)
     1907         frugs(i,nsrf) = MAX(frugs(i,nsrf),0.001)
     1908cccc        frugs(i,nsrf) = MAX(frugs(i,nsrf),0.000015)
    19051909      ENDDO
    19061910      ENDDO
     
    21102114     .        ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr
    21112115     .        ,clwcon0)
    2112 c          print*,'Apres conema3 '
     2116          print*,'Apres conema3 '
    21132117
    21142118c Calculer l'humidite relative pour diagnostique
     
    27702774      ENDDO
    27712775      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    2772       CALL histwrite(nid_day,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     2776      CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    27732777c     call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d,
    27742778c    .              'Precipitation','mm/day')
     
    27822786c
    27832787      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    2784       CALL histwrite(nid_day,"snow_cov",itau_w,zx_tmp_2d,iim*jjmp1,
     2788      CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    27852789     .               ndex2d)
    2786 c     call writephy(fid_day,prof2d_av,'snow_cov',zxsnow,
     2790c     call writephy(fid_day,prof2d_av,'snow_mass',zxsnow,
    27872791c    .              'Snow cover','mm')
    27882792c
     
    30253029      ENDDO
    30263030      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    3027       CALL histwrite(nid_mth,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     3031      CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    30283032c
    30293033      DO i = 1, klon
     
    30433047c
    30443048      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    3045       CALL histwrite(nid_mth,"snow_cov",itau_w,zx_tmp_2d,iim*jjmp1,
     3049      CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    30463050     .               ndex2d)
    30473051c
     
    32513255     .                                   iim*jjmp1*klev,ndex3d)
    32523256c
     3257      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d)
     3258      CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d,
     3259     .                                   iim*jjmp1*klev,ndex3d)
     3260c
    32533261      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
    32543262      CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d,
     
    32913299     .                                   iim*jjmp1*klev,ndex3d)
    32923300c
     3301      zpt_conv = 0.
     3302      where (ptconv) zpt_conv = 1.
    32933303      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d)
    32943304      CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d,
     
    34263436      ENDDO
    34273437      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
    3428       CALL histwrite(nid_ins,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
     3438      CALL histwrite(nid_ins,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
    34293439c
    34303440      DO i = 1, klon
     
    35423552c
    35433553      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
    3544       CALL histwrite(nid_ins,"snow_cov",itau_w,zx_tmp_2d,iim*jjmp1,
     3554      CALL histwrite(nid_ins,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
    35453555     .               ndex2d)
    35463556c
Note: See TracChangeset for help on using the changeset viewer.