Changeset 5075 for LMDZ6/trunk/libf/phylmd/readaerosol_mod.F90

Timestamp:

Jul 19, 2024, 10:05:57 AM (2 months ago)

Author:

abarral

Message:

[continued & end] replace netcdf by lmdz_netcdf.F90 wrapper
"use netcdf" is now only used in lmdz_netcdf.F90 (except ecrad and obsolete/)
<include "netcdf.inc"> is now likewise only used in lmdz_netcdf.F90.

systematically specify explicitely <USE lmdz_netcdf, ONLY:> (probably left some missing, to correct later on)

Further replacement of nf_put_* by nf90_put_* (same for _get_)

[minor] replace deprecated boolean operators along the way

File:

• LMDZ6/trunk/libf/phylmd/readaerosol_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/readaerosol_mod.F90

@@ -2 +2 @@
 !
 MODULE readaerosol_mod
+
+  USE lmdz_netcdf, ONLY: nf90_strerror,nf90_noerr,nf90_get_var,nf90_inq_varid,&
+          nf90_inquire_dimension,nf90_inq_dimid,nf90_open,nf90_nowrite,nf90_close
 
   REAL, SAVE :: not_valid=-333.
@@ -86 +89 @@
 ! Read data depending on actual year and interpolate if necessary
 !****************************************************************************************
-     IF (iyr_in .LT. 1850) THEN
+     IF (iyr_in < 1850) THEN
         cyear='.nat'
         WRITE(lunout,*) 'get_aero 1 iyr_in=', iyr_in,'   ',cyear
@@ -93 +96 @@
         CALL get_aero_fromfile(name_aero, cyear, filename, klev_src, pt_ap, pt_b, pt_out, psurf, load)
         
-     ELSE IF (iyr_in .GE. 2100) THEN
+     ELSE IF (iyr_in >= 2100) THEN
         cyear='2100'
         WRITE(lunout,*) 'get_aero 2 iyr_in=', iyr_in,'   ',cyear
@@ -103 +106 @@
         ! Read data from 2 decades and interpolate to actual year
         ! a) from actual 10-yr-period
-        IF (iyr_in.LT.1900) THEN
+        IF (iyr_in<1900) THEN
            iyr1 = 1850
            iyr2 = 1900
-        ELSE IF (iyr_in.GE.1900.AND.iyr_in.LT.1920) THEN
+        ELSE IF (iyr_in>=1900.AND.iyr_in<1920) THEN
            iyr1 = 1900
            iyr2 = 1920
@@ -174 +177 @@
 
 SUBROUTINE init_aero_fromfile(flag_aerosol, aerosol_couple)
-  USE netcdf
   USE mod_phys_lmdz_para
   USE mod_grid_phy_lmdz, ONLY: grid_type, unstructured
@@ -265 +267 @@
 !****************************************************************************************
 
-    USE netcdf
     USE dimphy
     USE mod_grid_phy_lmdz, ONLY: nbp_lon_=>nbp_lon, nbp_lat_=>nbp_lat, klon_glo, &
@@ -507 +508 @@
 !****************************************************************************************
           DO imth=1, 12
-             IF (imth.EQ.1) THEN
+             IF (imth==1) THEN
                 cvar=TRIM(varname)//'JAN'
-             ELSE IF (imth.EQ.2) THEN
+             ELSE IF (imth==2) THEN
                 cvar=TRIM(varname)//'FEB'
-             ELSE IF (imth.EQ.3) THEN
+             ELSE IF (imth==3) THEN
                 cvar=TRIM(varname)//'MAR'
-             ELSE IF (imth.EQ.4) THEN
+             ELSE IF (imth==4) THEN
                 cvar=TRIM(varname)//'APR'
-             ELSE IF (imth.EQ.5) THEN
+             ELSE IF (imth==5) THEN
                 cvar=TRIM(varname)//'MAY'
-             ELSE IF (imth.EQ.6) THEN
+             ELSE IF (imth==6) THEN
                 cvar=TRIM(varname)//'JUN'
-             ELSE IF (imth.EQ.7) THEN
+             ELSE IF (imth==7) THEN
                 cvar=TRIM(varname)//'JUL'
-             ELSE IF (imth.EQ.8) THEN
+             ELSE IF (imth==8) THEN
                 cvar=TRIM(varname)//'AUG'
-             ELSE IF (imth.EQ.9) THEN
+             ELSE IF (imth==9) THEN
                 cvar=TRIM(varname)//'SEP'
-             ELSE IF (imth.EQ.10) THEN
+             ELSE IF (imth==10) THEN
                 cvar=TRIM(varname)//'OCT'
-             ELSE IF (imth.EQ.11) THEN
+             ELSE IF (imth==11) THEN
                 cvar=TRIM(varname)//'NOV'
-             ELSE IF (imth.EQ.12) THEN
+             ELSE IF (imth==12) THEN
                 cvar=TRIM(varname)//'DEC'
              END IF
@@ -716 +717 @@
 
   SUBROUTINE check_err(status,text)
-    USE netcdf
     USE print_control_mod, ONLY: lunout
     IMPLICIT NONE