[WAM-IPE] WAM-IPE r96233: Commit to update the file .../NEMS/src/atmos/phys/i...

[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: 96233
Author:   weiyu.yang at noaa.gov
Date:     2017-08-02T20:08:29.137471Z
Message:
Commit to update the file .../NEMS/src/atmos/phys/ideaq_ion.f from Zhuxiao.


Differences:
Index: checkout/src/atmos/phys/idea_ion.f
===================================================================
--- checkout/src/atmos/phys/idea_ion.f	(revision 94277)
+++ checkout/src/atmos/phys/idea_ion.f	(revision 96233)
@@ -5,6 +5,9 @@
 ! Nov    2015   VAY  SUB-NE IDEA_ION_INIT(levs) & idea_ion_input.f 
 ! Nov    2016   VAY  Upgrades related to  idea_ion_input.f
 !                                          idea_ion_empirmodels.f   
+! Jul    2017   Zhuxiao Li, Rashid and Tim add the Joule heating factor(JH_fac) to 
+!               include the seasonal variation and semiannual variation of the
+!               Joule Heating.
 !====================================================================
 !=                    GetIonParams                      =
 !
@@ -139,7 +142,8 @@
       REAL  , INTENT(IN)   :: f107, f107d, kp  ! solar-geo inputs from different WAM applications
                                                ! CLIMATE-SWPC-RDATA, controlled in idea_phys
 !
-!      REAL, PARAMETER ::DTR=3.141592653/180.0
+!      REAL, PARAMETER :: DTR=3.141592653/180.0
+      REAL, PARAMETER :: pi = 3.141592653
       INTEGER, INTENT(IN)     :: ix !longitude dim size
       INTEGER, INTENT(IN)     :: im !number of logitude
       INTEGER, INTENT(IN)     :: levs ! number of pres grid
@@ -181,8 +185,10 @@
 !
 ! local
 !
-      real rlt(im),sza(im),jh(ix,levs)                                  
-     &,rinc(5)
+!      real rlt(im),sza(im),jh(ix,levs)                                  
+!     &,rinc(5)
+! define jh_fac, Zhuxiao Li 
+      real rlt(im),sza(im),jh(ix,levs),rinc(5),jh_fac                                 
       INTEGER   i,k
 ! get sza in rad
       sza=acos(cospass)
@@ -205,8 +211,14 @@
      &   dudt,dvdt,jh) 
 ! update to ........ K/sec
       do i=1,im
+!   Joule heating factor to consider the seasonal variation and
+!   semiannual variation, Zhuxiao.Li
+         jh_fac = 2.25+0.75*tanh(2.*rlat(i))*cos((dayno+9.)*2.*pi/365.)                    
+     &   + 0.5*(1.+cos(4.*pi*(dayno-80.)/365.))
       do k=1,levs
-         dtdt(i,k)=jh(i,k)/cp(i,k)
+
+           dtdt(i,k)=jh(i,k)*jh_fac/cp(i,k)
+
       enddo
       enddo
       return
-------------- next part --------------
Index: checkout/src/atmos/phys/idea_ion.f
===================================================================
--- checkout/src/atmos/phys/idea_ion.f	(revision 94277)
+++ checkout/src/atmos/phys/idea_ion.f	(revision 96233)
@@ -5,6 +5,9 @@
 ! Nov    2015   VAY  SUB-NE IDEA_ION_INIT(levs) & idea_ion_input.f 
 ! Nov    2016   VAY  Upgrades related to  idea_ion_input.f
 !                                          idea_ion_empirmodels.f   
+! Jul    2017   Zhuxiao Li, Rashid and Tim add the Joule heating factor(JH_fac) to 
+!               include the seasonal variation and semiannual variation of the
+!               Joule Heating.
 !====================================================================
 !=                    GetIonParams                      =
 !
@@ -139,7 +142,8 @@
       REAL  , INTENT(IN)   :: f107, f107d, kp  ! solar-geo inputs from different WAM applications
                                                ! CLIMATE-SWPC-RDATA, controlled in idea_phys
 !
-!      REAL, PARAMETER ::DTR=3.141592653/180.0
+!      REAL, PARAMETER :: DTR=3.141592653/180.0
+      REAL, PARAMETER :: pi = 3.141592653
       INTEGER, INTENT(IN)     :: ix !longitude dim size
       INTEGER, INTENT(IN)     :: im !number of logitude
       INTEGER, INTENT(IN)     :: levs ! number of pres grid
@@ -181,8 +185,10 @@
 !
 ! local
 !
-      real rlt(im),sza(im),jh(ix,levs)                                  
-     &,rinc(5)
+!      real rlt(im),sza(im),jh(ix,levs)                                  
+!     &,rinc(5)
+! define jh_fac, Zhuxiao Li 
+      real rlt(im),sza(im),jh(ix,levs),rinc(5),jh_fac                                 
       INTEGER   i,k
 ! get sza in rad
       sza=acos(cospass)
@@ -205,8 +211,14 @@
      &   dudt,dvdt,jh) 
 ! update to ........ K/sec
       do i=1,im
+!   Joule heating factor to consider the seasonal variation and
+!   semiannual variation, Zhuxiao.Li
+         jh_fac = 2.25+0.75*tanh(2.*rlat(i))*cos((dayno+9.)*2.*pi/365.)                    
+     &   + 0.5*(1.+cos(4.*pi*(dayno-80.)/365.))
       do k=1,levs
-         dtdt(i,k)=jh(i,k)/cp(i,k)
+
+           dtdt(i,k)=jh(i,k)*jh_fac/cp(i,k)
+
       enddo
       enddo
       return