Fixes #108

docs/source/conf.py

@@ -78,9 +78,9 @@
 # built documents.
 #
 # The short X.Y version.
-version = u'1.2.0'
+version = u'1.2.1'
 # The full version, including alpha/beta/rc tags.
-release = u'1.2.0'
+release = u'1.2.1'
 
 # The language for content autogenerated by Sphinx. Refer to documentation
 # for a list of supported languages.

falwa/__init__.py

@@ -4,7 +4,7 @@
 Author: Clare Huang, Christopher Polster
 """
 
-__version__ = "1.2.0"
+__version__ = "1.2.1"
 from .interpolate_fields import interpolate_fields
 from .interpolate_fields_direct_inv import interpolate_fields_direct_inv
 from .compute_qref_and_fawa_first import compute_qref_and_fawa_first

falwa/f90_modules/compute_flux_dirinv.f90

@@ -178,19 +178,17 @@ SUBROUTINE compute_flux_dirinv_nshem(pv,uu,vv,pt,tn0,qref,uref,tref,&
     ep2baro(:,:) = ep2baro(:,:)+ep2(:,:)*exp(-zk/h)*dc
     ep3baro(:,:) = ep3baro(:,:)+ep3(:,:)*exp(-zk/h)*dc
 
-    ! Bounds of y-indices in N/SHem
-    if (is_nhem) then  ! 5N and higher latitude
-      jstart = jb+1  ! 6
-      jend = nd
+    if (is_nhem) then
+      do j = jstart,jend
+        ubaro(:,j) = ubaro(:,j)+uu(:,nd-1+j,k)*exp(-zk/h)*dc
+        urefbaro(j) = urefbaro(j)+uref(j-jb,k)*exp(-zk/h)*dc
+      enddo
     else
-      jstart = 1
-      jend = nd-jb   ! nd - 5
+      do j = jstart,jend
+        ubaro(:,j) = ubaro(:,j)+uu(:,j,k)*exp(-zk/h)*dc
+        urefbaro(j) = urefbaro(j)+uref(j,k)*exp(-zk/h)*dc
+      enddo
     endif
-
-    do j = jstart,jend  ! ### yet to be multiplied by cosine
-      ubaro(:,j) = ubaro(:,j)+uu(:,j+nd-1,k)*exp(-zk/h)*dc
-      urefbaro(j) = urefbaro(j)+uref(j-jb,k)*exp(-zk/h)*dc
-    enddo
   enddo
 
 END SUBROUTINE