Changeset 5270 for LMDZ6/trunk/libf/dyn3d

Timestamp:

Oct 24, 2024, 1:55:38 PM (4 months ago)

Author:

abarral

Message:

Replace F77 netcdf library by F90 netcdf library

File:

• LMDZ6/trunk/libf/dyn3d/guide_mod.f90 (modified) (10 diffs)

LMDZ6/trunk/libf/dyn3d/guide_mod.f90

@@ -14 +14 @@
                     nf90_inq_dimid, nf90_inquire_dimension
   use pres2lev_mod, only: pres2lev
+  USE netcdf, ONLY: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, &
+          nf90_inq_dimid, nf90_inquire_dimension, nf90_float, nf90_def_var, &
+          nf90_create, nf90_def_dim, nf90_open, nf90_unlimited, nf90_write, nf90_enddef, nf90_redef, &
+          nf90_close, nf90_inq_varid, nf90_get_var, nf90_noerr, nf90_clobber, &
+          nf90_64bit_offset, nf90_inq_dimid, nf90_inquire_dimension, nf90_put_var
 
   IMPLICIT NONE
@@ -77 +82 @@
     INCLUDE "dimensions.h"
     INCLUDE "paramet.h"
-    INCLUDE "netcdf.inc"
 
     INTEGER                :: error,ncidpl,rid,rcod
@@ -224 +228 @@
 
     endif 
-    error=NF_INQ_DIMID(ncidpl,'LEVEL',rid)
-    IF (error.NE.NF90_NOERR) error=NF_INQ_DIMID(ncidpl,'PRESSURE',rid)
+    error=nf90_inq_dimid(ncidpl,'LEVEL',rid)
+    IF (error.NE.NF90_NOERR) error=nf90_inq_dimid(ncidpl,'PRESSURE',rid)
     IF (error.NE.NF90_NOERR) THEN
         CALL abort_gcm(modname,'Nudging: error reading pressure levels',1)
     ENDIF
-    error=NF_INQ_DIMLEN(ncidpl,rid,nlevnc)
+    error=nf90_inquire_dimension(ncidpl,rid,len=nlevnc)
     write(*,*)trim(modname)//' : number of vertical levels nlevnc', nlevnc 
     rcod = nf90_close(ncidpl)
@@ -1660 +1664 @@
     INCLUDE "dimensions.h"
     INCLUDE "paramet.h"
-    INCLUDE "netcdf.inc"
     INCLUDE "comgeom2.h"
     
@@ -1686 +1689 @@
 ! ----------------------------------------------
 ! Ouverture du fichier
-        ierr=NF_CREATE("guide_ins.nc",IOR(NF_CLOBBER,NF_64BIT_OFFSET),nid)
+        ierr=nf90_create("guide_ins.nc",IOR(nf90_clobber,nf90_64bit_offset),nid)
 ! Definition des dimensions
-        ierr=NF_DEF_DIM(nid,"LONU",iip1,id_lonu) 
-        ierr=NF_DEF_DIM(nid,"LONV",iip1,id_lonv) 
-        ierr=NF_DEF_DIM(nid,"LATU",jjp1,id_latu) 
-        ierr=NF_DEF_DIM(nid,"LATV",jjm,id_latv) 
-        ierr=NF_DEF_DIM(nid,"LEVEL",llm,id_lev)
-        ierr=NF_DEF_DIM(nid,"TIME",NF_UNLIMITED,id_tim)
+        ierr=nf90_def_dim(nid,"LONU",iip1,id_lonu)
+        ierr=nf90_def_dim(nid,"LONV",iip1,id_lonv)
+        ierr=nf90_def_dim(nid,"LATU",jjp1,id_latu)
+        ierr=nf90_def_dim(nid,"LATV",jjm,id_latv)
+        ierr=nf90_def_dim(nid,"LEVEL",llm,id_lev)
+        ierr=nf90_def_dim(nid,"TIME",nf90_unlimited,id_tim)
 
 ! Creation des variables dimensions
@@ -1710 +1713 @@
              varid_alpha_q)
         
-        ierr=NF_ENDDEF(nid)
+        ierr=nf90_enddef(nid)
 
 ! Enregistrement des variables dimensions
@@ -1730 +1733 @@
 ! Cr�ation des variables sauvegard�es
 ! --------------------------------------------------------------------
-        ierr = NF_REDEF(nid)
+        ierr = nf90_redef(nid)
 ! Pressure (GCM)
         dim4=(/id_lonv,id_latu,id_lev,id_tim/)
@@ -1764 +1767 @@
         ENDIF
         
-        ierr = NF_ENDDEF(nid)
-        ierr = NF_CLOSE(nid)
+        ierr = nf90_enddef(nid)
+        ierr = nf90_close(nid)
     ENDIF ! timestep=0
 
@@ -1771 +1774 @@
 ! Enregistrement du champ
 ! --------------------------------------------------------------------
-    ierr=NF_OPEN("guide_ins.nc",NF_WRITE,nid)
+    ierr=nf90_open("guide_ins.nc",nf90_write,nid)
 
     IF (varname=="SP") timestep=timestep+1
 
-    ierr = NF_INQ_VARID(nid,varname,varid)
+    ierr = nf90_inq_varid(nid,varname,varid)
     SELECT CASE (varname)
     CASE ("SP","ps")
@@ -1800 +1803 @@
 
     ierr = nf90_put_var(nid, varid, field2, start, count)
-    ierr = NF_CLOSE(nid)
+    ierr = nf90_close(nid)
 
   END SUBROUTINE guide_out