Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmd/regr_lat_time_coefoz_m.F90

Timestamp:

Jul 24, 2024, 4:23:34 PM (3 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/phylmd/regr_lat_time_coefoz_m.F90 (modified) (10 diffs)

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

@@ -33 +33 @@
     ! -- the latitude is in degrees and strictly monotonic (as all
     ! NetCDF coordinate variables should be);
-    ! -- time increases from January to December (even though we do
+    ! -- time increases from January to December (even though we DO
     ! not use values of the input time coordinate);
     ! -- pressure is in hPa and in strictly ascending order (even
@@ -40 +40 @@
     ! when we regrid in pressure).
 
-    use lmdz_grid_phy, ONLY: nbp_lat
-    use regr_conserv_m, ONLY: regr_conserv
-    use regr_lint_m, ONLY: regr_lint
-    use netcdf95, ONLY: nf95_open, nf95_close, nf95_inq_varid, nf95_get_var, &
+    USE lmdz_grid_phy, ONLY: nbp_lat
+    USE lmdz_regr_conserv, ONLY: regr_conserv
+    USE lmdz_regr_lint, ONLY: regr_lint
+    USE netcdf95, ONLY: nf95_open, nf95_close, nf95_inq_varid, nf95_get_var, &
          nf95_put_var, nf95_gw_var
-    use netcdf, ONLY: nf90_nowrite
-    use nrtype, ONLY: pi
-    use lmdz_regular_lonlat, ONLY: boundslat_reg, south
+    USE netcdf, ONLY: nf90_nowrite
+    USE lmdz_physical_constants, ONLY: pi
+    USE lmdz_regular_lonlat, ONLY: boundslat_reg, south
 
     ! Variables local to the procedure:
 
-    integer ncid_in, ncid_out ! NetCDF IDs for input and output files
-    integer n_plev ! number of pressure levels in the input data
-    integer n_lat! number of latitudes in the input data
-
-    real, allocatable:: latitude(:)
+    INTEGER ncid_in, ncid_out ! NetCDF IDs for input and output files
+    INTEGER n_plev ! number of pressure levels in the input data
+    INTEGER n_lat! number of latitudes in the input data
+
+    REAL, ALLOCATABLE:: latitude(:)
     ! (of input data, converted to rad, sorted in strictly ascending order)
 
-    real, allocatable:: lat_in_edg(:)
+    REAL, ALLOCATABLE:: lat_in_edg(:)
     ! (edges of latitude intervals for input data, in rad, in strictly
     ! ascending order)
 
-    real, allocatable:: plev(:) ! pressure level of input data
-    logical desc_lat ! latitude in descending order in the input file
-
-    real, allocatable:: o3_par_in(:, :, :) ! (n_lat, n_plev, 12)
+    REAL, ALLOCATABLE:: plev(:) ! pressure level of input data
+    LOGICAL desc_lat ! latitude in descending order in the input file
+
+    REAL, ALLOCATABLE:: o3_par_in(:, :, :) ! (n_lat, n_plev, 12)
     ! (ozone parameter from the input file)
     ! ("o3_par_in(j, l, month)" is at latitude "latitude(j)" and pressure
     ! level "plev(l)". "month" is between 1 and 12.)
 
-    real, allocatable:: v_regr_lat(:, :, :) ! (jjm + 1, n_plev, 0:13)
+    REAL, ALLOCATABLE:: v_regr_lat(:, :, :) ! (jjm + 1, n_plev, 0:13)
     ! (mean of a variable "v" over a latitude interval)
     ! (First dimension is latitude interval.
@@ -84 +84 @@
     ! Last dimension is month number.)
 
-    real, allocatable:: o3_par_out(:, :, :) ! (jjm + 1, n_plev, ndays)
+    REAL, ALLOCATABLE:: o3_par_out(:, :, :) ! (jjm + 1, n_plev, ndays)
     ! (regridded ozone parameter)
     ! ("o3_par_out(j, l, day)" is at latitude "rlatu(j)", pressure
@@ -90 +90 @@
     ! 360-day calendar.)
 
-    integer j
-    integer i_v ! index of ozone parameter
-    integer, parameter:: n_o3_param = 8 ! number of ozone parameters
+    INTEGER j
+    INTEGER i_v ! index of ozone parameter
+    INTEGER, parameter:: n_o3_param = 8 ! number of ozone parameters
 
     CHARACTER(LEN=11) name_in(n_o3_param)
@@ -100 +100 @@
     ! (name of NetCDF primary variable in the output file)
 
-    integer varid_in(n_o3_param), varid_out(n_o3_param), varid_plev, varid_time
-    integer ncerr, varid
+    INTEGER varid_in(n_o3_param), varid_out(n_o3_param), varid_plev, varid_time
+    INTEGER ncerr, varid
     ! (for NetCDF)
 
-    real, parameter:: tmidmonth(0:13) = (/(-15. + 30. * j, j = 0, 13)/)
+    REAL, parameter:: tmidmonth(0:13) = (/(-15. + 30. * j, j = 0, 13)/)
     ! (time to middle of month, in days since January 1st 0h, in a
     ! 360-day calendar)
@@ -110 +110 @@
     ! interpolated between the December and the January value.)
 
-    real, parameter:: tmidday(360) = (/(j + 0.5, j = 0, 359)/)
+    REAL, parameter:: tmidday(360) = (/(j + 0.5, j = 0, 359)/)
     ! (time to middle of day, in days since January 1st 0h, in a
     ! 360-day calendar)
@@ -165 +165 @@
     ! ascending order, so invert order if necessary:
     desc_lat = latitude(1) > latitude(n_lat)
-    if (desc_lat) latitude = latitude(n_lat:1:-1)
+    IF (desc_lat) latitude = latitude(n_lat:1:-1)
 
     ! Compute edges of latitude intervals:
@@ -202 +202 @@
 
        CALL nf95_get_var(ncid_in, varid_in(i_v), o3_par_in)
-       if (desc_lat) o3_par_in = o3_par_in(n_lat:1:-1, :, :)
+       IF (desc_lat) o3_par_in = o3_par_in(n_lat:1:-1, :, :)
 
        ! Regrid in latitude:
@@ -240 +240 @@
     ! dimensions and variables, and writes one of the coordinate variables.
 
-    use lmdz_grid_phy, ONLY: nbp_lat
-    use lmdz_assert_eq, ONLY: assert_eq
-
-    use netcdf95, ONLY: nf95_create, nf95_def_dim, nf95_def_var, &
+    USE lmdz_grid_phy, ONLY: nbp_lat
+    USE lmdz_assert_eq, ONLY: assert_eq
+
+    USE netcdf95, ONLY: nf95_create, nf95_def_dim, nf95_def_var, &
          nf95_put_att, nf95_enddef, nf95_copy_att, nf95_put_var
-    use netcdf, ONLY: nf90_clobber, nf90_float, nf90_copy_att, nf90_global
-    use nrtype, ONLY: pi
-    use lmdz_regular_lonlat, ONLY: lat_reg
-
-    integer, intent(in):: ncid_in, varid_in(:), n_plev
-    CHARACTER(LEN=*), intent(in):: name_out(:) ! of NetCDF variables
-    integer, intent(out):: ncid_out, varid_out(:), varid_plev, varid_time
+    USE netcdf, ONLY: nf90_clobber, nf90_float, nf90_copy_att, nf90_global
+    USE lmdz_physical_constants, ONLY: pi
+    USE lmdz_regular_lonlat, ONLY: lat_reg
+
+    INTEGER, INTENT(IN):: ncid_in, varid_in(:), n_plev
+    CHARACTER(LEN=*), INTENT(IN):: name_out(:) ! of NetCDF variables
+    INTEGER, INTENT(OUT):: ncid_out, varid_out(:), varid_plev, varid_time
 
     ! Variables local to the procedure:
 
-    integer ncerr
-    integer dimid_rlatu, dimid_plev, dimid_time
-    integer varid_rlatu
-    integer i, n_o3_param
+    INTEGER ncerr
+    INTEGER dimid_rlatu, dimid_plev, dimid_time
+    INTEGER varid_rlatu
+    INTEGER i, n_o3_param
 
     !---------------------------
@@ -328 +328 @@
     SUBROUTINE handle_err_copy_att(att_name)
 
-      use netcdf, ONLY: nf90_noerr, nf90_strerror
-
-      CHARACTER(LEN=*), intent(in):: att_name
+      USE netcdf, ONLY: nf90_noerr, nf90_strerror
+
+      CHARACTER(LEN=*), INTENT(IN):: att_name
 
       !----------------------------------------
 
-      if (ncerr /= nf90_noerr) THEN
+      IF (ncerr /= nf90_noerr) THEN
          print *, "prepare_out " // trim(name_out(i)) &
               // " nf90_copy_att " // att_name // " -- " &