Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.f90

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:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.f90 (modified) (21 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pentes_ini.f90

@@ -86 +86 @@
      PRINT*,'SCHEMA AMONT NOUVEAU'
      first=.FALSE.
-     do i=2,iip1
+     DO i=2,iip1
         coslon(i)=cos(rlonv(i))
         sinlon(i)=sin(rlonv(i))
@@ -187 +187 @@
   !CC
    IF(mode==2) THEN
-      do l=1,llm
+      DO l=1,llm
         s0s=0.
         s0n=0.
@@ -196 +196 @@
         smn=0.
         sms=0.
-        do i=1,iim
+        DO i=1,iim
            smn=smn+sm(i,1,l)
            sms=sms+sm(i,jjp1,l)
@@ -208 +208 @@
            dys2=dys2+coslondlon(i)*zz
         enddo
-        do i=1,iim
+        DO i=1,iim
            sy(i,1,l)=dyn1*sinlon(i)+dyn2*coslon(i)
            sy(i,jjp1,l)=dys1*sinlon(i)+dys2*coslon(i)
         enddo
-        do i=1,iim
+        DO i=1,iim
            s0(i,1,l)=s0n/smn+sy(i,1,l)
            s0(i,jjp1,l)=s0s/sms-sy(i,jjp1,l)
@@ -220 +220 @@
         s0(iip1,jjp1,l)=s0(1,jjp1,l)
 
-        do i=1,iim
+        DO i=1,iim
            sxn(i)=s0(i+1,1,l)-s0(i,1,l)
            sxs(i)=s0(i+1,jjp1,l)-s0(i,jjp1,l)
   !   on rerentre les masses
         enddo
-        do i=1,iim
+        DO i=1,iim
            sy(i,1,l)=sy(i,1,l)*sm(i,1,l)
            sy(i,jjp1,l)=sy(i,jjp1,l)*sm(i,jjp1,l)
@@ -233 +233 @@
         sxn(iip1)=sxn(1)
         sxs(iip1)=sxs(1)
-        do i=1,iim
+        DO i=1,iim
            sx(i+1,1,l)=0.25*(sxn(i)+sxn(i+1))*sm(i+1,1,l)
            sx(i+1,jjp1,l)=0.25*(sxs(i)+sxs(i+1))*sm(i+1,jjp1,l)
@@ -247 +247 @@
 
   IF (mode==4) THEN
-     do l=1,llm
-        do i=1,iip1
+     DO l=1,llm
+        DO i=1,iip1
            sx(i,1,l)=0.
            sx(i,jjp1,l)=0.
@@ -261 +261 @@
   ! CALL minmaxq(zq,1.e33,-1.e33,'avant advy     ')
   IF (mode==4) THEN
-     do l=1,llm
-        do i=1,iip1
+     DO l=1,llm
+        DO i=1,iip1
            sx(i,1,l)=0.
            sx(i,jjp1,l)=0.
@@ -273 +273 @@
    CALL advy( limit,.5*dtvr,pbarv,sm,s0,sx,sy,sz )
   ! CALL minmaxq(zq,1.e33,-1.e33,'avant advz     ')
-   do j=1,jjp1
-      do i=1,iip1
+   DO j=1,jjp1
+      DO i=1,iip1
          sz(i,j,1)=0.
          sz(i,j,llm)=0.
@@ -282 +282 @@
    CALL advz( limit,dtvr,w,sm,s0,sx,sy,sz )
   IF (mode==4) THEN
-     do l=1,llm
-        do i=1,iip1
+     DO l=1,llm
+        DO i=1,iip1
            sx(i,1,l)=0.
            sx(i,jjp1,l)=0.
@@ -293 +293 @@
     CALL limy(s0,sy,sm,pente_max)
    CALL advy( limit,.5*dtvr,pbarv,sm,s0,sx,sy,sz )
-   do l=1,llm
-      do j=1,jjp1
+   DO l=1,llm
+      DO j=1,jjp1
          sm(iip1,j,l)=sm(1,j,l)
          s0(iip1,j,l)=s0(1,j,l)
@@ -306 +306 @@
   ! CALL minmaxq(zq,1.e33,-1.e33,'avant advx     ')
   IF (mode==4) THEN
-     do l=1,llm
-        do i=1,iip1
+     DO l=1,llm
+        DO i=1,iip1
            sx(i,1,l)=0.
            sx(i,jjp1,l)=0.
@@ -354 +354 @@
      dqzpn=ssum(iim,sz(1,1,l),1)/masn
      dqzps=ssum(iim,sz(1,jjp1,l),1)/mass
-     do i=1,iip1
+     DO i=1,iip1
         q( i,1,llm+1-l,3)=dqzpn
         q( i,jjp1,llm+1-l,3)=dqzps
@@ -365 +365 @@
         dyn2=0.
         dys2=0.
-        do i=1,iim
+        DO i=1,iim
            dyn1=dyn1+sinlondlon(i)*sy(i,1,l)/sm(i,1,l)
            dyn2=dyn2+coslondlon(i)*sy(i,1,l)/sm(i,1,l)
@@ -371 +371 @@
            dys2=dys2+coslondlon(i)*sy(i,jjp1,l)/sm(i,jjp1,l)
         enddo
-        do i=1,iim
+        DO i=1,iim
            q(i,1,llm+1-l,2)= &
                  (sinlon(i)*dyn1+coslon(i)*dyn2)
@@ -387 +387 @@
         dyn2=0.
         dys2=0.
-        do i=1,iim
+        DO i=1,iim
            zz=s0(i,2,l)/sm(i,2,l)-q(i,1,llm+1-l,0)
            dyn1=dyn1+sinlondlon(i)*zz
@@ -395 +395 @@
            dys2=dys2+coslondlon(i)*zz
         enddo
-        do i=1,iim
+        DO i=1,iim
            q(i,1,llm+1-l,2)= &
                  (sinlon(i)*dyn1+coslon(i)*dyn2)/2.
@@ -407 +407 @@
         q(iip1,jjp1,llm+1-l,0)=q(1,jjp1,llm+1-l,0)
 
-        do i=1,iim
+        DO i=1,iim
            sxn(i)=q(i+1,1,llm+1-l,0)-q(i,1,llm+1-l,0)
            sxs(i)=q(i+1,jjp1,llm+1-l,0)-q(i,jjp1,llm+1-l,0)
@@ -413 +413 @@
         sxn(iip1)=sxn(1)
         sxs(iip1)=sxs(1)
-        do i=1,iim
+        DO i=1,iim
            q(i+1,1,llm+1-l,1)=0.25*(sxn(i)+sxn(i+1))
            q(i+1,jjp1,llm+1-l,1)=0.25*(sxs(i)+sxs(i+1))
@@ -426 +426 @@
 
   ! bouclage en longitude
-  do iq=0,3
-     do l=1,llm
-        do j=1,jjp1
+  DO iq=0,3
+     DO l=1,llm
+        DO j=1,jjp1
            q(iip1,j,l,iq)=q(1,j,l,iq)
         enddo
@@ -455 +455 @@
     ! PRINT*, '-------------------------------------------'
 
-   do l=1,llm
-      do j=1,jjp1
-       do i=1,iip1
+   DO l=1,llm
+      DO j=1,jjp1
+       DO i=1,iip1
          IF(q(i,j,l,0)<qmin) &
                PRINT*,'apres pentes, s0(',i,',',j,',',l,')=',q(i,j,l,0)