add BR to output arrays for Noah LSM snow albedo

git-svn-id: https://svn-wrf-model.cgd.ucar.edu/trunk@3794 b0b5d27b-6f0f-0410-a2a3-cb1e977edc3d
openwfm · Apr 7, 2009 · f20f927 · f20f927
1 parent f816c1b
commit f20f927
Showing 1 changed file with 7 additions and 4 deletions.
diff --git a/phys/module_sf_qnsesfc.F b/phys/module_sf_qnsesfc.F
@@ -70,6 +70,7 @@ SUBROUTINE QNSESFC(ITIMESTEP,HT,DZ                                &
      &                 ,LOWLYR,XLAND                                   &
      &                 ,USTAR,ZNT,Z0BASE,PBLH,MAVAIL,RMOL              &
      &                 ,AKHS,AKMS                                      &
+     &                 ,RIB                                            &
      &                 ,CHS,CHS2,CQS2,HFX,QFX,FLX_LH,FLHC,FLQC         &
      &                 ,QGH,CPM,CT                                     &
      &                 ,U10,V10,TSHLTR,TH10,QSHLTR,Q10,PSHLTR          &
@@ -105,6 +106,7 @@ SUBROUTINE QNSESFC(ITIMESTEP,HT,DZ                                &
 !
       REAL,DIMENSION(IMS:IME,JMS:JME),INTENT(INOUT) :: AKHS,AKMS       &
      &                                                ,PBLH,QSFC
+      REAL,DIMENSION(IMS:IME,JMS:JME),INTENT(OUT)   :: RIB
 !
       REAL,DIMENSION(IMS:IME,JMS:JME),INTENT(INOUT) :: QZ0,RMOL,THZ0   &
      &                                                ,USTAR,UZ0,VZ0   &
@@ -283,7 +285,7 @@ SUBROUTINE QNSESFC(ITIMESTEP,HT,DZ                                &
           CALL SFCDIF(NTSD,SEAMASK,THSK(I,J),QSFC(I,J),PSFC            &
      &               ,UZ0(I,J),VZ0(I,J),TZ0,THZ0(I,J),QZ0(I,J)         &
      &               ,USTAR(I,J),ZNT(I,J),Z0BASE(I,J),CT(I,J),RMOL(I,J) &
-     &               ,AKMS(I,J),AKHS(I,J),PBLH(I,J),MAVAIL(I,J)        &
+     &               ,AKMS(I,J),AKHS(I,J),RIB(I,J),PBLH(I,J),MAVAIL(I,J) &
      &               ,CHS(I,J),CHS2(I,J),CQS2(I,J)                     &
      &               ,HFX(I,J),QFX(I,J),FLX_LH(I,J)                    &
      &               ,FLHC(I,J),FLQC(I,J),QGH(I,J),CPM(I,J)            &
@@ -328,7 +330,7 @@ END SUBROUTINE QNSESFC
 !----------------------------------------------------------------------
       SUBROUTINE SFCDIF(NTSD,SEAMASK,THS,QS,PSFC                       &
      &                 ,UZ0,VZ0,TZ0,THZ0,QZ0                           &
-     &                 ,USTAR,Z0,Z0BASE,CT,RLMO,AKMS,AKHS,PBLH,WETM    &
+     &                 ,USTAR,Z0,Z0BASE,CT,RLMO,AKMS,AKHS,RIB,PBLH,WETM &
      &                 ,CHS,CHS2,CQS2,HFX,QFX,FLX_LH,FLHC,FLQC,QGH,CPM &
      &                 ,ULOW,VLOW,TLOW,THLOW,THELOW,QLOW,CWMLOW        &
      &                 ,ZSL,PLOW                                       &
@@ -357,7 +359,7 @@ SUBROUTINE SFCDIF(NTSD,SEAMASK,THS,QS,PSFC                       &
      &                  ,Z0BASE
 !
       REAL,INTENT(OUT) :: CHS,CHS2,CPM,CQS2,CT,FLHC,FLQC,FLX_LH,HFX    &
-     &                   ,PSHLTR,Q02,Q10,QFX,QGH,RLMO,TH02,TH10,U10,V10
+     &                   ,PSHLTR,Q02,Q10,QFX,QGH,RIB,RLMO,TH02,TH10,U10,V10
 !
       REAL,INTENT(INOUT) :: AKHS,AKMS,QZ0,THZ0,USTAR,UZ0,VZ0,Z0,QS
 !----------------------------------------------------------------------
@@ -368,7 +370,7 @@ SUBROUTINE SFCDIF(NTSD,SEAMASK,THS,QS,PSFC                       &
 !
       REAL :: A,B,BTGH,BTGX,CXCHL,CXCHS,DTHV,DU2,ELFC,FCT              &
      &       ,HLFLX,HSFLX,HV,PSH02,PSH10,PSHZ,PSHZL,PSM10,PSMZ,PSMZL   &
-     &       ,RDZ,RDZT,RIB,RLMA,RLMN,RLMP                              &
+     &       ,RDZ,RDZT,RLMA,RLMN,RLMP                                  &
      &       ,RLOGT,RLOGU,RWGH,RZ,RZST,RZSU,SIMH,SIMM,TEM,THM          &
      &       ,UMFLX,USTARK,VMFLX,WGHT,WGHTT,WGHTQ,WSTAR2               &
      &       ,X,XLT,XLT4,XLU,XLU4,XT,XT4,XU,XU4,ZETALT,ZETALU          &
@@ -944,6 +946,7 @@ SUBROUTINE SFCDIF(NTSD,SEAMASK,THS,QS,PSFC                       &
 !
       RLOW=PLOW/(R_D*TLOW)
       CHS=AKHS
+if(i.eq.27.and.j.eq.30)print *,' qnsesfc chs ',i,j,chs
       CHS2=AKHS02
       CQS2=AKHS02
       HFX=-RLOW*CP*HSFLX