Changeset 4009 for LMDZ6/trunk/libf/phylmdiso
- Timestamp:
- Nov 17, 2021, 11:23:39 PM (4 years ago)
- Location:
- LMDZ6/trunk/libf/phylmdiso
- Files:
-
- 1 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90
r4003 r4009 16 16 REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:) 17 17 !$OMP THREADPRIVATE(u_seri, v_seri) 18 REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:), l_mix(:,:,:), tke_dissip(:,:,:) 19 !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip )18 REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:), l_mix(:,:,:), tke_dissip(:,:,:), wprime(:,:,:) 19 !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip, wprime) 20 20 REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:) 21 21 !$OMP THREADPRIVATE(tr_seri) … … 693 693 ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev)) 694 694 ALLOCATE(u_seri(klon,klev),v_seri(klon,klev)) 695 ALLOCATE(l_mixmin(klon,klev+1,nbsrf), l_mix(klon,klev+1,nbsrf), tke_dissip(klon,klev+1,nbsrf) )696 l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ; tke_dissip(:,:,:)=0. ! doit etre initialse car pas toujours remplis695 ALLOCATE(l_mixmin(klon,klev+1,nbsrf), l_mix(klon,klev+1,nbsrf), tke_dissip(klon,klev+1,nbsrf), wprime(klon,klev+1,nbsrf)) 696 l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ; tke_dissip(:,:,:)=0. ; wprime(:,:,:)=0. ! doit etre initialse car pas toujours remplis 697 697 698 698 ALLOCATE(tr_seri(klon,klev,nbtr)) … … 1090 1090 DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri) 1091 1091 DEALLOCATE(u_seri,v_seri) 1092 DEALLOCATE(l_mixmin,l_mix, tke_dissip )1092 DEALLOCATE(l_mixmin,l_mix, tke_dissip, wprime) 1093 1093 1094 1094 DEALLOCATE(tr_seri) -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r4003 r4009 264 264 d_t_vdf_x, d_t_vdf_w, & 265 265 d_q_vdf_x, d_q_vdf_w, & 266 pbl_tke_input, &266 pbl_tke_input, tke_dissip, l_mix, wprime, & 267 267 t_therm, q_therm, u_therm, v_therm, & 268 268 cdragh_x, cdragh_w, & … … 1288 1288 integer iostat 1289 1289 1290 REAL, dimension(klon,klev+1) :: tke_dissip_ave, l_mix_ave, wprime_ave 1290 1291 REAL zzz 1291 1292 !albedo SB >>> … … 4468 4469 ! Computation of ratqs, the width (normalized) of the subrid scale 4469 4470 ! water distribution 4470 CALL calcratqs(klon,klev,prt_level,lunout, & 4471 4472 tke_dissip_ave(:,:)=0. 4473 l_mix_ave(:,:)=0. 4474 wprime_ave(:,:)=0. 4475 4476 DO nsrf = 1, nbsrf 4477 DO i = 1, klon 4478 tke_dissip_ave(i,:) = tke_dissip_ave(i,:) + tke_dissip(i,:,nsrf)*pctsrf(i,nsrf) 4479 l_mix_ave(i,:) = l_mix_ave(i,:) + l_mix(i,:,nsrf)*pctsrf(i,nsrf) 4480 wprime_ave(i,:) = wprime_ave(i,:) + wprime(i,:,nsrf)*pctsrf(i,nsrf) 4481 ENDDO 4482 ENDDO 4483 4484 CALL calcratqs(klon,klev,prt_level,lunout, & 4471 4485 iflag_ratqs,iflag_con,iflag_cld_th,pdtphys, & 4472 4486 ratqsbas,ratqshaut,ratqsp0, ratqsdp, & 4473 4487 tau_ratqs,fact_cldcon,wake_s, wake_deltaq, & 4474 4488 ptconv,ptconvth,clwcon0th, rnebcon0th, & 4475 paprs,pplay,q_seri,zqsat,fm_therm, & 4489 paprs,pplay,t_seri,q_seri, qtc_cv, sigt_cv, zqsat, & 4490 pbl_tke(:,:,is_ave),tke_dissip_ave,l_mix_ave,wprime_ave,t2m,q2m,fm_therm, & 4476 4491 ratqs,ratqsc,ratqs_inter) 4492 4477 4493 4478 4494 !
Note: See TracChangeset
for help on using the changeset viewer.