- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dtke.F90
r5119 r5158 40 40 41 41 ! calcul du detrainement 42 dok=1,nlay42 DO k=1,nlay 43 43 detr0(:,k)=fm0(:,k)-fm0(:,k+1)+entr0(:,k) 44 44 masse0(:,k)=(pplev(:,k)-pplev(:,k+1))/RG … … 51 51 detr(:,1)=0.5*detr0(:,1) 52 52 fm(:,1)=0. 53 dok=1,nlay-153 DO k=1,nlay-1 54 54 masse(:,k+1)=0.5*(masse0(:,k)+masse0(:,k+1)) 55 55 entr(:,k+1)=0.5*(entr0(:,k)+entr0(:,k+1)) … … 65 65 q(:,:)=tke(:,:,isrf) 66 66 ! calcul de la valeur dans les ascendances 67 doig=1,ngrid67 DO ig=1,ngrid 68 68 qa(ig,1)=q(ig,1) 69 69 enddo … … 71 71 72 72 IF (1==1) THEN 73 dok=2,nlay74 doig=1,ngrid73 DO k=2,nlay 74 DO ig=1,ngrid 75 75 IF ((fm(ig,k+1)+detr(ig,k))*ptimestep> & 76 76 1.e-5*masse(ig,k)) THEN … … 91 91 ! Calcul du flux subsident 92 92 93 dok=2,nlay94 doig=1,ngrid93 DO k=2,nlay 94 DO ig=1,ngrid 95 95 wqd(ig,k)=fm(ig,k)*q(ig,k) 96 96 IF (wqd(ig,k)<0.) THEN … … 99 99 enddo 100 100 enddo 101 doig=1,ngrid101 DO ig=1,ngrid 102 102 wqd(ig,1)=0. 103 103 wqd(ig,nlay+1)=0. … … 106 106 107 107 ! Calcul des tendances 108 dok=1,nlay109 doig=1,ngrid108 DO k=1,nlay 109 DO ig=1,ngrid 110 110 q(ig,k)=q(ig,k)+(detr(ig,k)*qa(ig,k)-entr(ig,k)*q(ig,k) & 111 111 -wqd(ig,k)+wqd(ig,k+1)) &
Note: See TracChangeset
for help on using the changeset viewer.