Skip to content

Commit

Permalink
Merge pull request #26 from fgeter/main
Browse files Browse the repository at this point in the history 
Added column units in the  header in the  file hru_cflux_stat.txt
  • Loading branch information
@odav
odav authored Nov 18, 2024
2 parents 1c25953 + 3ad364c commit 344efc6
Show file tree
Hide file tree
Showing 2 changed files with 52 additions and 0 deletions.
2 changes: 2 additions & 0 deletions src/output_landscape_init.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -513,11 +513,13 @@ subroutine output_landscape_init
if (pco%nb_hru%a == "y") then
write (4567,*) bsn%name, prog
write (4567,*) soil_org_flux_hdr
write (4567,*) soil_org_flux_hdr_units
write (9000,*) "HRU hru_cflux_stat.txt"
if (pco%csvout == "y") then
open (4568,file="hru_cflux_stat.csv", recl = 1500)
write (4568,*) bsn%name, prog
write (4568,'(*(G0.3,:,","))') soil_org_flux_hdr
write (4568,'(*(G0.3,:,","))') soil_org_flux_hdr_units
write (9000,*) "HRU hru_cflux_stat.csv"
end if
endif
Expand Down
50 changes: 50 additions & 0 deletions src/output_landscape_module.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -938,6 +938,56 @@ module output_landscape_module
! character(len=15) :: microb_c = " microb_c"
end type output_soil_org_flux_header
type (output_soil_org_flux_header) :: soil_org_flux_hdr

type output_soil_org_flux_header_units
character (len=11) :: day = " "
character (len=11) :: mo = " "
character (len=11) :: day_mo = " "
character (len=11) :: yrc = " "
character (len=16) :: isd = " "
character (len=22) :: id = " "
character (len=16) :: name = " "
character(len=15) :: cfmets1 = " kg_C/ha"
character(len=15) :: cfstrs1 = " kg_C/ha"
character(len=15) :: cfstrs2 = " kg_C/ha"
character(len=15) :: efmets1 = " kg_N/ha"
character(len=15) :: efstrs1 = " kg_N/ha"
character(len=15) :: efstrs2 = " kg_N/ha"
character(len=15) :: immmets1 = " kg_N/ha"
character(len=15) :: immstrs1 = " kg_N/ha"
character(len=15) :: immstrs2 = " kg_N/ha"
character(len=15) :: mnrmets1 = " kg_N/ha"
character(len=15) :: mnrstrs1 = " kg_N/ha"
character(len=15) :: mnrstrs2 = " kg_N/ha"
character(len=15) :: co2fmet = " kg_C/ha"
character(len=15) :: co2fstr = " kg_C/ha"
character(len=15) :: cfs1s2 = " kg_C/ha"
character(len=15) :: cfs1s3 = " kg_C/ha"
character(len=15) :: cfs2s1 = " kg_C/ha"
character(len=15) :: cfs2s3 = " kg_C/ha"
character(len=15) :: cfs3s1 = " kg_C/ha"
character(len=15) :: efs1s2 = " kg_N/ha"
character(len=15) :: efs1s3 = " kg_N/ha"
character(len=15) :: efs2s1 = " kg_N/ha"
character(len=15) :: efs2s3 = " kg_N/ha"
character(len=15) :: efs3s1 = " kg_N/ha"
character(len=15) :: imms1s2 = " kg_N/ha"
character(len=15) :: imms1s3 = " kg_N/ha"
character(len=15) :: imms2s1 = " kg_N/ha"
character(len=15) :: imms2s3 = " kg_N/ha"
character(len=15) :: imms3s1 = " kg_N/ha"
character(len=15) :: mnrs1s2 = " kg_N/ha"
character(len=15) :: mnrs1s3 = " kg_N/ha"
character(len=15) :: mnrs2s1 = " kg_N/ha"
character(len=15) :: mnrs2s3 = " kg_N/ha"
character(len=15) :: mnrs3s1 = " kg_N/ha"
character(len=15) :: co2fs1 = " kg_C/ha"
character(len=15) :: co2fs2 = " kg_C/ha"
character(len=15) :: co2fs3 = " kg_C/ha"
end type output_soil_org_flux_header_units
type (output_soil_org_flux_header_units) :: soil_org_flux_hdr_units



!!! NEW SOIL CARBON STAT OUTPUT

Expand Down

0 comments on commit 344efc6

Please sign in to comment.