From 2861005b97308669ddfe13aef89771b993f59a50 Mon Sep 17 00:00:00 2001
From: lucieb92 <lucie.bakels@univie.ac.at>
Date: Wed, 19 Jun 2024 10:03:19 +0200
Subject: [PATCH] bug fix writing out TROPOPAUSE when using mdomainfill
---
src/netcdf_output_mod.f90 | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/src/netcdf_output_mod.f90 b/src/netcdf_output_mod.f90
index 5d8655d6..585df93f 100644
--- a/src/netcdf_output_mod.f90
+++ b/src/netcdf_output_mod.f90
@@ -2357,10 +2357,14 @@ subroutine partoutput_netcdf(itime,field,np,imass,ncid)
call nf90_err(nf90_put_var(ncid,tempIDend,oro(0:nx-1,0:ny-1), (/ 1,1 /),(/ nx,ny /)))
topo_written=.true.
endif
- else !HM or TR
+ else if (partopt(np)%name.eq.'HM') then !HM
call nf90_err(nf90_inq_varid(ncid=ncid,name=trim(partopt(np)%short_name),varid=tempIDend))
call nf90_err(nf90_put_var(ncid,tempIDend,hmix(0:nx-1,0:ny-1,1,memind(1)), &
(/ tpointer_part,1,1 /),(/ 1,nx,ny /)))
+ else !TR
+ call nf90_err(nf90_inq_varid(ncid=ncid,name=trim(partopt(np)%short_name),varid=tempIDend))
+ call nf90_err(nf90_put_var(ncid,tempIDend,tropopause(0:nx-1,0:ny-1,1,memind(1)), &
+ (/ tpointer_part,1,1 /),(/ 1,nx,ny /)))
endif
else if (partopt(np)%name.eq.'MA') then
--
GitLab