[WAM-IPE] WAM-IPE r97098: Update the code change from Zhuxiao for src/atmos/p...

[Samuel.Trahan at noaa.gov]

Friendly WAM-IPE developers,

This is an automated email about a WAM-IPE commit.

Project: WAM-IPE
URL: https://svnemc.ncep.noaa.gov/projects/nems/branches/WAM-IPE/quasitrunk
Revision: 97098
Author:   weiyu.yang at noaa.gov
Date:     2017-09-02T19:54:10.561542Z
Message:
Update the code change from Zhuxiao for src/atmos/phys/idea_ion_empirmodels.f


Differences:
Index: checkout/src/atmos/phys/idea_ion_empirmodels.f
===================================================================
--- checkout/src/atmos/phys/idea_ion_empirmodels.f	(revision 96979)
+++ checkout/src/atmos/phys/idea_ion_empirmodels.f	(revision 97098)
@@ -2,8 +2,9 @@
 !r
 !r   chui_model.f       Chiu ionosphere, to return electron density.
 !r                      Converted to run F90, but not changed. mjh
+!r   tiros_ionize_data  added the scaling of qiont with large GW and high tiros_activity_level
+!r                      Zhuxiao Li, 8/31/2017 
 !r
-!r
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
       SUBROUTINE Earth_CHIU_MODEL(sda,sza,thmag,phimr,rlt,rlat,         
      &f107, dip, nday,ht1d,eden3d,ilon,lev1,ht_dim,lon_dim)
@@ -359,6 +360,12 @@
       l = tiros_activity_level - 2
       IF ( l.LT.1 ) l = 1
       IF ( l.GT.7 ) l = 7
+! define dfac to scale qiont later with large GW and tiros_activity_level
+! Added by Zhuxiao.Li
+      dfac = 1.0
+      IF (tiros_activity_level.gt.9 .and. GW.gt.96.0)
+     &   dfac = GW/96.0
+ 
 cc  **
       ri = essa1/18.0 + 11.
       i1 = ri                   ! i1 =int(ri) ?
@@ -511,6 +518,10 @@
      &*WIDTH_maxwell/RANGE_en/E0
   110 CONTINUE
       qiont(i) = qiont(i) + qion_maxwell(i)
+
+!  the qiont scaled by dfac with large GW and tiros_activity level
+!  added by Zhuxiao
+       qiont(i) = qiont(i)*dfac
 !
 !vay-2016, extra security >0 and non-zero
 !
@@ -620,6 +631,11 @@
 !
       IF ( l.LT.1 ) l = 1
       IF ( l.GT.7 ) l = 7
+! define dfac to scale qiont later with large GW and tiros_activity_level
+! Added by Zhuxiao.Li
+      dfac = 1.0
+      IF (tiros_activity_level.gt.9 .and. GW.gt.96.0)
+     &   dfac = GW/96.0
 !
 ! What's this VAY
 !
@@ -802,6 +818,10 @@
   110 CONTINUE   !  loop for300eV as Maxwellian with ch, and eflux
 !
       qiont(i) = qiont(i) + qion_maxwell(i)
+!  the qiont scaled by dfac with large GW and tiros_activity level
+!  added by Zhuxiao
+       qiont(i) = qiont(i)*dfac
+
 !vay-2016
       if (rr.gt.0.) then 
          eden_aurora1D(i)=sqrt(qiont(i)/rr)
-------------- next part --------------
Index: checkout/src/atmos/phys/idea_ion_empirmodels.f
===================================================================
--- checkout/src/atmos/phys/idea_ion_empirmodels.f	(revision 96979)
+++ checkout/src/atmos/phys/idea_ion_empirmodels.f	(revision 97098)
@@ -2,8 +2,9 @@
 !r
 !r   chui_model.f       Chiu ionosphere, to return electron density.
 !r                      Converted to run F90, but not changed. mjh
+!r   tiros_ionize_data  added the scaling of qiont with large GW and high tiros_activity_level
+!r                      Zhuxiao Li, 8/31/2017 
 !r
-!r
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
       SUBROUTINE Earth_CHIU_MODEL(sda,sza,thmag,phimr,rlt,rlat,         
      &f107, dip, nday,ht1d,eden3d,ilon,lev1,ht_dim,lon_dim)
@@ -359,6 +360,12 @@
       l = tiros_activity_level - 2
       IF ( l.LT.1 ) l = 1
       IF ( l.GT.7 ) l = 7
+! define dfac to scale qiont later with large GW and tiros_activity_level
+! Added by Zhuxiao.Li
+      dfac = 1.0
+      IF (tiros_activity_level.gt.9 .and. GW.gt.96.0)
+     &   dfac = GW/96.0
+ 
 cc  **
       ri = essa1/18.0 + 11.
       i1 = ri                   ! i1 =int(ri) ?
@@ -511,6 +518,10 @@
      &*WIDTH_maxwell/RANGE_en/E0
   110 CONTINUE
       qiont(i) = qiont(i) + qion_maxwell(i)
+
+!  the qiont scaled by dfac with large GW and tiros_activity level
+!  added by Zhuxiao
+       qiont(i) = qiont(i)*dfac
 !
 !vay-2016, extra security >0 and non-zero
 !
@@ -620,6 +631,11 @@
 !
       IF ( l.LT.1 ) l = 1
       IF ( l.GT.7 ) l = 7
+! define dfac to scale qiont later with large GW and tiros_activity_level
+! Added by Zhuxiao.Li
+      dfac = 1.0
+      IF (tiros_activity_level.gt.9 .and. GW.gt.96.0)
+     &   dfac = GW/96.0
 !
 ! What's this VAY
 !
@@ -802,6 +818,10 @@
   110 CONTINUE   !  loop for300eV as Maxwellian with ch, and eflux
 !
       qiont(i) = qiont(i) + qion_maxwell(i)
+!  the qiont scaled by dfac with large GW and tiros_activity level
+!  added by Zhuxiao
+       qiont(i) = qiont(i)*dfac
+
 !vay-2016
       if (rr.gt.0.) then 
          eden_aurora1D(i)=sqrt(qiont(i)/rr)