Make vertical remap working for 4D field

YM

Make vertical remap working for 4D field
YM
8a73a328 · ymeur · 71b9bed8 · 8a73a328
Commit 8a73a328 authored Jan 18, 2022 by ymeur
Hide whitespace changes
Inline Side-by-side

Showing with 8 additions and 5 deletions

src/vertical/vertical_remap.f90 src/vertical/vertical_remap.f90 +8 -5

No files found.
--- a/src/vertical/vertical_remap.f90
+++ b/src/vertical/vertical_remap.f90
@@ -20,7 +20,7 @@ CONTAINS
    REAL(rstd),POINTER :: out(:,:)
    REAL(rstd),POINTER :: p(:,:)
    
-    INTEGER :: ind
+    INTEGER :: ind,iq

    CALL allocate_field(f_p,field_t,type_real,llm+1)
    CALL pression(f_ps,f_p)
@@ -30,11 +30,14 @@ CONTAINS
      CALL swap_dimensions(ind)
      CALL swap_geometry(ind)
      p=f_p(ind)
-      in=field_in(ind)
-      out=field_out(ind)
-      CALL compute_vertical_remap(pressure_level,in,p,out)
+      DO iq=1,field_out(ind)%dim4
+        in=>field_in(ind)%rval4d(:,:,iq)
+        out=>field_out(ind)%rval4d(:,:,iq)
+        CALL compute_vertical_remap(pressure_level,in,p,out)
+      ENDDO
    ENDDO
-    
+    CALL deallocate_field(f_p)
+
  END SUBROUTINE  vertical_remap

  SUBROUTINE compute_vertical_remap(pressure_level,in,p,out)