Changeset 5270 for LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

Timestamp:

Oct 24, 2024, 1:55:38 PM (39 hours ago)

Author:

abarral

Message:

Replace F77 netcdf library by F90 netcdf library

File:

• LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

@@ -4568 +4568 @@
       USE mod_grid_phy_lmdz
       USE mod_phys_lmdz_para
+      USE netcdf, ONLY: nf90_open, nf90_close, nf90_inq_varid, nf90_nowrite, nf90_noerr, nf90_get_var
       IMPLICIT NONE
-
-      include "netcdf.inc"
 
       CHARACTER*800 filescaleparams
@@ -4581 +4580 @@
       integer start(4),count(4), status
 !      local
-      integer debutread,countread
       CHARACTER*104 varname
       CHARACTER*2 aux_2s
@@ -4587 +4585 @@
 !$OMP MASTER
       IF (is_mpi_root .AND. is_omp_root) THEN
-          !nci=NCOPN(trim(adjustl(filescaleparams)),NCNOWRIT,rcode)
-         ierr = NF_OPEN (trim(adjustl(filescaleparams)),NF_NOWRITE, nid)
-          if (ierr .EQ. NF_NOERR) THEN
-          debutread=step_sca
-          countread=1
+          !nci=nf90_open(trim(adjustl(filescaleparams)),nf90_nowrite,rcode)
+         ierr = nf90_open (trim(adjustl(filescaleparams)),nf90_nowrite, nid)
+          if (ierr .EQ. nf90_noerr) THEN
 
            do i=1,nbreg
@@ -4597 +4593 @@
             varname= trim(adjustl(paramname))//aux_2s
             print *,varname
-            ierr = NF_INQ_VARID (nid,trim(adjustl(varname)), nvarid)
-            ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debutread,          &
-                         countread, auxreal)
-            IF (ierr .NE. NF_NOERR) THEN
+            ierr = nf90_inq_varid(nid,trim(adjustl(varname)), nvarid)
+            ierr = nf90_get_var(nid, nvarid, auxreal, [step_sca])
+            IF (ierr .NE. nf90_noerr) THEN
              PRINT*, 'Pb de lecture pour modvalues'
        print *,'JE  scale_var, step_sca',trim(adjustl(varname)),step_sca
@@ -4613 +4608 @@
            enddo
 
-            ierr = NF_CLOSE(nid)
+            ierr = nf90_close(nid)
           else
            print *,'File '//trim(adjustl(filescaleparams))//' not found'