update to external_ic::remap_scale to skip remapping non-existent IC …

…tracers (#159)
NOAA-GFDL · Oct 14, 2021 · 4823307 · 4823307
1 parent 03e4164
commit 4823307
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/tools/external_ic.F90 b/tools/external_ic.F90
@@ -2853,7 +2853,7 @@ subroutine remap_scalar(Atm, km, npz, ncnst, ak0, bk0, psc, qa, zh, omga, t_in)
 
 ! map tracers
       do iq=1,ncnst
-         if (floor(qa(is,j,1,iq)) > -999) then !skip missing scalars
+        if (floor(qa(is,j,1,iq)) == -1000) cycle !skip missing scalars [floor(-999.99) is -1000]
          do k=1,km
             do i=is,ie
                qp(i,k) = qa(i,j,k,iq)
@@ -2871,7 +2871,6 @@ subroutine remap_scalar(Atm, km, npz, ncnst, ak0, bk0, psc, qa, zh, omga, t_in)
                Atm%q(i,j,k,iq) = qn1(i,k)
             enddo
          enddo
-         endif
       enddo
 
 !---------------------------------------------------