Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlspltqs_loc.f90

Timestamp:

Jul 24, 2024, 4:23:34 PM (2 months ago)

Author:

abarral

Message:

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

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlspltqs_loc.f90 (modified) (15 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlspltqs_loc.f90

@@ -59 +59 @@
   ije=ije_x
 
-  if (pole_nord.and.ijb==1) ijb=ijb+iip1
-  if (pole_sud.and.ije==ip1jmp1)  ije=ije-iip1
+  IF (pole_nord.AND.ijb==1) ijb=ijb+iip1
+  IF (pole_sud.AND.ije==ip1jmp1)  ije=ije-iip1
 
   IF (pente_max>-1.e-5) THEN
@@ -157 +157 @@
 !$OMP END DO NOWAIT
 
-  if (pole_nord) THEN
+  IF (pole_nord) THEN
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     DO l=1,llm
@@ -163 +163 @@
     ENDDO
 !$OMP END DO NOWAIT
-  endif
-
-  if (pole_sud)  THEN
+  ENDIF
+
+  IF (pole_sud)  THEN
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     DO l=1,llm
@@ -171 +171 @@
     ENDDO
 !$OMP END DO NOWAIT
-  endif
+  ENDIF
 
   !   calcul des flux a gauche et a droite
@@ -247 +247 @@
   !   indicage des mailles concernees par le traitement special
            DO ij=ijb,ije
-              IF(iadvplus(ij,l)==1.and.mod(ij,iip1)/=0) THEN
+              IF(iadvplus(ij,l)==1.AND.mod(ij,iip1)/=0) THEN
                  iju=iju+1
                  indu(iju)=ij
@@ -319 +319 @@
         !MVals: veiller a ce qu'on n'ait pas de denominateur nul
         masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),min_qMass)
-        if (q(ij,l,iq)>min_qParent) then ! modif 13 nov 2020
+        IF (q(ij,l,iq)>min_qParent) then ! modif 13 nov 2020
           Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
         else
@@ -448 +448 @@
   ijb=ij_begin-2*iip1
   ije=ij_end+2*iip1
-  if (pole_nord) ijb=ij_begin
-  if (pole_sud)  ije=ij_end
+  IF (pole_nord) ijb=ij_begin
+  IF (pole_sud)  ije=ij_end
   ij=3525
   l=3
-  if ((ij>=ijb).and.(ij<=ije)) THEN
+  IF ((ij>=ijb).AND.(ij<=ije)) THEN
     !WRITE(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
   ! &             ij,l,iq,ijb,q(ij,l,:)
-  endif
+  ENDIF
 
   IF(first) THEN
@@ -488 +488 @@
   !    le pole nord) qui sera utilisee pour evaluer les pentes au pole.
 
-  if (pole_nord) THEN
+  IF (pole_nord) THEN
     DO i = 1, iim
       airescb(i) = aire(i+ iip1) * q(i+ iip1,l,iq)
     ENDDO
     qpns   = SSUM( iim,  airescb ,1 ) / airej2
-  endif
-
-  if (pole_sud) THEN
+  ENDIF
+
+  IF (pole_sud) THEN
     DO i = 1, iim
       airesch(i) = aire(i+ ip1jm- iip1) * q(i+ ip1jm- iip1,l,iq)
     ENDDO
     qpsn   = SSUM( iim,  airesch ,1 ) / airejjm
-  endif
+  ENDIF
 
 
@@ -507 +507 @@
   ijb=ij_begin-2*iip1
   ije=ij_end+iip1
-  if (pole_nord) ijb=ij_begin
-  if (pole_sud)  ije=ij_end-iip1
+  IF (pole_nord) ijb=ij_begin
+  IF (pole_sud)  ije=ij_end-iip1
 
   DO ij=ijb,ije
@@ -520 +520 @@
   ijb=ij_begin-iip1
   ije=ij_end+iip1
-  if (pole_nord) ijb=ij_begin+iip1
-  if (pole_sud)  ije=ij_end-iip1
+  IF (pole_nord) ijb=ij_begin+iip1
+  IF (pole_sud)  ije=ij_end-iip1
 
   DO ij=ijb,ije
@@ -664 +664 @@
   ijb=ij_begin-iip1
   ije=ij_end+iip1
-  if (pole_nord) ijb=ij_begin+iip1
-  if (pole_sud)  ije=ij_end-iip1
+  IF (pole_nord) ijb=ij_begin+iip1
+  IF (pole_sud)  ije=ij_end-iip1
 
   DO ij=ijb,ije
@@ -680 +680 @@
   ijb=ij_begin-iip1
   ije=ij_end
-  if (pole_nord) ijb=ij_begin
-  if (pole_sud)  ije=ij_end-iip1
+  IF (pole_nord) ijb=ij_begin
+  IF (pole_sud)  ije=ij_end-iip1
 
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
@@ -708 +708 @@
   ijbm=ij_begin-iip1
   ijem=ij_end+iip1
-  if (pole_nord) ijb=ij_begin
-  if (pole_sud)  ije=ij_end
-  if (pole_nord) ijbm=ij_begin
-  if (pole_sud)  ijem=ij_end
+  IF (pole_nord) ijb=ij_begin
+  IF (pole_sud)  ije=ij_end
+  IF (pole_nord) ijbm=ij_begin
+  IF (pole_sud)  ijem=ij_end
 
   !WRITE(lunout,*) 'vlspltqs 737: iq,ijb,ije=',iq,ijb,ije
@@ -731 +731 @@
         !MVals: veiller a ce qu'on n'ait pas de denominateur nul
         !WRITE(lunout,*) 'ij,l,q(ij,l,iq)=',ij,l,q(ij,l,iq)
-        if (q(ij,l,iq)>min_qParent) then ! modif 13 nov 2020
+        IF (q(ij,l,iq)>min_qParent) then ! modif 13 nov 2020
           Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
         else
@@ -751 +751 @@
   ijb=ij_begin
   ije=ij_end
-  if (pole_nord) ijb=ij_begin+iip1
-  if (pole_sud)  ije=ij_end-iip1
+  IF (pole_nord) ijb=ij_begin+iip1
+  IF (pole_sud)  ije=ij_end-iip1
 
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)