From d88a4a7fb32a33b73183f49db5d45ba8082a4338 Mon Sep 17 00:00:00 2001
From: Joern Behrens <joern.behrens@uni-hamburg.de>
Date: Sat, 19 Mar 2022 22:31:47 +0100
Subject: [PATCH] fixed precision in netcdf driver
---
.../src/options-sphere/ADV_wind_stacked.f90 | 24 +++++++++----------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/flash2d/src/options-sphere/ADV_wind_stacked.f90 b/flash2d/src/options-sphere/ADV_wind_stacked.f90
index af2458c..8d7efcf 100644
--- a/flash2d/src/options-sphere/ADV_wind_stacked.f90
+++ b/flash2d/src/options-sphere/ADV_wind_stacked.f90
@@ -204,14 +204,14 @@
TYPE (control_struct), INTENT(in) :: p_control
!---------- local variables
- CHARACTER (len=io_fillen) :: c_uname, c_vname, c_wname
- INTEGER (KIND = GRID_SI) :: i_alct, i_ncstat
- INTEGER (KIND = GRID_SI) :: i_filexid, i_fileyid, i_filezid
- INTEGER (KIND = GRID_SI) :: i_dimid, i_varid
- INTEGER (KIND = GRID_SI) :: i_tm, i_ln, i_lt, i_ly
- REAL (KIND = GRID_SR), DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxx
- REAL (KIND = GRID_SR), DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxy
- REAL (KIND = GRID_SR), DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxz
+ CHARACTER (len=io_fillen) :: c_uname, c_vname, c_wname
+ INTEGER (KIND = GRID_SI) :: i_alct, i_ncstat
+ INTEGER (KIND = GRID_SI) :: i_filexid, i_fileyid, i_filezid
+ INTEGER (KIND = GRID_SI) :: i_dimid, i_varid
+ INTEGER (KIND = GRID_SI) :: i_tm, i_ln, i_lt, i_ly
+ REAL*4, DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxx
+ REAL*4, DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxy
+ REAL*4, DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxz
!---------- initialize file names
WRITE(GRID_parameters%ioout,*) &
@@ -572,11 +572,11 @@
IF(r_dy /= 0.0_GRID_SR) THEN
r_dyi= 1._GRID_SR/r_dy
- r_inter(1)= (r_hy* r_l1+ r_ly* r_h1)* r_dyi* r_scalx
- r_inter(2)= (r_hy* r_l2+ r_ly* r_h2)* r_dyi* r_scaly
+ r_inter(1)= (r_hy* r_l1+ r_ly* r_h1)* r_dyi* r_scaly
+ r_inter(2)= (r_hy* r_l2+ r_ly* r_h2)* r_dyi* r_scalx
ELSE
- r_inter(1)= r_l1* r_scalx
- r_inter(2)= r_l2* r_scaly
+ r_inter(1)= r_l1* r_scaly
+ r_inter(2)= r_l2* r_scalx
END IF
!---------- calculate wind vector in kartesian coordinates
--
GitLab