Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_interp.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_interp.F90 (modified) (7 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_interp.F90

@@ -8 +8 @@
   ! CALL to this routine only treats the aerosol "id_aero".
 
-  ! 1) Read in data for the whole year, only at first time step
-  ! 2) Interpolate to the actual day, only at new day
-  ! 3) Interpolate to the model vertical grid (target grid), only at new day
+  ! 1) Read in data for the whole year, ONLY at first time step
+  ! 2) Interpolate to the actual day, ONLY at new day
+  ! 3) Interpolate to the model vertical grid (target grid), ONLY at new day
   ! 4) Test for negative mass values
 
@@ -102 +102 @@
   ! Calculation to find if it is a new day
 
-  IF(mpi_rank == 0 .AND. debug)then
+  IF(mpi_rank == 0 .AND. debug)THEN
     PRINT*, 'CONTROL PANEL REGARDING TIME STEPING'
   ENDIF
@@ -125 +125 @@
   ENDIF
 
-  IF(mpi_rank == 0 .AND. debug)then
+  IF(mpi_rank == 0 .AND. debug)THEN
     ! 0.02 is about 0.5/24, namly less than half an hour
     OLDNEWDAY = (r_day - REAL(iday) < 0.02)
@@ -280 +280 @@
     month_mid(:) = month_start (:) + month_len(:) / 2.
 
-    if (debug) then
-      write(lunout, *)' month_len = ', month_len
-      write(lunout, *)' month_mid = ', month_mid
+    if (debug) THEN
+      WRITE(lunout, *)' month_len = ', month_len
+      WRITE(lunout, *)' month_mid = ', month_mid
     endif
 
@@ -300 +300 @@
     ! Find which months and days to use for time interpolation
     nbr_tsteps = 12
-    IF (nbr_tsteps == 12) then
+    IF (nbr_tsteps == 12) THEN
       IF (jDay < month_mid(im + 1)) THEN
         im2 = im - 1
@@ -319 +319 @@
         ENDIF
       ENDIF
-    ELSE IF (nbr_tsteps == 14) then
+    ELSE IF (nbr_tsteps == 14) THEN
       im = im + 1
       IF (jDay < month_mid(im)) THEN
@@ -335 +335 @@
       CALL abort_physic('readaerosol_interp', 'number of months undefined', 1)
     ENDIF
-    if (debug) then
-      write(lunout, *)' jDay, day1, day2, im, im2 = ', jDay, day1, day2, im, im2
+    if (debug) THEN
+      WRITE(lunout, *)' jDay, day1, day2, im, im2 = ', jDay, day1, day2, im, im2
     endif