Changeset 3965

Timestamp:

Jul 21, 2021, 4:05:15 PM (3 years ago)

Author:

acozic

Message:

update in LMDZ-INCA-Dyn branche function get_var2 in dyn3dmem/dynetat0_loc.F90 to allow 3d variables in start.nc to have OR not a fourth dimension with the time axis. For this copy what it was done for get_var1 function.
Commit not done on the trunk - waiting for information to know how do it in dyn3d/ directory

File:

• LMDZ6/branches/LMDZ-INCA-Dyn/libf/dyn3dmem/dynetat0_loc.F90 (modified) (1 diff)

LMDZ6/branches/LMDZ-INCA-Dyn/libf/dyn3dmem/dynetat0_loc.F90

@@ -229 +229 @@
   CHARACTER(LEN=*), INTENT(IN)  :: var
   REAL,             INTENT(OUT) :: v(:,:)
-  REAL,             ALLOCATABLE :: w4(:,:,:,:)
+  REAL,             ALLOCATABLE :: w4(:,:,:,:), w3(:,:,:)
   INTEGER :: nn(4), dids(4), k, nd
+
+
   CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
-  ierr=NF90_INQUIRE_VARIABLE(fID,vID,dimids=dids,ndims=nd)
+  ierr=NF90_INQUIRE_VARIABLE(fID,vID,ndims=nd)
+
+  IF(nd==1) THEN
+    CALL err(NF90_GET_VAR(fID,vID,v),"get",var); RETURN
+  END IF
+  ierr=NF90_INQUIRE_VARIABLE(fID,vID,dimids=dids)
+
   DO k=1,nd; ierr=NF90_INQUIRE_DIMENSION(fID,dids(k),len=nn(k)); END DO
-  ALLOCATE(w4(nn(1),nn(2),nn(3),nn(4)))
-  CALL err(NF90_GET_VAR(fID,vID,w4),"get",var)
-  v=RESHAPE(w4,[nn(1)*nn(2),nn(3)]); DEALLOCATE(w4)
+
+  SELECT CASE(nd)
+  CASE(3); ALLOCATE(w3(nn(1),nn(2),nn(3)))
+     CALL err(NF90_GET_VAR(fID,vID,w3),"get",var)
+     v=RESHAPE(w3,[nn(1)*nn(2),nn(3)]); DEALLOCATE(w3)
+  CASE(4);  ALLOCATE(w4(nn(1),nn(2),nn(3),nn(4)))
+     CALL err(NF90_GET_VAR(fID,vID,w4),"get",var)
+     v=RESHAPE(w4,[nn(1)*nn(2),nn(3)]); DEALLOCATE(w4)
+  END SELECT
 END SUBROUTINE get_var2