diff --git a/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/.nn_convection_flux.f90.swp b/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/.nn_convection_flux.f90.swp new file mode 100644 index 0000000..4ffa687 Binary files /dev/null and b/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/.nn_convection_flux.f90.swp differ diff --git a/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_convection_flux.f90 b/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_convection_flux.f90 index 98bf22a..c649e99 100644 --- a/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_convection_flux.f90 +++ b/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_convection_flux.f90 @@ -531,6 +531,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -547,7 +548,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_diffusion.f90 b/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_diffusion.f90 index 1b2f31a..ec2d81c 100644 --- a/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_diffusion.f90 +++ b/sam_code_NN/sam_cases/run_files_16_missing_bits_out_in_only/nn_diffusion.f90 @@ -421,7 +421,8 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. - do ib=0,nt + !round x down + do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -438,7 +439,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_convection_flux.f90 b/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_convection_flux.f90 index a4caf67..4085204 100644 --- a/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_convection_flux.f90 +++ b/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_convection_flux.f90 @@ -531,6 +531,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -547,7 +548,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_diffusion.f90 b/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_diffusion.f90 index 805b5f6..044988e 100644 --- a/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_diffusion.f90 +++ b/sam_code_NN/sam_cases/run_files_18_missing_bits_out_in_only/nn_diffusion.f90 @@ -421,6 +421,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -438,7 +439,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/.nn_diffusion.f90.swp b/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/.nn_diffusion.f90.swp deleted file mode 100644 index 1e07df5..0000000 Binary files a/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/.nn_diffusion.f90.swp and /dev/null differ diff --git a/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_convection_flux.f90 b/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_convection_flux.f90 index 8327071..da8bfe6 100644 --- a/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_convection_flux.f90 +++ b/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_convection_flux.f90 @@ -531,6 +531,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -547,7 +548,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_diffusion.f90 b/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_diffusion.f90 index 590ad58..8b1011a 100644 --- a/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_diffusion.f90 +++ b/sam_code_NN/sam_cases/run_files_20_missing_bits_out_in_only/nn_diffusion.f90 @@ -421,6 +421,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -438,7 +439,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_convection_flux.f90 b/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_convection_flux.f90 index 852d075..de4fb86 100644 --- a/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_convection_flux.f90 +++ b/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_convection_flux.f90 @@ -531,6 +531,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -547,7 +548,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp diff --git a/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_diffusion.f90 b/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_diffusion.f90 index dff0fc4..f57659d 100644 --- a/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_diffusion.f90 +++ b/sam_code_NN/sam_cases/run_files_22_missing_bits_out_in_only/nn_diffusion.f90 @@ -421,6 +421,7 @@ subroutine f_trunc_32( x ) ! bits 23-30 are exponent ! bits 0-22 are the mantissa, least significant bits are 0-7. +!round x down do ib=0,nt itemp = ibclr(itemp,ib) ! bit to zero enddo @@ -438,7 +439,7 @@ subroutine f_trunc_32( x ) endif enddo - if (abs(temp2 -x) < abs(temp-x)) then + if (abs(temp2 -x) < abs(temp-x)) then !tie breaking rule -> round toward zero x = temp2 else x = temp