Add cisun and cisha to history

mvdebolskiy · Dec 6, 2024 · 4ec97f2 · 4ec97f2
1 parent b9e4b4f
commit 4ec97f2
Showing 1 changed file with 24 additions and 0 deletions.
diff --git a/src/biogeophys/PhotosynthesisMod.F90 b/src/biogeophys/PhotosynthesisMod.F90
@@ -709,6 +709,30 @@ subroutine InitHistory(this, bounds)
 
     endif
 
+    ! megan variables
+    this%cisun_z_patch(begp:endp,:) = spval
+    this%cisha_z_patch(begp:endp,:) = spval
+    if (nlevcan>1) then 
+       call hist_addfld2d (fname='CISUN', units='(Pa)', type2d='nlevcan', &
+          avgflag='A', long_name='sunlit leaf intracellular CO2 pressure', &
+          ptr_patch=this%cisun_z_patch, set_lake=spval, set_urb=spval, default='inactive')
+
+       call hist_addfld2d (fname='CISHA', units='(Pa)', type2d='nlevcan', &
+          avgflag='A', long_name='shaded leaf intracellular CO2 pressure', &
+          ptr_patch=this%cisha_z_patch, set_lake=spval, set_urb=spval, default='inactive')
+    else
+       ptr_1d => this%cisun_z_patch(begp:endp,1)
+       call hist_addfld1d (fname='CISUN', units='(Pa)', &
+          avgflag='A', long_name='sunlit leaf intracellular CO2 pressure', &
+          ptr_patch=ptr_1d, default='inactive')
+
+       ptr_1d => this%cisha_z_patch(begp:endp,1)
+       call hist_addfld1d (fname='CISHA', units='(Pa)', &
+          avgflag='A', long_name='shaded leaf intracellular CO2 pressure', &
+          ptr_patch=ptr_1d, default='inactive')
+
+    endif
+
   end subroutine InitHistory
 
   !-----------------------------------------------------------------------