Update aeroFF

zhanglikate · Aug 29, 2024 · da5e6e0 · da5e6e0
1 parent 10e8817
commit da5e6e0
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/physics/MP/Thompson/mp_thompson.F90 b/physics/MP/Thompson/mp_thompson.F90
@@ -1197,13 +1197,20 @@ subroutine get_aero(aero3d,aeroFF, aerfld, nifa, nwfa, ncol, nlev,&
          aero3d(i,k,14)=aerfld(i,k,ntocb)/0.3232698*8*1.e6
          aero3d(i,k,15)=aerfld(i,k,ntocl)/0.3232698*8*1.e6
 
-            do nt =1, 15 
+            do nt =1, 5 
             if (aero3d(i,k,nt) >=1.e-15) then
             aeroFF(i,k,nt)=nifa(i,k)/aero3d(i,k,nt)
             else
             aero3d(i,k,nt) =100.
             endif
             enddo
+            do nt =6, 15
+            if (aero3d(i,k,nt) >=1.e-15) then
+            aeroFF(i,k,nt)=nwfa(i,k)/aero3d(i,k,nt)
+            else
+            aero3d(i,k,nt) =100.
+            endif
+            enddo
          aeroFF(i,k,12)=100.
          aeroFF(i,k,13)=100.
          aeroFF(i,k,14)=100.