Ignore:
Timestamp:
Jul 24, 2024, 4:23:34 PM (4 months ago)
Author:
abarral
Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90

    r5116 r5117  
    206206  ! @       do 200 k=2,nlp
    207207  ! @         do 190 i=1,len
    208   ! @          if((hm(i,k).lt.work(i)).and.
     208  ! @          if((hm(i,k).lt.work(i)).AND.
    209209  ! @      &      (hm(i,k).lt.hm(i,k-1)))THEN
    210210  ! @            work(i)=hm(i,k)
    … …  
    231231  ! @       do 240 k=minorig+1,nl
    232232  ! @         do 230 i=1,len
    233   ! @          if((hm(i,k).gt.work(i)).and.(k.le.ihmin(i)))THEN
     233  ! @          if((hm(i,k).gt.work(i)).AND.(k.le.ihmin(i)))THEN
    234234  ! @            work(i)=hm(i,k)
    235235  ! @            nk(i)=k
    … …  
    251251  ! -------------------------------------------------------------------
    252252  DO i = 1, len
    253     IF (((t(i, nk(i))<250.0) .OR. (q(i, nk(i))<=0.0)) & ! @      &       .or.(
     253    IF (((t(i, nk(i))<250.0) .OR. (q(i, nk(i))<=0.0)) & ! @      &       .OR.(
    254254            ! p(i,ihmin(i)).lt.400.0
    255255            ! )  )
    … …  
    296296  ! @      do 290 k=minorig,nl
    297297  ! @        do 280 i=1,len
    298   ! @          if((k.ge.(nk(i)+1)).and.(p(i,k).lt.plcl(i)))
     298  ! @          if((k.ge.(nk(i)+1)).AND.(p(i,k).lt.plcl(i)))
    299299  ! @     &    icb(i)=min(icb(i),k)
    300300  ! @ 280    continue
    … …  
    302302  ! @c
    303303  ! @      do 300 i=1,len
    304   ! @        if((icb(i).ge.nlm).and.(iflag(i).eq.0))iflag(i)=9
     304  ! @        if((icb(i).ge.nlm).AND.(iflag(i).EQ.0))iflag(i)=9
    305305  ! @ 300  continue
    306306
    … …  
    319319
    320320  DO i = 1, len
    321     ! @        if((icb(i).ge.nlm).and.(iflag(i).eq.0))iflag(i)=9
     321    ! @        if((icb(i).ge.nlm).AND.(iflag(i).EQ.0))iflag(i)=9
    322322    IF ((icb(i)==nlm) .AND. (iflag(i)==0)) iflag(i) = 9
    323323  END DO
    … …  
    683683  ! oct3       ath  = th(i,icb(i)-1) - dttrig
    684684  ! oct3
    685   ! oct3       if (tdif.lt.dtcrit .or. ath.gt.ath1) THEN
     685  ! oct3       if (tdif.lt.dtcrit .OR. ath.gt.ath1) THEN
    686686  ! oct3         do 60 k=1,nl
    687687  ! oct3            sig(i,k) = beta*sig(i,k) - 2.*alpha*tdif*tdif
    … …  
    794794  ! nn=0
    795795  ! do 101 i=1,len
    796   ! IF(iflag1(i).eq.0)THEN
     796  ! IF(iflag1(i).EQ.0)THEN
    797797  ! nn=nn+1
    798798  ! tra(nn,k,j)=tra1(i,k,j)
    799   ! endif
     799  ! END IF
    800800  ! 101  continue
    801801  ! 111  continue
    … …  
    11031103  ! inb(i)=k+1
    11041104  ! capem(i)=cape(i)
    1105   ! endif
    1106   ! endif
     1105  ! END IF
     1106  ! END IF
    11071107  ! 520    continue
    11081108  ! 530  continue
    … …  
    11301130  ! capem(i)=cape(i)
    11311131  ! byp(i)=(tvp(i,k+1)-tv(i,k+1))*dph(i,k+1)/p(i,k+1)
    1132   ! endif
    1133   ! endif
     1132  ! END IF
     1133  ! END IF
    11341134  ! 520    continue
    11351135  ! 530  continue
    … …  
    11531153  ! ori        do 520 i=1,ncum
    11541154  ! ori          IF(cape(i).lt.0.0)lcape(i)=.FALSE.
    1155   ! ori          if((k.ge.(icb(i)+1)).and.lcape(i))THEN
     1155  ! ori          if((k.ge.(icb(i)+1)).AND.lcape(i))THEN
    11561156  ! ori            by=(tvp(i,k)-tv(i,k))*dph(i,k)/p(i,k)
    11571157  ! ori            byp(i)=(tvp(i,k+1)-tv(i,k+1))*dph(i,k+1)/p(i,k+1)
    … …  
    15511551    ! do j=minorig,nl
    15521552    ! do il=1,ncum
    1553     ! IF( (i.ge.icb(il)).and.(i.le.inb(il)).and.
    1554     ! :       (j.ge.(icb(il)-1)).and.(j.le.inb(il)))THEN
     1553    ! IF( (i.ge.icb(il)).AND.(i.le.inb(il)).AND.
     1554    ! :       (j.ge.(icb(il)-1)).AND.(j.le.inb(il)))THEN
    15551555    ! traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
    15561556    ! :            +(1.-sij(il,i,j))*tra(il,nk(il),k)
    1557     ! endif
     1557    ! END IF
    15581558    ! enddo
    15591559    ! enddo
    … …  
    15711571    DO il = 1, ncum
    15721572      IF ((i>=icb(il)) .AND. (i<=inb(il)) .AND. (nent(il, i)==0)) THEN
    1573         ! @      IF(nent(il,i).eq.0)THEN
     1573        ! @      IF(nent(il,i).EQ.0)THEN
    15741574        ment(il, i, i) = m(il, i)
    15751575        qent(il, i, i) = rr(il, nk(il)) - ep(il, i) * clw(il, i)
    … …  
    15861586  ! do i=minorig+1,nl
    15871587  ! do il=1,ncum
    1588   ! if (i.ge.icb(il) .and. i.le.inb(il) .and. nent(il,i).eq.0) THEN
     1588  ! if (i.ge.icb(il) .AND. i.le.inb(il) .AND. nent(il,i).EQ.0) THEN
    15891589  ! traent(il,i,i,j)=tra(il,nk(il),j)
    1590   ! endif
     1590  ! END IF
    15911591  ! enddo
    15921592  ! enddo
    … …  
    17591759    ! do j=1,ntra
    17601760    ! do il=1,ncum
    1761     ! if ( i.ge.icb(il) .and. i.le.inb(il) .and. lwork(il)
    1762     ! :     .and. csum(il,i).lt.m(il,i) ) THEN
     1761    ! if ( i.ge.icb(il) .AND. i.le.inb(il) .AND. lwork(il)
     1762    ! :     .AND. csum(il,i).lt.m(il,i) ) THEN
    17631763    ! traent(il,i,i,j)=tra(il,nk(il),j)
    1764     ! endif
     1764    ! END IF
    17651765    ! enddo
    17661766    ! enddo
    … …  
    20942094          END IF
    20952095
    2096         END IF ! i.eq.1
     2096        END IF ! i.EQ.1
    20972097
    20982098        ! ***       find mixing ratio of precipitating downdraft     ***
    … …  
    23732373  ! :                     *(mp(il,2)*(trap(il,2,j)-tra(il,1,j))
    23742374  ! :             +am(il)*(tra(il,2,j)-tra(il,1,j)))
    2375   ! endif
     2375  ! END IF
    23762376  ! enddo
    23772377  ! enddo
    … …  
    24092409  ! ftra(il,1,k)=ftra(il,1,k)+0.1*work(il)*ment(il,j,1)
    24102410  ! :                *(traent(il,j,1,k)-tra(il,1,k))
    2411   ! endif
    2412 
    2413   ! endif
     2411  ! END IF
     2412
     2413  ! END IF
    24142414  ! enddo
    24152415  ! enddo
    … …  
    25312531    ! :         *(amp1(il)*(tra(il,i+1,k)-tra(il,i,k))
    25322532    ! :           -ad(il)*(tra(il,i,k)-tra(il,i-1,k)))
    2533     ! endif
    2534     ! endif
     2533    ! END IF
     2534    ! END IF
    25352535    ! enddo
    25362536    ! enddo
    … …  
    25802580    ! ftra(il,i,j)=ftra(il,i,j)+0.1*dpinv*ment(il,k,i)
    25812581    ! :        *(traent(il,k,i,j)-tra(il,i,j))
    2582     ! endif
    2583     ! endif
     2582    ! END IF
     2583    ! END IF
    25842584    ! enddo
    25852585    ! enddo
    … …  
    26142614    ! do k=i,nl+1
    26152615    ! do il=1,ncum
    2616     ! if (i.le.inb(il) .and. k.le.inb(il)) THEN
     2616    ! if (i.le.inb(il) .AND. k.le.inb(il)) THEN
    26172617    ! dpinv=1.0/(ph(il,i)-ph(il,i+1))
    26182618    ! cpinv=1.0/cpn(il,i)
    … …  
    26232623    ! ftra(il,i,j)=ftra(il,i,j)+0.1*dpinv*ment(il,k,i)
    26242624    ! :             *(traent(il,k,i,j)-tra(il,i,j))
    2625     ! endif
    2626     ! endif ! i and k
     2625    ! END IF
     2626    ! END IF ! i and k
    26272627    ! enddo
    26282628    ! enddo
    … …  
    26982698    ! :     *(mp(il,i+1)*(trap(il,i+1,j)-tra(il,i,j))
    26992699    ! :     -mp(il,i)*(trap(il,i,j)-tra(il,i-1,j)))
    2700     ! endif
    2701     ! endif ! i
     2700    ! END IF
     2701    ! END IF ! i
    27022702    ! enddo
    27032703    ! enddo
    … …  
    28442844    DO k = i, nl
    28452845      DO il = 1, ncum
    2846         ! test         if (i.ge.icb(il).and.i.le.inb(il).and.k.le.inb(il))
     2846        ! test         if (i.ge.icb(il).AND.i.le.inb(il).AND.k.le.inb(il))
    28472847        ! THEN
    28482848        IF (i<=inb(il) .AND. k<=inb(il)) THEN
    … …  
    32593259
    32603260  ! inputs:
    3261   integer ncum, nd, nloc
    3262   integer icb(nloc), inb(nloc)
    3263   real cape(nloc)
    3264   real clw(nloc, nd), lv(nloc, nd), t(nloc, nd), h(nloc, nd)
    3265   integer nk(nloc)
     3261  INTEGER ncum, nd, nloc
     3262  INTEGER icb(nloc), inb(nloc)
     3263  REAL cape(nloc)
     3264  REAL clw(nloc, nd), lv(nloc, nd), t(nloc, nd), h(nloc, nd)
     3265  INTEGER nk(nloc)
    32663266  ! inouts:
    3267   real ep(nloc, nd)
    3268   real hp(nloc, nd)
     3267  REAL ep(nloc, nd)
     3268  REAL hp(nloc, nd)
    32693269  ! outputs ou local
    3270   real epmax_diag(nloc)
     3270  REAL epmax_diag(nloc)
    32713271  ! locals
    3272   integer i, k
    3273   real hp_bak(nloc, nd)
     3272  INTEGER i, k
     3273  REAL hp_bak(nloc, nd)
    32743274  CHARACTER (LEN = 20) :: modname = 'cv30_epmax_fn_cape'
    32753275  CHARACTER (LEN = 80) :: abort_message
    … …  
    32773277  ! on recalcule ep et hp
    32783278
    3279   if (coef_epmax_cape>1e-12) THEN
     3279  IF (coef_epmax_cape>1e-12) THEN
    32803280    do i = 1, ncum
    32813281      epmax_diag(i) = epmax - coef_epmax_cape * sqrt(cape(i))
    … …  
    33003300    do k = minorig + 1, nl
    33013301      do i = 1, ncum
    3302         if((k>=icb(i)).and.(k<=inb(i)))THEN
     3302        IF((k>=icb(i)).AND.(k<=inb(i)))THEN
    33033303          hp(i, k) = h(i, nk(i)) + (lv(i, k) + (cpd - cpv) * t(i, k)) * ep(i, k) * clw(i, k)
    33043304        endif
    … …  
    33083308    do i = 1, ncum
    33093309      do k = 1, nl
    3310         if (abs(hp_bak(i, k) - hp(i, k))>0.01) THEN
     3310        IF (abs(hp_bak(i, k) - hp(i, k))>0.01) THEN
    33113311          WRITE(*, *) 'i,k=', i, k
    33123312          WRITE(*, *) 'coef_epmax_cape=', coef_epmax_cape
    … …  
    33263326      enddo !do k=1,nl
    33273327    enddo !do i=1,ncum
    3328   endif !if (coef_epmax_cape.gt.1e-12) THEN
     3328  ENDIF !if (coef_epmax_cape.gt.1e-12) THEN
    33293329END SUBROUTINE  cv30_epmax_fn_cape
    33303330
Note: See TracChangeset for help on using the changeset viewer.