Changeset 2904 for LMDZ5/trunk
- Timestamp:
- Jun 6, 2017, 4:31:18 PM (8 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd/dyn1d
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/dyn1d/1DUTILS.h
r2857 r2904 5025 5025 REAL tau 5026 5026 !c DATA tau /3600./ 5027 DATA tau /5400./ 5027 ! DATA tau /5400./ 5028 DATA tau /43200./ 5028 5029 ! 5029 5030 INTEGER k,i … … 5037 5038 DO k = 1,klev 5038 5039 DO i = 1,klon 5039 IF (paprs(i,1)-pplay(i,k) .GT. 10000.) THEN 5040 !CR: nudging everywhere 5041 ! IF (paprs(i,1)-pplay(i,k) .GT. 10000.) THEN 5040 5042 ! 5041 5043 d_u(i,k) = d_u(i,k) + 1./tau*(u_targ(i,k)-u(i,k)) … … 5044 5046 print *,' k,u,d_u,v,d_v ', & 5045 5047 k,u(i,k),d_u(i,k),v(i,k),d_v(i,k) 5046 ENDIF5048 ! ENDIF 5047 5049 ! 5048 5050 ENDDO -
LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90
r2899 r2904 1077 1077 endif 1078 1078 1079 !CRio:Attention modif spécifique cas de Caroline 1080 if (forcing_les) then 1081 fcoriolis=0. 1082 !Nudging 1083 1084 !on calcule dt_cooling 1085 do l=1,llm 1086 if (play(l).ge.20000.) then 1087 dt_cooling(l)=-1.5/86400. 1088 elseif ((play(l).ge.10000.).and.((play(l).lt.20000.))) then 1089 dt_cooling(l)=-1.5/86400.*(play(l)-10000.)/(10000.)-1./86400.*(20000.-play(l))/10000.*(temp(l)-200.) 1090 else 1091 dt_cooling(l)=-1.*(temp(l)-200.)/86400. 1092 endif 1093 enddo 1094 1095 endif 1096 !RC 1097 1079 1098 IF (prt_level >= 5) print*, 'fcoriolis, xlat,mxcalc ', & 1080 1099 fcoriolis, xlat,mxcalc
Note: See TracChangeset
for help on using the changeset viewer.