Fix CPU striping bugs due to not resetting data beforehand

geoschem · Jul 28, 2023 · 520fef1 · 520fef1
1 parent 596b28e
commit 520fef1
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/chemistry/mozart/mo_gas_phase_chemdr.F90 b/src/chemistry/mozart/mo_gas_phase_chemdr.F90
@@ -1050,7 +1050,9 @@ subroutine gas_phase_chemdr(lchnk, ncol, imozart, q, &
       State_Grid%AREA_M2(1,i) = real(col_area(i) * Re**2, fp)
       State_Met%AREA_M2(1,i)  = State_Grid%AREA_M2(1,i)
 
-      State_Met%TROPP(1,i)    = pint(i,troplev(i))
+      ! For some reason, setting this to hPa value causes artifacting
+      ! in the model output. Disabling this for now.
+      State_Met%TROPP(1,i)    = pmid(i,troplev(i))! * 0.01e+0_fp ! hPa
 
       do k = 1, pver
         State_Met%PMID(1,i,k) = pmid(i,pver+1-k) * 0.01e+0_fp ! hPa
@@ -1154,6 +1156,7 @@ subroutine gas_phase_chemdr(lchnk, ncol, imozart, q, &
     endif
 
     ! eventually run FAST_JX - photolysis rates will be available in ZPJ
+    ZPJ = 0.0_fp
     call FAST_JX(0, Input_Opt, State_Chm, State_Diag, State_Grid, State_Met, RC)
     if(RC /= GC_SUCCESS) then
         call endrun('gas_phase_chemdr: Fast-JX FAST_JX failed')