fix rhoa when lowest_temp is 0.0

* provide default value for rhoa when imported temp_height_lowest (Tair) is 0.0 * resolves seg fault when frac_grid=false and do_ca=true
DeniseWorthen · Sep 15, 2021 · 5cb78cd · 5cb78cd
1 parent cd021b5
commit 5cb78cd
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/cicecore/drivers/nuopc/cmeps/ice_import_export.F90 b/cicecore/drivers/nuopc/cmeps/ice_import_export.F90
@@ -582,7 +582,7 @@ subroutine ice_import( importState, rc )
                    rhoa(i,j,iblk) = inst_pres_height_lowest / &
                         (287.058_ESMF_KIND_R8*(1._ESMF_KIND_R8+0.608_ESMF_KIND_R8*Qa(i,j,iblk))*Tair(i,j,iblk))
                 else
-                   rhoa(i,j,iblk) = 0._ESMF_KIND_R8
+                   rhoa(i,j,iblk) = 1.2_ESMF_KIND_R8
                 endif
              end do !i
           end do !j