Changeset 3356 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/tend_to_tke.F90
- Timestamp:
- Jun 29, 2018, 12:31:11 PM (6 years ago)
- Location:
- LMDZ6/branches/DYNAMICO-conv
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv
- Property svn:mergeinfo changed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/tend_to_tke.F90
r2897 r3356 32 32 !************************************************************************************** 33 33 34 SUBROUTINE tend_to_tke(dt,plev,exner,temp,windu,windv,dt_a,du_a,dv_a, tke)34 SUBROUTINE tend_to_tke(dt,plev,exner,temp,windu,windv,dt_a,du_a,dv_a,pctsrf,tke) 35 35 36 36 USE dimphy, ONLY: klon, klev 37 37 USE indice_sol_mod, ONLY: nbsrf 38 39 IMPLICIT NONE 38 40 #include "YOMCST.h" 39 41 … … 53 55 REAL du_a(klon,klev) ! Zonal wind speed tendency [m/s], grid-cell average or for a one subsurface 54 56 REAL dv_a(klon,klev) ! Meridional wind speed tendency [m/s], grid-cell average or for a one subsurface 57 REAL pctsrf(klon,nbsrf+1) ! Turbulent Kinetic energy [m2/s2], grid-cell average or for a subsurface 55 58 56 59 ! Inputs/Outputs 57 60 !--------------- 58 REAL tke(klon,klev ,nbsrf) ! Turbulent Kinetic energy [m2/s2], grid-cell average or for a subsurface61 REAL tke(klon,klev+1,nbsrf+1) ! Turbulent Kinetic energy [m2/s2], grid-cell average or for a subsurface 59 62 60 63 … … 63 66 64 67 65 INTEGER i g,k,isrf ! indices68 INTEGER i,k,isrf ! indices 66 69 REAL masse(klon,klev) ! mass in the layers [kg/m2] 67 70 REAL unsmasse(klon,klev+1) ! linear mass in the layers [kg/m2] … … 121 124 DO isrf=1,nbsrf 122 125 DO k=1,klev 123 tke(:,k,isrf)= tke(:,k,isrf)+tendu(:,k)+tendv(:,k)+tendt(:,k) 124 tke(:,k,isrf)= max(tke(:,k,isrf),1.e-10) 126 DO i=1,klon 127 IF (pctsrf(i,isrf)>0.) THEN 128 tke(i,k,isrf)= tke(i,k,isrf)+tendu(i,k)+tendv(i,k)+tendt(i,k) 129 tke(i,k,isrf)= max(tke(i,k,isrf),1.e-10) 130 ENDIF 131 ENDDO 125 132 ENDDO 126 133 ENDDO 127 128 ! dtke_t(:,:)=tendt(:,:)129 ! dtke_u(:,:)=tendu(:,:)130 ! dtke_v(:,:)=tendv(:,:)131 134 132 135
Note: See TracChangeset
for help on using the changeset viewer.