Changeset 5270 for LMDZ6/trunk/libf/misc

Timestamp:

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

Author:

abarral

Message:

Replace F77 netcdf library by F90 netcdf library

File:

• LMDZ6/trunk/libf/misc/write_field.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/misc/write_field.f90

@@ -3 +3 @@
 !
 module write_field
+  USE netcdf, ONLY: nf90_sync, nf90_put_var, nf90_enddef, nf90_def_dim, nf90_unlimited, &
+          nf90_clobber, nf90_create, nf90_def_var, nf90_double
   USE strings_mod, ONLY: int2str
   IMPLICIT NONE; PRIVATE
@@ -75 +77 @@
     subroutine WriteField_gen(name,Field,dimx,dimy,dimz)
     implicit none
-    include 'netcdf.inc'
       character(len=*) :: name
       integer :: dimx,dimy,dimz
@@ -104 +105 @@
       count(4)=1
 
-      status = NF_PUT_VARA_DOUBLE(FieldId(Index),FieldVarId(Index),start,count,Field)
-      status = NF_SYNC(FieldId(Index))
+      status = nf90_put_var(FieldId(Index),FieldVarId(Index),Field,start,count)
+      status = nf90_sync(FieldId(Index))
       
     end subroutine WriteField_gen
@@ -111 +112 @@
     subroutine CreateNewField(name,dimx,dimy,dimz)
     implicit none
-    include 'netcdf.inc'  
       character(len=*) :: name
       integer :: dimx,dimy,dimz
@@ -123 +123 @@
       
       
-      status = NF_CREATE(TRIM(ADJUSTL(name))//'.nc', NF_CLOBBER, FieldId(NbField))
-      status = NF_DEF_DIM(FieldId(NbField),'X',dimx,TabDim(1))
-      status = NF_DEF_DIM(FieldId(NbField),'Y',dimy,TabDim(2))
-      status = NF_DEF_DIM(FieldId(NbField),'Z',dimz,TabDim(3))
-      status = NF_DEF_DIM(FieldId(NbField),'iter',NF_UNLIMITED,TabDim(4))
-      status = NF_DEF_VAR(FieldId(NbField),FieldName(NbField),NF_DOUBLE,4,TabDim,FieldVarId(NbField))
-      status = NF_ENDDEF(FieldId(NbField))
+      status = nf90_create(TRIM(ADJUSTL(name))//'.nc', nf90_clobber, FieldId(NbField))
+      status = nf90_def_dim(FieldId(NbField),'X',dimx,TabDim(1))
+      status = nf90_def_dim(FieldId(NbField),'Y',dimy,TabDim(2))
+      status = nf90_def_dim(FieldId(NbField),'Z',dimz,TabDim(3))
+      status = nf90_def_dim(FieldId(NbField),'iter',nf90_unlimited,TabDim(4))
+      status = nf90_def_var(FieldId(NbField),FieldName(NbField),nf90_double,TabDim,FieldVarId(NbField))
+      status = nf90_enddef(FieldId(NbField))
 
     end subroutine CreateNewField