Changeset 4287

Timestamp:

Sep 29, 2022, 9:19:32 PM (2 years ago)

Author:

lguez

Message:

Bug fix: add count argument to nf90_get_var

Because, in the rico case, NetCDF variables which are read have
dimensions (..., lat, lon), in NetCDF order, with lon = lat =

1. Without the count argument, for apbp for example, nf90_get_var

tried to read nlevel + 1 subscript values in the dimension lon, 1
subscript value in the lat dimension, and 1 subscript value in the
nlevp1 dimension: not what we want. I have checked that the cases
ARMCU/REF, amma and rico work with this revision.

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read2.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read2.F90

@@ -711 +711 @@
         !-----------------------------------------------------------------------
         if(i.LE.4) then     ! Lecture des coord pression en (nlevelp1,lat,lon)
-           ierr = NF90_GET_VAR(nid,var3didin(i),apbp)
+           ierr = NF90_GET_VAR(nid,var3didin(i),apbp, &
+                count = [1, 1, nlevel + 1])
            print *,'read2_cas(apbp), on a lu ',i,name_var(i)
            if(ierr/=NF_NOERR) then
@@ -719 +720 @@
            !-----------------------------------------------------------------------
         else if(i.gt.4.and.i.LE.45) then   ! Lecture des variables en (time,nlevel,lat,lon)
-           ierr = NF90_GET_VAR(nid,var3didin(i),resul)
+           ierr = NF90_GET_VAR(nid,var3didin(i),resul, &
+                count = [1, 1, nlevel, ntime])
            print *,'read2_cas(resul), on a lu ',i,name_var(i)
            if(ierr/=NF_NOERR) then
@@ -727 +729 @@
            !-----------------------------------------------------------------------
         else if (i.gt.45.and.i.LE.51) then   ! Lecture des variables en (time,lat,lon)
-           ierr = NF90_GET_VAR(nid,var3didin(i),resul2)
+           ierr = NF90_GET_VAR(nid,var3didin(i),resul2, count = [1, 1, ntime])
            print *,'read2_cas(resul2), on a lu ',i,name_var(i)
            if(ierr/=NF_NOERR) then