Changeset 4288 for LMDZ6

Timestamp:

Oct 1, 2022, 2:48:32 PM (2 years ago)

Author:

fhourdin

Message:

Une derniere correction pour le 1D

File:

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

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

@@ -576 +576 @@
        'radT','uw','vw','q1','q2','sens','flat','ts','ustar'/
 
+
   do i=1,nbvar3d
      print *,'Dans read_cas2, on va lire ',nid,i,name_var(i)
@@ -589 +590 @@
      print *,'Dans read_cas2, on va lire ',var3didin(i),name_var(i)
      if(i.LE.35) then
-        ierr = NF90_GET_VAR(nid,var3didin(i),resul)
+        ierr = NF90_GET_VAR(nid,var3didin(i),resul, count = [1, 1, nlevel, ntime])
         print *,'Dans read_cas2, on a lu ',ierr,var3didin(i),name_var(i)
         if(ierr/=NF_NOERR) then
@@ -596 +597 @@
         endif
      else
-        ierr = NF90_GET_VAR(nid,var3didin(i),resul1)
         print *,'Dans read_cas2, on a lu ',ierr,var3didin(i),name_var(i)
+        ierr = NF90_GET_VAR(nid,var3didin(i),resul1, count = [1, 1, ntime])
         if(ierr/=NF_NOERR) then
            print *,'Pb a la lecture de cas.nc: ',name_var(i)
@@ -711 +712 @@
         !-----------------------------------------------------------------------
         if(i.LE.4) then     ! Lecture des coord pression en (nlevelp1,lat,lon)
-           ierr = NF90_GET_VAR(nid,var3didin(i),apbp, &
-                count = [1, 1, nlevel + 1])
+           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
@@ -720 +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, &
-                count = [1, 1, nlevel, ntime])
+           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