Changeset 4884
- Timestamp:
- Mar 28, 2024, 10:15:51 PM (10 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90
r4881 r4884 348 348 DO igrid=1,ngrid 349 349 tke(igrid,nlay+1)=0. 350 eps(igrid,nlay+1)=0. 350 351 END DO 351 352 … … 357 358 ustar=sqrt(cdrag_uv(igrid)*(wind_u(igrid,1)**2+wind_v(igrid,1)**2)) 358 359 tke(igrid,1)=ctkes*(ustar**2) 360 eps(igrid,1)=0. ! arbitrary as TKE is not properly defined at the surface 359 361 END DO 360 362 -
LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90
r4882 r4884 708 708 * ( 1329.31 + 0.0074615 * ( pplay(i) - 0.85e5 ) & 709 709 - 40.637 * ( temp(i) - 275. ) ) 710 air_thermal_conduct=(5.69+0.017*(temp(i)-RTT))*1.e-3*4.184 ! thermal conductivity of the air, SI 710 711 711 712 !--In clear air … … 715 716 nb_snowflake_clr = snowclr(i) / precipfracclr(i) / snow_fallspeed_clr & 716 717 / ( 4. / 3. * RPI * r_snow**3. * rho_snow ) 717 air_thermal_conduct=(5.69+0.017*(temp(i)-RTT))*1.e-3*4.184 ! thermal conductivity of the air, SI718 718 dqsclrmelt = - nb_snowflake_clr * 4. * RPI * air_thermal_conduct & 719 719 * capa_snowflake / RLMLT * coef_ventil & -
LMDZ6/trunk/libf/phylmd/yamada4.F90
r4881 r4884 688 688 689 689 ! Diagnostics 690 691 eps(:,:)= dissip(:,:)690 691 eps(:,:)=0. 692 692 wprime(1:ngrid,:,nsrf)=0. 693 693 DO k=2,klev 694 694 DO ig=1,ngrid 695 eps(ig,k)=dissip(ig,k) 695 696 jg=ni(ig) 696 697 wprime(jg,k,nsrf)=sqrt(MAX(1./3*q2(ig,k),0.))
Note: See TracChangeset
for help on using the changeset viewer.