Ignore:
Timestamp:
Mar 4, 2015 5:44:27 PM (6 years ago)
Author:
maronga
Message:

LSM output of r_a and r_s added

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r1552 r1555  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! Added output of r_a and r_s.
    2424!
    2525! Former revisions:
     
    160160    USE land_surface_model_mod,                                                &
    161161        ONLY:   dots_soil, ghf_eb, land_surface, m_soil, nzb_soil, nzt_soil,   &
    162                 qsws_eb, qsws_liq_eb, qsws_soil_eb, qsws_veg_eb, shf_eb,       &
    163                 t_soil
     162                qsws_eb, qsws_liq_eb, qsws_soil_eb, qsws_veg_eb, r_a, r_s,     &
     163                shf_eb, t_soil
    164164
    165165    USE pegrid
     
    717717
    718718             IF ( land_surface )  THEN
    719                 sums_l(nzb,93,tn) = sums_l(nzb,93,tn) + ghf_eb(j,i)
    720                 sums_l(nzb,94,tn) = sums_l(nzb,94,tn) + shf_eb(j,i)
    721                 sums_l(nzb,95,tn) = sums_l(nzb,95,tn) + qsws_eb(j,i)
    722                 sums_l(nzb,96,tn) = sums_l(nzb,96,tn) + qsws_liq_eb(j,i)
    723                 sums_l(nzb,97,tn) = sums_l(nzb,97,tn) + qsws_soil_eb(j,i)
    724                 sums_l(nzb,98,tn) = sums_l(nzb,98,tn) + qsws_veg_eb(j,i)
     719                sums_l(nzb,93,tn)  = sums_l(nzb,93,tn) + ghf_eb(j,i)
     720                sums_l(nzb,94,tn)  = sums_l(nzb,94,tn) + shf_eb(j,i)
     721                sums_l(nzb,95,tn)  = sums_l(nzb,95,tn) + qsws_eb(j,i)
     722                sums_l(nzb,96,tn)  = sums_l(nzb,96,tn) + qsws_liq_eb(j,i)
     723                sums_l(nzb,97,tn)  = sums_l(nzb,97,tn) + qsws_soil_eb(j,i)
     724                sums_l(nzb,98,tn)  = sums_l(nzb,98,tn) + qsws_veg_eb(j,i)
     725                sums_l(nzb,99,tn)  = sums_l(nzb,99,tn) + r_a(j,i)
     726                sums_l(nzb,100,tn) = sums_l(nzb,100,tn)+ r_s(j,i)
    725727             ENDIF
    726728
    727729             IF ( radiation )  THEN
    728                 sums_l(nzb,99,tn)  = sums_l(nzb,99,tn)  + rad_net(j,i)
    729                 sums_l(nzb,100,tn) = sums_l(nzb,100,tn) + rad_sw_in(j,i)
     730                sums_l(nzb,101,tn)  = sums_l(nzb,101,tn)  + rad_net(j,i)
     731                sums_l(nzb,102,tn) = sums_l(nzb,102,tn) + rad_sw_in(j,i)
    730732             ENDIF
    731733
     
    12931295          hom(:,1,91,sr) = sums(:,91)              ! m_soil
    12941296                                                   ! 92 is initial m_soil profile
    1295           hom(:,1,93,sr) = sums(:,93)              ! ghf_eb
    1296           hom(:,1,94,sr) = sums(:,94)              ! shf_eb
    1297           hom(:,1,95,sr) = sums(:,95)              ! qsws_eb
    1298           hom(:,1,96,sr) = sums(:,96)              ! qsws_liq_eb
    1299           hom(:,1,97,sr) = sums(:,97)              ! qsws_soil_eb
    1300           hom(:,1,98,sr) = sums(:,98)              ! qsws_veg_eb
     1297          hom(:,1,93,sr)  = sums(:,93)             ! ghf_eb
     1298          hom(:,1,94,sr)  = sums(:,94)             ! shf_eb
     1299          hom(:,1,95,sr)  = sums(:,95)             ! qsws_eb
     1300          hom(:,1,96,sr)  = sums(:,96)             ! qsws_liq_eb
     1301          hom(:,1,97,sr)  = sums(:,97)             ! qsws_soil_eb
     1302          hom(:,1,98,sr)  = sums(:,98)             ! qsws_veg_eb
     1303          hom(:,1,99,sr)  = sums(:,99)             ! r_a
     1304          hom(:,1,100,sr) = sums(:,100)            ! r_s
     1305
    13011306       ENDIF
    13021307
    13031308       IF ( radiation )  THEN
    1304           hom(:,1,99 ,sr) = sums(:,99)             ! rad_net
    1305           hom(:,1,100,sr) = sums(:,100)            ! rad_sw_in
     1309          hom(:,1,101 ,sr) = sums(:,101)            ! rad_net
     1310          hom(:,1,102,sr)  = sums(:,102)            ! rad_sw_in
    13061311       ENDIF
    13071312
     
    14601465          ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr)           ! qsws_soil_eb
    14611466          ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr)           ! qsws_veg_eb
     1467          ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr)           ! r_a
     1468          ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr)          ! r_s
    14621469       ENDIF
    14631470!
     
    32633270          ts_value(dots_soil+4,sr) = hom(nzb,1,97,sr)           ! qsws_soil_eb
    32643271          ts_value(dots_soil+5,sr) = hom(nzb,1,98,sr)           ! qsws_veg_eb
     3272          ts_value(dots_soil+6,sr) = hom(nzb,1,99,sr)           ! r_a
     3273          ts_value(dots_soil+7,sr) = hom(nzb,1,100,sr)          ! r_s
    32653274       ENDIF
    32663275!
Note: See TracChangeset for help on using the changeset viewer.