From da5e6e01e684d1d51b493dc408bc7dafb0f04bd5 Mon Sep 17 00:00:00 2001 From: "Kate.Zhang" Date: Thu, 29 Aug 2024 18:01:17 +0000 Subject: [PATCH] Update aeroFF --- physics/MP/Thompson/mp_thompson.F90 | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/physics/MP/Thompson/mp_thompson.F90 b/physics/MP/Thompson/mp_thompson.F90 index 88feb2c57..fcb6ce21d 100644 --- 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.