Changeset 585 for codes/icosagcm/devel/src/diagnostics/diagflux.f90

Timestamp:

10/16/17 13:07:30 (7 years ago)

Author:

dubos

Message:

devel : reconstruct fluxes at cell centers

File:

• codes/icosagcm/devel/src/diagnostics/diagflux.f90 (modified) (1 diff)

codes/icosagcm/devel/src/diagnostics/diagflux.f90

@@ -37 +37 @@
     END DO
   END SUBROUTINE zero_qfluxt
+
+  SUBROUTINE flux_centered_lonlat(f_flux, f_flux_lon, f_flux_lat)
+    TYPE(t_field),POINTER :: f_flux(:), f_flux_lon(:), f_flux_lat(:)
+    REAL(rstd), POINTER :: flux(:,:,:), flux_lon(:,:,:), flux_lat(:,:,:)
+    INTEGER :: ind, itrac
+    DO ind=1,ndomain
+       IF (.NOT. assigned_domain(ind)) CYCLE
+       CALL swap_dimensions(ind)
+       CALL swap_geometry(ind)
+       flux=f_flux(ind)
+       DO itrac=1,nqtot
+          CALL compute_flux_centered_lonlat(flux(:,:,itrac), flux_lon(:,:,itrac), flux_lat(:,:,itrac))
+       END DO
+    END DO
+  END SUBROUTINE flux_centered_lonlat
   
+  SUBROUTINE compute_flux_centered_lonlat(flux, flux_lon, flux_lat)
+    REAL(rstd), INTENT(IN) :: flux(3*iim*jjm,llm)
+    REAL(rstd), INTENT(OUT) :: flux_lon(iim*jjm,llm), flux_lat(iim*jjm,llm)
+    REAL(rstd) :: flux_3d(iim*jjm,llm,3)
+    CALL compute_flux_centered(flux, flux_3d)
+    CALL compute_wind_centered_lonlat_compound(flux_3d, flux_lon, flux_lat)
+  END SUBROUTINE compute_flux_centered_lonlat
+
 END MODULE diagflux_mod