diff --git a/src/netcdf_output_mod.f90 b/src/netcdf_output_mod.f90
index 5d8655d6cf8506079384f10027b1d73f8eab3b17..585df93fba86ffe0991f18cce93bfcb26f314352 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