Changeset 4412
- Timestamp:
- Feb 1, 2023, 5:35:29 PM (22 months ago)
- Location:
- LMDZ6/trunk
- Files:
-
- 9 edited
-
DefLists/field_def_lmdz.xml (modified) (1 diff)
-
libf/phylmd/lscp_ini_mod.F90 (modified) (3 diffs)
-
libf/phylmd/lscp_mod.F90 (modified) (21 diffs)
-
libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
-
libf/phylmd/physiq_mod.F90 (modified) (13 diffs)
-
libf/phylmd/phytrac_mod.F90 (modified) (1 diff)
-
libf/phylmd/tracinca_mod.F90 (modified) (1 diff)
-
libf/phylmd/write_bilKP_ave.h (modified) (1 diff)
-
libf/phylmd/write_bilKP_ins.h (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/DefLists/field_def_lmdz.xml
r4248 r4412 583 583 <field id="concdust" long_name="Concentration of Dust " unit="kg/m3" /> 584 584 <field id="ec550aer" long_name="Extinction at 550nm" unit="m^-1" /> 585 <field id="lwcon" long_name="Cloud liquid water content " unit="kg/kg" />586 <field id="iwcon" long_name="Cloud ice water content " unit="kg/kg" />585 <field id="lwcon" long_name="Cloud liquid water content seen by radiation" unit="kg/kg" /> 586 <field id="iwcon" long_name="Cloud ice water content seen by radiation" unit="kg/kg" /> 587 587 <field id="temp" long_name="Air temperature" unit="K" /> 588 588 <field id="theta" long_name="Potential air temperature" unit="K" /> -
LMDZ6/trunk/libf/phylmd/lscp_ini_mod.F90
r4380 r4412 35 35 !$OMP THREADPRIVATE(iflag_mpc_bl) 36 36 37 LOGICAL, SAVE :: ok_rad liq_snow=.false. ! take into account the mass of ice precip in the cloud ice content seen by radiation38 !$OMP THREADPRIVATE(ok_rad liq_snow)37 LOGICAL, SAVE :: ok_radocond_snow=.false. ! take into account the mass of ice precip in the cloud ice content seen by radiation 38 !$OMP THREADPRIVATE(ok_radocond_snow) 39 39 40 40 … … 56 56 CALL getin_p('rain_int_min',rain_int_min) 57 57 CALL getin_p('iflag_mpc_bl',iflag_mpc_bl) 58 CALL getin_p('ok_rad liq_snow',ok_radliq_snow)58 CALL getin_p('ok_radocond_snow',ok_radocond_snow) 59 59 WRITE(lunout,*) 'lscp, ninter:', ninter 60 60 WRITE(lunout,*) 'lscp, iflag_evap_prec:', iflag_evap_prec … … 62 62 WRITE(lunout,*) 'lscp, rain_int_min:', rain_int_min 63 63 WRITE(lunout,*) 'lscp, iflag_mpc_bl:', iflag_mpc_bl 64 WRITE(lunout,*) 'lscp, ok_rad liq_snow:', ok_radliq_snow64 WRITE(lunout,*) 'lscp, ok_radocond_snow:', ok_radocond_snow 65 65 66 66 ! check for precipitation sub-time steps -
LMDZ6/trunk/libf/phylmd/lscp_mod.F90
r4397 r4412 9 9 paprs,pplay,t,q,ptconv,ratqs, & 10 10 d_t, d_q, d_ql, d_qi, rneb, rneblsvol, rneb_seri, & 11 rad liq, radicefrac, rain, snow, &11 radocond, radicefrac, rain, snow, & 12 12 frac_impa, frac_nucl, beta, & 13 13 prfl, psfl, rhcl, zqta, fraca, & … … 41 41 ! References: 42 42 ! 43 ! - Bony, S., & Emanuel, K. A. 2001, JAS, doi: 10.1175/1520-0469(2001)058<3158:APOTCA>2.0.CO;2 43 44 ! - Hourdin et al. 2013, Clim Dyn, doi:10.1007/s00382-012-1343-y 44 45 ! - Jam et al. 2013, Boundary-Layer Meteorol, doi:10.1007/s10546-012-9789-3 46 ! - Jouhaud, et al. 2018. JAMES, doi:10.1029/2018MS001379 45 47 ! - Madeleine et al. 2020, JAMES, doi:10.1029/2020MS002046 48 ! - Touzze-Peifert Ludo, PhD thesis, p117-124 46 49 ! ------------------------------------------------------------------------------- 47 50 ! Code structure: … … 50 53 ! P1> Evaporation of the precipitation (falling from the k+1 level) 51 54 ! P2> Cloud formation (at the k level) 52 ! P2.A.0> Cloud properties calculation from a rectangular pdf 53 ! P2.A.1> With the new PDFs, calculation of cloud properties using the inital 55 ! P2.A.1> With the PDFs, calculation of cloud properties using the inital 54 56 ! values of T and Q 55 57 ! P2.A.2> Coupling between condensed water and temperature 56 58 ! P2.A.3> Calculation of final quantities associated with cloud formation 57 ! P2.B> 'All or nothing' cloud 58 ! P2.C> Release of Latent heat after cloud formation 59 ! P2.B> Release of Latent heat after cloud formation 59 60 ! P3> Autoconversion to precipitation (k-level) 60 61 ! P4> Wet scavenging … … 65 66 ! water used in the physics and the dynamics 66 67 ! 67 ! During the autoconversion to precipitation (P3 step), rad liq(cloud water used68 ! During the autoconversion to precipitation (P3 step), radocond (cloud water used 68 69 ! by the radiation scheme) is calculated as an average of the water that remains 69 70 ! in the cloud during the precipitation and not the water remaining at the end … … 75 76 ! Radiation: 76 77 ! xflwc(newmicro)+xfiwc(newmicro) = 77 ! cldliq(physiq)=radliq(fisrt)=lwcon(nc)+iwcon(nc)78 ! radocond=lwcon(nc)+iwcon(nc) 78 79 ! 79 80 ! Notetheless, be aware of: 80 81 ! 81 ! rad liq(fisrt).NE. ocond(nc)82 ! radocond .NE. ocond(nc) 82 83 ! i.e.: 83 84 ! lwcon(nc)+iwcon(nc) .NE. ocond(nc) 85 ! but oliq+(ocond-oliq) .EQ. ocond 84 86 ! (which is not trivial) 85 87 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ … … 91 93 92 94 USE lscp_ini_mod, ONLY : seuil_neb, ninter, iflag_evap_prec, t_coup, DDT0, ztfondue, rain_int_min 93 USE lscp_ini_mod, ONLY : iflag_mpc_bl, ok_rad liq_snow, a_tr_sca95 USE lscp_ini_mod, ONLY : iflag_mpc_bl, ok_radocond_snow, a_tr_sca 94 96 95 97 … … 115 117 REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay ! mid-layer pressure [Pa] 116 118 REAL, DIMENSION(klon,klev), INTENT(IN) :: t ! temperature (K) 117 REAL, DIMENSION(klon,klev), INTENT(IN) :: q ! specific humidity[kg/kg]119 REAL, DIMENSION(klon,klev), INTENT(IN) :: q ! total specific humidity (= vapor phase) [kg/kg] 118 120 INTEGER, INTENT(IN) :: iflag_cld_th ! flag that determines the distribution of convective clouds 119 121 INTEGER, INTENT(IN) :: iflag_ice_thermo! flag to activate the ice thermodynamics … … 152 154 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rneb ! cloud fraction [-] 153 155 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rneblsvol ! cloud fraction per unit volume [-] 154 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rad liq! condensed water used in the radiation scheme [kg/kg]156 REAL, DIMENSION(klon,klev), INTENT(OUT) :: radocond ! condensed water used in the radiation scheme [kg/kg] 155 157 REAL, DIMENSION(klon,klev), INTENT(OUT) :: radicefrac ! ice fraction of condensed water for radiation scheme 156 158 REAL, DIMENSION(klon,klev), INTENT(OUT) :: rhcl ! clear-sky relative humidity [-] 157 REAL, DIMENSION(klon), INTENT(OUT) :: rain ! large-scale rainfall [kg/s/m2]158 REAL, DIMENSION(klon), INTENT(OUT) :: snow ! large-scale snowfall [kg/s/m2]159 REAL, DIMENSION(klon), INTENT(OUT) :: rain ! surface large-scale rainfall [kg/s/m2] 160 REAL, DIMENSION(klon), INTENT(OUT) :: snow ! surface large-scale snowfall [kg/s/m2] 159 161 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qsatl ! saturation specific humidity wrt liquid [kg/kg] 160 162 REAL, DIMENSION(klon,klev), INTENT(OUT) :: qsats ! saturation specific humidity wrt ice [kg/kg] … … 231 233 REAL velo(klon,klev), vr 232 234 REAL qlmpc, qimpc, rnebmpc 233 REAL rad liqi(klon,klev), radliql(klon,klev)235 REAL radocondi(klon,klev), radocondl(klon,klev) 234 236 235 237 INTEGER i, k, n, kk … … 295 297 d_qi(:,:) = 0.0 296 298 rneb(:,:) = 0.0 297 rad liq(:,:) = 0.0299 radocond(:,:) = 0.0 298 300 radicefrac(:,:) = 0.0 299 301 frac_nucl(:,:) = 1.0 … … 350 352 ! P0> Thermalization of precipitation falling from the overlying layer 351 353 ! -------------------------------------------------------------------- 352 ! Computes air temperature variation due to latent heattransported by354 ! Computes air temperature variation due to enthalpy transported by 353 355 ! precipitation. Precipitation is then thermalized with the air in the 354 356 ! layer. 355 357 ! The precipitation should remain thermalized throughout the different 356 ! thermodynamical transformations. The corresponding water mass should 358 ! thermodynamical transformations. 359 ! The corresponding water mass should 357 360 ! be added when calculating the layer's enthalpy change with 358 361 ! temperature 362 ! See lmdzpedia page todoan 363 ! todoan: check consistency with ice phase 364 ! todoan: understand why several steps 359 365 ! --------------------------------------------------------------------- 360 366 … … 409 415 IF (zrfl(i)+zifl(i).GT.0.) THEN 410 416 411 ! L TP: we only account for precipitation evaporation in the clear-sky (iflag_evap_prec=4).417 ! LudoTP: we only account for precipitation evaporation in the clear-sky (iflag_evap_prec=4). 412 418 ! c_iso: likely important to distinguish cs from neb isotope precipitation 413 419 … … 475 481 ! redistribute zqev0 conserving the ratio liquid/ice 476 482 483 ! todoan: check the consistency of this formula 477 484 IF (zqevt+zqevti.GT.zqev0) THEN 478 485 zqev=zqev0*zqevt/(zqevt+zqevti) … … 494 501 zq(i) = zq(i) - (zrfln(i)+zifln(i)-zrfl(i)-zifl(i)) & 495 502 * (RG/(paprs(i,k)-paprs(i,k+1)))*dtime 503 496 504 ! precip thermalization 497 505 zmqc(i) = zmqc(i) + (zrfln(i)+zifln(i)-zrfl(i)-zifl(i)) & … … 525 533 526 534 ! Melting: 527 zmelt = ((zt(i)- 273.15)/(ztfondue-273.15)) ! JYG535 zmelt = ((zt(i)-RTT)/(ztfondue-RTT)) ! JYG 528 536 ! precip fraction that is melted 529 537 zmelt = MIN(MAX(zmelt,0.),1.) 530 538 531 ! Icemelting539 ! update of rainfall and snowfall due to melting 532 540 IF (iflag_evap_prec.EQ.4) THEN 533 541 zrflclr(i)=zrflclr(i)+zmelt*ziflclr(i) 534 542 zrflcld(i)=zrflcld(i)+zmelt*ziflcld(i) 535 543 zrfl(i)=zrflclr(i)+zrflcld(i) 544 545 ziflclr(i)=ziflclr(i)*(1.-zmelt) 546 ziflcld(i)=ziflcld(i)*(1.-zmelt) 547 zifl(i)=ziflclr(i)+ziflcld(i) 548 536 549 ELSE 537 550 zrfl(i)=zrfl(i)+zmelt*zifl(i) 538 ENDIF 551 552 zifl(i)=zifl(i)*(1.-zmelt) 553 ENDIF 554 555 539 556 ! c_iso: melting of isotopic precipi with zmelt (no fractionation) 540 557 … … 543 560 *RLMLT/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i))) 544 561 545 IF (iflag_evap_prec.EQ.4) THEN546 ziflclr(i)=ziflclr(i)*(1.-zmelt)547 ziflcld(i)=ziflcld(i)*(1.-zmelt)548 zifl(i)=ziflclr(i)+ziflcld(i)549 ELSE550 zifl(i)=zifl(i)*(1.-zmelt)551 ENDIF552 !c_iso: same for isotopic precip.553 562 554 563 ELSE … … 908 917 ENDIF 909 918 910 ! Initialisation of zoliq and rad liqvariables919 ! Initialisation of zoliq and radocond variables 911 920 912 921 DO i = 1, klon … … 919 928 iwc(i) = 0. 920 929 zneb(i) = MAX(rneb(i,k),seuil_neb) 921 rad liq(i,k) = zoliq(i)/REAL(ninter+1)930 radocond(i,k) = zoliq(i)/REAL(ninter+1) 922 931 radicefrac(i,k) = zfice(i) 923 rad liqi(i,k)=zoliq(i)*zfice(i)/REAL(ninter+1)924 rad liql(i,k)=zoliq(i)*(1.-zfice(i))/REAL(ninter+1)932 radocondi(i,k)=zoliq(i)*zfice(i)/REAL(ninter+1) 933 radocondl(i,k)=zoliq(i)*(1.-zfice(i))/REAL(ninter+1) 925 934 ENDDO 926 935 … … 992 1001 ! c_iso: call isotope_conversion (for readibility) or duplicate 993 1002 994 rad liq(i,k) = radliq(i,k) + zoliq(i)/REAL(ninter+1)995 rad liql(i,k) = radliql(i,k) + zoliql(i)/REAL(ninter+1)996 rad liqi(i,k) = radliqi(i,k) + zoliqi(i)/REAL(ninter+1)1003 radocond(i,k) = radocond(i,k) + zoliq(i)/REAL(ninter+1) 1004 radocondl(i,k) = radocondl(i,k) + zoliql(i)/REAL(ninter+1) 1005 radocondi(i,k) = radocondi(i,k) + zoliqi(i)/REAL(ninter+1) 997 1006 998 1007 ENDIF ! rneb >0 … … 1109 1118 1110 1119 ! Calculation of the concentration of condensates seen by the radiation scheme 1111 ! for liquid phase, we take rad liql1112 ! for ice phase, we take rad liqi if we neglect snowfall, otherwise (ok_radliq_snow=true)1113 ! we recaulate rad liqi to account for contributions from the precipitation flux1114 1115 IF ((ok_rad liq_snow) .AND. (k .LT. klev)) THEN1120 ! for liquid phase, we take radocondl 1121 ! for ice phase, we take radocondi if we neglect snowfall, otherwise (ok_radocond_snow=true) 1122 ! we recaulate radocondi to account for contributions from the precipitation flux 1123 1124 IF ((ok_radocond_snow) .AND. (k .LT. klev)) THEN 1116 1125 ! for the solid phase (crystals + snowflakes) 1117 ! we recalculate rad liqi by summing1126 ! we recalculate radocondi by summing 1118 1127 ! the ice content calculated in the mesh 1119 1128 ! + the contribution of the non-evaporated snowfall … … 1121 1130 DO i=1,klon 1122 1131 IF (ziflprev(i) .NE. 0.0) THEN 1123 rad liqi(i,k)=zoliq(i)*zfice(i)+zqpreci(i)+ziflprev(i)/zrho(i,k+1)/velo(i,k+1)1132 radocondi(i,k)=zoliq(i)*zfice(i)+zqpreci(i)+ziflprev(i)/zrho(i,k+1)/velo(i,k+1) 1124 1133 ELSE 1125 rad liqi(i,k)=zoliq(i)*zfice(i)+zqpreci(i)1134 radocondi(i,k)=zoliq(i)*zfice(i)+zqpreci(i) 1126 1135 ENDIF 1127 rad liq(i,k)=radliql(i,k)+radliqi(i,k)1136 radocond(i,k)=radocondl(i,k)+radocondi(i,k) 1128 1137 ENDDO 1129 1138 ENDIF 1130 1139 1131 ! caculate the percentage of ice in "rad liq" so cloud+precip seen by the radiation scheme1140 ! caculate the percentage of ice in "radocond" so cloud+precip seen by the radiation scheme 1132 1141 DO i=1,klon 1133 IF (rad liq(i,k) .GT. 0.) THEN1134 radicefrac(i,k)=MIN(MAX(rad liqi(i,k)/radliq(i,k),0.),1.)1142 IF (radocond(i,k) .GT. 0.) THEN 1143 radicefrac(i,k)=MIN(MAX(radocondi(i,k)/radocond(i,k),0.),1.) 1135 1144 ENDIF 1136 1145 ENDDO -
LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90
r4059 r4412 1413 1413 'ec550aer', 'Extinction at 550nm', 'm^-1', (/ ('', i=1, 10) /)) 1414 1414 TYPE(ctrl_out), SAVE :: o_lwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1415 'lwcon', 'Cloud liquid water content ', 'kg/kg', (/ ('', i=1, 10) /))1415 'lwcon', 'Cloud liquid water content seen by radiation', 'kg/kg', (/ ('', i=1, 10) /)) 1416 1416 TYPE(ctrl_out), SAVE :: o_iwcon = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), & 1417 'iwcon', 'Cloud ice water content ', 'kg/kg', (/ ('', i=1, 10) /))1417 'iwcon', 'Cloud ice water content seen by radiation', 'kg/kg', (/ ('', i=1, 10) /)) 1418 1418 TYPE(ctrl_out), SAVE :: o_temp = ctrl_out((/ 2, 3, 4, 10, 10, 10, 11, 11, 11, 11/), & 1419 1419 'temp', 'Air temperature', 'K', (/ ('', i=1, 10) /)) -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4389 r4412 827 827 REAL dialiq(klon,klev) ! eau liquide nuageuse 828 828 REAL diafra(klon,klev) ! fraction nuageuse 829 REAL cldliq(klon,klev) ! eau liquide nuageuse829 REAL radocond(klon,klev) ! eau condensee nuageuse 830 830 ! 831 831 !XXX PB … … 3637 3637 t_seri, q_seri,ptconv,ratqs, & 3638 3638 d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, & 3639 cldliq, picefra, rain_lsc, snow_lsc, &3639 radocond, picefra, rain_lsc, snow_lsc, & 3640 3640 frac_impa, frac_nucl, beta_prec_fisrt, & 3641 3641 prfl, psfl, rhcl, & … … 3649 3649 CALL fisrtilp(phys_tstep,paprs,pplay, & 3650 3650 t_seri, q_seri,ptconv,ratqs, & 3651 d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, cldliq, &3651 d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, radocond, & 3652 3652 rain_lsc, snow_lsc, & 3653 3653 pfrac_impa, pfrac_nucl, pfrac_1nucl, & … … 3696 3696 cldfra(i,k) = rneb(i,k) 3697 3697 !CR: a quoi ca sert? Faut-il ajouter qs_seri? 3698 IF (.NOT.new_oliq) cldliq(i,k) = ql_seri(i,k)3698 IF (.NOT.new_oliq) radocond(i,k) = ql_seri(i,k) 3699 3699 ENDDO 3700 3700 ENDDO … … 3760 3760 DO i = 1, klon 3761 3761 IF (diafra(i,k).GT.cldfra(i,k)) THEN 3762 cldliq(i,k) = dialiq(i,k)3762 radocond(i,k) = dialiq(i,k) 3763 3763 cldfra(i,k) = diafra(i,k) 3764 3764 ENDIF … … 3797 3797 DO i=1,klon 3798 3798 IF (ptconv(i,k).AND.ptconvth(i,k)) THEN 3799 cldliq(i,k)=cldliq(i,k)+rnebcon(i,k)*clwcon(i,k)3799 radocond(i,k)=radocond(i,k)+rnebcon(i,k)*clwcon(i,k) 3800 3800 cldfra(i,k)=min(cldfra(i,k)+rnebcon(i,k),1.) 3801 3801 ELSE IF (ptconv(i,k)) THEN 3802 3802 cldfra(i,k)=rnebcon(i,k) 3803 cldliq(i,k)=rnebcon(i,k)*clwcon(i,k)3803 radocond(i,k)=rnebcon(i,k)*clwcon(i,k) 3804 3804 ENDIF 3805 3805 ENDDO … … 3810 3810 DO i=1,klon 3811 3811 cldfra(i,k)=min(cldfra(i,k)+rnebcon(i,k),1.) 3812 cldliq(i,k)=cldliq(i,k)+rnebcon(i,k)*clwcon(i,k)3812 radocond(i,k)=radocond(i,k)+rnebcon(i,k)*clwcon(i,k) 3813 3813 ENDDO 3814 3814 ENDDO … … 3828 3828 IF (ptconv(i,k).AND. .NOT.ptconvth(i,k)) THEN 3829 3829 cldfra(i,k)=rnebcon(i,k) 3830 cldliq(i,k)=rnebcon(i,k)*clwcon(i,k)3830 radocond(i,k)=rnebcon(i,k)*clwcon(i,k) 3831 3831 ENDIF 3832 3832 ENDDO … … 3839 3839 ! Ancienne version 3840 3840 cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.) 3841 cldliq(:,:)=cldliq(:,:)+rnebcon(:,:)*clwcon(:,:)3841 radocond(:,:)=radocond(:,:)+rnebcon(:,:)*clwcon(:,:) 3842 3842 ENDIF 3843 3843 … … 3859 3859 DO i = 1, klon 3860 3860 IF (diafra(i,k).GT.cldfra(i,k)) THEN 3861 cldliq(i,k) = dialiq(i,k)3861 radocond(i,k) = dialiq(i,k) 3862 3862 cldfra(i,k) = diafra(i,k) 3863 3863 ENDIF … … 4168 4168 ENDIF 4169 4169 CALL newmicro (flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, & 4170 paprs, pplay, t_seri, cldliq, picefra, cldfra, &4170 paprs, pplay, t_seri, radocond, picefra, cldfra, & 4171 4171 cldtau, cldemi, cldh, cldl, cldm, cldt, cldq, & 4172 4172 flwp, fiwp, flwc, fiwc, & … … 4176 4176 ELSE 4177 4177 CALL nuage (paprs, pplay, & 4178 t_seri, cldliq, picefra, cldfra, cldtau, cldemi, &4178 t_seri, radocond, picefra, cldfra, cldtau, cldemi, & 4179 4179 cldh, cldl, cldm, cldt, cldq, & 4180 4180 ok_aie, & … … 5003 5003 presnivs, pphis, pphi, albsol1, & 5004 5004 sh_in, ch_in, rhcl, cldfra, rneb, & 5005 diafra, cldliq, itop_con, ibas_con, &5005 diafra, radocond, itop_con, ibas_con, & 5006 5006 pmflxr, pmflxs, prfl, psfl, & 5007 5007 da, phi, mp, upwd, & -
LMDZ6/trunk/libf/phylmd/phytrac_mod.F90
r4389 r4412 184 184 REAL,DIMENSION(klon),INTENT(IN) :: pphis 185 185 REAL,DIMENSION(klev),INTENT(IN) :: presnivs 186 REAL,DIMENSION(klon,klev),INTENT(IN) :: cldliq ! eau liquide nuageuse186 REAL,DIMENSION(klon,klev),INTENT(IN) :: cldliq ! eau condensee totale 187 187 REAL,DIMENSION(klon,klev),INTENT(IN) :: cldfra ! fraction nuageuse (tous les nuages) 188 188 REAL,DIMENSION(klon,klev),INTENT(IN) :: diafra ! fraction nuageuse (convection ou stratus artificiels) -
LMDZ6/trunk/libf/phylmd/tracinca_mod.F90
r4389 r4412 89 89 REAL,DIMENSION(klon,klev),INTENT(IN) :: pphi ! geopotentiel 90 90 REAL,DIMENSION(klon),INTENT(IN) :: pphis 91 REAL,DIMENSION(klon,klev),INTENT(IN) :: cldliq ! eau liquide nuageuse91 REAL,DIMENSION(klon,klev),INTENT(IN) :: cldliq ! eau condensee pour le radiatif 92 92 REAL,DIMENSION(klon,klev),INTENT(IN) :: cldfra ! fraction nuageuse (tous les nuages) 93 93 REAL,DIMENSION(klon,klev),INTENT(IN) :: diafra ! fraction nuageuse (convection ou stratus artificiels) -
LMDZ6/trunk/libf/phylmd/write_bilKP_ave.h
r3102 r4412 49 49 CALL histwrite_phy(nid_bilKPave,"play",itau_w,paprs) 50 50 c 51 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)51 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, radocond, zx_tmp_3d) 52 52 CALL histwrite_phy(nid_bilKPave,"oliq",itau_w,cldliq) 53 53 c -
LMDZ6/trunk/libf/phylmd/write_bilKP_ins.h
r3102 r4412 49 49 CALL histwrite_phy(nid_bilKPins,"play",itau_w,paprs) 50 50 c 51 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d)51 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, radocond, zx_tmp_3d) 52 52 CALL histwrite_phy(nid_bilKPins,"oliq",itau_w,cldliq) 53 53 c
Note: See TracChangeset
for help on using the changeset viewer.
