Changeset 3773 for palm/trunk/SOURCE


Ignore:
Timestamp:
Mar 1, 2019 8:56:57 AM (3 years ago)
Author:
maronga
Message:

added output of theta_2m_av and minor revisions in palm_csd

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r3655 r3773  
    2525! -----------------
    2626! $Id$
     27! Added output of theta_2m*_xy_av
     28!
     29! 3655 2019-01-07 16:51:22Z knoop
    2730! Implementation of the PALM module interface
    2831!
     
    436439             ENDIF
    437440
     441         CASE ( 'theta_2m*' )
     442             IF ( ALLOCATED( tsurf_av ) ) THEN
     443                DO  i = nxlg, nxrg
     444                   DO  j = nysg, nyng
     445                      pt_2m_av(j,i) = pt_2m_av(j,i) / REAL( average_count_3d, KIND=wp )
     446                   ENDDO
     447                ENDDO
     448                CALL exchange_horiz_2d( pt_2m_av, nbgp )
     449             ENDIF
     450             
    438451          CASE ( 't*' )
    439452             IF ( ALLOCATED( ts_av ) ) THEN
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r3761 r3773  
    2525! -----------------
    2626! $Id$
     27! Added output of theta_2m*_xy_av
     28!
     29! 3761 2019-02-25 15:31:42Z raasch
    2730! unused variables removed
    2831!
     
    920923             ENDIF
    921924
     925          CASE ( 'theta_2m*' )
     926             IF ( ALLOCATED( pt_2m_av ) ) THEN   
     927                DO  i = nxl, nxr
     928                   DO  j = nys, nyn
     929                      match_def = surf_def_h(0)%start_index(j,i) <=            &
     930                                  surf_def_h(0)%end_index(j,i)
     931                      match_lsm = surf_lsm_h%start_index(j,i) <=               &
     932                                  surf_lsm_h%end_index(j,i)
     933                      match_usm = surf_usm_h%start_index(j,i) <=               &
     934                                  surf_usm_h%end_index(j,i)
     935
     936                      IF ( match_def )  THEN
     937                         m = surf_def_h(0)%end_index(j,i)
     938                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
     939                                         surf_def_h(0)%pt_2m(m)
     940                      ELSEIF ( match_lsm  .AND.  .NOT. match_usm )  THEN
     941                         m = surf_lsm_h%end_index(j,i)
     942                        pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
     943                                         surf_lsm_h%pt_2m(m)
     944                      ELSEIF ( match_usm )  THEN
     945                         m = surf_usm_h%end_index(j,i)
     946                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
     947                                         surf_usm_h%pt_2m(m)
     948                      ENDIF
     949                   ENDDO
     950                ENDDO
     951             ENDIF
     952             
     953             
    922954          CASE ( 't*' )
    923955             IF ( ALLOCATED( ts_av ) ) THEN
Note: See TracChangeset for help on using the changeset viewer.