From c9987a0230f3b29152dbea362ca2f81b85af7108 Mon Sep 17 00:00:00 2001 From: Michael Barry Date: Thu, 16 Feb 2023 21:34:45 +1000 Subject: [PATCH] Suggests unit change for diags to per day --- src/aed_pathogens.F90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/aed_pathogens.F90 b/src/aed_pathogens.F90 index def536e..858725e 100644 --- a/src/aed_pathogens.F90 +++ b/src/aed_pathogens.F90 @@ -617,11 +617,11 @@ SUBROUTINE aed_calculate_pathogens(data,column,layer_idx) !----------------------------------------------------------------- ! SET DIAGNOSTICS - _DIAG_VAR_(data%id_total(pth_i)) = pth_f + pth_a + pth_d ! orgs/m3/s + _DIAG_VAR_(data%id_total(pth_i)) = pth_f + pth_a + pth_d ! orgs/m3 IF ( diag_level >= 10 ) THEN - _DIAG_VAR_(data%id_growth(pth_i)) = growth*(pth_f + pth_a) ! orgs/m3/s - _DIAG_VAR_(data%id_sunlight(pth_i)) = light*pth_f + (light/2.)*pth_a ! orgs/m3/s - _DIAG_VAR_(data%id_mortality(pth_i)) = mortality*(pth_f + pth_a) ! orgs/m3/s + _DIAG_VAR_(data%id_growth(pth_i)) = growth*(pth_f + pth_a) * secs_per_day ! orgs/m3/d + _DIAG_VAR_(data%id_sunlight(pth_i)) = (light*pth_f + (light/2.)*pth_a) * secs_per_day ! orgs/m3/d + _DIAG_VAR_(data%id_mortality(pth_i)) = mortality*(pth_f + pth_a) * secs_per_day ! orgs/m3/d ENDIF ENDDO END SUBROUTINE aed_calculate_pathogens