Changeset 3940 for TOOLS/MOSAIX/src/MOSAIX/interpol.f90

Timestamp:

06/25/18 09:56:20 (6 years ago)

Author:

lbedidi

Message:

12. Bedidi: fixing 0D case

File:

• TOOLS/MOSAIX/src/MOSAIX/interpol.f90 (modified) (1 diff)

TOOLS/MOSAIX/src/MOSAIX/interpol.f90

@@ -101 +101 @@
       lat_src (:,:) = SPREAD ( lat_src(1,:), DIM=1, ncopies=nj_src)
    CASE default
-      IF ( nj_src == 1 ) THEN
+      IF ( nj_src == 1 .AND. ni_src==1) THEN
+         WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- lecture lon lat src 2D ', 6I9)") rank
+         CALL xios_get_domain_attr ("domain_src", lonvalue_2d=lon_src(:,:), latvalue_2d=lat_src(:,:) )
+      ELSE IF ( nj_src == 1 ) THEN
          WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- lecture lon lat src 1D ', 6I9)") rank
          CALL xios_get_domain_attr ("domain_src", lonvalue_1d=lon_src(:,1), latvalue_1d=lat_src(:,1) )