Ignore:
Timestamp:
Aug 12, 2016 3:14:59 PM (5 years ago)
Author:
suehring
Message:

Prescribing scalar flux at model top; several bugfixes concering data output of scalars and output of flight data

File:
1 edited

Legend:

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

    r1961 r1992  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Adapted for top_scalarflux
    2222!
    2323! Former revisions:
     
    10171017
    10181018       ELSE
    1019           r_lower = 'q(0)     = q(1)'
     1019          r_lower = 'q(0)      = q(1)'
    10201020       ENDIF
    10211021       IF ( ibc_q_t == 0 )  THEN
    1022           r_upper =  'q(nzt = q_top'
     1022          r_upper =  'q(nzt+1) = q_top'
    10231023       ELSE
    1024           r_upper =  'q(nzt)   = q(nzt-1) + dq/dz'
     1024          r_upper =  'q(nzt+1) = q(nzt) + dq/dz'
    10251025       ENDIF
    10261026       WRITE ( io, 301 ) 'q', r_lower, r_upper
     
    10291029    IF ( passive_scalar )  THEN
    10301030       IF ( ibc_s_b == 0 )  THEN
    1031           r_lower = 's(0)     = s_surface'
     1031          r_lower = 's(0)      = s_surface'
    10321032       ELSE
    1033           r_lower = 's(0)     = s(1)'
     1033          r_lower = 's(0)      = s(1)'
    10341034       ENDIF
    10351035       IF ( ibc_s_t == 0 )  THEN
    1036           r_upper =  's(nzt)   = s_top'
    1037        ELSE
    1038           r_upper =  's(nzt)   = s(nzt-1) + ds/dz'
     1036          r_upper =  's(nzt+1) = s_top'
     1037       ELSEIF ( ibc_s_t == 1 )  THEN
     1038          r_upper =  's(nzt+1) = s(nzt)'
     1039       ELSEIF ( ibc_s_t == 2 )  THEN
     1040          r_upper =  's(nzt+1) = s(nzt) + ds/dz'
    10391041       ENDIF
    10401042       WRITE ( io, 301 ) 's', r_lower, r_upper
     
    10731075          WRITE ( io, 316 )
    10741076       ENDIF
    1075        IF ( ocean  .AND.  constant_top_salinityflux )  THEN
     1077       IF ( ocean  .AND.  constant_top_salinityflux )                          &
    10761078          WRITE ( io, 309 )  top_salinityflux
    1077        ENDIF
    10781079       IF ( humidity       )  WRITE ( io, 315 )
    1079        IF ( passive_scalar )  WRITE ( io, 315 )
     1080       IF ( passive_scalar .AND.  constant_top_scalarflux )                    &
     1081          WRITE ( io, 302 ) top_scalarflux
    10801082    ENDIF
    10811083
     
    20332035313 FORMAT ('       Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
    20342036314 FORMAT ('       Predefined scalar value at the surface')
    2035 315 FORMAT ('       Humidity / scalar flux at top surface is 0.0')
     2037302 FORMAT ('       Predefined constant scalarflux:   ',F9.6,' kg/(m**2 s)')
     2038315 FORMAT ('       Humidity flux at top surface is 0.0')
    20362039316 FORMAT ('       Sensible heatflux and momentum flux from coupled ', &
    20372040                    'atmosphere model')
Note: See TracChangeset for help on using the changeset viewer.