Changeset 5084 for LMDZ6/trunk/libf/phylmd/readaerosolstrato.F90

Timestamp:

Jul 19, 2024, 6:40:44 PM (4 months ago)

Author:

Laurent Fairhead

Message:

Reverting to r4065. Updating fortran standard broke too much stuff. Will do it by smaller chunks
AB, LF

File:

• LMDZ6/trunk/libf/phylmd/readaerosolstrato.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/readaerosolstrato.F90

@@ -3 +3 @@
     use netcdf95, only: nf95_close, nf95_gw_var, nf95_inq_dimid, & 
                         nf95_inq_varid, nf95_open
-    use lmdz_netcdf, only: nf90_get_var, nf90_noerr, nf90_nowrite
+    use netcdf, only: nf90_get_var, nf90_noerr, nf90_nowrite
 
     USE phys_cal_mod, ONLY : mth_cur
@@ -68 +68 @@
 
 !--only read file if beginning of run or start of new month
-    IF (debut.OR.mth_cur/=mth_pre) THEN
+    IF (debut.OR.mth_cur.NE.mth_pre) THEN
 
 !--only root reads
     IF (is_mpi_root.AND.is_omp_root) THEN
 
-    IF (nbands/=2) THEN
+    IF (nbands.NE.2) THEN 
         abort_message='nbands doit etre egal a 2 dans readaerosolstrat'
         CALL abort_physic(modname,abort_message,1)
@@ -83 +83 @@
     CALL nf95_gw_var(ncid_in, varid, lev)
     n_lev = size(lev)
-    IF (n_lev/=klev) THEN
+    IF (n_lev.NE.klev) THEN 
        abort_message='Le nombre de niveaux n est pas egal a klev'
        CALL abort_physic(modname,abort_message,1)
@@ -93 +93 @@
     WRITE(lunout,*) 'LAT aerosol strato=', n_lat, latitude
     IF (grid_type/=unstructured) THEN
-      IF (n_lat/=nbp_lat) THEN
+      IF (n_lat.NE.nbp_lat) THEN 
          abort_message='Le nombre de lat n est pas egal a nbp_lat'
          CALL abort_physic(modname,abort_message,1)
@@ -104 +104 @@
     IF (grid_type/=unstructured) THEN
       WRITE(lunout,*) 'LON aerosol strato=', n_lon, longitude
-      IF (n_lon/=nbp_lon) THEN
+      IF (n_lon.NE.nbp_lon) THEN 
          abort_message='Le nombre de lon n est pas egal a nbp_lon'
          CALL abort_physic(modname,abort_message,1)
@@ -114 +114 @@
     n_month = size(time)
     WRITE(lunout,*) 'TIME aerosol strato=', n_month, time
-    IF (n_month/=12) THEN
+    IF (n_month.NE.12) THEN 
        abort_message='Le nombre de month n est pas egal a 12'
        CALL abort_physic(modname,abort_message,1)
@@ -131 +131 @@
 
 !---select the correct month
-    IF (mth_cur<1.OR.mth_cur>12) THEN
+    IF (mth_cur.LT.1.OR.mth_cur.GT.12) THEN
      WRITE(lunout,*) 'probleme avec le mois dans readaerosolstrat =', mth_cur
     ENDIF