Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advn.F

Timestamp:

Jul 23, 2024, 3:29:36 PM (3 months ago)

Author:

abarral

Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advn.F (modified) (11 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/advn.F

@@ -300 +300 @@
             zqv(ij,l)=q(ij,l)
             zqv(ip1jm-iip1+ij,l)=q(ip1jm+ij,l)
-            extremum(ij)=.true.
-            extremum(ip1jmp1-iip1+ij)=.true.
+            extremum(ij)=.TRUE.
+            extremum(ip1jmp1-iip1+ij)=.TRUE.
          enddo
 
@@ -404 +404 @@
             zqw(ij,2)=q(ij,1)
             zqw(ij,llm)=q(ij,llm)
-            extremum(ij,1)=.true.
-            extremum(ij,llm)=.true.
+            extremum(ij,1)=.TRUE.
+            extremum(ij,llm)=.TRUE.
          enddo
 
@@ -502 +502 @@
                zdq=qd(ij,l)-qg(ij,l)
 c              if((qd(ij,l)-q(ij,l))*(q(ij,l)-qg(ij,l)).lt.0.) then
-c                 print*,'probleme au point ij=',ij,'  l=',l
-c                 print*,qd(ij,l),q(ij,l),qg(ij,l)
+c                 PRINT*,'probleme au point ij=',ij,'  l=',l
+c                 PRINT*,qd(ij,l),q(ij,l),qg(ij,l)
 c                 qd(ij,l)=q(ij,l)
 c                 qg(ij,l)=q(ij,l)
@@ -512 +512 @@
 c                 if(.not.(zsigd(ij,l).ge.0..and.zsigd(ij,l).le.1. .and.
 c    s               zsigg(ij,l).ge.0..or.zsigg(ij,l).le.1.) ) then
-c                    print*,'probleme au point ij=',ij,'  l=',l
-c                    print*,'sigg=',zsigg(ij,l),'  sigd=',zsigd(ij,l)
-c                    print*,'q d,c,g ',qd(ij,l),q(ij,l),qg(ij,l),zdq
+c                    PRINT*,'probleme au point ij=',ij,'  l=',l
+c                    PRINT*,'sigg=',zsigg(ij,l),'  sigd=',zsigd(ij,l)
+c                    PRINT*,'q d,c,g ',qd(ij,l),q(ij,l),qg(ij,l),zdq
 c                    stop
 c                 endif
@@ -566 +566 @@
           endif
 c         if(zsig.lt.0.) then
-c            print*,'au point ij=',ij,'  l=',l,'  sig=',zsig
+c            PRINT*,'au point ij=',ij,'  l=',l,'  sig=',zsig
 c            stop
 c         endif
@@ -611 +611 @@
                enddo
                niju=iju
-c              print*,'niju,nl',niju,nl(l)
+c              PRINT*,'niju,nl',niju,nl(l)
 
 c  traitement des mailles
@@ -759 +759 @@
                zdq=qn(ij,l)-qs(ij,l)
 c              if((qn(ij,l)-q(ij,l))*(q(ij,l)-qs(ij,l)).lt.0.) then
-c                 print*,'probleme au point ij=',ij,'  l=',l,'  advnqx'
-c                 print*,qn(ij,l),q(ij,l),qs(ij,l)
+c                 PRINT*,'probleme au point ij=',ij,'  l=',l,'  advnqx'
+c                 PRINT*,qn(ij,l),q(ij,l),qs(ij,l)
 c                 qn(ij,l)=q(ij,l)
 c                 qs(ij,l)=q(ij,l)
@@ -769 +769 @@
 c                 if(.not.(zsign(ij).ge.0..and.zsign(ij).le.1. .and.
 c    s               zsigs(ij).ge.0..or.zsigs(ij).le.1.) ) then
-c                    print*,'probleme au point ij=',ij,'  l=',l
-c                    print*,'sigs=',zsigs(ij),'  sign=',zsign(ij)
+c                    PRINT*,'probleme au point ij=',ij,'  l=',l
+c                    PRINT*,'sigs=',zsigs(ij),'  sign=',zsign(ij)
 c                    stop
 c                 endif
@@ -886 +886 @@
                zdq=qb(ij,l)-qh(ij,l)
 c              if((qh(ij,l)-q(ij,l))*(q(ij,l)-qb(ij,l)).lt.0.) then
-c                 print*,'probleme au point ij=',ij,'  l=',l
-c                 print*,qh(ij,l),q(ij,l),qb(ij,l)
+c                 PRINT*,'probleme au point ij=',ij,'  l=',l
+c                 PRINT*,qh(ij,l),q(ij,l),qb(ij,l)
 c                 qh(ij,l)=q(ij,l)
 c                 qb(ij,l)=q(ij,l)
@@ -903 +903 @@
        enddo
 
-c      print*,'ok1'
+c      PRINT*,'ok1'
 c   calcul de la pente maximum dans la maille en valeur absolue
        do l=2,llm
@@ -947 +947 @@
          enddo
       enddo
-c     print*,'ok3'
+c     PRINT*,'ok3'
       RETURN
       END