<p><b>laura@ucar.edu</b> 2011-11-01 14:59:05 -0600 (Tue, 01 Nov 2011)</p><p>corrected typo in the definition of the season as a function of the Julian day<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_physics_landuse.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_landuse.F        2011-11-01 20:33:46 UTC (rev 1160)
+++ branches/atmos_physics/src/core_physics/module_physics_landuse.F        2011-11-01 20:59:05 UTC (rev 1161)
@@ -5,7 +5,8 @@
 !=============================================================================================
  module module_physics_landuse
  use configure,only: input_landuse_data, &amp;
-                     config_sfc_albedo
+                     config_sfc_albedo,  &amp;
+                     config_frac_seaice
  use dmpar
  use grid_types
 
@@ -83,6 +84,8 @@
 !---------------------------------------------------------------------------------------------
  write(0,*)
  write(0,*) '--- enter subroutine landuse_init_forMPAS: julian day=', julday
+ write(0,*) '--- config_frac_seaice = ',config_frac_seaice
+ write(0,*) '--- xice_threshold     = ',xice_threshold
 
  nCells = mesh % nCells
  latCell  =&gt; mesh % latCell % array
@@ -179,20 +182,16 @@
  DM_BCAST_MACRO(therin)
  DM_BCAST_MACRO(sfhc)
  DM_BCAST_MACRO(scfx)
+ write(0,*) '--- isice   =',isice
+ write(0,*) '--- iswater =',iswater
 
-!defines sea-ice threshold:
- if(frac_seaice == 0) then
-    xice_threshold = 0.5
- elseif(frac_seaice == 1) then
-    xice_threshold = 0.02
- endif
-
 !defines the surface properties over the entire domain:
  do iCell = 1, nCells
 
-    !finds the season as function of julian day (summer=1, winter=2):
+    !finds the season as function of julian day (summer=1, winter=2): summer in the Northern
+    !Hemisphere is defined between March 15th and September 15th (winter otherwise).
     isn = 1
-    if(julday.lt.0.05 .or. julday.gt.288) isn=2
+    if(julday.lt.105 .or. julday.ge.288) isn=2
     if(latCell(iCell) .lt. 0.) isn=3-isn
 
 !   is = nint(ivgtyp(iCell))
@@ -235,7 +234,8 @@
        xland(iCell)  = 1.0
        albbck(iCell) = albd(isice,isn) / 100.
        embck(iCell)  = sfem(isice,isn)
-       if(frac_seaice .eq. 1) then
+       
+       if(config_frac_seaice) then
           !0.08 is the albedo over open water.
           !0.98 is the emissivity over open water.
           sfc_albedo(iCell) = xice(iCell)*albbck(iCell) + (1-xice(iCell))*0.08
@@ -255,7 +255,7 @@
  write(0,*) '--- end subroutine landuse_init_forMPAS'
 
 !formats:
- 101 format(i6,8e15.8)
+ 101 format(i6,8(1x,e15.8))
  102 format(i6,1x,8e15.8)
 
  end subroutine landuse_init_forMPAS

</font>
</pre>