Ignore:
Timestamp:
Feb 22, 2016 3:49:32 PM (6 years ago)
Author:
maronga
Message:

some changes in land surface model, radiation model, nudging and some minor updates

File:
1 edited

Legend:

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

    r1750 r1757  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! Minor fixes.
    2222!
    2323! Former revisions:
     
    358358
    359359!
    360 !--       Invert array and switch sign, then calculate Obukhov length and bulk
    361 !--       Richardson number
    362           z0h_min = MINVAL(z0h)
    363           z0_min  = MINVAL(z0)
     360!--       Save roughness lengths to temporary variables
     361          z0h_min = z0h(nys,nxl)
     362          z0_min  = z0(nys,nxl)
    364363         
    365364!
     
    641640             DO  j = nys, nyn
    642641
    643                 k   = nzb_s_inner(j,i)
    644 
    645642!
    646643!--             If the bulk Richardson number is outside the range of the lookup
     
    10551052             !$acc loop independent
    10561053             DO  j = nysg, nyng
    1057                 e(nzb_s_inner(j,i)+1,j,i) = ( us(j,i) / 0.1_wp )**2
     1054                k = nzb_s_inner(j,i)
     1055                e(k+1,j,i) = ( us(j,i) / 0.1_wp )**2
    10581056!
    10591057!--             As a test: cm = 0.4
    1060 !               e(nzb_s_inner(j,i)+1,j,i) = ( us(j,i) / 0.4_wp )**2
    1061                 e(nzb_s_inner(j,i),j,i)   = e(nzb_s_inner(j,i)+1,j,i)
     1058!               e(k+1,j,i) = ( us(j,i) / 0.4_wp )**2
     1059                e(k,j,i)   = e(k+1,j,i)
    10621060             ENDDO
    10631061          ENDDO
Note: See TracChangeset for help on using the changeset viewer.