Changeset 4013 for LMDZ6/branches/Ocean_skin/libf/phylmd/yamada4.F90
- Timestamp:
- Nov 19, 2021, 4:58:59 PM (3 years ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
- Property svn:mergeinfo changed
-
LMDZ6/branches/Ocean_skin/libf/phylmd/yamada4.F90
r3798 r4013 6 6 USE dimphy 7 7 USE ioipsl_getin_p_mod, ONLY : getin_p 8 USE phys_local_var_mod, only: tke_dissip 8 USE phys_local_var_mod, only: tke_dissip,wprime 9 9 10 10 IMPLICIT NONE … … 726 726 lyam(1:ngrid, 2:klev)*5.17*kn(1:ngrid, 2:klev)*n2(1:ngrid, 2:klev)/ & 727 727 sqrt(q2(1:ngrid,2:klev)) 728 728 729 729 t2yam(1:ngrid, 2:klev) = 9.1*kn(1:ngrid, 2:klev)* & 730 730 dtetadz(1:ngrid, 2:klev)**2/sqrt(q2(1:ngrid,2:klev))* & … … 750 750 751 751 !============================================================================ 752 ! Diagnostique de la dissipation 752 ! Diagnostique de la dissipation et vitesse verticale 753 753 !============================================================================ 754 754 755 755 ! Diagnostics 756 756 tke_dissip(1:ngrid,:,nsrf)=0. 757 ! DO k=2,klev 758 ! DO ig=1,ngrid 759 ! jg=ni(ig) 760 ! tke_dissip(jg,k,nsrf)=dissip(ig,k) 761 ! ENDDO 762 ! ENDDO 757 wprime(1:ngrid,:,nsrf)=0. 758 DO k=2,klev 759 DO ig=1,ngrid 760 jg=ni(ig) 761 wprime(jg,k,nsrf)=sqrt(MAX(1./3*q2(ig,k),0.)) 762 tke_dissip(jg,k,nsrf)=dissip(ig,k) 763 ENDDO 764 ENDDO 763 765 764 766 !=============================================================================
Note: See TracChangeset
for help on using the changeset viewer.