Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/top_bound_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/top_bound_loc.f90 (modified) (11 diffs)

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

@@ -60 +60 @@
   !   ----------
 
-  real,intent(inout) :: ucov(iip1,jjb_u:jje_u,llm) ! covariant zonal wind
-  real,intent(inout) :: vcov(iip1,jjb_v:jje_v,llm) ! covariant meridional wind
-  real,intent(inout) :: teta(iip1,jjb_u:jje_u,llm) ! potential temperature
-  real,intent(in) :: masse(iip1,jjb_u:jje_u,llm) ! mass of atmosphere
-  real,intent(in) :: dt ! time step (s) of sponge model
+  REAL,INTENT(INOUT) :: ucov(iip1,jjb_u:jje_u,llm) ! covariant zonal wind
+  REAL,INTENT(INOUT) :: vcov(iip1,jjb_v:jje_v,llm) ! covariant meridional wind
+  REAL,INTENT(INOUT) :: teta(iip1,jjb_u:jje_u,llm) ! potential temperature
+  REAL,INTENT(IN) :: masse(iip1,jjb_u:jje_u,llm) ! mass of atmosphere
+  REAL,INTENT(IN) :: dt ! time step (s) of sponge model
 
    ! REAL dv(iip1,jjb_v:jje_v,llm),du(iip1,jjb_u:jje_u,llm)
@@ -78 +78 @@
   INTEGER :: i
   REAL,SAVE :: rdamp(llm)
-  real,save :: lambda(llm) ! inverse or quenching time scale (Hz)
+  REAL,save :: lambda(llm) ! inverse or quenching time scale (Hz)
   LOGICAL,SAVE :: first=.TRUE.
   INTEGER :: j,l,jjb,jje
 
 
-  if (iflag_top_bound == 0) return
-
-  if (first) THEN
+  IF (iflag_top_bound == 0) return
+
+  IF (first) THEN
 !$OMP BARRIER
 !$OMP MASTER
-     if (iflag_top_bound == 1) THEN
+     IF (iflag_top_bound == 1) THEN
   ! sponge quenching over the topmost 4 atmospheric layers
          lambda(:)=0.
@@ -95 +95 @@
          lambda(llm-2)=tau_top_bound/4.
          lambda(llm-3)=tau_top_bound/8.
-     else if (iflag_top_bound == 2) THEN
+     ELSE IF (iflag_top_bound == 2) THEN
   ! sponge quenching over topmost layers down to pressures which are
   ! higher than 100 times the topmost layer pressure
@@ -110 +110 @@
      WRITE(lunout,*)'p (Pa)  z(km)  tau(s)   1./tau (Hz)'
      do l=1,llm
-       if (rdamp(l)/=0.) THEN
+       IF (rdamp(l)/=0.) THEN
          WRITE(lunout,'(6(1pe12.4,1x))') &
                presnivs(l),log(preff/presnivs(l))*scaleheight, &
@@ -119 +119 @@
 !$OMP END MASTER
 !$OMP BARRIER
-  endif ! of if (first)
+  ENDIF ! of if (first)
 
 
@@ -125 +125 @@
 
   ! compute zonal average of vcov (or set it to zero)
-  if (mode_top_bound>=2) THEN
+  IF (mode_top_bound>=2) THEN
    jjb=jj_begin
    jje=jj_end
@@ -150 +150 @@
    enddo
 !$OMP END DO NOWAIT
-  endif ! of if (mode_top_bound.ge.2)
+  ENDIF ! of if (mode_top_bound.ge.2)
 
   ! compute zonal average of u (or set it to zero)
-  if (mode_top_bound>=2) THEN
+  IF (mode_top_bound>=2) THEN
    jjb=jj_begin
    jje=jj_end
@@ -177 +177 @@
    enddo
 !$OMP END DO NOWAIT
-  endif ! of if (mode_top_bound.ge.2)
+  ENDIF ! of if (mode_top_bound.ge.2)
 
   ! compute zonal average of potential temperature, if necessary
-  if (mode_top_bound>=3) THEN
+  IF (mode_top_bound>=3) THEN
    jjb=jj_begin
    jje=jj_end
@@ -198 +198 @@
    enddo
 !$OMP END DO NOWAIT
-  endif ! of if (mode_top_bound.ge.3)
-
-  if (mode_top_bound>=1) THEN
+  ENDIF ! of if (mode_top_bound.ge.3)
+
+  IF (mode_top_bound>=1) THEN
    ! Apply sponge quenching on vcov:
    jjb=jj_begin
@@ -233 +233 @@
    enddo
 !$OMP END DO NOWAIT
-  endif ! of if (mode_top_bound.ge.1)
-
-  if (mode_top_bound>=3) THEN
+  ENDIF ! of if (mode_top_bound.ge.1)
+
+  IF (mode_top_bound>=3) THEN
    ! Apply sponge quenching on teta:
    jjb=jj_begin
@@ -252 +252 @@
    enddo
 !$OMP END DO NOWAIT
-  endif ! of if (mode_top_bond.ge.3)
+  ENDIF ! of if (mode_top_bond.ge.3)
 
 END SUBROUTINE top_bound_loc