Changeset 5270 for LMDZ6/trunk/libf/phylmd/Dust/read_surface.f90

Timestamp:

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

Author:

abarral

Message:

Replace F77 netcdf library by F90 netcdf library

File:

• LMDZ6/trunk/libf/phylmd/Dust/read_surface.f90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/Dust/read_surface.f90

@@ -10 +10 @@
        USE mod_phys_lmdz_para
        USE iophy
-       USE netcdf, ONLY: nf90_get_var
+       USE netcdf, ONLY: nf90_inq_varid,nf90_noerr,nf90_get_var,nf90_nowrite,nf90_inq_varid,nf90_open
        IMPLICIT NONE
-
-       INCLUDE "netcdf.inc"
        INCLUDE "dimensions.h"
        INCLUDE "paramet.h"
@@ -31 +29 @@
        real surfa_glo(klon_glo,5)
 !
-       integer ncid
-       integer varid
-       real rcode
+       integer ncid, varid, rcode
        integer start(2),count(2),status
        integer i,j,l,ig
@@ -43 +39 @@
       real, dimension(jjp1) :: lats
       real, dimension(nbp_lat) :: lats_glo
-      real :: rcode2
       integer, dimension(1) :: startj,endj
 !JE20140526>>
@@ -50 +45 @@
 
        print*,'Lecture du fichier donnees_lisa.nc' 
-       ncid=NCOPN('donnees_lisa.nc',NCNOWRIT,rcode)
+       ncid=nf90_open('donnees_lisa.nc',nf90_nowrite,rcode)
 
 !JE20140526<<: check if are inversed or not the latitude grid in donnes_lisa
@@ -62 +57 @@
        if (i==4) aux4s='Latu'
        if (i==5) aux4s='latU'
-       status = NF_INQ_VARID (ncid, aux4s, rcode)
+       status = nf90_inq_varid(ncid, aux4s, rcode)
 !       print *,'stat,i',status,i,outcycle,aux4s
-!       print *,'ifclause',status.NE. NF_NOERR ,outcycle == .false.
-       IF ((.not.(status.NE. NF_NOERR) ).and.( .not. outcycle )) THEN
+!       print *,'ifclause',status.NE. nf90_noerr ,outcycle == .false.
+       IF ((.not.(status.NE. nf90_noerr) ).and.( .not. outcycle )) THEN
          outcycle=.true.
          latstr=aux4s
@@ -73 +68 @@
 !      endj(1)=jjp1
       endj(1)=nbp_lat
-      varid=NCVID(ncid,latstr,rcode)
+      varid=nf90_inq_varid(ncid,latstr,rcode)
 
           status = nf90_get_var(ncid, varid, lats_glo, startj, endj)
 !      print *,latstr,varid,status,jjp1,rcode
-!      IF (status .NE. NF_NOERR) print*,'NOOOOOOO'
+!      IF (status .NE. nf90_noerr) print*,'NOOOOOOO'
 !      print *,lats
 !stop
@@ -90 +85 @@
           varname=trim(name)//str1
        print*,'lecture variable:',varname
-          varid=NCVID(ncid,trim(varname),rcode)
-!          varid=NCVID(ncid,varname,rcode)
+          varid=nf90_inq_varid(ncid,trim(varname),rcode)
+!          varid=nf90_inq_varid(ncid,varname,rcode)
 
 !  dimensions pour les champs scalaires et le vent zonal