Ignore:
Timestamp:
Jul 24, 2024, 4:23:34 PM (3 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/dyn3d_common/advn.f90

    r5116 r5117  
    148148  REAL :: dxqu(ip1jmp1),zqu(ip1jmp1)
    149149  REAL :: zqmax(ip1jmp1),zqmin(ip1jmp1)
    150   logical :: extremum(ip1jmp1)
     150  LOGICAL :: extremum(ip1jmp1)
    151151
    152152  INTEGER :: mode
    … …  
    156156  !   calcul des pentes en u:
    157157  !   -----------------------
    158   if (mode==0) THEN
     158  IF (mode==0) THEN
    159159     do l=1,llm
    160160        do ij=1,ip1jm
    … …  
    230230     enddo
    231231     do ij=iip2,ip1jm-1
    232         if (extremum(ij).and..not.extremum(ij+1)) &
     232        IF (extremum(ij).and..not.extremum(ij+1)) &
    233233              qg(ij+1,l)=q(ij,l)
    234234     enddo
    … …  
    2392398888   continue
    240240  enddo
    241   endif
     241  ENDIF
    242242  RETURN
    243243END SUBROUTINE advnqx
    … …  
    265265  REAL :: dyqv(ip1jm),zqv(ip1jm,llm)
    266266  REAL :: zqmax(ip1jm),zqmin(ip1jm)
    267   logical :: extremum(ip1jmp1)
     267  LOGICAL :: extremum(ip1jmp1)
    268268
    269269  INTEGER :: mode
    … …  
    271271  data mode/1/
    272272
    273   if (mode==0) THEN
     273  IF (mode==0) THEN
    274274     do l=1,llm
    275275        do ij=1,ip1jmp1
    … …  
    318318           qs(ij,l)=q(ij,l)
    319319           qn(ij,l)=q(ij,l)
    320            ! if (.not.extremum(ij-iip1)) qs(ij-iip1,l)=q(ij,l)
    321            ! if (.not.extremum(ij+iip1)) qn(ij+iip1,l)=q(ij,l)
     320           ! if (.NOT.extremum(ij-iip1)) qs(ij-iip1,l)=q(ij,l)
     321           ! if (.NOT.extremum(ij+iip1)) qn(ij+iip1,l)=q(ij,l)
    322322        else
    323323           qs(ij,l)=zqv(ij,l)
    … …  
    334334
    335335  enddo
    336   endif
     336  ENDIF
    337337  RETURN
    338338END SUBROUTINE advnqy
    … …  
    361361  REAL :: dzqw(ip1jmp1,llm+1),zqw(ip1jmp1,llm+1)
    362362  REAL :: zqmax(ip1jmp1,llm),zqmin(ip1jmp1,llm)
    363   logical :: extremum(ip1jmp1,llm)
     363  LOGICAL :: extremum(ip1jmp1,llm)
    364364
    365365  INTEGER :: mode
    … …  
    371371  !   -----------------------
    372372
    373   if (mode==0) THEN
     373  IF (mode==0) THEN
    374374     do l=1,llm
    375375        do ij=1,ip1jmp1
    … …  
    436436  !    do ij=1,ip1jmp1
    437437  !       IF(extremum(ij,l)) THEN
    438   !          if (.not.extremum(ij,l-1)) qh(ij,l-1)=q(ij,l)
    439   !          if (.not.extremum(ij,l+1)) qb(ij,l+1)=q(ij,l)
     438  !          if (.NOT.extremum(ij,l-1)) qh(ij,l-1)=q(ij,l)
     439  !          if (.NOT.extremum(ij,l+1)) qb(ij,l+1)=q(ij,l)
    440440  !       endif
    441441  !    enddo
    … …  
    449449  enddo
    450450
    451   endif
     451  ENDIF
    452452
    453453  RETURN
    … …  
    491491  REAL :: zm,zq,zsigm,zsigp,zqm,zqp,zu
    492492
    493   logical :: ladvplus(ip1jmp1,llm)
     493  LOGICAL :: ladvplus(ip1jmp1,llm)
    494494
    495495  REAL :: prec
    … …  
    506506           !    qd(ij,l)=q(ij,l)
    507507           !    qg(ij,l)=q(ij,l)
    508            ! endif
     508           ! END IF
    509509           IF(abs(zdq)>prec) THEN
    510510              zsigd(ij,l)=(q(ij,l)-qg(ij,l))/zdq
    511511              zsigg(ij,l)=1.-zsigd(ij,l)
    512               ! IF(.not.(zsigd(ij,l).ge.0..and.zsigd(ij,l).le.1. .and.
     512              ! IF(.NOT.(zsigd(ij,l).ge.0..and.zsigd(ij,l).le.1. .AND.
    513513  !    s               zsigg(ij,l).ge.0..or.zsigg(ij,l).le.1.) ) THEN
    514514              !    PRINT*,'probleme au point ij=',ij,'  l=',l
    … …  
    516516              !    PRINT*,'q d,c,g ',qd(ij,l),q(ij,l),qg(ij,l),zdq
    517517              !    stop
    518               ! endif
     518              ! END IF
    519519           else
    520520              zsigd(ij,l)=0.5
    … …  
    530530   do l=1,llm
    531531   do ij=iip2,ip1jm-1
    532       if (u_m(ij,l)>=0.) THEN
     532      IF (u_m(ij,l)>=0.) THEN
    533533         zsigp=zsigd(ij,l)
    534534         zsigm=zsigg(ij,l)
    … …  
    549549      zsig=zu/zm
    550550      IF(zsig==0.) zsigp=0.1
    551       if (mode==1) THEN
    552          if (zsig<=zsigp) THEN
     551      IF (mode==1) THEN
     552         IF (zsig<=zsigp) THEN
    553553             u_mq(ij,l)=u_m(ij,l)*zqp
    554          else if (mode==1) THEN
     554         ELSE IF (mode==1) THEN
    555555             u_mq(ij,l)= &
    556556                   sign(zm,u_m(ij,l))*(zsigp*zqp+(zsig-zsigp)*zqm)
    557557         endif
    558558      else
    559          if (zsig<=zsigp) THEN
     559         IF (zsig<=zsigp) THEN
    560560             u_mq(ij,l)=u_m(ij,l)*(zqp-0.5*zsig/zsigp*(zqp-zq))
    561561         else
    … …  
    568568      !    PRINT*,'au point ij=',ij,'  l=',l,'  sig=',zsig
    569569      !    stop
    570       ! endif
     570      ! END IF
    571571  enddo
    572572  enddo
    … …  
    605605  !   indicage des mailles concernees par le traitement special
    606606           do ij=iip2,ip1jm
    607               IF(ladvplus(ij,l).and.mod(ij,iip1)/=0) THEN
     607              IF(ladvplus(ij,l).AND.mod(ij,iip1)/=0) THEN
    608608                 iju=iju+1
    609609                 indu(iju)=ij
    … …  
    639639      ! goto 8888
    640640            zz=0.5*(zsig-zsigd(ijq,l))/zsigg(ijq,l)
    641             IF(.not.(zz>0..and.zz<=0.5)) THEN
     641            IF(.NOT.(zz>0..and.zz<=0.5)) THEN
    642642                 WRITE(lunout,*)'probleme2 au point ij=',ij, &
    643643                       '  l=',l
    … …  
    669669        ! goto 9999
    670670            zz=0.5*(zsig-zsigg(ijq,l))/zsigd(ijq,l)
    671             IF(.not.(zz>0..and.zz<=0.5)) THEN
     671            IF(.NOT.(zz>0..and.zz<=0.5)) THEN
    672672                 WRITE(lunout,*)'probleme22 au point ij=',ij &
    673673                       ,'  l=',l
    … …  
    685685        endif
    686686     enddo
    687   endif  ! n0.gt.0
     687  ENDIF  ! n0.gt.0
    688688
    689689  !   bouclage en latitude
    … …  
    763763           !    qn(ij,l)=q(ij,l)
    764764           !    qs(ij,l)=q(ij,l)
    765            ! endif
     765           ! END IF
    766766           IF(abs(zdq)>prec) THEN
    767767              zsign(ij)=(q(ij,l)-qs(ij,l))/zdq
    768768              zsigs(ij)=1.-zsign(ij)
    769               ! IF(.not.(zsign(ij).ge.0..and.zsign(ij).le.1. .and.
     769              ! IF(.NOT.(zsign(ij).ge.0..and.zsign(ij).le.1. .AND.
    770770  !    s               zsigs(ij).ge.0..or.zsigs(ij).le.1.) ) THEN
    771771              !    PRINT*,'probleme au point ij=',ij,'  l=',l
    772772              !    PRINT*,'sigs=',zsigs(ij),'  sign=',zsign(ij)
    773773              !    stop
    774               ! endif
     774              ! END IF
    775775           else
    776776              zsign(ij)=0.5
    … …  
    782782
    783783   do ij=1,ip1jm
    784       if (v_m(ij,l)>=0.) THEN
     784      IF (v_m(ij,l)>=0.) THEN
    785785         zsigp=zsign(ij+iip1)
    786786         zsigm=zsigs(ij+iip1)
    … …  
    799799      zsig=abs(v_m(ij,l))/zm
    800800      IF(zsig==0.) zsigp=0.1
    801       if (zsig<=zsigp) THEN
     801      IF (zsig<=zsigp) THEN
    802802          v_mq(ij,l)=v_m(ij,l)*(zqp-0.5*zsig/zsigp*(zqp-zq))
    803803      else
    … …  
    890890           !    qh(ij,l)=q(ij,l)
    891891           !    qb(ij,l)=q(ij,l)
    892            ! endif
     892           ! END IF
    893893
    894894           IF(abs(zdq)>prec) THEN
    … …  
    907907   do l=2,llm
    908908   do ij=1,ip1jmp1
    909       if (w_m(ij,l)>=0.) THEN
     909      IF (w_m(ij,l)>=0.) THEN
    910910         zsigp=zsigb(ij,l)
    911911         zsigm=zsigh(ij,l)
    … …  
    924924      zsig=abs(w_m(ij,l))/zm
    925925      IF(zsig==0.) zsigp=0.1
    926       if (zsig<=zsigp) THEN
     926      IF (zsig<=zsigp) THEN
    927927          w_mq(ij,l)=w_m(ij,l)*(zqp-0.5*zsig/zsigp*(zqp-zq))
    928928      else
Note: See TracChangeset for help on using the changeset viewer.