Ignore:
Timestamp:
Aug 2, 2024, 2:12:03 PM (3 months ago)
Author:
abarral
Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dtke.F90

    r5119 r5158  
    4040
    4141!   calcul du detrainement
    42       do k=1,nlay
     42      DO k=1,nlay
    4343         detr0(:,k)=fm0(:,k)-fm0(:,k+1)+entr0(:,k)
    4444         masse0(:,k)=(pplev(:,k)-pplev(:,k+1))/RG
     
    5151      detr(:,1)=0.5*detr0(:,1)
    5252      fm(:,1)=0.
    53       do k=1,nlay-1
     53      DO k=1,nlay-1
    5454         masse(:,k+1)=0.5*(masse0(:,k)+masse0(:,k+1))
    5555         entr(:,k+1)=0.5*(entr0(:,k)+entr0(:,k+1))
     
    6565   q(:,:)=tke(:,:,isrf)
    6666!   calcul de la valeur dans les ascendances
    67       do ig=1,ngrid
     67      DO ig=1,ngrid
    6868         qa(ig,1)=q(ig,1)
    6969      enddo
     
    7171
    7272    IF (1==1) THEN
    73       do k=2,nlay
    74          do ig=1,ngrid
     73      DO k=2,nlay
     74         DO ig=1,ngrid
    7575            IF ((fm(ig,k+1)+detr(ig,k))*ptimestep>  &
    7676           1.e-5*masse(ig,k)) THEN
     
    9191! Calcul du flux subsident
    9292
    93       do k=2,nlay
    94          do ig=1,ngrid
     93      DO k=2,nlay
     94         DO ig=1,ngrid
    9595            wqd(ig,k)=fm(ig,k)*q(ig,k)
    9696            IF (wqd(ig,k)<0.) THEN
     
    9999         enddo
    100100      enddo
    101       do ig=1,ngrid
     101      DO ig=1,ngrid
    102102         wqd(ig,1)=0.
    103103         wqd(ig,nlay+1)=0.
     
    106106
    107107! Calcul des tendances
    108       do k=1,nlay
    109          do ig=1,ngrid
     108      DO k=1,nlay
     109         DO ig=1,ngrid
    110110            q(ig,k)=q(ig,k)+(detr(ig,k)*qa(ig,k)-entr(ig,k)*q(ig,k)  &
    111111                 -wqd(ig,k)+wqd(ig,k+1))  &
Note: See TracChangeset for help on using the changeset viewer.