added comments in phys/module_fr_sfire_atm.F

openwfm · Feb 21, 2011 · 89a2e77 · 89a2e77
1 parent 627fa64
commit 89a2e77
Showing 1 changed file with 34 additions and 0 deletions.
diff --git a/wrfv2_fire/phys/module_fr_sfire_atm.F b/wrfv2_fire/phys/module_fr_sfire_atm.F
@@ -316,6 +316,13 @@ subroutine interpolate_atm2fire(id,               & ! for debug output, <= 0 no
 endif
 7001 format(a,' dimensions ',a4,':',i6,' to ',i6,' by ',i6,' to ',i6)
 
+!**********************************************************
+!*                                                        *
+!*           find the altitude of the w points            *
+!*                                                        *
+!**********************************************************
+!! in future, replace by z8w & test if the same
+
     kdmax=kde-1   ! max layer to interpolate from, can be less
 
     do j = jtst,jtet
@@ -339,6 +346,12 @@ subroutine interpolate_atm2fire(id,               & ! for debug output, <= 0 no
 !$OMP END CRITICAL(SFIRE_DRIVER_CRIT)
 endif
 
+!**********************************************************
+!*                                                        *
+!*  interpolate the altitude from w points to u points    *
+!*                                                        *
+!**********************************************************
+
     do j = jtsu,jteu          
       do k=kds,kdmax+1
         do i = itsu,iteu       
@@ -364,6 +377,13 @@ subroutine interpolate_atm2fire(id,               & ! for debug output, <= 0 no
   call message(msg)
 !$OMP END CRITICAL(SFIRE_DRIVER_CRIT)
 endif
+
+!**********************************************************
+!*                                                        *
+!*  interpolate the altitude from w points to v points    *
+!*                                                        *
+!**********************************************************
+
     do j = jtsv,jtev          
       do k=kds,kdmax+1
         do i = itsv,itev       
@@ -385,6 +405,13 @@ subroutine interpolate_atm2fire(id,               & ! for debug output, <= 0 no
 #endif
     logfwh = log(fire_wind_height)
 
+
+!**********************************************************
+!*                                                        *
+!*  interpolate u to the fire mesh at fire_wind_height    *
+!*                                                        *
+!**********************************************************
+
     ! interpolate u, staggered in X
 
     do j = jtsu,jteu            ! compute on domain by 1 smaller, otherwise z0 and ph not available
@@ -415,6 +442,13 @@ subroutine interpolate_atm2fire(id,               & ! for debug output, <= 0 no
       enddo
     enddo
 
+
+!**********************************************************
+!*                                                        *
+!*  interpolate v to the fire mesh at fire_wind_height    *
+!*                                                        *
+!**********************************************************
+
     ! interpolate v, staggered in Y
 
     do j = jtsv,jtev