update module_sf_noahmp_glacier.F90 and module_sf_noahmplsm.F90 accor…

…ding to review

physics/SFC_Models/Land/Noahmp/module_sf_noahmp_glacier.F90

@@ -842,12 +842,14 @@ subroutine radiation_glacier (dt      ,tg      ,sneqvo  ,sneqv   ,cosz   , & !in
 
 ! snow albedos: age even when sun is not present
 
+  if(cosz > 0) then
   if(opt_alb == 1) &
      call snowalb_bats_glacier (nband,cosz,fage,albsnd,albsni)
   if(opt_alb == 2) then
      call snowalb_class_glacier(nband,qsnow,dt,alb,albold,albsnd,albsni)
      albold = alb
   end if
+  end if 
 
 ! zero summed solar fluxes
 

physics/SFC_Models/Land/Noahmp/module_sf_noahmplsm.F90

@@ -3024,7 +3024,7 @@ subroutine albedo (parameters,vegtyp ,ist    ,ice    ,nsoil  , & !in
 
   call snow_age (parameters,dt,tg,sneqvo,sneqv,tauss,fage) 
 
-  if(cosz <= 0) goto 100
+  if(cosz <= 0) 
 
 ! weight reflectance/transmittance by lai and sai
 
@@ -3045,6 +3045,8 @@ subroutine albedo (parameters,vegtyp ,ist    ,ice    ,nsoil  , & !in
      albold = alb
   end if
 
+100 continue
+
 ! ground surface albedo
 
   call groundalb (parameters,nsoil   ,nband   ,ice     ,ist     , & !in
@@ -9143,10 +9145,6 @@ subroutine groundwater(parameters,nsnow  ,nsoil  ,dt     ,sice   ,zsoil  , & !in
 
 ! recharge rate qin to groundwater
 
-!      ka  = hk(iwt)
-! harmonic average, c.he changed based on gy niu's update
-!      ka  = 2.0*(hk(iwt)*parameters%dksat(iwt)*1.0e3) / (hk(iwt)+parameters%dksat(iwt)*1.0e3)
-!     tried one suggesteed by gy niu
       ka =  0.5*(hk(iwt)+parameters%dksat(iwt)*1.0e3)      
 
       wh_zwt  = - zwt * 1.e3                          !(mm)